[R] censored counts and glmer/glmmADMB

Bert Gunter gunter.berton at gene.com
Thu Dec 12 16:58:05 CET 2013


This post should go the r-sig-mixed-models list, where you are much
more likely to get useful help than here, which is a general R
programming help list.

Cheers,
Bert

On Thu, Dec 12, 2013 at 7:41 AM, Eric Elguero <Eric.Elguero at ird.fr> wrote:
> dear R-users,
>
> I have to model counts where all counts above some threshold
> have been censored. In the same dataset I have too many zeroes for
> a Poisson or even a negative binomial distribution to make
> sense, so I would need a zero-inflated-censored negative binomial
> family for use in glmer (or glmmADMB?). That seems not to exist.
>
> my question is :
> how could I add a custom-built family of distributions that
> I could call in glmer/glmmADMM ?
>
> if it's not possible, I am considering imputing fake values
> to replace the censored ones, but I am unsure whether this
> is bad or very bad...
>
> Eric Elguero
> MIVEGEC (UM1- UM2 -CNRS 5290-IRD 224)
> Maladies infectieuses et vecteurs :
> écologie, génétique, évolution et contrôle
> Centre IRD de Montpellier
> 911 Av Agropolis - BP 64501
> 34394 Montpellier Cedex
>
> ______________________________________________
> 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.



-- 

Bert Gunter
Genentech Nonclinical Biostatistics

(650) 467-7374



More information about the R-help mailing list