• R语言绘制环状柱状堆积图+分组+显著性

    无叠加、显著性的代码：

    1. #设置工作环境
    2. rm(list=ls())
    3. setwd("D:/Desktop/0000/code-main/条形图")
    4.  
    5. #加载R包
    6. library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics
    7. library(tidyverse) # Easily Install and Load the 'Tidyverse'
    8.  
    9. #加载数据（随机编写，无实际意义）
    10. df <- read.table("无叠加.csv", header = T, sep = ",")
    11. head(df)
    12. df$group <- factor(df$group,levels = c("M","N","G","F")) # 这几个字母表示不同的组，对应数据表格的第一个大组：group
    13. head(df)
    14.  
    15. # ----------------------------------------
    16. ##计算标签角度
    17. df2 <- as.data.frame(df[c(1:10,11:22,23:34,35:46,47:48),])  # 重点：分别为不同组的对应的行。注意对应数据理解
    18. rownames(df2) <- 1:48
    19. df2$group <- factor(df2$group,levels = c("M","N","G","F"))
    20. df2$ID <- as.numeric(rownames(df2))
    21. number_of_bar <- nrow(df2)
    22. angle <-  90 - 360 * (df2$ID-0.5) /number_of_bar
    23. df2$hjust<-ifelse(angle < -90, 1, 0)
    24. df2$angle<-ifelse(angle < -90, angle+180, angle)
    25. #创建标签数据及位置
    26. df3 <- df2 %>% 
    27.   group_by(group) %>% 
    28.   summarize(start=min(ID), end=max(ID) - 4) %>% 
    29.   rowwise() %>% 
    30.   mutate(title=mean(c(start, end)))
    31. df3$group <- factor(df3$group,levels = c("M","N","G","F"))
    32.  
    33. # ----------------------------------------
    34. #绘图
    35. ggplot()+
    36.   #手动添加辅助线
    37.   geom_hline(yintercept = 0,lty="solid", color = "black",linewidth=0.6)+
    38.   geom_hline(yintercept = 20,lty="solid", color = "grey80")+
    39.   geom_hline(yintercept = 60,lty="solid", color = "grey80")+
    40.   geom_hline(yintercept = 40,lty="solid", color = "grey80")+
    41.   geom_hline(yintercept = 80,lty="solid", color = "grey80")+
    42.   geom_hline(yintercept = 100,lty="solid", color = "grey80")+
    43.   #柱状堆积图绘制
    44.   geom_col(df, mapping=aes(group3, value, fill = group), color = "grey20", linewidth = 0.25, width = 0.8)+
    45.   # y轴范围确定
    46.   scale_y_continuous(limits = c(-25,150))+
    47.   # 颜色
    48.   scale_fill_manual(values = c("#4fbb98","#f46024","#dd6ab0","#7c8ebe"))+
    49.   # 主题
    50.   theme_void()+
    51.   theme(legend.position = 'none'
    52.   )+
    53.   # 手动添加显著性
    54.   geom_text(data=df2, aes(x=ID, y=103, label=c("Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  "," ",
    55.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  ",
    56.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  ",
    57.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  "),
    58.                           hjust=hjust,color=group),
    59.             fontface="bold", size=3,
    60.             angle= df2$angle, inherit.aes = F)+
    61.   
    62.   #手动添加坐标
    63.   geom_text(data=df2, x=-0.3,y=5,label="0",color="black",size=3)+
    64.   geom_text(data=df2, x=-0.3,y=25,label="20",color="black",size=3)+
    65.   geom_text(data=df2, x=-0.3,y=45,label="40",color="black",size=3)+
    66.   geom_text(data=df2, x=-0.3,y=65,label="60",color="black",size=3)+
    67.   geom_text(data=df2, x=-0.3,y=85,label="80",color="black",size=3)+
    68.   geom_text(data=df2, x=-0.3,y=105,label="100",color="black",size=3)+
    69.   #极坐标转换
    70.   coord_polar(direction=1)+
    71.   ##分组标签
    72.   geom_text(data=df3, aes(x = title, y = 140, label=group,color=group), 
    73.             hjust=c(1,1,0,0), angle=c(335,250,135,60), size=5, 
    74.             fontface="bold", inherit.aes = F)+
    75.   #颜色
    76.   scale_color_manual(values = c("#4fbb98","#f46024","#dd6ab0","#7c8ebe"))

    带叠加、显著性的代码：

    1. #设置工作环境
    2. rm(list=ls())
    3. setwd("D:/Desktop/0000/code-main/条形图")
    4.  
    5. #加载R包
    6. library(ggplot2) # Create Elegant Data Visualisations Using the Grammar of Graphics
    7. library(tidyverse) # Easily Install and Load the 'Tidyverse'
    8.  
    9. #加载数据（随机编写，无实际意义）
    10. df <- read.table("叠加.csv", header = T, sep = ",")
    11. head(df)
    12. df$group <- factor(df$group,levels = c("M","N","G","F")) # 这几个字母表示不同的组，对应数据表格的第一个大组：group
    13. head(df)
    14.  
    15. # ----------------------------------------
    16. ##计算标签角度
    17. df2 <- as.data.frame(df[c(1:10,51:62,103:114,155:166,207:208),]) 
    18. rownames(df2) <- 1:48
    19. df2$group <- factor(df2$group,levels = c("M","N","G","F"))
    20. df2$ID <- as.numeric(rownames(df2))
    21. number_of_bar <- nrow(df2)
    22. angle <-  90 - 360 * (df2$ID-0.5) /number_of_bar
    23. df2$hjust<-ifelse(angle < -90, 1, 0)
    24. df2$angle<-ifelse(angle < -90, angle+180, angle)
    25. #确定显著性标签位置
    26. result <- aggregate(value ~ group3, data = df, sum)
    27. #创建标签数据及位置
    28. df3 <- df2 %>% 
    29.   group_by(group) %>% 
    30.   summarize(start=min(ID), end=max(ID) - 4) %>% 
    31.   rowwise() %>% 
    32.   mutate(title=mean(c(start, end)))
    33. df3$group <- factor(df3$group,levels = c("M","N","G","F"))
    34.  
    35. # ----------------------------------------
    36. #绘图
    37. ggplot()+
    38.   #手动添加辅助线
    39.   geom_hline(yintercept = 0,lty="solid", color = "black",linewidth=0.6)+
    40.   geom_hline(yintercept = 20,lty="solid", color = "grey80")+
    41.   geom_hline(yintercept = 60,lty="solid", color = "grey80")+
    42.   geom_hline(yintercept = 40,lty="solid", color = "grey80")+
    43.   geom_hline(yintercept = 80,lty="solid", color = "grey80")+
    44.   geom_hline(yintercept = 100,lty="solid", color = "grey80")+
    45.   #柱状堆积图绘制
    46.   geom_col(df, mapping=aes(group3, value, fill = group), color = "grey20", linewidth = 0.5, width = 0.8)+
    47.   #y轴范围确定
    48.   scale_y_continuous(limits = c(-25,150))+
    49.   #颜色
    50.   scale_fill_manual(values = c("#4fbb98","#f46024","#dd6ab0","#7c8ebe"))+
    51.   #主题
    52.   theme_void()+
    53.   theme(legend.position = 'none'
    54.   )+
    55.   #手动添加显著性
    56.   geom_text(data=df2, aes(x=ID, y=103, label=c("Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  "," ",
    57.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  ",
    58.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  ",
    59.                                                "Contral", "ADOM", "LHy1","LHA", "LFA","SHy1","SHA","SFA","RHy1","RHA","  ","  "),
    60.                           hjust=hjust,color=group), 
    61.             fontface="bold", size=3, 
    62.             angle= df2$angle, inherit.aes = F)+
    63.   #手动添加标签
    64.   geom_text(data=df2, aes(x=ID, y=result$value+4, 
    65.                           label=c("", "", "***","", "","","**","","","","","",
    66.                                   "", "", "**","", "","","","***","","","","",
    67.                                   "", "***", "***","***", "**","***","**","***","**","***","","",
    68.                                   "", "***", "","", "","","","**","**","","",""),
    69.                           color=group), 
    70.             fontface="bold", size=4, 
    71.             angle= df2$angle+90, inherit.aes = F)+ #标签
    72.   #手动添加坐标及标题
    73.   geom_text(data=df2,x=11,y=30, label="Biodegradation rate(%)",color="black",size=3.5)+
    74.   geom_text(data=df2, x=-0.2,y=5,label="0",color="black",size=3)+
    75.   geom_text(data=df2, x=-0.2,y=25,label="20",color="black",size=3)+
    76.   geom_text(data=df2, x=-0.2,y=45,label="40",color="black",size=3)+
    77.   geom_text(data=df2, x=-0.2,y=65,label="60",color="black",size=3)+
    78.   geom_text(data=df2, x=-0.2,y=85,label="80",color="black",size=3)+
    79.   geom_text(data=df2, x=-0.2,y=105,label="100",color="black",size=3)+
    80.   #极坐标转换
    81.   coord_polar(direction=1)+
    82.   ##分组标签
    83.   geom_text(data=df3, aes(x = title, y = 140, label=group,color=group), 
    84.             hjust=c(1,1,0,0), angle=c(335,250,135,60), size=5, 
    85.             fontface="bold", inherit.aes = F)+
    86.   #颜色
    87.   scale_color_manual(values = c("#4fbb98","#f46024","#dd6ab0","#7c8ebe"))

  • 相关阅读:
    智能汽车能否真正实现无人驾驶,为什么?
    bind服务编译安装及压力测试
    5分钟掌握接口自动化测试，4个知识点简单易学！
    【C语言】详解 memset() 函数用法
    面试背诵版—操作系统
    Vue2和Vue3的响应式原理及区别
    SSO单点登录
    linux基础：4：gdb的使用
    design compiler中的drc规则详解
    计算机毕业设计（附源码）python疫情综合管控系统平台
  • 原文地址：https://blog.csdn.net/qq_45100200/article/details/133198468