单纯绘制分组的boxplot很简单,用fill指定组内的类别即可
your_data形状:
|group |subgroup |cor|
|--|--|--|
| ASD | IN |0.1|
|--|--|--|
| TD | ON |0.3|
ggplot() +
geom_boxplot(data = your_data, aes(x = group,fill = subgroup, y = cor),
color = "black", outlier.size = 0.2, width = 0.4, size = 0.4) +
scale_fill_manual(values = c('white', 'white','white', 'white'))
scale_fill_manual是指定四个Box都是白色。绘制结果为:
但是要加上errorbar和point就比较麻烦,errorbar和point用单独的数据框指定了group和subgroup
<- data.frame(
# subgroup = c("in","on","in","on"),
# group = c('ASD','ASD','TD','TD'),
COR = c(0.3083, 0.3086,0.2614, 0.2000),
upper = c(0.3587,0.3531,0.2980,0.2326),
lower = c(0.2579,0.2640,0.2247,0.1675),
)
但是这样在绘制errorbar和point时不能使用同样的分组,只能指定x = group,此时组内两个errorbar的x无法区分:
忘了因为什么还会提示不存在’group’
position_dodge也不太管用
最后解决办法:对于errorbar和point,指定位置而不是x轴名称:
cor_es_combined$positions = c(0.9,1.1,1.9,2.1)
geom_errorbar(data = cor_es_combined, aes(x = positions, ymin = lower, ymax = upper),
color = "black", width = 0.03, size = 0.6) +
geom_point(data = cor_es_combined, aes(x = positions, y = COR),
color = '#9F0900', size = 1)
这个position就是差0.1,多了就偏了
最后效果
这是目前想到的比较好的办法了,要不然只能根据Group和subgroup再打一列label将数据分为四类,然后做四个等间距的Box了,这个的errorbar和point倒是可以指定x,但间距无法调节