lmerTest在控制台中的结果，但在针织PDF中没有显示。

Question

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

我只是上了R的统计学课程，以帮助我在实验室做行为实验，所以我对此非常陌生。使用results时，我可以在挂载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

lmerTest包需要在lmer模型的规范之前附加。下列措施应能发挥作用：

图书馆(LmerTest) lmsocial <- lmer(social~ grps + stim + grps*stim +(1连队列)+(1区subj)) anova(lmsocial)