回答編集履歴
1
記述例の改善
test
CHANGED
|
@@ -4,8 +4,7 @@
|
|
|
4
4
|
|
|
5
5
|
[multcomp のマニュアル](https://cran.r-project.org/web/packages/multcomp/multcomp.pdf) の `cld` の部分に目を通しましたが値の大小関係と連動してアルファベットを設定する機能は読み取れませんでした。現在その機能には対応していないのかもしれません。
|
|
6
6
|
|
|
7
|
-
なお，下記の記述例の最後に記載のように値の大小関係に合わせて `x` のレベルを設定し多重比較すれば値の大小関係とアルファベットを対応させることはできるようです
|
|
7
|
+
なお，下記の記述例の最後に記載のように値の大小関係に合わせて `x` のレベルを設定し多重比較すれば値の大小関係とアルファベットを対応させることはできるようです
|
|
8
|
-
|
|
9
8
|
|
|
10
9
|
```R
|
|
11
10
|
library(multcomp, quietly=TRUE)
|
|
@@ -14,8 +13,8 @@
|
|
|
14
13
|
N <- 10
|
|
15
14
|
tuk1 <- data.frame(
|
|
16
15
|
x = factor(rep(c("apple", "peach", "melon", "mango"), rep(N, 4))),
|
|
17
|
-
y = c(s
|
|
16
|
+
y = c(as.integer(rnorm(N, 25, 2.5)), as.integer(rnorm(N, 20, 2.5)),
|
|
18
|
-
s
|
|
17
|
+
as.integer(rnorm(N, 15, 2.5)), as.integer(rnorm(N, 30, 2.5))))
|
|
19
18
|
|
|
20
19
|
tuk1$x <- factor(tuk1$x, levels = c("apple", "peach", "melon", "mango"))
|
|
21
20
|
tuk2 <- aov(y ~ x, data = tuk1)
|
|
@@ -31,12 +30,12 @@
|
|
|
31
30
|
##
|
|
32
31
|
## Linear Hypotheses:
|
|
33
32
|
## Estimate Std. Error t value Pr(>|t|)
|
|
34
|
-
## peach - apple == 0 -4
|
|
33
|
+
## peach - apple == 0 -3.400 1.148 -2.962 0.0263 *
|
|
35
|
-
## melon - apple == 0 -9.
|
|
34
|
+
## melon - apple == 0 -9.400 1.148 -8.189 <0.001 ***
|
|
36
|
-
## mango - apple == 0
|
|
35
|
+
## mango - apple == 0 4.400 1.148 3.833 0.0026 **
|
|
37
|
-
## melon - peach == 0 -
|
|
36
|
+
## melon - peach == 0 -6.000 1.148 -5.227 <0.001 ***
|
|
38
|
-
## mango - peach == 0
|
|
37
|
+
## mango - peach == 0 7.800 1.148 6.795 <0.001 ***
|
|
39
|
-
## mango - melon == 0 1
|
|
38
|
+
## mango - melon == 0 13.800 1.148 12.021 <0.001 ***
|
|
40
39
|
## ---
|
|
41
40
|
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
|
|
42
41
|
## (Adjusted p values reported -- single-step method)
|