ML.k {asbio} | R Documentation |
The function uses the maximum likelihood method described by Bliss and R. A. Fisher (1953) to determine maximum likelihood estimates for the binomial parameters m (the mean) and k (a parameter describing aggregation/dispersion).
ML.k(f, x, res = 1e-06)
f |
A vector of frequencies for objects in |
x |
A vector of counts, must be sequential integers. |
res |
Resolution for the ML estimator. |
Returns a list with two items
k |
The negative binomial dispersion parameter, k |
m |
The negative binomial distribution mean, m |
The program is slow at the current resolution. Later iterations will use linear interpolation, or Fortran loops, or both.
Ken Aho
Bliss, C. I., and R. A. Fisher (1953) Fitting the negative binomial distribution to biological data. Biometrics 9: 176-200.
mites <- seq(0, 8)
freq <- c(70, 38, 17, 10, 9, 3, 2, 1, 0)
ML.k(freq, mites)