因子レベルが 10 のデータに対してポストホック Tukey HSD を実行しています。この表は巨大であり、p 値をペアごとの表で読者に提示することを望んでいました。45 行の表は付録に残しておきます。
データセットの例を次に示します。
set.seed(42)
x <- rnorm(100,1,2)
category <- letters[1:10]
data <- cbind.data.frame(x, category)
summary(data.aov <- aov(x~category, data = data))
data.hsd<-TukeyHSD(data.aov)
data.hsd.result<-data.frame(data.hsd$category)
data.hsd.result
結果は、45 行のテーブルです。代わりに、行と列の名前として因子レベルを持ち、セルに p 値があり、2 つが有意に異なるかどうかを示すテーブルが必要です。X やアンダースコアなど、繰り返しまたは不必要な比較を表す可能性のあるもの。このようなもの：
a b c d e f ... j
a X 0.97 1 0.99 0.89 0.99 ... 0.99
b X X 0.99 0.89 0.94 0.92 ... 0.97
c X X X 0.85 0.93 0.96 ... 0.98
| ... ... ... ... ... ... ... ...
i X X X X X X ... 0.84
等々。
このような表を自動的に作成する方法はありますか?