> tomato.fit <- aov(mfwlog~tg430+yr+yr:tg430+Error(entry/tg430),tomatos)
> tomato.fit

Call:
aov(formula = mfwlog ~ tg430 + yr + yr:tg430 + Error(entry/tg430), data = 
        tomatos)

Grand Mean: 0.7813939 

Stratum 1: entry

Terms:
                   tg430       yr yr:tg430 Residuals 
 Sum of Squares  2.40849  0.00001  0.00029  12.45864
Deg. of Freedom        2        1        1        91

Residual standard error: 0.3700109 
1 out of 5 effects not estimable
Estimated effects may be unbalanced

Stratum 2: Within

Terms:
                      yr yr:tg430 Residuals 
 Sum of Squares 0.000383 0.075300  1.647066
Deg. of Freedom        1        2        66

Residual standard error: 0.1579733 
Estimated effects may be unbalanced
>
> summary(tomato.fit)
Error: entry 
          Df Sum of Sq  Mean Sq  F Value     Pr(F) 
    tg430  2   2.40849 1.204245 8.796009 0.0003218
       yr  1   0.00001 0.000013 0.000095 0.9922260
 yr:tg430  1   0.00029 0.000290 0.002121 0.9633697
Residuals 91  12.45864 0.136908                   

Error: Within 
          Df Sum of Sq    Mean Sq  F Value     Pr(F) 
       yr  1  0.000383 0.00038333 0.015361 0.9017410
 yr:tg430  2  0.075300 0.03765023 1.508692 0.2287282
Residuals 66  1.647066 0.02495555