anova_joint function

Joint analysis of variance