lmer從anova獲取p值

lmerTest被設計爲允許使用Satterthwaite分母自由度（ddf）估計從lmer混合模型分析估計p-values的包裝。但lmerTest現在似乎被打破。它現在返回一條消息，指出存在內部計算錯誤，並只返回lmer結果（沒有p-values）。我可以用summary（）函數計算p-values，使用Dan Mirman的優秀代碼計算Kenward-Rogers估計的ddf。但是我找不到相同的代碼來計算lmer模型的anova調用中的p-values。我懷疑只需要餵食anova() ddf，但我無法弄清楚如何做到這一點。lmer從anova獲取p值

在此先感謝任何能爲此問題提出解決方案的人。每當在得到薩特斯韋特的近似（諸如例如在計算漸近方差協方差矩陣）發生某些計算誤差

拉里Hunsicker

來源

2016-07-11 Larry Hunsicker

lmerTest返回lme4封裝的ANOVA輸出。 lmerTest沒有被破壞，只是當Satterthwaite的逼近不能被計算時可能會有例子。根據我的經驗，這種情況並不常見。

來源

2016-07-12 09:50:49 Alexandra

我的歉意，亞歷山德拉。我現在明白，問題在於計算Satterthwaite的近似值，而不是使用lmerTest。 –

我現在也找到了我上面的問題的答案。 lmerTest包中的anova函數提供了請求Kenwood-Roger近似ddf（ddf ='Ken'）的選項，並且您的anova調用給出了我正在尋找的內容。我還沒有找到一種方法來請求lmer函數中的K-R ddf近似值。在你的lmer中有沒有辦法使用K-R？如果不是，這可能是你的軟件包的一個很好的補充。感謝lmerTest。 –

要使用Kenward-Roger對anova的近似值，可以編寫anova（model，ddf =「kenw」），並在彙總函數中使用它，類似地可以這樣寫：summary（model，ddf =「kenw」）。 lmerTest包使用pbkrtest包實現了Kenward-Roger的近似值 – Alexandra

lmer從anova獲取p值

回答

相關問題