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
censored counts and glmer/glmmADMB
2 messages · Eric Elguero, Bert Gunter
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