lmer输出非常低的P值

时间:2016-12-21 10:46:23

标签: lme4

我正在尝试获取使用lmer构建的混合效果模型(m)的P值。

我一直在使用以下命令来提取P值,但不幸的是,对于那些< 10-16,我得到P = 0.0000。请参阅以下示例:

wicket-glue

如何计算那些的提取物P值?

1 个答案:

答案 0 :(得分:1)

您应该能够使用move *OFF *IN30 *IN30 DOWEQ *OFF READ PFILE 30 *IN30 IFEQ *ON LEAVE ENDIF fILEid IFNE COMPAREid ITER ENDIF DELET PFILEREC LEAVE ENDDO 参数更准确地计算lower.tail=FALSE

1-p

如果你有真正的小p值(即< 1e-300而不是< 1e-16),你仍然可以在log-p范围内计算它们:

coefs$p.z <- 2 * pnorm(abs(coefs$t.value),lower.tail=FALSE)