我有几个适合lmer
( lme4 包)的随机效果模型。例如:
a.1=lmer(songTrait ~ 1 + (1|Year) + (1|ID), REML=F)
a.2=lmer(songTrait ~ Date + (1|Year) + (1|ID), REML=F)
我想打印这些模型的表格摘要,理想情况下,类似于你得到的打印
a.1
或
show(a.1)
此代码......
pdf("c:/temp/x.pdf")
textplot(anova(a.1,a.2,method="ML"),cex=.8,halign="center",cmar=2)
dev.off()
可以很好地打印模型比较,但我还想打印每个模型的详细摘要(因为我有18个songTrait
变量我正在拟合模型)。不幸的是,textplot
( gplots 包)无法处理S4输出,我还没有找到一个好的解决方法。任何帮助将不胜感激!