Because of overdispersion, I try to use beta-binomial method fit the data. However, it seems emmeans and glht did not support the model by betabin. I test by the example data or