2

我有一个关于 lmerTest 的问题，用于近似线性混合模型的自由度和 p 值。

我刚刚在 R 中上过统计课来帮助我在实验室进行行为实验，所以我对此很陌生。使用 R Studio 时，我可以在安装 lmerTest 库后运行 anova()，我可以在控制台中看到结果，如下图所示：

lmsocial <- lmer(social~ grps + stim + grps*stim + (1|cohort) +(1|subj))
library(lmerTest)
anova(lmsocial)

Analysis of Variance Table of type 3  with  Satterthwaite 
approximation for degrees of freedom
      Sum Sq Mean Sq NumDF DenDF F.value   Pr(>F)   
grps       22471   22471     1    21  5.5922 0.027747 * 
stim       54289   54289     1    22 13.5107 0.001326 **
grps:stim  40423   40423     1    22 10.0599 0.004416 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

然而，下面是我在 R Studio 中编织时在我的 PDF 上得到的读数（与 HTML 中的编织相同）。编织时没有错误，只是缺少信息：

lmsocial <- lmer(social~ grps + stim + grps*stim + (1|cohort) +(1|subj))
library(lmerTest)
anova(lmsocial)

## Analysis of Variance Table
##           Df Sum Sq Mean Sq F value
## grps       1  22471   22471  5.5922
## stim       1  54289   54289 13.5107
## grps:stim  1  40423   40423 10.0599

我错过了什么吗？在尝试生成报告时，最好以 lmerTest 的形式显示生成的 ANOVA 表。

我怎样才能让它正确编织？

4

2 回答 2

1

lmerTest正在做一些偷偷摸摸的事情，我还没弄清楚是什么。同时，只要您添加明确的print()语句，它似乎就可以工作：

library("lmerTest")
fm1 <- lmer(Reaction~Days + (Days|Subject),sleepstudy)
print(a1 <- anova(fm1))

（或print(anova(fm1))，或a1 <- anova(fm1); print(a)）

于 2015-03-19T20:31:10.420 回答
0

lmerTest 包需要在 lmer 模型的规范之前附加。以下应该有效：

图书馆(lmerTest) lmsocial <- lmer(social~ grps + stim + grps*stim + (1|cohort) +(1|subj)) 方差分析(lmsocial)

于 2015-03-19T22:12:49.290 回答