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