使用plotRCS绘制生存样条图时出现rms::datadist常量列警告的问题求助

我来帮你分析下这个问题的原因，以及对应的解决办法：

为什么会出现这些警告？ #

你虽然没直接调用这些列，但plotRCS包的rcsplot函数底层依赖了rms包的建模逻辑，其中datadist()函数会扫描整个传入的数据集，而不是只扫描你指定的outcome、exposure等参数用到的列。

另外你提到这些列“只包含0或1”，但datadist提示“constant”说明这些列实际上是单值常量——要么全是0，要么全是1（可能你之前误解了，以为有两种值，但实际某一类没有样本）。你可以先验证这一点：

# 检查警告中列的唯一值数量
check_constant <- c("HAG_13", "HAG_13_21", "HAG_21", "AF_Flutter", "HF", "MI")
sapply(check_constant, function(col) {
  cat(col, "的唯一值数量：", length(unique(heart_data[[col]])), "\n")
})

如果输出都是1，就说明这些列确实是全同值的常量列。

怎么解决这个问题？ #

有两种简单的处理方式：

1. 提前移除数据集中的常量列（推荐）

既然这些列对模型没有任何贡献（全同值无法提供信息），直接从数据里删掉再运行代码就不会有警告了：

# 过滤掉所有常量列
heart_data_clean <- heart_data[, !sapply(heart_data, function(x) length(unique(x)) == 1)]

# 用清理后的数据集绘图
plot_object <- rcsplot(
  data = heart_data_clean,
  outcome = 'MORTALITY',
  time = 'Echo_to_censure_date',
  exposure = 'HAG',
  covariates = 'Age',
  positive = NULL,
  group = NULL,
  knots = c(0.25, 0.5, 0.75),
  knots.line = FALSE,
  ref.value = 0,
  ref.line = T,
  conf.int = T,
  conf.level = 0.95,
  conf.type = c("shape", "line"),
  pvalue = F,
  pvalue.digits = 3,
  pvalue.position = c(1, 1),
  pvalue.label.overall = "P for overall",
  pvalue.label.nonlinear = "P for nonlinear",
  fontsize = 10,
  fontfamily = "serif",
  linesize = 0.25,
  linecolor = "#0072B5FF",
  alpha = 0.1,
  xbreaks = c(-5, 0, 5, 10, 15, 20, 25, 30, 35),
  ybreaks = NULL,
  xlab = "Heart Age Gap (years)",
  ylab = "Mortality Hazard Ratio (95% CI)",
  explain = TRUE,
  xlim = c(-5, 30),
  ylim = c(0, 8),
  main = "Cardiovascular Mortality Hazard with increasing Heart Age Gap"
)

2. 手动指定datadist的变量范围（进阶）

如果你不想删除这些列，可以手动创建只包含你用到的变量的datadist对象，再传递给模型（不过plotRCS可能需要额外配置，不如第一种方法直接）：

# 只保留需要用到的变量
used_vars <- c("MORTALITY", "Echo_to_censure_date", "HAG", "Age")
dd <- datadist(heart_data[, used_vars])
options(datadist = "dd")

# 再运行rcsplot
plot_object <- rcsplot(...)

额外说明 #

这些警告本身不会影响你的生存样条图结果，只是datadist在提示你数据里存在无信息的常量列。清理数据不仅能消除警告，还能让建模过程更高效~

备注：内容来源于stack exchange，提问作者joeke