This function is used for extracting the asymptotic variance-covariance matrix of the variance parameters.
devfun_vp(varpar, devfun, reml)
the REML or ML deviance.
variance parameters; varpar = c(theta, sigma)
.
deviance function as a function of theta only.
if TRUE
the REML deviance is computed;
if FALSE
, the ML deviance is computed.
Rune Haubo B. Christensen. Adapted to pbkrtest by Søren Højsgaard.