Block spatiotemporal kriging with separable covariance model

Dear list,

is there a way to have block-kriging predictions when using spatio-temporal
with separable covariance model, in particular by using the gstat/spacetime
suite. ?
I saw examples of block kriging when using a 3D metric kriging, with
rescaled time axis, so that a parallelepiped can be set in the krige()
function.

I'd call it a rectangle.

Setting a block in the krigeST does not work, seemingly.

That is correct -- I've now added this to the docs, too (in svn on r-forge).

Thanks for any help,
Piero

	[[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
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
http://www.52north.org/geostatistics      e.pebesma at wwu.de

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130113/12af17cb/attachment.pl>

Dear Edzer,
thanks for the help.
Is block-kriging possible with this kind of covariance model, or is it just
not implemented in the `spacetime` package?

it is just not implemented.

On 13 January 2013 18:31, Edzer Pebesma <edzer.pebesma at uni-muenster.de>wrote:

I saw examples of block kriging when using a 3D metric kriging, with
rescaled time axis, so that a parallelepiped can be set in the krige()
function.

I'd call it a rectangle.

..why would you? Wouldn't you have a 3D xyt block ?

yes, but with right angles.

I don't know what you saw, with the 3D metric model. Specifying a block
in gstat with block = c(dx, dy, dz), you get a rectangular block.

You can estimate block means for arbitrary shapes when you specify a
polygon (only 2D), or a set of points (1, 2 or 3D) that discretizes your
shape, but I've never actually seen people trying to do this for
parallelepipeds.

Bests,

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

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20130113/1e70ce1f/attachment.pl>

Dear,

How can I set the parameter horizontal band and vertical band of 
experimental variogram in gstat package? I can only transforming 
coordinates.

Tiago Vieira.

Dear,

How can I set the parameter horizontal band and vertical band of 
experimental variogram in gstat package? I can only transforming 
coordinates.

Tiago Vieira.