[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