簡単なエラーバーの描き方

plot(
    c(1, 3, 9, 10),
    c(5, 6, 2, 8),
    xlim=c(0.5, 10),
    ylim=c(0.5, 10)
)

arrows(
    c(1, 3), c(5, 6),
    c(9, 10), c(2, 8),
    col=c("black", "red")
)

Aの平均 = 2
Bの平均 = 3
Aの標準偏差 = 0.4  
Bの標準偏差 = 0.7

fig <- function(){
    plot(
    c(1,2), c(Aの平均, Bの平均),
    ylim = c(0,5), xlim = c(0.5, 2.5)
    )
}

fig()

fig()
arrows(
    c(1, 2), c(Aの平均-Aの標準偏差, Bの平均-Bの標準偏差),
    c(1, 2), c(Aの平均+Aの標準偏差, Bの平均+Bの標準偏差),
)

fig()
arrows(
    c(1, 2), c(Aの平均-Aの標準偏差, Bの平均-Bの標準偏差),
    c(1, 2), c(Aの平均+Aの標準偏差, Bの平均+Bの標準偏差),
    code=3
)

fig()
arrows(
    c(1, 2), c(Aの平均-Aの標準偏差, Bの平均-Bの標準偏差),
    c(1, 2), c(Aの平均+Aの標準偏差, Bの平均+Bの標準偏差),
    code=3, angle=90
)

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

# データ点の定義
x <- c(1, 2)
y <- c(2, 3)

# 平均値の計算
mean_A <- mean(y[1])
mean_B <- mean(y[2])

# 標準偏差の指定
sd_A <- 0.4
sd_B <- 0.7

# グラフを描画
plot(x, y, xlim = c(0.5, 2.5), ylim = c(1, 4), main = "散布図", xlab = "", ylab = "y", pch = 16, col = "blue", xaxt = "n", yaxt = "s")

# データ点をAとBで表示
points(x[1], y[1], pch = 16, col = "red", cex = 1.5, lwd = 2)
points(x[2], y[2], pch = 16, col = "green", cex = 1.5, lwd = 2)

# エラーバーを追加
arrows(x[1], mean_A - sd_A, x[1], mean_A + sd_A, angle = 90, code = 3, length = 0.1, col = "red")
arrows(x[2], mean_B - sd_B, x[2], mean_B + sd_B, angle = 90, code = 3, length = 0.1, col = "green")

# X軸の目盛りとラベルを設定
axis(1, at = c(1, 2), labels = c("A", "B"))

# 凡例を表示
legend("topright", legend = c("A", "B"), pch = 16, col = c("red", "green"), bty = "n")

# エラーバーのテキストを消す（空のテキストを追加）
text(x[1], mean_A, "", col = "red")
text(x[2], mean_B, "", col = "green")