Is it a recent R/gstat version? If yes, could you provide me with a
reproducable example?
--
Edzer
Sven Burbeck wrote:
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
_______________________________________________
R-sig-Geo mailing list
[email protected]
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/
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
