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.dbin=TRUE, mesh.dbout=TRUE, mesh.list=NA, flag.Q=FALSE, flag.G=TRUE, flag.std=FALSE, flag.blin=FALSE, flag.sparse=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.dbin |
When TRUE, the location corresponding to the Input Data are systematically included in the meshing |
mesh.dbout |
When TRUE, the location corresponding to the Output Targets are systematically included in the meshing |
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. |
flag.sparse |
When TRUE (if the package 'Matrix' is loaded), the resulting matrices are returned as a sparse matrix. |
verbose |
Verbose flag. |
Several vectors or matrices used in the SPDE framework are provided such as:
Q The precision matrix
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
The matrices 'Q', 'G' and 'B'are sparse. They are returned in the 'triplet' form if 'flag.sparse' is FALSE. Each non-zero element of the sparse matrix is specified by its row, column and the item value. If 'flag.sparse'=TRUE, the matrics are returned as sparse matrices (of the package 'Matrix').