データがあり、ANOVAとpost-hocテストを行う場合、Rの外でFigureを編集するのではなく、post-hoc分類を自動的に追加するboxplotを作成するには?anova boxplot上のポストホックラベルR
例えば、ここに始めるためにいくつかのデータは以下のとおりです。
install.packages("reshape", dependencies=T)
library(reshape)
x <- rnorm(30)
y <- rnorm(30)+1
z <- rnorm(30)+0.5
data.1 <- data.frame(x, y, z)
data.2 <- melt(data.1)
ここでは、単純な一方向ANOVAおよび計画外の比較です事後テストの実行中のすべてのコードです:で
linear.model <- lm(value~variable, data=data.2)
anova(linear.model)
# Analysis of Variance Table
# Response: value
# Df Sum Sq Mean Sq F value Pr(>F)
# variable 2 10.942 5.4710 5.8628 0.004087 **
# Residuals 87 81.185 0.9332
TukeyHSD(aov(linear.model))
# Tukey multiple comparisons of means
# 95% family-wise confidence level
# Fit: aov(formula = linear.model)
# $variable
# diff lwr upr p adj
# y-x 0.8344105 0.2396705 1.42915051 0.0034468
# z-x 0.2593612 -0.3353788 0.85410126 0.5539050
# z-y -0.5750493 -1.1697893 0.01969078 0.0602975
は、この時点で、グループ「a」のx、グループ「b」のy、グループ「a、b」のzを分類したいと考えています。私はボックスプロットを作ることができますが、どのようにそれを文字で注釈しますか?
boxplot(value~variable, data=data.2)
これは機能します。ベースのRグラフィックスで同じアノテーションを行う方法はありますか? – user1006644
はい、例えば、 'boxplot'関数の後に' text(x = 1:3、y = 3、c( "a"、 "b"、 "b、c")) ' – rcs
ありがとうJames&RCS! – user1006644