[gstat] simple reboxing test fails
On 11/07/2012 04:39 AM, Tom Roche wrote:
...
But all my predicted box values=NA; how to fix? details:
...
# Here I fail to get the values for the output grid :-( data3D.out <- krige(formula=v ~ 1, data3D.in, grid3D.out) summary(data3D.out) # gets #> Object of class SpatialGridDataFrame #> Coordinates: #> min max #> [1,] 0 8 #> [2,] 0 8 #> [3,] 0 8 #> Is projected: NA #> proj4string : [NA] #> Grid attributes: #> cellcentre.offset cellsize cells.dim #> 1 0.5 1 8 #> 2 0.5 1 8 #> 3 4.0 8 1 #> Data attributes: #> var1.pred var1.var #> Min. :1 Min. : NA #> 1st Qu.:1 1st Qu.: NA #> Median :1 Median : NA #> Mean :1 Mean :NaN #> 3rd Qu.:1 3rd Qu.: NA #> Max. :1 Max. : NA #> NA's :64 What should I call to get the values for the output grid boxes?
Did you notice the predictions in var1.pred values are all 1? (inverse distance always returns NA values for the prediction variance var1.var)
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