error in dbeta (PR#970)
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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._