Displaying 1 result from an estimated 1 matches for "nonug".
Did you mean:
nobug
2009 Jan 29
0
assign point values as labels in spplot
...X+Y
pred.dat<-read.table("C:/temp/Predict_Location_XY.txt", header = FALSE)
names(pred.dat)<-c("x","y","p")
my.pred.dat<-data.frame(pred.dat)
coordinates(my.pred.dat)<-~x+y
gridded(my.pred.dat) = TRUE
#Set up 2 exponential variograms
lzm.vgm.exp.noNug<-vgm(1.1,"Exp",2000,0)
lzm.vgm.exp.Nug<-vgm(0.7,"Exp",3000,0.4)
#Krige the exponential variograms
lzm.krig.exp.noNug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.noNug)
lzm.krig.exp.Nug<-krige(Kh_Log10~1,my.dat,my.pred.dat,model=lzm.vgm.exp.Nug)
lzm.krig....