From tlm21@cam.ac.uk  Wed May 16 12:49:35 2001
From: tlm21@cam.ac.uk (tlm21@cam.ac.uk)
Date: Wed, 16 May 2001 13:49:35 +0200 (MET DST)
Subject: [Rd] Mistake in qbeta.c ? (PR#941)
Message-ID: <200105161149.NAA27502@pubhealth.ku.dk>

Full_Name: Tim Massingham
Version: 1.2.2
OS: Debian/Linux
Submission from: (NULL) (131.111.8.68)


In 1.2.2 sources (also in 0.90.1. I haven't been able to check other versions)
Line 103 in qbeta.c should read:
w = y * sqrt(h + r) / h - (t - s) * (r + 5. / 6. - 2 / (3 * h));

since otherwise the 5 / 6 will evaluate to zero (I think).
The Statlib fortran code uses five / six instead.

Cheers, thanks for great software.

	  TimM.

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


[Rd] Mistake in qbeta.c ? (PR#941)
tlm21@cam.ac.uk