rr_em: Estimate the parameters of a spatial mixed effects model

In jelsema/RRSM: Methods for reduced rank spatial models

Description

Computes estimates of parameters for spatial mixed effects models using an EM algorithm.

Usage

rr_em(Y, X, S, epsilon = 1e-05, max_iter = 500, ...)

rr_em_constrained(Y, X, X2 = NULL, S, cone = NULL, epsilon = 1e-05,
  max_iter = 500, mySolver = "LS", ...)

Arguments

Y	the observed data.
X	the fixed-effects design matrix.
S	the matrix of basis functions.
epsilon	the stopping condition for EM algorithms.
max_iter	for EM algorithms, the maximum number of iterations.
...	space for additional arguments to be passed to the estimation function.
X2	for constrained EM algorithm, optional design matrix for unconstrained fixed-effects.
cone	for constrained EM algorithm, matrix of -1, 0, and 1's defining the order restriction.
mySolver	the solver to be used in the order restriction (passed to isotone::activeSet.

Note

Some of the estimation functions do not estimate the fixed effects parameter (e.g., rr_em_nobeta). In this case, bhat in the return object is NULL. In the estimation functions of RRSM, if bhat is NULL, then a GLS estimate will be computed automatically.

See Also

rr_est mom_bec mom_fit rr_universal_krige

jelsema/RRSM documentation built on May 19, 2019, 4:02 a.m.