Question

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

我刚刚在R中使用了stats类来帮助我在实验室中进行行为实验，所以我非常这个新手。使用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

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

如何才能正确编织？

Answer 1

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

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

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

Answer 2

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

库（lmerTest） lmsocial＆lt; - lmer（社交~grps +刺激+ grps *刺激+（1 |群组）+（1 | subj）） ANOVA（lmsocial）

lmerTest导致控制台但不显示在针织PDF中

2 个答案: