[R-sig-Geo] Trouble with universal kriging
Dominik.Schneider at colorado.edu
Fri Jun 21 22:03:48 CEST 2013
Thanks for the heads up on +proj. I had not noticed that difference.
I read chapter 8, twice now. It is not clear to me what the error message
I'm getting means:
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 5742 rows but variable(s) found have 200 rows
2: 'newdata' had 5742 rows but variable(s) found have 200 rows
I thought newdata was supposed to be the new grid on which to predict.
Here is code for working data.
recon[1:50]=snotel[1:50] # i'm not sure how to introduce spatial correlation
in this random data.
coordinates(swe.grid) <- ~ x+y
gridded(swe.grid) <- TRUE
R version 2.15.2 (2012-10-26)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
attached base packages:
 grid stats graphics grDevices utils datasets methods
other attached packages:
 raster_2.1-25 maptools_0.8-23 lattice_0.20-13 foreign_0.8-52
loaded via a namespace (and not attached):
 intervals_0.14.0 spacetime_1.0-4 tools_2.15.2 xts_0.9-3
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Trouble-with-universal-kriging-tp7583848p7583875.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
More information about the R-sig-Geo