Hi. I have some problems in prediction of cokriging 1)- I don't know how make a grid for prediction, I have a polygone and I select 3000 points into this polygone. I want make grid like "meuse.grid"
nc.sids <- readShapePoly("City Boundary.shp") pts <- spsample(nc.sids, n=3000, type="regular") grd<-data.frame(pts) grided(grd)~c("x","y") 2)- I have problem for predict and plot results of cokriging: coordinates(dd) = ~x+y m<-vgm( 0.22493303,"Gau",78006.24,0,cutoff=100607.2, width = 100607.2/12) g <- gstat(NULL, id = " B", form = B ~1,data=dd,fill.cross = F) g <- gstat(g, id = "A", form = A~1,data=dd,fill.cross = F) g <- gstat(g, id = "T", form = T~ 1,data=dd,fill.cross = F) vm <- variogram(g,covariance=F) plot(vm) vm.fit <- fit.lmc(vm, g, model=m,fit.ranges =F) plot(vm, vm.fit) z<-predict(vm.fit,grd) R output: Error in function (classes, fdef, mtable) : unable to find an inherited method for function "proj4string", for signature "data.frame" when I use pts <- spsample(nc.sids, n=3000, type="regular") grd <- as(pts, "SpatialPixels") grided(grd)~c("x1","x2") . . . z<-predict(vm.fit,grd) and when I use spplot(z,contour=T) it is plot not all of predictions and thier variances only plot some of them. when I want plot predictions or covariances image(z["B.var"],col = "beige") R output is: Error in `[.data.frame`(x@data, i, j, ..., drop = FALSE) : undefined columns selected how can I plot all three prediction variables and thier variances???? _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo