Re: [R] Negative Binomial Regression

From: Ben Bolker <bolker_at_ufl.edu>
Date: Tue, 29 Jul 2008 12:53:35 +0000 (UTC)

jcarmichael <jcarmichael314 <at> gmail.com> writes:

>
>
> Hello.
>
> I am attempting to duplicate a negative binomial regression in R. SAS uses
> generalized estimating equations for model fitting in the GENMOD procedure.
>
> proc genmod data=mydata (where=(gender='F'));
> by agegroup;
> class id gender type;
> model count = var1 var2 var3 /dist=NB link=log offset=lregtm;
> repeated subject=id /type=exch;
> run;
>
> Since my dataset has several observations for each subject, I need the
> REPEATED statement in order to indicate dependence among observations with
> the same subject ID and independence amongst those with distinct subject
> IDs. The TYPE statement goes on to specify the structure of the correlation
> matrix to be used (exchangeable in this case).

  I would try glmmPQL in the MASS package. I don't think you can *quite* get negative binomial regression this way, but you can definitely get a quasipoisson model. I think exchangeable correlation corresponds to correlation=corCompSymm() in your glmmPQL command.

  Ben Bolker



R-help_at_r-project.org mailing list
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. Received on Tue 29 Jul 2008 - 13:00:14 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 29 Jul 2008 - 15:32:46 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive