Re: [R] Modelling poisson distribution with variance structure

[Ben Bolker]

Karen Moore  tcd.ie> writes:

> 
> I'm dealing with count data that's nested and has spatial dependence.
> I ran a glmm in lmer with a random factor for nestedness. Spatial dependence
> seems to have been accommodated by model. However I can't add a variance
> strcuture to this model (to accommodate heterogeneity).
> 
> Is there a model that can have a poisson distribution *AND*  a variance
> structure *AND* have AIC in output (for model comparison and selection)?
> Some we've looked at that can't:
> 
>- glmmPQL  - can add structures BUT can't have AIC (you can calculate it
>but it doesn't give correct AIC with this model)
>- glmm in lme4 (lmer)  - won't allow variance structure
>- gls -  can add variance but can't have Poisson


  [Any further discussion should probably go to 
r-sig-mixed-mod...@r-project.org ...]

  I'm not sure I know what you mean by Poisson + variance structure --
if the data are really Poisson (not overdispersed in some way), then
the variance structure is completely defined.  If you want to deal
with overdispersion, and have a well-defined AIC, you may be able
to add a per-observation random effect in lme4.  Alternatively,
you could just use a weights= argument in glmmPQL to set some sensible
mean-variance relationship, overlooking the fact that the data
are discrete and positive rather than being normally distributed
with an equivalent variance structure.  
may also be useful.

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

[R] Modelling poisson distribution with variance structure

[Karen Moore]

I'm dealing with count data that's nested and has spatial dependence.
I ran a glmm in lmer with a random factor for nestedness. Spatial dependence
seems to have been accommodated by model. However I can't add a variance
strcuture to this model (to accommodate heterogeneity).

Is there a model that can have a poisson distribution *AND*  a variance
structure *AND* have AIC in output (for model comparison and selection)?
Some we've looked at that can't:

   - glmmPQL  - can add structures BUT can't have AIC (you can calculate it
   but it doesn't give correct AIC with this model)
   - glmm in lme4 (lmer)  - won't allow variance structure
   - gls -  can add variance but can't have Poisson

Thanks so much,

Karen Moore
PhD Researcher,
FORESTBIO,
Department of Botany,
Trinity College Dublin
Ireland

[[alternative HTML version deleted]]

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