Displaying 1 result from an estimated 1 matches for "snfact".
Did you mean:
infact
2013 Mar 26
2
NaNS Error Message
...s 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;
}
}...