2

線形混合モデルの自由度と p 値を近似するための lmerTest について質問があります。

ラボでの行動実験を支援するために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

何か不足していますか？レポートを生成しようとするとき、結果の ANOVA テーブル フォーム lmerTest を表示すると便利です。

どうしたら綺麗に編めるようになりますか？

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 に答える