Richard and Barbara Males
2008-Jul-25 22:43 UTC
[R] fit.dist gnlm question, NaN and Inf results
I am attempting to fit discrete data (daily counts of arrivals of recreational vessels at locks on a river) using the fit.dist package. Some distributions return values of NaN and Inf for certain situations, an example with Inf values is shown below. # of vessels: 1 2 3 4 5 6 7 8 9 10 11 # of days with # of vessels: 35 20 10 5 6 3 1 3 1 0 1 (stored in rTemp$counts) Can anyone tell me under what conditions I will get these Inf/NaN? Thanks in advance Richard Males Cincinnati, Ohio, USA my function call is as follows, for each of the distributions. fit.dist(c(1:length(rTemp$counts)),rTemp$counts,"binomial") output from fit.dist for above input situation binomial distribution, n = 85 mean variance nu.hat 2.6352941 4.6552249 0.2635294 -log likelihood AIC Inf Inf beta binomial distribution, n = 85 mean variance nu.hat rho.hat 2.6352941 4.6552249 0.2947923 0.1403921 -log likelihood AIC Inf Inf Poisson distribution, n = 85 mean variance mu.hat 2.635294 4.655225 2.635294 -log likelihood AIC 28.86352 29.86352