Extraction functions for Lmekin
Extract the fixed effects, random effects, variance of the fixed effects, or variance of the random effects from a linear mixed effects model fit with lmekin.
## S3 method for class 'lmekin' fixef(object, ...) ## S3 method for class 'lmekin' ranef(object, ...) ## S3 method for class 'lmekin' vcov(object, ...) ## S3 method for class 'lmekin' VarCorr(x, ...) ## S3 method for class 'lmekin' logLik(object, ...)
object |
an object inheriting from class |
x |
an object inheriting from class |
... |
some methods for this generic require additional arguments. None are used in this method. |
For the random effects model y = X beta
+ Zb + e, let s^2 be the variance of the error term
e.
Let A = sP be the variance of the random effects
b. There is a computational advantage to solving the problem
in terms of P instead of A, and that is what is
stored in the returned lmekin object.
The VarCorr function returns elements of P; the print
and summary functions report values of A.
Pinhiero and Bates call P the precision factor.
the fixed effects are a vector and vcov returns their variance/covariance matrix. The random effects are a list with one element for each random effect. The ranef component contains the coefficients and VarCorr the estimated variance/covariance matrix. The logLik method returns the loglikelihood along with its degrees of freedom.
Terry Therneau
J Pinheiro and D Bates, Mixed-effects models in S and S-Plus. Springer, 2000.
lmekin, random.effects,
fixed.effects, link{vcov}, VarCorr
data(ergoStool, package="nlme") # use a data set from nlme efit <- lmekin(effort ~ Type + (1|Subject), ergoStool) ranef(efit)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.