比較した箱ひげにアスタリスクを入れる(segments)

備忘録

2023.04.122024.10.26

head(dat)

box <- boxplot(dat$BP変化量~dat$治療, ylim=c(-25,5)) 
#四分位点
box$stats

segments関数：(x0, y1)から(x1, y1)に直線を引く

A == 1
B == 2
Y == 2

#横棒の描き方
xx <- segments(x0=A, y0=Y, x1=B, y1=Y) 

#縦棒の描き方
segments(
    x0=A, y0=Y, 
    x1=A, y1=box$stats[5,1] + 2
)

segments(
    x0=B, y0=Y,
    x1=B, y1=box$stats[5,2] + 2
)

#アスタリスクを書く位置
text(1.5, 3, paste("**"), cex = 1.5)