Extracting (approximate) variance-covariance matrices for fixed-effect parameters, variance components, ratios of variance components to error variance, or the response variable.
# S3 method for varComp
vcov(object, what = c("fixed", "beta", "random", "varComp",
"var.ratio", "tau", "response", "Y"), drop = TRUE, beta.correction=TRUE, ...)An object of class varComp.
A character vector (only the first element will be used) specifying what variance-covariance matrices are requested. "fixed" or "beta" request approximate variance-covariance of fixed-effect parameters (see details). "random" or "varComp" request the approximate variance-covariance matrix of variance components computed from the expected information matrix. "var.ratio" or "tau" requests approximate variance-covariance matrix of ratio of variance components to the error variance computed from the observed information matrix. "response" or "Y" request the marginal variance of the response variable computed from the plug-in estimate.
A logical scalar, indicating whether zero variance components should be dropped from the results.
A logical scalar, only applicable when what='beta', indicating whether the variance-covariance matrix for fixed effect estimates is corrected according to Kackar and Harville (1984).
Place holder.
A numeric matrix of the requested variance-covariance.
For fixed-effect parameters, the results is the plug-in estimate variance of generalized least squares estimates when beta.correction=FALSE; Otherwise, the Kackar and Harville (1984) correction will be used (default). For ratios of variance components to error variance, the result is the Hessian matrix. For response variable, the result is the plug-in estimate of the marginal variance. For variance components, the result is the plug-in estimate of inverse expected information matrix from the restricted likelihood.
Raghu N. Kackar and David A. Harville (1984) Approximations for standard errors of estimators of fixed and random effect in mixed linear models. Journal of the American Statistical Association 79, 853--862
varComp for the varComp object;
KR.varComp for testing fixed effect parameters accounting for uncertainty in variance parameter estimates.