Is there any package in R for fitting random effects model arising from mixture of two Poisson models?

Question

That is, the model contains two random effects (independent/correlated) at each level of the hierarchy.