Hi John,
The Gamlss.mx <http://gamlss.mx/> package can accommodate variables that follow
negative binomial (and other) distributions in multilevel models.
Mark
> On Jun 25, 2016, at 11:00 PM, John Sorkin <jsor...@grecc.umaryland.edu> wrote:
>
> Is there a function that will run a model appropriate for over dispersed data
> (such as a negative binomial or quasipoisson)
> with a random effects (or mixed effects) model in R? GLMER will not accept:
> family=quasipoisson(link="log") or
> family=negbinomial(link="log")
>
> I want to run something like the following:
> fit0 <- glmer(Fall ~
> Group+(1|PID)+offset(log(TimeYrs)),family=quasipoisson(link="log"),data=data)
> Thank you
> John
> John David Sorkin M.D., Ph.D.
> Professor of Medicine
> Chief, Biostatistics and Informatics
> University of Maryland School of Medicine Division of Gerontology and
> Geriatric Medicine
> Baltimore VA Medical Center
> 10 North Greene Street
> GRECC (BT/18/GR)
> Baltimore, MD 21201-1524
> (Phone) 410-605-7119
> (Fax) 410-605-7913 (Please call phone number above prior to faxing)
>
> Confidentiality Statement:
> This email message, including any attachments, is for ...{{dropped:16}}
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
