R语言绘图

多组火山图

数据准备:

将CSV文件同一在一个路径下,用代码合并

确保文件列名正确

R 复制代码
library(fs)
library(dplyr)
library(tidyr)
library(stringr)
library(ggplot2)
library(ggfun)
library(ggrepel)

# 获取文件列表
file_paths <- dir_ls(path = "E:\\proteomic\\fuxian_bacteria\\multivolcano", pattern = "_DEG\\.csv$", full.names = T)
# 检查文件列表
print(file_paths)

# 读取和处理数据
all_deg <- map(file_paths, function(x){
  # 检查文件扩展名,忽略大小写
  if (str_ends(str_to_lower(x), ".csv")) {
    cat("Processing file:", x, "\n") # 打印正在处理的文件名
    data <- read_delim(file = x, col_names = T, delim = ",")
    if (is.data.frame(data)) {
      group <- str_remove(basename(x), pattern = "_DEG.*")
      return(dplyr::mutate(data, group = group))
    } else {
      cat("Failed to read file:", x, "\n")
      return(NULL)
    }
  } else {
    cat("Skipping non-CSV file:", x, "\n")
    return(NULL)
  }
}) %>%
  compact() %>% # 移除列表中的 NULL 元素
  do.call(rbind, .)

ggplot(data = all_deg) +
  geom_jitter(data = all_deg %>% dplyr::filter(change == "Normal"),
              aes(x = group, y = logFC, color = change, size = abs(logFC), alpha = abs(logFC)),
              width = 0.4) +
  geom_jitter(data = all_deg %>% dplyr::filter(change != "Normal"),
              aes(x = group, y = logFC, color = change, size = abs(logFC), alpha = abs(logFC)),
              width = 0.4) +
  geom_jitter(data = all_deg %>% dplyr::group_by(group) %>%
                dplyr::arrange(desc(abs(logFC))) %>%
                dplyr::slice_head(n = 15) %>%
                dplyr::ungroup() %>%
                na.omit(),
              aes(x = group, y = logFC, size = abs(logFC)),
              width = 0.4, shape = 21, fill = "#e77381") +
  geom_text_repel(data = all_deg %>% dplyr::group_by(group) %>%
                    dplyr::arrange(desc(abs(logFC))) %>%
                    dplyr::slice_head(n = 15) %>%
                    dplyr::ungroup() %>%
                    na.omit(),
                  aes(x = group, y = logFC, label = SYMBOL)) +
  geom_tile(aes(x = group, y = 0, fill = group), height = 0.4) +
  geom_text(data = all_deg %>% dplyr::select(group) %>% dplyr::distinct(group, .keep_all = T),
            aes(x = group, y = 0, label = group), size = 6) +
  geom_hline(yintercept = c(-log2(1.5), log2(1.5))) +
  scale_y_continuous(limits = c(-5, 5)) +
  scale_size(range = c(1, 10)) +
  scale_alpha(range = c(0.1, 1)) +
  scale_color_manual(values = c("Up" = "#f46d43", "Normal" = "#bdbdbd", "Down" = "#3288bd")) +
  scale_fill_manual(values = c('#8dd3c7','#ffffb3','#bebada','#fb8072','#80b1d3')) +
  theme_bw() +
  theme(axis.text = element_text(color = "#000000", size = 12),
        axis.title = element_text(color = "#000000", size = 15),
        panel.grid = element_blank(),
        legend.background = element_roundrect(color = "#969696")) +
  labs(y = "log2(Fold Change)")
相关推荐
侃侃_天下2 天前
最终的信号类
开发语言·c++·算法
echoarts2 天前
Rayon Rust中的数据并行库入门教程
开发语言·其他·算法·rust
Aomnitrix2 天前
知识管理新范式——cpolar+Wiki.js打造企业级分布式知识库
开发语言·javascript·分布式
每天回答3个问题2 天前
UE5C++编译遇到MSB3073
开发语言·c++·ue5
伍哥的传说2 天前
Vite Plugin PWA – 零配置构建现代渐进式Web应用
开发语言·前端·javascript·web app·pwa·service worker·workbox
小莞尔2 天前
【51单片机】【protues仿真】 基于51单片机八路抢答器系统
c语言·开发语言·单片机·嵌入式硬件·51单片机
我是菜鸟0713号2 天前
Qt 中 OPC UA 通讯实战
开发语言·qt
JCBP_2 天前
QT(4)
开发语言·汇编·c++·qt·算法
Brookty2 天前
【JavaEE】线程安全-内存可见性、指令全排序
java·开发语言·后端·java-ee·线程安全·内存可见性·指令重排序
百锦再2 天前
[特殊字符] Python在CentOS系统执行深度指南
开发语言·python·plotly·django·centos·virtualenv·pygame