[R-sig-ME] need help in glm

[hossain zamani]

Dear all,
I am a R using. I have a mixed Poisson regression. 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

 and I've attached in a pdf file also.This is a poisson mixed with a weighted exponential distribution in the case that E(theta)=1(theta=random.effect). in fact this model 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 then compare with Poisson regression(glm).Which package in R can be used.
 with the best wishes
 Zamani.H




      
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mixed poisson.pdf
Type: application/pdf
Size: 55072 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20090929/d6644c83/attachment.pdf>