我正在尝试检查 p 上的声明。23使用回归和多级/分层模型进行数据分析
例如,考虑两个独立研究,其效应估计和标准误差分别为 25 ± 10 和 10 ± 10。第一项研究在 1% 的水平上具有统计显着性
在 R 中
> estimate <- 25
> se <- 10
> z <- estimate / se
> z
[1] 2.5
随着样本量的增加,t分布收敛于正态(Cassella Berger 2ed 练习 5.18)
然后我插入z分布函数以获得观察任何等于|z|或大于(p值)ISL p 的值的概率。67
> pnorm(2.5, lower.tail=FALSE)
[1] 0.006209665
而如果我使用方法 [3] 中的公式
> exp(-0.717 * z - 0.416 * z * z)
[1] 0.01236977
0.0062 是因为近似差吗?或者它不应该那样计算?
[3]:林俊泰。在袖珍计算器上使用正态尾概率及其倒数的近似值。应用统计 1989;38:69-70。