From fa85a3a9f5c3071b9f62e1fb844a313ff77bb7e9 Mon Sep 17 00:00:00 2001 From: agebhard <> Date: Sun, 19 May 2002 15:11:20 +0000 Subject: [PATCH] mehr experimente --- devel/fresh/.Rhistory | 40 ++++++++++++++++++------------------ devel/fresh/leman-example.R | 4 ++-- devel/fresh/leman-example.R~ | 20 ++++++++++++++++++ src/bkgrid.f | 4 ++-- 4 files changed, 44 insertions(+), 24 deletions(-) diff --git a/devel/fresh/.Rhistory b/devel/fresh/.Rhistory index 53d665e..5466b4f 100644 --- a/devel/fresh/.Rhistory +++ b/devel/fresh/.Rhistory @@ -1,23 +1,3 @@ -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() diff --git a/devel/fresh/leman-example.R b/devel/fresh/leman-example.R index 8098b74..a4f8d51 100644 --- a/devel/fresh/leman-example.R +++ b/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: diff --git a/devel/fresh/leman-example.R~ b/devel/fresh/leman-example.R~ index 9a02785..8098b74 100644 --- a/devel/fresh/leman-example.R~ +++ b/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") + + diff --git a/src/bkgrid.f b/src/bkgrid.f index f2c4647..02a9e68 100644 --- a/src/bkgrid.f +++ b/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 -- GitLab