[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-devel
Subject:    [Rd] besselI inaccurate for negative integer order (PR#13556)
From:       Jerry.Lewis () biogenidec ! com
Date:       2009-02-26 15:15:11
Message-ID: 20090226151511.C1D272834311 () mail ! pubhealth ! ku ! dk
[Download RAW message or body]

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


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

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

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

______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
[prev in list] [next in list] [prev in thread] [next in thread]