library(reshape2) CW = read.table('../../Data/cloth_waste.txt', header=FALSE, skip=20, na.strings='*') colnames(CW) = c('A', 'B', 'C', 'D', 'E') ##print(CW) ## Convert from wide to long: ## CW_long = melt(CW, id.vars=c(), variable.name='supplier', value.name='run_up') boxplot(run_up ~ supplier, data=CW_long, xlab='supplier', ylab='run_up') print(anova(lm(run_up ~ supplier, data=CW_long)))