In this post I see how they calculate the confidence interval on the theta parameter of a negative binomial GLM: Confidence Interval for the Dispersion parameter of negbin distribution. I typically take into account overdispersion either using a quasibinomial GLM or using a binomial GLMM by incorporating an observation-level random effect (from which one could calculate the intraclass correlation as a measure of overdispersion, cf. https://royalsocietypublishing.org/doi/10.1098/rsif.2017.0213) though.
Would anyone happen to know how to calculate the confidence intervals on the dispersion coefficients if overdispersion is taken into account in those ways? Via nonparametric bootstrapping, or is there other faster approaches?