Почему я получаю ошибку с гамма-дистрибутивом с GLMM?

Question

Я задал этот вопрос более подробно на Stack Exchange. Здесь я приведу более короткую версию. Приведенный ниже код не работает для гамма-дистрибутива, почему?

df <- data.frame(month=sample(1:12, 3000, replace = T),
                 location=sample(1:93, 3000, replace = T),
                 river=sample(1:60, 3000, replace = T),
                 pco2=rgamma(3000, shape = 2, rate = 1),
                 chla=rlnorm(3000, 15, 1))

mod1 <- glmmTMB::glmmTMB(chla ~ pco2 + (1|river/location/month), family = "Gamma", data = df)
mod2 <- lme4::glmer(chla ~ pco2 + (1|river/location/month), family = "Gamma", data = df)