在R中绘制对数累积风险图：如何设置X轴为对数刻度？

你在验证比例风险假设（PH假设）时，用survival包生成的对数累积风险图X轴不是对数刻度，这里给你两种实用解决方法：

方法1：直接在绘图时启用对数X轴 #

这是最简便的方式，给plot()函数添加log="x"参数即可，同时优化代码写法让结构更清晰：

library(survival)

# 加载肺癌生存数据
lung <- lung

# 创建生存对象与拟合模型
surv_object <- Surv(lung$time, lung$status)
survfitobj <- survfit(surv_object ~ sex, data = lung)

# 绘制对数累积风险图，设置X轴为对数刻度
plot(survfitobj, 
     col = c("black", "red"),
     fun = "cloglog",
     log = "x")  # 关键参数：将X轴设为对数刻度

方法2：手动自定义对数X轴（适合精细控制刻度场景） #

如果需要指定特定的刻度值，比如聚焦关键时间点，可以先隐藏默认轴，再手动添加对数刻度：

library(survival)

lung <- lung
surv_object <- Surv(lung$time, lung$status)
survfitobj <- survfit(surv_object ~ sex, data = lung)

# 先绘制图表，关闭默认坐标轴
plot(survfitobj, 
     col = c("black", "red"),
     fun = "cloglog",
     axes = FALSE)

# 添加Y轴（保留默认样式）
axis(2)
# 添加自定义对数刻度的X轴，示例指定了几个关键时间点
axis(1, at = c(10, 30, 100, 300, 1000), log = TRUE)
# 补画图表边框
box()
# 可选：添加轴标签让图表更清晰
title(xlab = "生存时间（对数刻度）", ylab = "对数累积风险")

补充说明：当PH假设成立时，对数累积风险图中不同组的曲线应该大致平行，调整后的图表样式更便于你观察验证结果。

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