[R-sig-Geo] error from universal kriging with krige()

Sven Burbeck sburbeck at bfs.de
Wed Jul 23 18:36:53 CEST 2008


Hello,

when I try to apply universal kriging by krige() I produce the following
error:
>myVg_2 <- variogram(log10(M_BL) ~  GEO_UNIT - 1, data = myBlt_utm)
>myVgm <- vgm(psill = 0.05,model = "Sph",range = 25000,nugget = 0.1)
>myVgmf <- fit.variogram(object = myVg_2,model = myVgm,fit.sills = TRUE,
fit.ranges = TRUE,fit.method = 7)
>myBlt_grid_UK <- krige(log10(M_BL) ~  GEO_UNIT - 1, locations =
myBlt_utm, newdata = myPointGrid_utm, model = myVgmf)
Error in gstat.formula.predict(d$formula, newdata, na.action = na.action,  :
  NROW(locs) != NROW(X): this should not occur
In addition: Warning messages:
1: 'newdata' had 2701 rows but variable(s) found have 911 rows
2: 'newdata' had 2701 rows but variable(s) found have 911 rows


I know there was a posting about this error. There it was a coordinate
problem but still I don't know how to solve the problem.

Can anybody give me a hint?

Thanks in advance,
Sven




More information about the R-sig-Geo mailing list