[Rd] error in dbeta (PR#970)

[ripley@stats.ox.ac.uk]

On 6 Jun 2001, Peter Dalgaard BSA wrote:

> pwolf@wiwi.uni-bielefeld.de writes:
>
> > Full_Name: Hans Peter Wolf
> > Version: 1.2.1
> > OS: hpux10.20
> > Submission from: (NULL) (129.70.84.25)
> >
> >
> > dbeta computes a wrong result with parameters (1.3,1)
> >
> > > version
> > platform hppa2.0-hp-hpux10.20
> > arch     hppa2.0
> > os       hpux10.20
> > system   hppa2.0, hpux10.20
> > status
> > major    1
> > minor    2.1
> > year     2001
> > month    01
> > day      15
> > language R
> > > dbeta((0:5)/5,1.3,1.0)
> > [1] 0.0 0.0 0.0 0.0 0.0 1.3
> > > dbeta((0:5)/5,1.30001,1.0)
> > [1] 0.0000000 0.8021373 0.9875537 1.1152952 1.2158296 1.3000100
> > > dbeta((0:5)/5,1.29999,1.0)
> > [1] 0.0000000 0.8021508 0.9875566 1.1152895 1.2158163 1.2999900
>
> Hm. Doesn't happen with 1.2.3 on Linux and dbeta.c hasn't changed
> inbetween. Any chance of getting you to step through the computations
> with a debugger on you machine?

It does happen on Solaris 2.7, though, which might give Peter a means
to look at it on an earlier Solaris.

-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272860 (secr)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595


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