Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

getGB

Giesbrecht & Burns Approximation of the Variance-Covariance Matrix of Variance Components


Description

Compute variance covariance matrix of variance components of a linear mixed model via the method stated in Giesbrecht and Burns (1985).

Usage

getGB(obj, tol = 1e-12)

Arguments

obj

(object) with list-type structure, e.g. VCA object fitted by ANOVA or a premature VCA object fitted by REML

tol

(numeric) values < 'tol' will be considered being equal to zero

Details

This function is not intended to be called by users and therefore not exported.

Value

(matrix) corresponding to the Giesbrecht & Burns approximation of the variance-covariance matrix of variance components

Author(s)

References

Searle, S.R, Casella, G., McCulloch, C.E. (1992), Variance Components, Wiley New York

Giesbrecht, F.G. and Burns, J.C. (1985), Two-Stage Analysis Based on a Mixed Model: Large-Sample Asymptotic Theory and Small-Sample Simulation Results, Biometrics 41, p. 477-486

See Also

Examples

## Not run: 
data(dataEP05A2_3)
fit <- anovaVCA(y~day/run, dataEP05A2_3)
fit <- solveMME(fit)		# some additional matrices required
getGB(fit)

## End(Not run)

VCA

Variance Component Analysis

v1.4.3
GPL (>= 3)
Authors
Andre Schuetzenmeister [aut, cre], Florian Dufey [aut]
Initial release
2020-06-17

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.