Dear all,
I have a mixed Poisson regression model . my function is:
## P(N=k)=[(a+2)/a
][(lambda^k)]*[((a+2/a+1)+lambda)^(-k-1)-(a+2+lambda)^(-k-1)];
lambda=a0+a1*x1+a2*x2
This is same as a Poisson-inverse
gaussian or Poisson-lognormal regression.I want to estimate the
parameters by fitting this model on my data and also I need loglike value. Since
some of the explanatory variables (x1,x2,...) may be nominal or ordinal; which
package in R can be used?
I am looking for any answer.
with the best wishes
Zamani.H
[[alternative HTML version deleted]]