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
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.