[Rd] Non central chi squared bug (PR#9406)

From: <alan.bain_at_gmail.com>
Date: Fri 01 Dec 2006 - 10:44:34 GMT

Full_Name: Alan Bain
Version: 2.4.0
Submission from: (NULL) (

Code for pnchisq contains following

if (tSml) {

if (x> f+ theta+ 5*sqrt(2*(f+2*theta) ))){

     return 1.; /* better than 0 --- but definately FIXME */


This needs to check which tail has been requested; it is only correct if the default lower_tail=1 has been requested; for upper tail should return 0 if more than 5 std devs away from mean under these circumstances.

R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Sat Dec 02 08:33:10 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 01 Dec 2006 - 23:30:47 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.