Hi All.
I need a bit of help extracting the residual error variance from the VarCorr
structure from lmer.
#Here's a 2-way random effects model
lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat)
#Get the structure
vc.fit <- VarCorr(lmer.1)
#results in.....
$person
1 x 1 Matrix of class "dpoMatrix"
(Intercept)
(Intercept) 0.7755392
$rater
1 x 1 Matrix of class "dpoMatrix"
(Intercept)
(Intercept) 0.2054469
attr(,"sc")
[1] 0.5051518
#I can pull out the person and rater variance components easy enough. For
example...
vc.person <- [EMAIL PROTECTED]
I'm sure it's simple but I have not been able to grab the residual variance
in the last matrix. I simply wish to asign the residual to a scalar
variable. Any suggestions would be appreciated!
Thanks!
Rick DeShon
--
Rick DeShon
306 Psychology Building
Department of Psychology
Michigan State University
East Lansing, MI 48824-1116
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
