[R] ANOVA table look

Keith Jones keithlj at keithljelp.com
Wed Jan 26 00:14:42 CET 2011


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



More information about the R-help mailing list