我正在尝试获取使用lmer构建的混合效果模型(m)的P值。
我一直在使用以下命令来提取P值,但不幸的是,对于那些< 10-16,我得到P = 0.0000。请参阅以下示例:
wicket-glue
如何计算那些的提取物P值?
答案 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)