[Rd] qbeta function (FYI, compiler bug)

[apjaworski@mmm.com]

I am not sure if this is still true.

Hre is a quick check:

(1) R-1.4.0 (Prof. Ripley build), Win2000

     > qbeta(0.025,4,0.5)
     [1] 0.5137795

(2) R-1.4.0 (compiled from source), RH-7.2 (gcc 2.96-98)

     > qbeta(0.025,4,0.5)
     [1] 0.5137795

Andy

__________________________________
Andy Jaworski
Engineering Systems Technology Center
3M Center, 518-1-01
St. Paul, MN 55144-1000
-----
E-mail: apjaworski@mmm.com
Tel:  (651) 733-6092
Fax:  (651) 736-3122


                                                                                                                                               
                    Morten Welinder                                                                                                            
                    <terra@diku.dk>      To:     r-devel@hypatia.math.ethz.ch                                                                  
                                         cc:     (bcc: Andrzej P. Jaworski/US-Corporate/3M/US)                                                 
                    01/07/2002           Subject:     [Rd] qbeta function (FYI, compiler bug)                                                  
                    09:19 AM                                                                                                                   
                                                                                                                                               
                                                                                                                                               






Hi there,

this is just to let you know that the qbeta function, which was
copied from R into Gnumeric, has been confirmed to be miscompiled
by gcc 2.96 on Linux.  (That's Red Hat's compiler.)

This shows by qbeta(0.025,4,0.5) ending up taking the wrong
branch of "if (alpha <= 0.5)".

We compile things in a different context, so this may or may not
affect you.  The qbeta function was taken from R 0.64.

Morten

http://www.gnome.org/gnumeric/
http://bugzilla.gnome.org/show_bug.cgi?id=68057
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
_._._




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