For more accuracy you can check the stats/distribs addon. http://www.jsoftware.com/jwiki/Addons/stats/distribs
load 'stats/distribs' (,: pnorm) - i. _10 _9 _8 _7 _6 _5 _4 _3 _2 _1 0 1.14519e_19 6.28471e_16 1.27998e_12 9.86588e_10 2.86652e_7 3.16712e_5 0.0013499 0.0227501 0.158655 0.5 On Thu, Mar 12, 2015 at 3:42 PM, Ryan <[email protected]> wrote: > 'normalcdf' gives different a result than R and Python for y < -7 > > In J: > > load'stats' > ([,:normalcdf) - i._10 > _9 _8 _7 _6 _5 _4 > _3 _2 _1 0 > _2.22045e_16 1.38778e_15 1.27998e_12 9.86588e_10 2.86652e_7 3.16712e_5 > 0.0013499 0.0227501 0.158655 0.5 > > In R: > > cat(pnorm(-9:0)) >> > 1.128588e-19 6.220961e-16 1.279813e-12 9.865876e-10 2.866516e-07 > 3.167124e-05 0.001349898 0.02275013 0.1586553 0.5 > > Does anyone know why? > > thanks, > Ryan > > > The information in this e-mail is intended only for the person to whom it > is > addressed. If you believe this e-mail was sent to you in error and the > e-mail > contains patient information, please contact the Partners Compliance > HelpLine at > http://www.partners.org/complianceline . If the e-mail was sent to you in > error > but does not contain patient information, please contact the sender and > properly > dispose of the e-mail. > > ---------------------------------------------------------------------- > For information about J forums see http://www.jsoftware.com/forums.htm > ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
