aov(), asking if case 
and education and their interaction have any association 
with age.
> anova(aov(age~case*education,infert))
Analysis of Variance Table
 
Response: age
                Df Sum Sq Mean Sq F value    Pr(>F)
case             1    0.1     0.1  0.0034    0.9536
education        2  753.3   376.7 15.0456 6.938e-07
case:education   2    0.2     0.1  0.0031    0.9969
Residuals      242 6058.4    25.0 
We're getting pretty much the same answers, except that the variance has been
partitioned appropriately for looking at the more complex question. 
anova() produces the ANOVA table from the output of
aov().