[R] fit.dist gnlm question, NaN and Inf results

From: Richard and Barbara Males <rbmales_at_gmail.com>
Date: Fri, 25 Jul 2008 18:43:09 -0400


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



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Fri 25 Jul 2008 - 22:46:51 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 25 Jul 2008 - 23:32:34 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive