Prior to the steps below, enter the data as appropriate for the analyses (described elsewhere).
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)