Когда я смотрю на ANOVA для своего набора данных, я получаю значение P 0,018, которое является значимым при альфа = 0,05. Однако, когда я провожу тест HSD Тьюки на этом ANOVA, все попадает в группу A. Может ли кто-нибудь объяснить, что здесь происходит?
Вот мои данные:
dataHodag <- read.table(text="
Plot Rate Rep HollowHeartPercent
1 2155 1 1 0.0
2 2252 1 2 0.2
3 2343 1 3 0.2
4 2433 1 4 0.2
5 2211 2 2 0.0
6 2326 2 3 0.0
7 2412 2 4 0.0
8 2136 2 1 0.2
9 2124 3 1 0.0
10 2334 3 3 0.0
11 2421 3 4 0.0
12 2236 3 2 0.2
13 2246 4 2 0.0
14 2354 4 3 0.0
15 2446 4 4 0.0
16 2116 4 1 0.1
17 2145 5 1 0.0
18 2225 5 2 0.0
19 2314 5 3 0.0
20 2451 5 4 0.0", header=TRUE)
Я использую пакет Agricolae для своего анализа. Это код, который дает мне ANOVA:
TukeyAOV <- aov(HollowHeartPercent ~ Rep + Rate, data=dataHodag)
summary(TukeyAOV)
# Df Sum Sq Mean Sq F value Pr(>F)
# Rep 1 0.00250 0.00250 0.406 0.533
# Rate 1 0.04225 0.04225 6.857 0.018 *
# Residuals 17 0.10475 0.00616
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
А вот код моего теста HSD Тьюки:
HSD <- agricolae::HSD.test(TukeyAOV, "Rate", alpha = 0.05, group=TRUE)
print(HSD$groups)
# HollowHeartPercent groups
# 1 0.150 a
# 2 0.050 a
# 3 0.050 a
# 4 0.025 a
# 5 0.000 a
Любая помощь приветствуется. Благодаря!!
1 ответ
Вы заметили несоответствие степеней свободы? Вы читаете Rate и Rep как числовые переменные, а не как факторы. Сравнить с:
> dataHodag$Plot <- as.factor(dataHodag$Plot)
> dataHodag$Rate <- as.factor(dataHodag$Rate)
> TukeyAOV <- aov(HollowHeartPercent ~ Rep + Rate, data=dataHodag)
> summary(TukeyAOV)
Df Sum Sq Mean Sq F value Pr(>F)
Rep 1 0.0025 0.002500 0.368 0.554
Rate 4 0.0520 0.013000 1.916 0.164
Residuals 14 0.0950 0.006786
Надеюсь, это поможет.
Похожие вопросы
Новые вопросы
r
R - это бесплатный язык программирования с открытым исходным кодом и программная среда для статистических вычислений, биоинформатики, визуализации и общих вычислений. Пожалуйста, предоставьте минимальные и воспроизводимые примеры вместе с желаемым результатом. Используйте dput () для данных и укажите все небазовые пакеты с вызовами library (). Не вставляйте изображения для данных или кода, вместо этого используйте блоки кода с отступом. По вопросам, связанным со статистикой, используйте https://stats.stackexchange.com.