Compute Deviance of an LMM as a Function of Variance Parameters
This function is used for extracting the asymptotic variance-covariance matrix of the variance parameters.
devfun_vp(varpar, devfun, reml)
varpar
: variance parameters; varpar = c(theta, sigma)
.devfun
: deviance function as a function of theta only.reml
: if TRUE
the REML deviance is computed; if FALSE
, the ML deviance is computed.the REML or ML deviance.
Rune Haubo B. Christensen