一、环境
二、绘制bar图
2.1 什么是errorbar
首先你得知道这个errorbar是什么意思,简单来说,像一组10个样本的数据,该组数据计算后的均值±标准差,就是能够绘制出bar图。具体可以参考下面的代码和图片。l另外补充一点,既可以用标准差(SD)也可以用标准误(SE)来做计算bar值,她两的区别在于标准误会使bar值变的短些,看着更美观。
data <- data.frame(mean = c(10, 15), sd = c(12, 17))
rownames(data) <- c("case", "control")
par(lwd = 2)
b <- barplot(data\$mean, names.arg = rownames(data), col = c("red", "blue"), ylim = c(0, 25), axes = F, font = 2)
arrows(b[1], data\$mean[1], b[1], data\$sd[1], angle = 90)
arrows(b[2], data\$mean[2], b[2], data\$sd[2], angle = 90)
lines( x = c(b[1], b[1], b[2], b[2]), y = c( data\$sd[1] * 1.05 , data\$sd[2] * 1.1, data\$sd[2] * 1.1, data\$sd[2] * 1.05), lty = 2)
text( x = b[1] + (b[2] - b[1]) / 2, y = data$sd[2] * 1.1, label = "****", cex = 2, adj = c(0.5, 0))
axis(side = 2, lwd = 2, font = 2, cex = 1.5)
2.2 软件推荐
你如果是用R做某种分析,可以顺带上统计,如果你单纯想要统计,那我不推荐你用R,因为有更“香”的软件-----Graphpad(有免费资源的尽量别买🐶)。因为当你想用R做统计图时遇到种种问题,还要考虑代码报错的问题;而用Graphpad你只用输入数据,简直是太香了。关于Graphpad的教程,可以戳这里👉B站教程👈
以下是我用Graphpad做的图