From maechler at stat.math.ethz.ch  Thu Mar  5 09:45:36 2009
From: maechler at stat.math.ethz.ch (maechler at stat.math.ethz.ch)
Date: Thu,  5 Mar 2009 09:45:36 +0100 (CET)
Subject: [Rd] besselI inaccurate for negative integer order (PR#13556)
Message-ID: <20090305084536.8D05C2832191@mail.pubhealth.ku.dk>

>>>>> "JL" == Jerry Lewis <Jerry.Lewis at biogenidec.com>
>>>>>     on Thu, 26 Feb 2009 16:15:11 +0100 (CET) writes:

    JL> Full_Name: Jerry W. Lewis
    JL> Version: 2.8.1
    JL> OS: Windows XP Professional
    JL> Submission from: (NULL) (198.180.131.16)


    JL> It should be the case that
    JL> besselI(x,-nu) == besselI(x,nu) == besselI(x,abs(nu))
    JL> for integer nu, yet R currently can return ridiculous values when nu is a
    JL> negative integer.

    JL> For instance, besselI(9.6,-44) returns -234626490 instead of the correct value
    JL> of 5.9041042646307223e-25, while besselI(9.6,44) gives essentially machine
    JL> accuracy.

Yes.  I have committed a bug-fix to both R-patched and R-devel.

Thank you for the report!
Martin Maechler, ETH Zurich


    JL> This is more than an idle mathematical curiosity, since
    JL> one consequence is that dskellam in the VGAM package can
    JL> return values <0 or >1.



[Rd] besselI inaccurate for negative integer order (PR#13556)
maechler at stat.math.ethz.ch