Description Usage Arguments Value
Fits a spatio-temporal marked point process model using INLA coupled with the SPDE approch
1 2 3 4 5 6 7 8 | fit.marked.lgcp(mesh = NULL, locs = NULL, t.index = NULL, mark = NULL,
covariates = NULL, mark.family = "gaussian", verbose = FALSE,
prior.rho = list(theta = list(prior = "pccor1", param = c(0, 0.9))),
prior.range = c(400, 0.5), prior.sigma = c(1, 0.05), hyper = list(theta
= list(prior = "normal", param = c(0, 10))), control.compute = list(dic =
TRUE, waic = TRUE, cpo = TRUE, config = TRUE), control.inla = list(strategy
= "gaussian", int.strategy = "eb"), link = NULL, pp.int = FALSE,
mark.int = FALSE, ...)
|
mesh |
delauney triangulation of area, an object returned by make.mesh is suitable. |
locs |
a matrix of |
t.index |
a vector of length |
mark |
a vector of length |
covariates |
a matrix of covariates (for the mark) with named columns. |
mark.family |
assumed likelihood for mark, by defalt "gaussian". Ignored if |
verbose |
Logical, if |
prior.range |
pc prior for the range of the latent field (rnage0,Prange) (i.e., P(range < range0) = Prange NOTE should be changed to reflect range of the domain by default this is 400, see Martins, Thiago G., et al. "Penalising model component complexity: A principled, practical approach to constructing priors." (2014) for more details. |
prior.sigma |
pc prior for the sd of the latent field (sigma0,Psigma) by default c(1,0.05) i.e., prob sigma > 1 = 0.05 |
hyper |
prior for the copy parameter by default is a N(0,10) i.e., list(theta=list(prior='normal', param=c(0,10))) |
control.compute |
a list of fit statistics the user wants INLA to return. By default this
is |
control.inla |
a list which controls the fitting procedures INLA uses see Blangiardo, Marta, et al.
"Spatial and spatio-temporal models with R-INLA." Spatial and spatio-temporal epidemiology 4 (2013)
by default this is |
pp.int |
Logical, should an intercept term for the lgcp be included. By default FALSE. If TURE alpha0 estimated. |
mark.int |
Logical, should an interpect term for the mark be included. By default FALSE. If TRUE beta0 estimated. |
... |
other |
An inla model fit object,
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.