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