[Rd] error in dbeta (PR#970)

Peter Dalgaard BSA p.dalgaard@biostat.ku.dk
06 Jun 2001 18:48:08 +0200


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?

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