R可视化——基于ggplot2绘制环状柱状图

代码展示

1、加载R包

#设置工作环境rm(list=ls())setwd("D:\\环状柱形图")#加载R包library(tidyverse)library(reshape2)library(ggplot2)library(ggprism)

2、数据

#加载数据df <- read.table("data.txt",header = T, check.names = F)#转换数据data=melt(df)data$G<-rep(c("T","F","H"), each = 24)data_label <- datadata_label$ID <- as.numeric(rownames(data_label))

3、绘图——无分组情况

#计算标签角度number_of_bar <- nrow(data_label)angle <- 90 - 360 * (data_label$ID-0.5) /number_of_bardata_label$hjust<-ifelse(angle < -90, 1, 0)data_label$angle<-ifelse(angle < -90, angle+180, angle)#绘图p1 <- ggplot(data_label, aes(x=ID, y=value))+ geom_bar(stat="identity", fill="blue", alpha=0.7) + ylim(-75,75) +#y轴范围，控制内圆大小与条形大小 theme_minimal() +#主题 theme(axis.text = element_blank(), axis.title = element_blank(), panel.grid = element_blank(), plot.margin = unit(rep(-1,4), "cm")) +#调整边缘以使得标签不会被截断 coord_polar(start = 0) +#极坐标 geom_text(data=data_label, aes(x=ID, y=value+10, label=variable, hjust=hjust), color="black", fontface="bold",alpha=0.6, size=2.5, angle= data_label$angle, inherit.aes = F) #标签

p1

4、绘图——添加分组并增加分组间隔

#调整柱子显示顺序data_label = data_label %>% arrange(G, value)#设置分组间的空白间隔data_label$G<-as.factor(data_label$G)number_empty_bar <- 3to_add <- data.frame(matrix(NA, number_empty_bar*nlevels(as.factor(data_label$G)), ncol(data_label)) )colnames(to_add) <- colnames(data_label)to_add$G <- rep(levels(data_label$G), each=number_empty_bar)data_label <- rbind(data_label, to_add)data_label <- data_label %>% arrange(G)data_label$ID <- seq(1, nrow(data_label))##标签设置number_of_bar <- nrow(data_label)angle <- 90 - 360 * (data_label$ID-0.5) /number_of_bardata_label$hjust<-ifelse(angle < -90, 1, 0)data_label$angle<-ifelse(angle < -90, angle+180, angle)#绘图p2 <- ggplot(data_label, aes(x=ID, y=value, fill=G)) + geom_bar(stat="identity", alpha=0.5) + ylim(-75,75) + theme_minimal() + theme(legend.position = "none", axis.text = element_blank(), axis.title = element_blank(), panel.grid = element_blank(), plot.margin = unit(rep(-1,4), "cm")) + coord_polar() + #极坐标 geom_text(data=data_label, aes(x=ID, y=value+10, label=variable, hjust=hjust), color="black", fontface="bold",alpha=0.6, size=2.5, angle= angle, inherit.aes = F)+#标签 scale_fill_prism(palette = "candy_bright")#使用ggprism包修改颜色
p2

5、为分组添加标签

#创建标签数据及位置base_data <- data_label %>% group_by(G) %>% summarize(start=min(ID), end=max(ID) - number_empty_bar) %>% rowwise() %>% mutate(title=mean(c(start, end)))#绘图p2+geom_segment(data=base_data, aes(x = start, y = -5, xend = end, yend = -5), color = "red", alpha=0.8, size=0.8 , inherit.aes = F)+#添加分组线 geom_text(data=base_data, aes(x = title, y = -18, label=G), hjust=c(1,1,0), colour = "black", alpha=0.8, size=4, fontface="bold", inherit.aes = F)#分组标签

参考：https://r-graph-gallery.com/297-circular-barplot-with-groups.html

我就知道你“在看”