# 安装包
if (!requireNamespace("grafify", quietly = TRUE)) {
install.packages("grafify")
}if (!requireNamespace("dplyr", quietly = TRUE)) {
install.packages("dplyr")
}
# 加载包
library(grafify)
library(dplyr)
点图(标准差)
注记
Hiplot 网站
本页面为 Hiplot Point (SD)
插件的源码版本教程,您也可以使用 Hiplot 网站实现无代码绘图,更多信息请查看以下链接:
用于展示多组数据点标准差。
环境配置
系统: Cross-platform (Linux/MacOS/Windows)
编程语言: R
依赖包:
grafify
;dplyr
数据准备
# 加载数据
<- read.delim("files/Hiplot/142-point-sd-data.txt", header = T)
data
# 整理数据格式
<- "Doubling_time"
y <- "Student"
group <- factor(data[, group], levels = unique(data[, group]))
data[, group] <- data %>%
data mutate(median = median(get(y), na.rm = TRUE),
mean = mean(get(y), na.rm = TRUE))
# 查看数据
head(data)
Experiment Student Doubling_time facet median mean
1 Exp1 A 17.36765 F1 20.18114 19.91642
2 Exp1 B 18.04119 F1 20.18114 19.91642
3 Exp1 C 18.70120 F1 20.18114 19.91642
4 Exp1 D 20.06762 F1 20.18114 19.91642
5 Exp1 E 20.19807 F2 20.18114 19.91642
6 Exp1 F 22.11908 F2 20.18114 19.91642
可视化
# 点图(标准差)
<- plot_point_sd(data = data, Student, Doubling_time, symsize = 5,
p symthick = 0.5, s_alpha = 1, ewid = 0, symshape = 21,
all_alpha = 0) +
geom_hline(aes(yintercept = median), colour = 'black', linetype = 2,
size = 0.5) +
xlab(group) + ylab(y) +
guides(fill = guide_legend(title = group)) +
ggtitle("Point-SD") +
theme(text = element_text(family = "Arial"),
plot.title = element_text(size = 12,hjust = 0.5),
axis.title = element_text(size = 12),
axis.text = element_text(size = 10),
axis.text.x = element_text(angle = 0, hjust = 0.5,vjust = 1),
legend.position = "right",
legend.direction = "vertical",
legend.title = element_text(size = 10),
legend.text = element_text(size = 10))
p
