At the optimum the covariance is available as vcov(lmer-model)
. This function
computes cov(beta)
at non (RE)ML estimates of varpar
.
get_covbeta(varpar, devfun)
The covariances matrix of the fixed effects at supplied varpar
values.
variance parameters; varpar = c(theta, sigma)
.
deviance function as a function of theta only.
Rune Haubo B. Christensen. Adapted to pbkrtest by Søren Højsgaard.