[R] neg-bin clustered analysis in R?

[Lutz Ph. Breitling]

Dear All,

I'm analysing a negative binomial dataset from a population-based
study. Many covariates were determined on household level, so all
members of a household have the same value for those covariates.

In STATA, there seems to be an option for 'clustered analysis' for
neg-bin regression. Does an equivalent exist for R(MASS)'s glm.nb or a
comparable function?

Many thanks for all help!
Lutz
-- 
Lutz Ph. Breitling, Marie Curie EST Research Fellow
Integrative & Comparative Biology (Irene Manton 8.93e)
University of Leeds, LS2 9JT Leeds/UK