2016-05-20 2 views
1

summary()の固定効果部分をdata.frameとして抽出したいと思います。私は、次のモデルを実行するためにlme4を使用しています:lme4から要約の一部を抽出する方法は?

SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy) 
summary(SleepStudy) 

私は次のように使用してsummaryのランダムなエフェクト・セクションを抽出できることを知っている:

SleepStudy_RE <- as.data.frame(VarCorr(SleepStudy)) 

はのためのコードのような行があります推定値、標準誤差、自由度、正確なp値を含む固定効果?

ありがとうございます。

答えて

4

する必要があります。

library(lmerTest) 
SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy) 
coef(summary(SleepStudy)) 
##    Estimate Std. Error  df t value Pr(>|t|) 
## (Intercept) 251.40510 9.7467163 22.8102 25.79383  0 
## Days   10.46729 0.8042214 161.0036 13.01543  0 

私はp値が正確ゼロである理由を知りません:あなたはp値が必要な場合

library(lme4) 
SleepStudy <- lmer(Reaction ~ Days + (1|Subject), data = sleepstudy) 
coef(summary(SleepStudy)) 
##    Estimate Std. Error t value 
## (Intercept) 251.40510 9.7467163 25.79383 
## Days   10.46729 0.8042214 13.01543 

あなたは(あなたがモデルを再フィットする必要があります)lmerTest必要この場合;メンターがlmerTestといっしょに取るべきことかもしれません。

また、broom packageに興味があります。

+0

丸めのため「_正確に」ゼロですか? 'dt(26、22)'は '2.135502e-18'を生成しますか?私はそれを前に 'lmer'と見ました。小さなP値は見積もりが難しいことは分かっています。 –

+0

ありがとう、私は正確なp値が必要なので、信頼限界と効果の大きさを計算できます。うわー、私は前に 'broom'パッケージを見たことがなかった!素晴らしい! – user2716568

+0

2.0e-18は '.Machine $ double.eps'(最小値は1.0 + x!= 1.0)より下ですが、' .Machine $ double.xmin'の上にあります(最小値はx!= 0です)。 .. –

関連する問題