data2.aov<-aov(data2$y~factor(data2$x))
print(summary(data2.aov))
boxplot(data2$y~data2$x)