datafilename="http://personality-project.org/r/datasets/R.appendix2.data" #read the data into a table data.ex2=read.table(datafilename,header=T) #show the data data.ex2 #do the analysis of variance aov.ex2 = aov(Alertness~1, data=data.ex2) summary(aov.ex2) #report the means and the number of subjects/cell print(model.tables(aov.ex2,"means"),digits=3) #graphical summary of means of the 4 cells boxplot(Alertness~Dosage*Gender,data=data.ex2)