I’m working on a model that relies on the
normal_lcdf(x|mu, sigma) function, but I appear to be running into numerical precision problems with x >> mu. R’s
pnorm(x,mu,sigma,log.p=T) doesn’t exhibit the same problem (at least at the scale of x I’m looking at):
I found this thread where the similar case x << mu was discussed and Bob mentions concerns with the lcdf functions.
Has there been any movement with the implementation of *_lcdf functions? Or does anyone know of an approximation of normal_lcdf that works with large x?