Prior to the steps below, enter the data as appropriate for the analyses (described elsewhere). As always, the following commands should be typed directly in the R console window.
Get the mean and standard deviation of the dependent variable for each of the levels.
by(Outcome,Factor,mean)
by(Outcome,Factor,sd)
Results <- aov(Outcome~Factor)
model.tables(Results,"means")
Test each pairwise comparison for statistical significance.
TukeyHSD(Results)