shunyaowu 为什么我计算的pvalue和论文计算的不一样？ 论文pvalue的计算值网址：http://www.pnas.org/content/suppl/2007/05/03/0701361104.DC1/01361Table5.pdf 论文计算z=8.4的pvalue是1.6x10-17 我使用 2*pnorm(-abs(8.4)) 得到的是 4.464786e-17。 其他的值也都不同
havocor 回复 第1楼 的 shunyaowu：求 Φ(z) 一直都需要用數值方法. 有多種數值方法可以用, 如最簡單的 梯形法, Simpson 法, 或 Gaussian quadrature. 當 z 值小時, 也可以 用 e^x 的 Taylor's expansion 做近似. 因此，可能用的是不同的計算方法，以致有不同的計算結果。
havocor 回复 第1楼 的 shunyaowu：應該用pnorm(-abs(8.4)) 可以了吧， > pnorm(-abs(8.4)) [1] 2.232393e-17 數值還是不一樣，應該是求 Φ(z) 一直都需要用數值方法. 有多種數值方法可以用, 如最簡單的 梯形法, Simpson 法, 或 Gaussian quadrature.
