理学療法士がまとめたノート

統計学備忘録(R言語のメモ)

since2016 ときどきTEXのメモ

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

アスタリスク 比較したラインを箱ひげの上限から2メモリ離した位置まで伸ばす・・・

治療 <- c("A","A","A","A","A","A","A","A","A","A","B","B","B","B","B","B","B","B","B","B")
治療前BP <- c(140,135,145,141,145,159,141,150,149,160,170,158,165,170,168,159,168,168,178,175)
治療後BP <- c(135,132,140,130,138,145,140,139,142,145,155,150,155,150,160,148,150,157,160,158)
BP変化量 <- c(治療後BP - 治療前BP)

bp <- boxplot(data$BP変化量~data$治療, ylim=c(-25,5))
#アスタリスクの下のライン
x1 <- 1 # x=Aの位置
x2 <- 2 # x=Bの位置
y1 <- 2 # yの位置
yh <- y1 - bp$stats[5,1] - 2 # x=A,y=y1から下に伸ばす距離(上限から2メモリ離す)
xx <- segments(x0=x1, y0=y1, x1=x2) 
text(1.5, 3, paste("**"),cex = 1.5)  #アスタリスクを書く位置
segments(x1, y1, x1, y1-yh)
segments(x2, y1, x2, y1-(yh + bp$stats[5,1]-bp$stats[5,2]))

[:plai]