From maechler at stat.math.ethz.ch  Tue Feb  7 11:04:16 2006
From: maechler at stat.math.ethz.ch (maechler@stat.math.ethz.ch)
Date: Tue,  7 Feb 2006 11:04:16 +0100 (CET)
Subject: [Rd] Problems in d-p-q-r functions (was PR#8528, but unrelated)
Message-ID: <20060207100416.88F423F13C@slim.kubism.ku.dk>

>>>>> "BDR" == Prof Brian Ripley <ripley at stats.ox.ac.uk>
>>>>>     on Tue, 7 Feb 2006 08:52:43 +0000 (GMT) writes:

    BDR> For the record, some of these claims are untrue:
    >> df(0, 2, 2)
    BDR> [1] 1
    >> df(0, 1.3, 2)
    BDR> [1] Inf

Well, these first two I had fixed in the mean time.
So Ian was right about reporting them.

    >> x <- 1e-170
    >> pbeta(x, x, x)
    BDR> [1] 0.5

( but Ian only explicitly mentioned dbeta(x,x,x)  for x = 1e-162 
  and that *did* need to be fixed -- and I see you did fix it; 
  thanks a lot! )

Martin

 BDR> ............................
 BDR> ............................



[Rd] Problems in d-p-q-r functions (was PR#8528, but unrelated)
maechler@stat.math.ethz.ch