Linear part of RMmodel
RFlinearpart returns the linear part of a model
RFlinearpart(model, x, y = NULL, z = NULL, T = NULL, grid=NULL,
data, params, distances, dim, set=0, ...)model,params |
object of class |
x |
vector of x coordinates, or object of class |
y,z |
optional vectors of y (z) coordinates, which should not be given if |
T |
optional vector of time coordinates, |
grid |
logical; the function finds itself the correct value in nearly all cases, so that usually |
distances,dim |
another alternative for the argument |
data |
matrix, data.frame or object of class |
set |
integer. See section Value for details. |
... |
for advanced use: further options and control arguments for the simulation that are passed to and processed by |
RFlinearpart returns a list
of three components, Y, X, vdim returning
the deterministic trend, the design matrix, and the multivariability,
respectively.
If set is positive, Y and X contain
the values for the set-th set of coordinates.
Else, Y and X are both lists containing
the values for all the sets.
In the linear part of the model specification the parameters
that are NA must be the first model part. I.e.
NA * sin(R.p(new="isotropic")) + NA + R.p(new="isotropic")
is OK, but not
sin(R.p(new="isotropic")) * NA + NA + R.p(new="isotropic")
Martin Schlather, schlather@math.uni-mannheim.de, https://www.wim.uni-mannheim.de/schlather/
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set ## RFoptions(seed=NA) to make them all random again x <- seq(0, pi, len=10) trend <- 2 * sin(R.p(new="isotropic")) + 3 model <- RMexp(var=2, scale=1) + trend print(RFlinearpart(model, x=x)) ## only a deterministic part trend <- NA * sin(R.p(new="isotropic")) + NA + R.p(new="isotropic") / pi model <- RMexp(var=NA, scale=NA) + trend print(RFlinearpart(model, x=x))
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.