fit.marked.lgcp: Fitting a marked point process model

Description Usage Arguments Value

View source: R/m_pp_fit.r

Description

Fits a spatio-temporal marked point process model using INLA coupled with the SPDE approch

Usage

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, ...)

Arguments

mesh

delauney triangulation of area, an object returned by make.mesh is suitable.

locs

a matrix of nrow locations in ncol dimesions.

t.index

a vector of length nrow of time index units refering to each point location given in locs.

mark

a vector of length nrow of marks refering to each point location

covariates

a matrix of covariates (for the mark) with named columns.

mark.family

assumed likelihood for mark, by defalt "gaussian". Ignored if mark is NULL.

verbose

Logical, if TRUE, model fitting is output the console.

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 list(dic = TRUE, waic = TRUE,cpo = TRUE, config = TRUE).

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 list(strategy='gaussian',int.strategy = 'eb') for quick and dirty fitting.

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 inla arguments

Value

An inla model fit object,


cmjt/lgcpSPDE documentation built on July 25, 2019, 3:05 p.m.