mehr experimente

devel/fresh/.Rhistory

-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=10, ny=10, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean")
-leman.bk1 <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior.1,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=10, ny=10, trend=0, rsearch = 10, extrap = F,border=leman.bank, duplicate="mean",trend=1)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-gc()
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-library(baykrig)
-q()
 library(baykrig)
 library(baykrig)
 library(baykrig)
@@ -30,3 +10,23 @@ leman.bk1$z
 leman.bk1$snb
 leman.bk1$
 q()
+source("leman-example.R")
+leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean",name.x="x",name.y="y")
+leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean",namx="x",namy="y")
+leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior,duplicate="mean",namx="x",namy="y"))
+leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior,duplicate="mean",namx="x",namy="y")
+save.image()
+leman.bk.0 <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior.0,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")
+leman.prior.0
+rm(leman.prior.0)
+leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean",namx="x",namy="y")
+leman.prior.0
+leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior.0,duplicate="mean",namx="x",namy="y")
+save.image()
+leman.prior.0
+leman.bk.0 <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior.0,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")
+str(leman.bk.0)
+save.image()
+leman.bk.0$z
+gc()
+q()

devel/fresh/leman-example.R

@@ -46,9 +46,9 @@ leman.prior<-empirical.prior(leman.78,cadpbm~x+y,leman.88.vmsph,prior=leman.prio
 # 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=1, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean")
 
-leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean")
+leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean",namx="x",namy="y")
 
-leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior,duplicate="mean")
+leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior.0,duplicate="mean",namx="x",namy="y")
 
 
 # crash:

devel/fresh/leman-example.R~

@@ -35,3 +35,23 @@ plot(leman.88.ev30,var.mod.obj=leman.88.vmsph)
 # leman.krg<-krige.grid(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)
 
 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")
+
+
+# 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=1, rsearch = 10, extrap = T,border=leman.bank, duplicate="mean")
+
+leman.prior.0<-empirical.prior(leman.83,cadpbm~1,leman.88.vmsph,duplicate="mean")
+
+leman.prior.0<-empirical.prior(leman.78,cadpbm~1,leman.88.vmsph,prior=leman.prior,duplicate="mean")
+
+
+# crash:
+#leman.bk.0 <- bk.grid(point = leman.88.pt , at = "cadpbm", prior=leman.prior.0,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")
+
+

src/bkgrid.f

@@ -533,11 +533,11 @@ c           extract results for this tile
                   DOG(I,J)=0
                   ZG(I,J)=0
                   VARG(I,J)=0
-c                  write(*,*)"x"
+                  write(*,*)"x"
                ELSE
                   ZG(I,J)=Z0
                   VARG(I,J)=VAR0
-c                  write(*,*)"o"
+                  write(*,*)"o"
 c                  VARG(I,J)=IWORK(1)*1.0D0
 c                  VARG(I,J)=lambd0
                END IF