Question

我在R：

中完成了这个公式

> 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

但是，我不确定stat，ndf，ddf，F.scaling，p.value是什么。我可以问一下如何报告此输出的统计结果？为什么在输出的右侧，列的头部末端有“.U”，数字与左边相似？

Answer 1

这有点灰色，但它可能介于统计问题和编程问题之间。

Effect：效果名称
stat：计算出F statistic
ndf：分子自由度（效果使用的参数数量）
ddf：分母自由度（测试效果的有效残余自由度），使用pbkrtest::KRmodcomp
F.scaling：从K-R
p.value：K-R估计p值

Kenward-Roger校正做了两件事：（1）它计算分母df的有效数; （2）它按统计量计算统计量。我认为（但不确定）这些列的.U变体（我认为）是在没有缩放的情况下计算的等效指标。

您会注意到.U版本几乎与第一组值相同;这在很大程度上是正确的，因为分母的自由度都很大（例如> 50） - 当ddf仅为27.4（＆＃34; vowel3＆＃34;效果）时，注意到最大的差异。有了这么大的ddf，你使用什么价值几乎不重要 - 你可以通过似然比检验获得相同的答案。

有关更多信息，我需要深入挖掘 - 这似乎没有记录。原始计算在pbkrtest::.KR_adjust;据我所知，ndf.U和ddf.U参数实际上与非.U对应参数相同，而F.scaling.U始终为NA;它只是统计数据和p值不同（它们与F.scaling接近1的程度非常相似......）

我可以通过电子邮件发送维护者建议对输出和文档进行一些小的调整（即ndf.U，ddf.U和F.scaling.U列应该被抑制;其他列可以记录下来更好）...

如何在afex函数中报告混合（MODEL）的结果？

1 个答案: