Commit b890d4c0 authored by agebhard's avatar agebhard
Browse files

T -> TRUE

parent 6297df40
......@@ -45,7 +45,7 @@
n <- dim(x)[1]
}
else{
ord <- (hist(idup,plot=F,freq=T,breaks=seq(0.5,max(idup)+0.5,1))$counts==1)
ord <- (hist(idup,plot=FALSE,freq=TRUE,breaks=seq(0.5,max(idup)+0.5,1))$counts==1)
x <- x[ord,]
n <- dim(x)[1]
}
......@@ -57,7 +57,7 @@
if(!is.null(var.mod)){
c0 <- var.mod$parameters[1]+var.mod$parameters[2]
xx.dist<-as.matrix(dist(cbind(x[,namx], x[,namy]),
diag=T,upper=T))
diag=TRUE,upper=TRUE))
covmat <- c0-var.mod$model(xx.dist,var.mod$parameters)
}
else
......
......@@ -82,7 +82,7 @@ require(rgeostat)
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)
rsearch = 7.5,extrap = F,border=leman.bank, get.lm=TRUE)
plot(leman.krg,show.snb=T)
}
\keyword{ baykrig }
......@@ -57,7 +57,7 @@ empirical.prior(x, formula=NULL, var.mod=NULL, prior=NULL, method="gqr", namx=NU
leman.88.ev30<-est.variogram(leman.88.pt,leman.88.pr.30,"cadpbm")
# fit a sperical variogram
if(length(ls(pat="leman.88.vmsph"))==0)
leman.88.vmsph<-fit.variogram("spherical",leman.88.ev30,0.1,0.2,25,plot.it=T)
leman.88.vmsph<-fit.variogram("spherical",leman.88.ev30,0.1,0.2,25,plot.it=TRUE)
plot(leman.88.ev30,var.mod.obj=leman.88.vmsph)
......
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