求助：在R中计算并可视化转移时间(TTM)的分组条形图

1. 数据读取与筛选 #

先完成数据加载和转移患者的筛选，需要用到readxl读取Excel文件，dplyr做数据过滤：

# 首次运行请先安装依赖包
install.packages(c("readxl", "ggplot2", "dplyr"))

# 加载所需包
library(readxl)
library(ggplot2)
library(dplyr)

# 替换成你的Excel文件实际路径
df <- read_excel("your_patient_data.xlsx")

# 筛选出发生转移的患者（Metastasis=1）
M1 <- df %>% filter(Metastasis == 1)

注意：R语言对变量名大小写敏感，你的数据中转移时间列是time，但原代码写的是Time，必须统一变量名才能正常运行。

2. 绘制目标条形图 #

你需要的是2个条形的分组图，通常用于展示接受/未接受免疫治疗的转移患者的转移时间汇总值（比如均值），代码如下：

# 将Immunotherapy转为因子类型，让x轴标签更直观
M1 <- M1 %>% 
  mutate(Immunotherapy = factor(Immunotherapy, 
                                levels = c(0, 1), 
                                labels = c("未接受免疫治疗", "接受免疫治疗")))

# 绘制带均值的条形图，可选添加误差棒展示数据波动
ggplot(M1, aes(x = Immunotherapy, y = time)) +
  stat_summary(fun = mean, geom = "col", fill = "#4285F4", width = 0.6) +
  stat_summary(fun.data = mean_se, geom = "errorbar", width = 0.2) +
  labs(x = "免疫治疗状态", y = "转移时间（TTM）", title = "转移患者转移时间与免疫治疗的关系") +
  theme_minimal()

如果你的需求是展示每个转移患者的个体转移时间条形，按免疫治疗分组排列，可通过position_dodge实现分组效果：

ggplot(M1, aes(x = Immunotherapy, y = time, fill = Immunotherapy)) +
  geom_col(position = position_dodge(width = 0.8), width = 0.7) +
  labs(x = "免疫治疗状态", y = "转移时间（TTM）", title = "转移患者转移时间分布") +
  scale_fill_manual(values = c("#4285F4", "#EA4335")) +
  theme_minimal()

关键注意事项 #

• 变量名必须和数据列名完全匹配（比如time和Time在R中是两个不同变量）
• 若Excel文件路径包含中文或空格，务必用引号将路径包裹
• 可通过修改fill参数调整条形颜色，theme系列参数自定义图表整体样式

内容的提问来源于stack exchange，提问作者JMD