[R-sig-ME] glmmTMB negbinom not working with spatial autocorrelation
ud|t@@b@n@@|17 @end|ng |rom |mper|@|@@c@uk
Wed Oct 16 19:33:23 CEST 2019
I had been running a mixed model with poisson distribution of the following type, with a spatial autocorrelation term, which works fine:
Y(count data) ~ x1 + square(x1) + x2 + square(x2) + exp( ) + (1|population/species)
I realized that my dataset has a lot of small values (mostly 1 and 2) and some large values, so that the data is highly skewed and over dispersed. So I tried to run the following negbinom1 model:
Y(count data) ~ x1 + square(x1) + x2 + square(x2) + exp( ) + (1|population/species) + ziformula = ~.
This time the model doesn’t run and says it cannot find one of the independent variables in the dataset. If I remove that variable from the model then it says so for another variable and so on. If I remove the factor for spatial autocorrelation, the model seems to work fine. Can anyone tell me what’s happening and if what I am doing is appropriate for a highly skewed and over dispersed dataset?
Centre for Ecological Sciences
Indian Institute of Science
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models