library(glmmTMB) library(dplyr) library(DHARMa) ml <- glmmTMB(menthlth ~ sex + medcost + ..., ziformula = ~ sex + medcost + X_educag + internet + ..., data = ...

I have relational data, i.e. observations for pairs of objects. More specifically these are migration rates between plant populations, which I would like to explain by a predictor. The migration rates ...

I fitted a complex betabinomial model with very large values for both variables comprised in the response variable (cbin(Methylated_C, unMethylated_C)). model <- glmmTMB( ...

Is there a way to obtain estimated expected values from a glmmTMB object? Specifically, if Y is the response vector and X,Z are design matrices for fixed and random effects respectively, I want an ...

I would like to obtain the estimated covariance matrix of the (random) linear predictors from a glmmTMB object in R, or at least the diagonal of that matrix, that is, the variances. In lme4 the full ...

I'm a newb using glmmTMB, especially with smooths. The first model (m1) gives me what I expect: parallel smooths for each site since there are random intercepts for site. My question is: is the second ...

Suppose I have the following model: library(glmmTMB) Owls <- transform( Owls, Nest = reorder(Nest, NegPerChick), NCalls = SiblingNegotiation, FT = FoodTreatment ) fit_zipoisson <- ...

I am using buildglmmTMB to build a model for several different taxa. I have 700 ish taxa and ~10 effects. I have a script that loops through all the taxa and applies buildglmmTMB to each taxa ...

I used glmmTMB to run GLMMs with a beta distribution to evaluate the proportion of time an animal spends doing a behavior. My observations were of different lengths, and I weighted the observations ...

Consider the following estimate: library(glmmTMB) library(mice) Owls <- Owls Owls[sample(1:nrow(Owls),10), "ArrivalTime"] <- NA predictorMatrix <- quickpred( Owls, include = ...

I have a dataset with 125 animals across 3 sites and 100500 observations. Both show up properly when looking at the structure of the data but when I run the model with an updated data frame (I added a ...

I have a dataset of bird counts, on a spatial grid, with environmental and prey data. I have applied a correction factor for distance from observers, so I have semi-continuous values, with lots of 0's ...

I am working on a grant proposal evaluating the effectiveness of an intervention on several different groups of hospitals and after some statistical consulting have landed on using a generalized ...

I'm treating a data with one continuous variable and one categorical variable. I built a model as following using glmmTMB: m=glmmTMB(y~x+(1+x|group), family=Gamma(link="log"), data=dat) '...

I'm using the emmeans package with a negative-binomial model implemented using the glmmTMB package. I'm trying to bias adjust my backtransformed emmeans per the workflow illustrated here: https://cran....