[R] lmer and negative binomial family

Ben Bolker bolker at ufl.edu
Thu Oct 29 16:38:21 CET 2009




John Sorkin wrote:
> 
> Patrick,
> I am not sure about this but 
> (1) I think you need to load the lme4 package, not MASS
> (2) I don't think neg.bin is allowable. I suggest you load lme4 and then
> type ?family
> 
> John
> 
> 

This conversation would probably go better on
r-sig-mixed-models at r-project.org ...
In the meanwhile --
   for a long time (I thought always, but I guess that's not true) the
allowable
families have been hard-coded within lme4, for purposes of efficiency.
The bottom line is that you can only use the following families:

famNms <- c("binomial", "gaussian", "Gamma", "inverse.gaussian",
            "poisson", "quasibinomial", "quasipoisson", "quasi")

'quasipoisson' might be good, if you're willing to settle for a quasi-
answer,
although you should be careful with random effects variances -- in earlier
versions there are some problems with the scaling there

   If you really need NB I would suggest glmmADMB ...


-- 
View this message in context: http://www.nabble.com/lmer-and-negative-binomial-family-tp26112313p26115443.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list