Commit 662c81c4 authored by agebhard's avatar agebhard

??

parent 2c15e551
......@@ -367,7 +367,8 @@ bk.grid <- function(point.obj,
npr<-prior$n
typpr<-prior$info
typpr[prior$type=="subjective"]<-typpr[prior$type=="subjective"]*(-1)
browser()
ans<-.C("bk_grid",
xsw = as.double(xsw),
ysw = as.double(ysw),
......
No preview for this file type
leman.prior.ord
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=5, ny=5, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
plot.krige.map(leman.bk.ord)
str(leman.bk.ord$var)
str(leman.bk.ord$z)
leman.bk.ord$z
leman.bk.ord$var
image(leman.bk.ord$z)
contour(leman.bk.ord$z)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=50, ny=50, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean"); save.image()
save.image()
plot.krige.map(leman.bk.ord)
contour(leman.bk.ord$z)
image(leman.bk.ord$z)
image(leman.bk.ord$var)
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
system("echo $DISPLAY")
x11(display="zeta.home.algepop.net:10.0")
image(leman.bk.ord$z)
lines(leman.bank)
plot.krige.map(leman.bk.ord$z)
plot.krige.map(leman.bk.ord)
library(rgeostat);library(baykrig);
str(ans)
library(rgeostat);library(baykrig);
str(ans)
str(leman.bk.ord)
str(ans)
str(leman.bk.ord)
str(ans)
library(rgeostat);library(baykrig);
str(ans)
str(leman.bk.ord)
str(leman.bk.ord)
rm(leman.bk.ord)
rm(leman.bk.ord)
str(leman.bk.ord)
library(rgeostat);library(baykrig);leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean"); save.image()
str(ans)
str(leman.bk.ord)
q()
library(rgeostat);library(baykrig);
library(rgeostat);library(baykrig);
library(rgeostat);library(baykrig);
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=25, ny=25, trend=0, rsearch = 15, extrap = F,border=leman.bank, duplicate="mean");
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=20, ny=20, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean");
ls(ans)
ls()
library(rgeostat);library(baykrig);
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=5, ny=5, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean");
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=5, ny=5, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean");
library(rgeostat);library(baykrig);
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=15, ny=15, trend=0, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean");
leman.bk.ord <- bk.grid(point = leman.pt, at = "cadpbm", prior=leman.prior.ord,var.mod.obj = leman.sph, xsw=min(leman.bank$x),ysw=min(leman.bank$y), xne=max(leman.bank$x), yne=max(leman.bank$y), nx=5, ny=5, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean");
#include <R.h>
#include <Rinternals.h>
// #define TRANSIENT
void bk_grid(double *xsw,
double *ysw,
double *xne,
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment