[R] ANOVA table look

[Keith Jones]

Y'all,

I need to get the look of a "standard" fixed effect ANOVA table:

anova(aov(meas~op*part,data=fs))
Analysis of Variance Table

Response: meas
           Df  Sum Sq Mean Sq F value Pr(>F)
op         2    2.62   1.308  1.3193 0.2750
part      19 1185.43  62.391 62.9151 <2e-16 ***
op:part   38   27.05   0.712  0.7178 0.8614
Residuals 60   59.50   0.992

when I perform a two-factor factorial with random factors analysis  
using:

fs.anova=lmer(meas~(1|op)+(1|part),data=fs)
 > print(fs.anova)
Linear mixed model fit by REML
Formula: meas ~ (1 | op) + (1 | part)
    Data: fs
    AIC   BIC logLik deviance REMLdev
  417.6 428.7 -204.8    410.7   409.6
Random effects:
  Groups   Name        Variance   Std.Dev.
  part     (Intercept) 1.0250e+01 3.2015e+00
  op       (Intercept) 1.8497e-16 1.3600e-08
  Residual             8.9167e-01 9.4428e-01
Number of obs: 120, groups: part, 20; op, 3

Fixed effects:
             Estimate Std. Error t value
(Intercept)  22.3917     0.7211   31.05

The data is twenty parts, three operators with two replications of  
measurements on the twenty parts.

Any help would be appreciated.

Thanks,

Keith Jones