> 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/ ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
