應各位老師的需求復現一篇文章的中的某個圖
具體復現圖5的整個思路圖，這里沒有原始數據，所以我使用虛擬生產的metadata進行畫圖
不廢話直接上代碼，先上python的代碼的結果圖
import matplotlib.pyplot as plt
import numpy as np# 數據（這里使用示例數據，您需要替換為實際數據）
data = {'IGHD': [20, 15, 25, 20, 15, 25, 10, 5, 10, 5, 2, 5, 2, 1, 5],'IGHM': [65, 70, 45, 65, 70, 45, 75, 80, 65, 70, 55, 25, 10, 5, 10],'IGHA1': [5, 5, 10, 5, 5, 10, 5, 5, 10, 10, 15, 30, 30, 35, 25],'IGHA2': [2, 2, 5, 2, 2, 5, 2, 2, 5, 5, 10, 10, 10, 10, 15],'IGHG1': [3, 3, 5, 3, 3, 5, 3, 3, 5, 5, 8, 15, 30, 30, 25],'IGHG2': [2, 2, 3, 2, 2, 3, 2, 2, 3, 2, 5, 10, 10, 10, 10],'IGHG3': [2, 2, 1, 2, 2, 1, 2, 2, 1, 2, 3, 3, 5, 5, 5],'IGHG4': [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5]
}# 設置顏色
colors = ['#ADD8E6', "#4169E1", '#90EE90', '#32CD32', '#FFC0CB', '#FF0000', '#800000', '#FFA500']# 創建圖表
fig, ax = plt.subplots(figsize=(12, 6))# 繪制堆疊柱狀圖
bottom = np.zeros(15)
for cell_type, percentage in data.items():ax.bar(range(15), percentage, bottom=bottom, width=0.8, label=cell_type, color=colors[list(data.keys()).index(cell_type)])bottom += percentage# 設置x軸標簽
x_labels = ['HD', 'EBnor', 'EBpro'] * 5
ax.set_xticks(range(15))
ax.set_xticklabels(x_labels, rotation=45)# 添加組標簽
group_labels = ['Btr', 'Bn', 'Bm', 'AtM', 'ASC']
for i, label in enumerate(group_labels):ax.text(i*3 + 1, 105, label, ha='center', fontweight='bold')# 設置y軸
ax.set_ylim(0, 110)
ax.set_ylabel('Fraction of cells (%)')# 添加圖例
ax.legend(bbox_to_anchor=(1.05, 1), loc='upper left')# 調整布局
plt.tight_layout()# 顯示圖表
plt.show()
貌似R畫的更好
library(ggplot2)
library(dplyr)
library(tidyr)# 創建示例數據
data <- data.frame(Group = rep(c("Btr", "Bn", "Bm", "AtM", "ASC"), each = 3),Condition = rep(c("HD", "EBnor", "EBpro"), 5),IGHD = c(20,15,25, 20,15,25, 10,5,10, 5,2,5, 2,1,5),IGHM = c(65,70,45, 65,70,45, 75,80,65, 70,55,25, 10,5,10),IGHA1 = c(5,5,10, 5,5,10, 5,5,10, 10,15,30, 30,35,25),IGHA2 = c(2,2,5, 2,2,5, 2,2,5, 5,10,10, 10,10,15),IGHG1 = c(3,3,5, 3,3,5, 3,3,5, 5,8,15, 30,30,25),IGHG2 = c(2,2,3, 2,2,3, 2,2,3, 2,5,10, 10,10,10),IGHG3 = c(2,2,1, 2,2,1, 2,2,1, 2,3,3, 5,5,5),IGHG4 = c(1,1,1, 1,1,1, 1,1,1, 1,2,2, 3,4,5)
)# 轉換數據格式并歸一化
data_long <- data %>%pivot_longer(cols = IGHD:IGHG4, names_to = "CellType", values_to = "Percentage") %>%group_by(Group, Condition) %>%mutate(Percentage = Percentage / sum(Percentage) * 100) %>%ungroup()# 設置細胞類型的順序
cell_type_order <- c("IGHG4", "IGHG3", "IGHG2", "IGHG1", "IGHA2", "IGHA1", "IGHM", "IGHD")
data_long$CellType <- factor(data_long$CellType, levels = cell_type_order)# 創建圖表
ggplot(data_long, aes(x = Condition, y = Percentage, fill = CellType)) +geom_bar(stat = "identity", position = "stack") +facet_grid(~ Group, scales = "free_x", space = "free_x") +scale_fill_manual(values = c("#FFA500", "#8B4513", "#FF0000", "#FF69B4", "#32CD32", "#90EE90", "#4169E1", "#ADD8E6")) +labs(y = "Fraction of cells (%)", x = NULL) +theme_minimal() +theme(axis.text.x = element_text(angle = 45, hjust = 1),panel.grid.major.x = element_blank(),panel.grid.minor.x = element_blank(),strip.background = element_blank(),strip.text = element_text(size = 12, face = "bold")) +scale_y_continuous(expand = c(0, 0), limits = c(0, 100))
)
