too large alpha or beta in dbeta ? (PR#643)

Dear friends.

Is this as expected ? Is alpha and beta too large simply ?

 > dbeta(.1,534,646)

[1] NaN
Warning message:
NaNs produced in: dbeta(x, shape1, shape2, log)

well, it should work, but the correct answer is effectively zero.
  pbeta(.1,534,646) gives 3.6e-213

Perhaps more worrying is

dbeta(.25,534,646)

[1] Inf


	-thomas

Thomas Lumley
Assistant Professor, Biostatistics
University of Washington, Seattle


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