Yes, it is possible with the mmc function in the HH package.
install.packages(HH) ## if you don't have it yet.
library(HH)
?MMC
Look at the maiz example, the long last example in ?MMC.
When you are ready for followup questions, reply to this message to
R-help (not directly to me and
not to nabble). Please turn off the html mail option. The html
version of R code is garbled
to illegibility on this plain text email list.
Please read the posting guide linked at the bottom of all R-help mail messages.
Please read also read
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
Rich
On Tue, May 12, 2015 at 10:38 AM, Md Newaz msne...@lakeheadu.ca wrote:
Dear R-help,
can you please post the following message to the r-nabble forum.
Thanks!
-
Dear R users,
I have been attempting to carry out post-hoc tests on a split-plot design:
Model:
Yijkl = µ + Ci + ω(i)j + δ(ij) + Tk + CTik + ωT(i)jk + Pl + CPil + ωP(i)jl
+ TPkl + CTPikl + ωTP(i)jkl + Ɛ(ijkl)
I have successfully matched the appropriate degrees of freedom and mean
squares presented in the table below using aov().
*EMS Table:*
2
2
2
3
F
R
F
F
i
j
k
l
EMS
df
F(1,2)
Ci
0
2
2
3
δ2 + 6δ2δ + 6δ2ω + 12Φ(C)
1
ω(i)j
1
1
2
3
δ2 + 6δ2δ + 6δ2ω
2
δ(ij)
1
1
2
3
δ2 + 6δ2δ
0
F(1,2)
Tk
2
2
0
3
δ2 + 3δ2ωT + 12Φ(T)
1
F(1,2)
CTik
0
2
0
3
δ2 + 3δ2ωT + 6Φ(CT)
1
ωT(i)jk
1
1
0
3
δ2 + 3δ2ωT
2
F(2,4)
Pl
2
2
2
0
δ2 + 2δ2ωP + 8Φ(P)
2
F(2,4)
CPil
0
2
2
0
δ2 + 2δ2ωP + 4Φ(CP)
2
ωP(i)jl
1
1
2
0
δ2 + 2δ2ωP
4
F(2,4)
TPkl
2
2
0
0
δ2 + δ2ωTP + 4Φ(TP)
2
F(2,4)
CTPikl
0
2
0
0
δ2 + δ2ωTP + 2Φ(CTP)
2
ωTP(i)jkl
1
1
0
0
δ2 + δ2ωTP
4
Ɛ(ijkl)
1
1
1
1
δ2
0
Total
23
mod - aov(Budburst ~ CO2*SoilTemp*Photoperiod +
Error(Greenhouse/(SoilTemp*Photoperiod)), data = data)
summary(mod)
Error: Greenhouse
DfSum Sq Mean Sq F value Pr(F)
CO2 1 1465.21465.2 38.81 0.0248
Residuals2 75.5 37.8
Error: Greenhouse:SoilTemp
DfSum Sq Mean Sq F value Pr(F)
SoilTemp 1 238.00238.00 80.570.0122
CO2:SoilTemp 1 145.70145.70 49.320.0197
Residuals 2 5.912.95
Error: Greenhouse:Photoperiod
DfSum Sq Mean Sq F value Pr(F)
Photoperiod 2986.9 493.46.965 0.0498
CO2:Photoperiod20.2 0.10.001 0.9989
Residuals 4283.4 70.8
Error: Greenhouse:SoilTemp:Photoperiod
Df Sum Sq Mean Sq
F value Pr(F)
SoilTemp:Photoperiod 214.56 7.28
0.514 0.6330
CO2:SoilTemp:Photoperiod2186.31 93.15 6.576
0.0544
Residuals 456.67
14.17
Error: Within
Df Sum Sq Mean Sq F value Pr(F)
Residuals 216 288713.37
However, as neither TukeyHSD() nor glht() accept objects of class
“aovlist”, I cannot carry out the post-hoc tests. Is there any way to run a
post-hoc test on an object of class aovlist?
Alternatively, I tried modelling the data using lme() and lmer(), but the
problem is that I cannot match the appropriate degrees of freedom and mean
squares obtained from the above included expected mean squares table using
lme() or lmer().
Has anyone else encountered and overcome this issue?
Thanks in advance,
Md. Shah Newaz
Faculty of Natural Resources Management
Lakehead University
Thunder Bay, Ontario, Canada
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.