[R] censored counts and glmer/glmmADMB

2013-12-12 Thread Eric Elguero

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@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.


Re: [R] censored counts and glmer/glmmADMB

2013-12-12 Thread 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.elgu...@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@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

__
R-help@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.