Данные, которые мне нужно проанализировать, следующие, и это код, который я использовал:
gw <-read.table ("C: / Пользователи / Dell / Documents / DRY .txt ", header = T, sep =", ") </p>
гВт
день производства
1 A 0.9916
2 A 0.9978
3 A 1.0123
4 B 2.1468
5 B 1.8116
6 B 1.9617
7 C 5.2277
8 C 5.3413
9 C 5.1211
10 D 7.3134
11 D 7.0067
12 D 7,6201
13 E 8,9215
14 E 8,1246
15 E 8,5231
anova1 <- aov (производство ~ день, данные = gw) </p>
сводка (anova1) Df Сум Sq Среднее Sq F значение Pr (> F)
день 4 128,37 32,09 547,2 1,16e-11 ***
Остатки 10 0,59 0,06
Signif. коды: 0 '' 0,001 '' 0,01 '' 0,05 '.' 0,1 '' 1
TukeyHSD (anova1) Множественное сравнение средних значений по Тьюки
95% семейного уровня достоверности
Fit: aov (формула = производство ~ день, данные = gw)
$ разность дней день вверх п ад
BA 0,972800 0,3220216 1,623578 0,0042593
C -A 4,229467 3,55786883 4,880245 0,0000000
DA 6,312833 5,6662050 6.963612 0.0000000
EA 7.522500 6.8717216 8.173278 0.0000000
C -B 3.256667 2.6058883 3.907445 0.0000001
DB 5.340033 4.6892550 5.990812 0.0000000
10.700 10.7000000000000000000000
D- C 2.083367 1.4325883 2.734145 0.0000076
E- C 3.293033 2.6422550 3.943812 0.0000001
ED 1.209667 0.5588883 1.860445 0.0008304
* 1083 Пока здесь все в порядке, проблема в том, что когда я хочу использовать HSD.tukey для того, чтобы иметь графику c (boxplot) с разными буквами, я не знаю, как это сделать ... Если кто-нибудь может помочь ??