In R, I have used the glmmTMB package to run a zero-inflated negative binomial model. This model had random and nested factors in the model. I originally did a quick and dirty posteriori power analysis using pwr.f2.test from the pwr package. Obviously that analysis was not appropriate because this function is for general linear models. I was planning on using the SIMR package to calculate power for my model since SIMR can be used for generalized linear mixed models. It was built for LME4 models, but I think it can be used for models from other packages. Can I use a model from glmmTMB with the powerSim function from the SIMR package? If not, what other statistical power testing functions can I use that will complement my glmmTMB model?

Documentation about the SIMR package: https://besjournals.onlinelibrary.wiley.com/doi/10.1111/2041-210X.12504

Thank you!