Displaying 1 result from an estimated 1 matches for "regrexpdat".
2012 Jun 26
1
How to estimate variance components with lmer for models with random effects and compare them with lme results
...mily, method="REML")
However, in some cases, lme does not converge, hence I used lmer from the
lme4 package:
lmer(Trait~1+(Treatment*Source|Family),data=DATA)
Where I extract the variances from the model using the summary function:
model<-lmer(Trait~1+(Treatment*Source|Family),data=regrexpdat)
results<-model at REmat
variances<-results[,3]
I get the same values as with the VarCorr function. I use then these values
to calculate the actual percentage of variation taking the sum as the total
variation.
Where I am struggling is with the interpretation of the results from the
initial...