Hello, At present, can generalized linear mixed models with negative binomial distribution and estimating the shape parameter be fit using R? I am aware of glm.nb but am wondering about incorporation of mixed effects. Thanks in advance, Brian Aukema