[Rd] Mistake in qbeta.c ? (PR#941)

[Peter Dalgaard BSA]

tlm21@cam.ac.uk writes:

> 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.

Whoops! Fortunately it sits in code that calculates an initial
approximation, so one would hope it only affects speed of convergence. 
Thanks for pointing that out.

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._