I am using iteratively weighted least squares algorithm to obtain ML estimates in Poisson regression. Due to the factorial term in the log-likelihood equation, I can not compute the log-likelihood value which becomes -Inf. I think it is because of some big values of the dependent variable computed from the Poisson regression.
Do you have any solution?