[Rd] lambda==0 in dpois() (PR#459)

Ben Bolker bolker@zoo.ufl.edu
Fri, 25 Feb 2000 14:50:52 -0500 (EST)


On 25 Feb 2000, Peter Dalgaard BSA wrote:

> bolker@zoo.ufl.edu writes:

> >   Here's a patch to src/nmath/dpois.c that appears to work (although it could ? be more
> > elegant).
> 
> Thanks. We may just be able to squeeze that one in before 1.0.0. But
> shouldn't it return -Inf  (ML_NEGINF) for the log case rather than NaN?
> 
> 

  I think you're right.
  By the way, it looks like qpois() may be subject to a similar problem,
although I haven't dug in.  Don't know about any of the other
distributions (dbinom() looks OK).

  Ben

-- 
Ben Bolker                                  bolker@zoo.ufl.edu
Zoology Department, University of Florida   http://www.zoo.ufl.edu/bolker
318 Carr Hall/Box 118525                    tel: (352) 392-5697
Gainesville, FL 32611-8525                  fax: (352) 392-3704

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