This function computes all the matrices needed to build the precision matrix for spatio-temporal model, as in Lindgren et. al. (2023)
stModel.matrices(smesh, tmesh, model, constr = FALSE)
a list containing needed objects for model definition.
'manifold' to spedify the a string with the model identification
a length three vector with the constants c1
, c2
and c3
the vector d
the matrix T
the model matrices M_1
, ..., M_m
a mesh object over the spatial domain.
a mesh object over the time domain.
a string identifying the model. So far we have the following models: '102', '121', '202' and '220' models.
logical to indicate if the integral of the field over the domain is to be constrained to zero. Default value is FALSE.
See the paper for details.