From thomas@biostat.washington.edu  Thu Aug 24 21:45:54 2000
From: thomas@biostat.washington.edu (thomas@biostat.washington.edu)
Date: Thu, 24 Aug 2000 22:45:54 +0200 (MET DST)
Subject: [Rd] Re: [R] too large alpha or beta in dbeta ? (PR#643)
Message-ID: <200008242045.WAA25703@pubhealth.ku.dk>

On Thu, 24 Aug 2000, Troels Ring wrote:

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


[Rd] Re: [R] too large alpha or beta in dbeta ? (PR#643)
thomas@biostat.washington.edu