Hi, I need to fit a general generalized linear model, for observations on a response Y which is Gamma-distributed and observed in clusters. So, if E(Y) = mu, then for a suitable link function f, f(mu) is a (linear) function of both fixed and random effects. Are there R packages or functions which fit these models, using either MLE or estimating-equation approaches? Many thanks, Sam -- Professor Samuel D. Oman Department of Statistics Hebrew University Mount Scopus Jerusalem 91905 Israel tel: + 972 2 5883 442 fax: + 972 2 5883 549