Displaying 1 result from an estimated 1 matches for "trunpoisml".
Did you mean:
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 &...