devfun_vp function

Compute deviance of a linear mixed model 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)

Arguments

• 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.

Returns

the REML or ML deviance.

Author(s)

Rune Haubo B. Christensen. Adapted to pbkrtest by Søren Højsgaard.