R语言生存分析报错：二元运算符的非数值参数

问题情况 #

• 执行base R生存曲线绘图代码时触发报错，提示non-numeric argument to binary operator，但使用ggsurvplot的代码却能正常运行。
• 报错代码：

plot(survfit(Surv(time,DEATH_EVENT) ~ hypertension, data=HF), main = "Hypertension Survival Distributions", xlab = "Length of Survival",ylab="Probability of Survival",col=c("blue","red")) +
  legend("topright", legend=c("Absent", "Present"),fill=c("blue","red"),bty="n")

• 报错信息：

Error in plot(survfit(Surv(time, DEATH_EVENT) ~ hypertension, data = HF),  : 
  non-numeric argument to binary operator

• 正常运行代码：

ggsurvplot(survfit(Surv(time,DEATH_EVENT) ~ hypertension, data=HF),
           data = HF, 
           censor.shape="|",  
           conf.int = FALSE, 
           ggtheme = theme_bw())

报错原因 #

核心是语法规则混淆：

• plot()是base R的原生绘图函数，执行后不会返回可链式操作的对象，它的作用是直接在绘图设备上生成图像，不支持ggplot2的+拼接语法。
• 你试图用+把base R的legend()函数和plot()的执行结果绑定，这完全违背了base R的语法逻辑，因此触发报错。
• ggsurvplot是survminer包中基于ggplot2构建的函数，天然支持+语法来添加图层元素，所以代码能正常运行。

解决方法 #

将legend()单独作为一行代码执行，不要用+和plot()绑定，修改后的代码如下：

# 先绘制生存曲线
plot(survfit(Surv(time,DEATH_EVENT) ~ hypertension, data=HF), 
     main = "Hypertension Survival Distributions", 
     xlab = "Length of Survival",
     ylab="Probability of Survival",
     col=c("blue","red"))
# 单独添加图例
legend("topright", legend=c("Absent", "Present"),fill=c("blue","red"),bty="n")

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