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値を含む固定効果?
ありがとうございます。
丸めのため「_正確に」ゼロですか? 'dt(26、22)'は '2.135502e-18'を生成しますか?私はそれを前に 'lmer'と見ました。小さなP値は見積もりが難しいことは分かっています。 –
ありがとう、私は正確なp値が必要なので、信頼限界と効果の大きさを計算できます。うわー、私は前に 'broom'パッケージを見たことがなかった!素晴らしい! – user2716568
2.0e-18は '.Machine $ double.eps'(最小値は1.0 + x!= 1.0)より下ですが、' .Machine $ double.xmin'の上にあります(最小値はx!= 0です)。 .. –