anova <- aov(values ~ ind, data = df_stacked) summary(anova)