如何在afex函數中報告混合（MODEL）的結果？

我在R來此公式：如何在afex函數中報告混合（MODEL）的結果？

> mixed3 <- mixed(peak_Mid ~ (1|item) + (1+vowel3|speaker) + sex*vowel3*Language, data=data1.frame, na.action=na.omit) 
Fitting 9 (g)lmer() models: 
[.........] 
Obtaining 8 p-values: 
[Note: method with signature ‘sparseMatrix#ANY’ chosen for function ‘kronecker’, 
target signature ‘dgCMatrix#ngCMatrix’. 
"ANY#sparseMatrix" would also be valid 
........] 
> summary(mixed3) 
       Effect  stat ndf  ddf F.scaling  p.value  stat.U ndf.U  ddf.U F.scaling.U p.value.U 
1   (Intercept) 9500.922104 1 70.40672 1.0000000 7.529698e-77 9500.922104  1 70.40672   NA 7.529698e-77 
2     sex 15.980281 1 71.52842 1.0000000 1.538529e-04 15.980281  1 71.52842   NA 1.538529e-04 
3    vowel3 8.596702 2 27.40531 0.9916348 1.264905e-03 8.669222  2 27.40531   NA 1.209863e-03 
4   Language 3.996819 2 70.74337 0.9909675 2.267036e-02 4.033250  2 70.74337   NA 2.194066e-02 
5   sex:vowel3 1.746398 2 75.92257 0.9870432 1.813334e-01 1.769323  2 75.92257   NA 1.774036e-01 
6  sex:Language 4.136050 2 170.78334 0.9964821 1.761500e-02 4.150652  2 170.78334   NA 1.737140e-02 
7  vowel3:Language 1.573332 4 66.15951 0.9799283 1.917146e-01 1.605559  4 66.15951   NA 1.832701e-01 
8 sex:vowel3:Language 1.239002 4 195.29430 0.9894859 2.956981e-01 1.252168  4 195.29430   NA 2.903144e-01

但是，我不知道什麼是統計，NDF，DDF，F.scaling，p.value。我可以問我如何報告這個輸出的統計結果嗎？爲什麼在輸出的右邊，在列頭的最後有'.U'，數字與左邊相似？

來源

2014-09-01 user3288202

這有點灰色地帶，但它可能介於統計問題和編程問題之間。

Effect：的效果
stat：計算F統計
ndf：分子自由度
ddf（的用於效果參數數）：分母自由度（有效剩餘度的測試效果的自由度），根據Kenward-Roger校正計算使用pbkrtest::KRmodcomp
F.scaling：F-統計計算從KR
p.value：K-R估計p值

的Kenward-羅傑校正做兩個獨立的事情：（1）它計算分母DF的有效號碼; （2）按照計算出的量對統計量進行縮放。我認爲（但不確定）這些列的.U變體是（我認爲是）在沒有縮放的情況下計算的等效指標。您會注意到.U版本與第一組值相差無幾;這在很大程度上是正確的，因爲分母的自由度都很大（例如> 50） - 當ddf僅爲27.4（「元音3」效應）時會注意到最大差異。有了這麼大的ddf，它幾乎不重要你使用什麼價值 - 你會從似然比測試中得到相同的答案。

欲瞭解更多信息，我需要深入挖掘 - 這似乎沒有記錄。原始計算在pbkrtest::.KR_adjust;據我所知，參數ndf.U和ddf.U實際上與它們的非.U相同，並且F.scaling.U總是NA;它只是統計值和p值不同（它們將非常類似於F.scaling接近1的程度）。

我可能會通過電子郵件向維護人員發送電子郵件，建議對輸出進行一些微小調整和文件（即ndf.U,ddf.U和F.scaling.U列應該可能被壓制;其他人可能被證明更好）...

來源

2014-09-01 21:11:33

非常感謝本。你的回答啓發了我。是的，我會通過電子郵件詢問維護人員。也許這是新的，因爲afex文件（cran.r-project.org/web/packages/afex/afex.pdf）在2014年8月發佈，這就是爲什麼很多人還沒有使用它。 – user3288202 2014-09-02 06:25:39

Ben Bolker，Henrik回覆我說'afex是基於pbkrtest包進行計算併產生結果'，他建議我諮詢'pbkrtest的維護者'。不過，他會根據你的建議做一些改變。 – user3288202 2014-09-02 13:29:48

這裏有兩點值得注意。（1）當只調用'mixed3'而不是'summary（mixed3）'時，會給出一個很好的格式化輸出，它已經省略了不必要的列（當調用'summary'或'anova'方法時，新版本的afex也會產生這個結果）。（2）在當前版本的afex中，「stat」已被稱爲「F」（僅適用於LMM）。 – Henrik 2014-09-02 13:39:08

如何在afex函數中報告混合（MODEL）的結果？

回答

相關問題