Я выполняю следующие ezANOVA:
RMANOVAGHB1 <- ezANOVA(GHB1, dv=DIF.SCORE.STARTLE, wid=RAT.ID, within=TRIAL.TYPE, between=GROUP, detailed = TRUE, return_aov = TRUE)
Мой набор данных выглядит так:
RAT.ID DIF.SCORE.STARTLE GROUP TRIAL.TYPE
1 1 170.73 SAL TONO
2 1 80.07 SAL NOAL
3 2 456.40 PROP TONO
4 2 290.40 PROP NOAL
5 3 507.20 SAL TONO
6 3 261.60 SAL NOAL
7 4 208.67 PROP TONO
8 4 137.60 PROP NOAL
9 5 500.50 SAL TONO
10 5 445.73 SAL NOAL
до rat.id 16.
Мои руководители не работают с R, поэтому они не могут мне помочь. Мне нужен код, который даст мне все постфактум контрасты, но поиск его только сбивает меня с толку все больше и больше. Я уже пытался выполнить TukeyHSD на выводе aov ezANOVA, а затем попробовал pairwise.t.test (поскольку я обнаружил, что bonferroni является более подходящим исправлением в этом случае), но, похоже, ни один из них не работает. Я также нашел информацию об использовании линейной модели, а затем multcomp, но я не знаю, будет ли это хорошим решением в этом случае. Я чувствую, что проблема со всем, что я пробовал, заключается либо в том, что у меня есть переменные между и внутри, либо в том, что мой набор данных настроен неправильно. Еще одним усложняющим фактором является то, что я только начинаю работать с R, и мои статистические знания все еще довольно базовые, поскольку это один из моих первых практических опытов проведения анализа.
Если это важно, вот результат анова:
$ANOVA
Effect DFn DFd SSn SSd F p p<.05 ges
1 (Intercept) 1 14 1233568.9 1076460.9 16.043280 0.001302172 * 0.508451750
2 GROUP 1 14 212967.9 1076460.9 2.769771 0.118272657 0.151521743
3 TRIAL.TYPE 1 14 137480.6 116097.9 16.578499 0.001143728 * 0.103365833
4 GROUP:TRIAL.TYPE 1 14 11007.2 116097.9 1.327335 0.268574391 0.009145489
$aov
Call:
aov(formula = formula(aov_formula), data = data)
Grand Mean: 196.3391
Stratum 1: RAT.ID
Terms:
GROUP Residuals
Sum of Squares 212967.9 1076460.9
Deg. of Freedom 1 14
Residual standard error: 277.2906
1 out of 2 effects not estimable
Estimated effects are balanced
Stratum 2: RAT.ID:TRIAL.TYPE
Terms:
TRIAL.TYPE GROUP:TRIAL.TYPE Residuals
Sum of Squares 137480.6 11007.2 116097.9
Deg. of Freedom 1 1 14
Residual standard error: 91.0643
Estimated effects may be unbalanced