简化正态累积函数的除法

Question

我正在努力研究如何在R中简化两个正态概率函数的商。实际上，我在计算一个条件斜率正态密度，这是两个函数的除法：

pnorm(alpha0+t(alpha2)%*%chol2inv(chol(omega2))%*%t(y2-xi2.1))/pnorm(tau2.1)

其中alpha0+t(alpha2)%*%chol2inv(chol(omega2))%*%t(y2-xi2.1)和tau2.1产生实数。例如，有时我会遇到pnorm(-50)/pnorm(-40)，例如不一致的0/0。但这些值不是零，R只是近似值。我尝试使用erf函数，但我得到了相同的问题(0/0)。

有什么关于如何克服这个问题的提示吗？

Answer 1

pnorm有一个log参数，这使得它返回log(p)。将公式更改为exp(log(p1) - log(p2))：

exp(pnorm(-50, log = TRUE) - pnorm(-40, log = TRUE))
#[1] 2.95577e-196