1

我在 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'，并且数字与左侧相似？

4

1 回答 1

4

这是一个灰色地带，但它可能介于统计问题和编程问题之间。

  • Effect: 效果名称
  • stat: 计​​算的 F 统计量
  • ndf：分子自由度（用于效果的参数数量）
  • ddf: 分母自由度（用于测试效果的有效剩余自由度），使用 Kenward-Roger 校正计算得出pbkrtest::KRmodcomp
  • F.scaling：从 KR 扩展 F 统计计算
  • p.value: KR 估计的 p 值

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

您会注意到.U版本几乎与第一组值相同。这在很大程度上是正确的，因为分母的自由度都很大（例如 >50）——当 ddf 仅为 27.4（“元音 3”效应）时，会发现最大的差异。使用如此大的 ddf，您使用什么值几乎无关紧要——您会从似然比测试中得到相同的答案。

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

我可能会通过电子邮件向维护人员发送电子邮件，建议对输出和文档进行一些小的调整（即，ndf.U可能应该取消列；其他列可以更好地记录）...ddf.UF.scaling.U

于 2014-09-01T21:11:33.397 回答