library(glmmTMB)
library(glmmADMB)
library(lme4)

# same problems on larger data
data = read.csv("/home/psinger/WikiLinks/data/r_sample_1000.csv", header=TRUE)
data$id = factor(data$id)

data$nz <- as.numeric(data$y>0)

head(data)

nrow(data)

m1 <- glmer(nz~1+(1|id), data=data, family="binomial")

summary(m1)

m2 <- glmer(nz~1+a+(1|id), data=s, family="binomial")

m3 <- glmmTMB(nz~1+a+(1|id), data=data, family="binomial")

# this works
m4 <- glmer(nz~1+a+(a|id), data=data, family="binomial")

summary(m4)

# this as well
m5 <- glmer(nz~1+b+(1|id), data=data, family="binomial")

summary(m5)

m6 <- glmmTMB(y~1+a+(1|id), data=subset(data,y>0), family="nbinom2")

summary(m6)

# when I use a larger dataset, I receive thousands of warnings like:
# 50: In f(par, order = order, ...) : value out of range in 'lgamma'