-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
project.org] On Behalf Of Sahana Srinivasan
Sent: Tuesday, March 26, 2013 4:52 PM
To: Berend Hasselman
Cc: r-help at r-project.org
Subject: Re: [R] NaNS Error Message
Hi, sorry I can't provide the reproducible code - the formula is
currently being worked on for a research paper so I can't really give
out details in that regard.
That being said, k is initialized to 1 at the beginning of the loop. n
has various values (whole numbers, 0 or greater, no upper limit).
On Tue, Mar 26, 2013 at 3:41 PM, Berend Hasselman <bhh at xs4all.nl>
wrote:
On 26-03-2013, at 16:25, Sahana Srinivasan
<sahanasrinivasan.91 at gmail.com>
wrote:
Hi,
I'm using R to do a series of calculation and I have gotten several
warnings that say "NaNS produced". Whatever I could read on line
gives me an idea that this warning is produced when the number is
use is a
log or otherwise mathematically problematic.
I'm getting this error while using factorial() and gamma () on
strictly positive numbers (always greater than zero).
Here is a snippet of my code:
k goes from 1 to a positive limit.
S is always 100.
Only positive values of n are allowed into the loop.
if(n>0)
{
while(k<=lim)
{
sn<-(S-n);snfact<-gamma(sn);
sn2<-(2-n+S);sn2gam<-gamma(sn2);
num<-(ngam*sn2gam);
nk2<-(2+k-(2*n)+S);
nk2gam<-gamma(nk2);
den<-(k*nk2gam);
prob<-(num/den);
sum<-(as.numeric((k*prob))+sum);
k<-k+1;
}
}
The error message received for every instance of this loop is :
In gamma(sn) : NaNs produced
In gamma(sn2) : NaNs produced
In gamma(nk2) : NaNs produced
This not reproducible code.
Where is ngam?
Where is k set to 1?
What value does n have?
Berend
[[alternative HTML version deleted]]