Message-ID: <509A1160.4010101@uni-muenster.de>
Date: 2012-11-07T07:44:32Z
From: Edzer Pebesma
Subject: [gstat] simple reboxing test fails
In-Reply-To: <87vcdi3wuv.fsf@pobox.com>
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