[R-sig-Geo] [gstat] simple reboxing test fails

[Edzer Pebesma]

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