Assume we parametrize the variance in a negative binomial regression model as:
variance = mu(1+mu/theta).
Fitting the model with glmmTMB, I check theta from:
fixef(fit)$disp
In mgcv, I look at:
fit$family$getTheta()
Are those equivalent in meaning theta as in the variance formula?