diff --git a/devel/fresh/leman-example.R b/devel/fresh/leman-example.R index bad0374738e85e906f3023a10b43e211879ee40b..fe3738620c214eab3601146801802d1dd7e22ae3 100644 --- a/devel/fresh/leman-example.R +++ b/devel/fresh/leman-example.R @@ -36,8 +36,13 @@ plot(leman.88.ev30,var.mod.obj=leman.88.vmsph) leman.krg<-krige.grid.c(min(leman.bank$x),min(leman.bank$y), max(leman.bank$x), max(leman.bank$y), nx=100, ny=100, point = leman.88.pt, at = "cadpbm", var.mod.obj = leman.88.vmsph, rsearch = 7.5, extrap = F,border=leman.bank) +plot(leman.krg,show.snb=T) + leman.prior<-empirical.prior(leman.83,cadpbm~x+y,leman.88.vmsph,duplicate="mean") leman.prior<-empirical.prior(leman.78,cadpbm~x+y,leman.88.vmsph,prior=leman.prior,duplicate="mean") -leman.bk <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior,var.mod.obj = leman.88.vmsph, 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") + +# crash: +#leman.bk <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior,var.mod.obj = leman.88.vmsph, 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") +