Well, I'm going to reply to my own thread with a solution here, turns out
one attempt we made last week nearly had it, slight adjustment made it work.
For anyone that is interested / in the future wants to achieve the same
thing >
*varcomp <- matrix(nrow=0, ncol=3)
for (i in 1:nlevels(narrow$gene)) {
x<-lme4::VarCorr(rg.lmer[[i]])
varcomp <- rbind(varcomp, c(as.numeric(x$"line:sex"),
as.numeric(x$"line"), attr(x, "sc")^2))}
varcomp<-data.frame(varcomp, row.names=levels(narrow$gene))
colnames(varcomp)<-c("sex.line", "line", "residual")*
--
View this message in context:
http://r.789695.n4.nabble.com/Extract-Variance-Components-tp4629895p4629932.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
