[R-sig-ME] Question about zero-inflated Poisson glmer

David Duffy David.Duffy at qimrberghofer.edu.au
Tue Jun 28 10:51:26 CEST 2016


On Tue, 28 Jun 2016, Philipp Singer wrote:

> You can find a sample of the data here:
> https://www.dropbox.com/s/kqqxmc3wp225lug/r_sample.csv.gz?dl=1
>
> You can think of the setting as popularity of items "y" inside stores 
> "id" explained by two features "a" and "b" whereas "a" is more of a 
> control covariate and I am interested in whether "b" has a positive impact.
>

Probably not very helpful (what a horrible distribution y has), but
library(ordinal)
x$cats <- cut(x$y, c(-1,0,10,100,1000,20000))
c1 <- clmm2(cats ~ a + (1|id), data=x)
c2 <- clmm2(cats ~ a + b + (1|id), data=x)

does run...

Cheers, David Duffy.



More information about the R-sig-mixed-models mailing list