anova_fn function

Compares differences between baseline means using Carlisle's montecarlo anova method