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

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

since2016 ときどきTEXのメモ

stripchartの横にエラーバー

サンプル

ID <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
治療 <- c("A","B","A","B","B","A","A","A","B","B","A","B","A","A","A","B","A","B","B","B")
BP変化量 <- c(治療後BP - 治療前BP)
dat1 <- data.frame(id, 治療, 治療前BP, 治療後BP, BP変化量)

標準偏差を挿入

stripchart(BP変化量~治療, data = dat1, method="jitter", 
    jitter = 0.1, vertical=TRUE, pch=20, cex=1.5, xlim=c(0.5,2.5), ylim=c(-20, -0),
    xlab="治療", ylab="BP変化量")

par(new=T) #図を追加

Aの平均 = meanAB[[1]] #ここに平均を代入します
Bの平均 = meanAB[[2]]

Aの標準偏差 = sdAB[[1]]
Bの標準偏差= sdAB[[2]]

dplot <- plot(c(1.15,2.15), c(Aの平均, Bの平均), ylim=c(-20, -0), xlim = c(0.5, 2.5), xaxt="n",xlab="", ylab="")

arrows(1.15:2.15, c(Aの平均, Bの平均) - c(Aの標準偏差, Bの標準偏差), 1.15:2.15,  c(Aの平均, Bの平均) + c(Aの標準偏差, Bの標準偏差), lwd=1.0,angle=90,length=0.1,code=3)

f:id:yoshida931:20201218162617p:plain:w400

標準誤差を挿入

stripchart(BP変化量~治療, data = dat1, method="jitter", 
    jitter = 0.1, vertical=TRUE, pch=20, cex=1.5, xlim=c(0.5,2.5), ylim=c(-20, -0),
    xlab="治療", ylab="BP変化量")

par(new=T) #図を追加

Aの平均 = meanAB[[1]] #ここに平均を代入します
Bの平均 = meanAB[[2]]

dA <- dim(subset(dat1, dat1$治療=="A"))
dB <- dim(subset(dat1, dat1$治療=="B"))
Aの標準誤差 = sdAB[[1]]/sqrt(dA[1])
Bの標準誤差 = sdAB[[2]]/sqrt(dB[1])

dplot <- plot(c(1.15,2.15), c(Aの平均, Bの平均), ylim=c(-20, -0), xlim = c(0.5, 2.5), xaxt="n",xlab="", ylab="")

arrows(1.15:2.15, c(Aの平均, Bの平均) - c(Aの標準誤差, Bの標準誤差), 1.15:2.15,  c(Aの平均, Bの平均) + c(Aの標準誤差, Bの標準誤差), lwd=1.0,angle=90,length=0.1,code=3)

f:id:yoshida931:20201218163524p:plain:w400

参考
yoshida931.hatenablog.com