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.
Results <- aov(Outcome~FactorA*FactorB)
model.tables(Results,"means")
Get the analysis of variance source table and a test of statistical significance.
summary(Results)