spatialfusion-package | R Documentation |
Multivariate modelling of geostatistical (point), lattice (areal) and point pattern data in a unifying spatial fusion framework. Details are given in Wang and Furrer (2021) <doi: 10.1016/j.csda.2021.107240>. Model inference is done using either 'Stan' <https://mc-stan.org/> or 'INLA' <https://www.r-inla.org/>.
Package: | spatialfusion |
Type: | Package |
Version: | 0.6-6 |
Date: | 2022-08-14 |
License: | GPL (>= 3) |
LazyLoad: | yes |
fusionData()
is used to set up the data structure needed for spatial fusion modelling. Depending on the chosen 'method'
, either a dstan
or a dinla
object is returned. This object is then suppied to the 'data'
argument in fusion()
for fitting a spatial fusion model.
In terms of the 'method'
, Stan provides Hamiltonian Monte Carlo-based full Bayesian inference, while INLA provides approximate Bayesian inferece at a much faster computation speed. Their results are very similar in our simulation studies (Wang and Furrer, 2021 <doi: 10.1016/j.csda.2021.107240>).
IMPORTANT: Users should be familiar with either rstan or INLA packages themselves. For Stan, users should know how to choose priors appropriately. For INLA, users should know how to set up an appropriate mesh.
fusion()
is used to fit a spatial fusion model. The most related publication is by Wang and Furrer (2021) <doi: 10.1016/j.csda.2021.107240>, which introduced the framework.
We suggest users to test their model on smaller sub-sampled dataset first, to check model fitting issues such as convergence, identifiability etc. It also helps to get an idea of the computation time required. Afterwards, users can fit the model to their full dataset. The output has a class fusionModel
.
Common generic functions such as fitted()
, predict()
, summary()
and plot()
are available for fusionModel
objects. Diagnostics of spatial fusion models should be done in the same way as for a Stan or a INLA model, depending on the chosen method.
Craig Wang craigwang247@gmail.com
## Not run: ## Citations citation('spatialfusion') ## Vignette: short demo vignette("spatialfusion_vignette", package = "spatialfusion") ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.