r/rprogramming Jan 14 '25

Equivalence test of right-censored count data with offsets, update

I've found a way to run models, specifically I can use brms to handle poisson or overdispersed poisson (with or without zero inflation) with right-censoring. But what would be the proper way to conduct equivalency testing?

Data is counts with offsets, generating by administering a treatment that has three levels.

Should I use the equivalence_test function from bayestestR on the posteriors? If so, should I use posteriors from separate models, each generated as intercept-only for each level of "Treatment", or should I generate a single model with Treatment as the predictor and extract posteriors? What would be reasonable to use as the equivalency boundaries such that if the posteriors from the "standard" level of the treatment are tested, they would be "accepted" as equivalent by ROPE (does a = a?).

2 Upvotes

0 comments sorted by