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 at data, i, j, ..., drop = FALSE) :
undefined columns selected
how can I plot all three prediction variables and thier variances????
problem in cokriging plot
2 messages · Saman Monfared, Edzer Pebesma
Saman, in communicating these things it would be helpful if you provided
examples that we can reproduce, with data that are available in the
packages.
For instance, I'm sure your command
grided(grd)~c("x","y")
will have caused an error message, because it is plain wrong. Whether
this is a mistake in your email, or whether you simply ignored the error
message (never a good idea) is now unclear.
On 01/11/2012 11:30 AM, Saman Monfared wrote:
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 at 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 at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de