data <- matrix(scan(),ncol=2,byrow=TRUE)
  1 6.08   1 22.29  1 7.51   1 34.36  1 23.68 
  2 30.45  2 22.71  2 44.52  2 31.47  2 36.81 
  3 32.04  3 28.03  3 32.74  3 23.84  3 29.64 
  
data
tab312 <- data.frame(data)
colnames(tab312) <- c("treatment","y")

tab312$treatment <- as.factor(tab312$treatment)    # to make sure that treatment is a factor, not numeric

rm(data)
tab312

boxplot(y ~ treatment, data=tab312, main="Table 3.1.2 data")

anova(lm(y ~ treatment, data=tab312))             # standard ANOVA

kruskal.test(y ~ treatment, data=tab312)          # Kruskal-Wallis test

rm(tab312)                                        # remove objects when finished