search for: trunpoismle

Displaying 1 result from an estimated 1 matches for "trunpoismle".

2008 Jun 25
0
Goodness-of-fit for zero-truncated poisson distribution
...nce, 3 counted twice and 2 counted 4 times. Applying table() here does not return the frequency of individuals counted 3 times (which is 0 in this case). How can I achieve such a contingency table where I could end up with something like Fi 1 2 3 4 3 3 0 2 2) function to estimate lambda with MLE trunpoismle <- function(xbar, tol = .Machine$double.eps^0.25, nmax = 10000) { lambda <- xbar i <- 1 diff <- 1 while (abs(diff) > tol && i < nmax) { diff <- (lambda - xbar * (1 - exp(-lambda))) / (1 - xbar * exp(-lambda)) lambda &l...