title: Is there an R function to fit a GLMM for count data with range 0-15? (Possible right censoring needed?)
body: <p>We are due to collect some survey data that has a range 0-15 score, potentially skewed, and has a multilevel structure (repeated measures and clustering). I'm anticipating that fitting a linear mixed model in R lmer will be problematic given the outcome distribution.</p>
  <p>I am considering whether some sort of right-censored generalized linear mixed model (Poisson) may be a solution but I'm struggling to find something to fit this model.</p>
  <p>I think the closest that I can find is the <code>VGAM::vglm</code> with <code>family = cens.poisson</code> but, as far as I can tell, it cannot include multilevel structure?</p>
  <p>Does anyone know any R functions that would permit this model? If so, is there an equivalent power calc function or would this be written as a simulation?</p>