> pnorm(37:39,lower.tail=FALSE) > [1] 5.725571e-300 0.000000e+00 0.000000e+00 > > This is just a limitation of double precision floating-point arithmetic > ... > > curve(pnorm(x,lower.tail=FALSE),from=30,to=40,log="y") > .Machine$double.xmin But note curve(pnorm(x,lower.tail=FALSE, log=T),from=30,to=1000) Hadley -- http://had.co.nz/