A function to create a regular grid over an observation window in order to model the spatial randome effects as a Gaussian
Markov random field.
Usage
getgrd(shape, cellwidth)
Arguments
shape
an object of class SpatialPolygons or SpatialPolygonsDataFrame
cellwidth
a scalar, the width of the grid cells
Value
a SpatialPolygons object: the grid on which prediction of the spatial effects will occur
References
Benjamin M. Taylor. Auxiliary Variable Markov Chain Monte Carlo for Spatial Survival and Geostatistical Models. Benjamin M. Taylor. Submitted.http://arxiv.org/abs/1501.01665
Finn Lindgren, Havard Rue, Johan Lindstrom. An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach. Journal of the Royal Statistical Society: Series B 73(4)