> suntan.aov <- aov(clrchng ~ brand + form + brand:form)
> summary.lm(suntan.aov)

Call:
aov(formula = clrchng ~ brand + form + brand:form)

Residuals:
    Min      1Q  Median      3Q     Max 
-3.9875 -1.1025 -0.0938  0.5219  5.6825 

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)    13.402      1.176  11.394 1.16e-09 ***
brand2         18.995      1.663  11.419 1.12e-09 ***
form2          -3.903      1.663  -2.346  0.03063 *  
form3          -4.303      1.663  -2.586  0.01862 *  
brand2:form2   -7.195      2.352  -3.058  0.00677 ** 
brand2:form3   -3.597      2.352  -1.529  0.14359    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 2.352 on 18 degrees of freedom
Multiple R-squared: 0.9455,     Adjusted R-squared: 0.9304 
F-statistic: 62.48 on 5 and 18 DF,  p-value: 9.676e-11 

> anova(suntan.aov)
Analysis of Variance Table

Response: clrchng
           Df  Sum Sq Mean Sq  F value    Pr(>F)    
brand       1 1422.50 1422.50 257.0423 4.210e-12 ***
form        2  254.48  127.24  22.9924 1.103e-05 ***
brand:form  2   51.77   25.88   4.6772   0.02313 *  
Residuals  18   99.61    5.53                       
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 
> 
>