spde.matrices {RGeostats} | R Documentation |
Build and return the Precision matrix
spde.matrices(dbin=NA, dbout=NA, model, triswitch="nqQ", nostat=NA, gext=NA, mesh.list=NA, flag.Q=FALSE, flag.G=TRUE, flag.std=FALSE, flag.blin=FALSE, verbose = FALSE)
dbin |
The |
dbout |
The |
model |
The |
triswitch |
Command line for the internal triangulation step. For more information see
|
nostat |
List of non-stationary parameters.
For details see |
gext |
The 'dbout' may be dilated by gext. This argument designates an array, with its dimension equal to the dimension of the space and which contains the extension in each direction. If not defined, the 'dbout' is not dilated and the simulated results may suffer some edge effect problems. |
mesh.list |
The list containing the meshing information. Check the output of
|
flag.Q |
When TRUE, the Q sparse matrix is output (in triplet format). |
flag.G |
When TRUE, the B and G sparse matrices (in triplet format), as well as the TildeC and Lambda vectors are output. |
flag.std |
When TRUE, the Vars vector is output. |
flag.blin |
When TRUE, the blin vector is output. |
verbose |
Verbose flag. |
Several vectors or matrices used in the SPDE framework are provided such as:
Q The precision matrix (in triplet format)
G The G array
B The B array: C^(-1/2) * G * C^(-1/2)
TildeC The vector C vector (used for anisotropy and range)
Lambda The L vector (used for sill): L = (Lb * C)^(1/2)
Vars The diagonal of the inverse Q matrix
blin The array of coefficients used to construct Q
indCo The list of coarse(1) or fine(0) samples
L The L array (connectivity sparse matrix)
dupl.data Indices of the data samples close to a grid node
dupl.grid Indices of the target points close to a data sample
order Meshing order: 1 (Input,Output,Steiner); 2 (Output,Input,Steiner)
inddat which gives the indices of the data samples in Q
nbnodup which gives the number of samples of Dbin which do not coincide with samples of Dbout.
Note that:
Q is produced if 'flag.Q' is TRUE
G, B, TildeC and Lambda are produced if 'flag.G' is TRUE
Stdv is produced if 'flag.std' is TRUE
blin is produced if 'flag.blin' is TRUE
All matrices are sparse. They are returned in the 'triplet' form. Each non-zero element of the sparse matrix is specified by its row, column and the item value.