R中coxph公式使用survival::strata与strata结果差异的技术咨询

原本认为以下两种coxph模型调用等价：
首先安装并加载survival包：

install.packages("survival")
library(survival)

模型1（预期分层效果） #

coxph(Surv(futime, fustat) ~ age + strata(rx), data=ovarian)

# 输出结果
Call:
coxph(formula = Surv(futime, fustat) ~ age + strata(rx), data = ovarian)

       coef exp(coef) se(coef)     z       p
age 0.13735   1.14723  0.04741 2.897 0.00376

Likelihood ratio test=12.69  on 1 df, p=0.0003678
n= 26, number of events= 12 

模型2（非预期的固定效应处理） #

coxph(Surv(futime, fustat) ~ age + survival::strata(rx), data=ovarian)

# 输出结果
Call:
coxph(formula = Surv(futime, fustat) ~ age + survival::strata(rx), 
    data = ovarian)

                             coef exp(coef) se(coef)      z       p
age                       0.14733   1.15873  0.04615  3.193 0.00141
survival::strata(rx)rx=2 -0.80397   0.44755  0.63205 -1.272 0.20337

Likelihood ratio test=15.89  on 2 df, p=0.0003551
n= 26, number of events= 12 

实际结果并不等价：使用survival::strata的模型2将rx视为固定效应变量，而非在不同分层水平内分别拟合模型。不确定这是：

• A) R语言公式内使用package::function调用的问题
• B) survival包的命名空间/依赖问题

需求：开发包时不想加载整个体积较大的survival包，希望找到正确实现分层模型的方法。

核心原因：R公式的解析特性 #

这是R公式解析的机制问题，而非survival包的bug。当在公式中使用package::function()形式的调用时，R的公式解析器不会将其识别为特殊的模型构造函数（比如strata这种用于声明分层的专用函数），而是将其当作普通函数调用处理，最终把生成的结果作为模型的协变量。

survival包的strata在library(survival)后能实现分层效果，是因为coxph的公式解析逻辑会主动匹配函数名是否为strata；但加上survival::前缀后，函数名变成了survival::strata，coxph无法识别这是需要处理的分层声明，因此将其当作普通变量处理。

可行解决方案 #

方案1：包开发标准做法——导入strata函数 #

在你开发的包的NAMESPACE文件中添加：

importFrom(survival, strata)

这样在包代码中直接使用strata(rx)即可，无需加survival::前缀，同时不会将整个survival包加载到用户环境中，符合R包开发的规范。

方案2：手动指定公式解析环境 #

如果不想修改NAMESPACE，可以通过构造公式字符串并指定解析环境的方式实现：

# 构造公式字符串
formula_str <- "Surv(futime, fustat) ~ age + strata(rx)"
# 转换为公式对象并指定解析环境为survival命名空间
formula_obj <- as.formula(formula_str, envir = asNamespace("survival"))
# 调用coxph
survival::coxph(formula_obj, data = ovarian)

这种方式让公式在survival的命名空间内解析，确保strata被正确识别为分层函数。

方案3：局部环境临时绑定strata #

在局部环境中临时将strata绑定到survival::strata，避免污染全局环境：

local({
  strata <- survival::strata
  survival::coxph(Surv(futime, fustat) ~ age + strata(rx), data = ovarian)
})

局部环境的绑定会让公式解析时优先找到survival::strata，从而实现正确的分层效果。

以上三种方案都能得到和模型1完全一致的分层模型结果，其中方案1是包开发的首选方式，代码可读性和规范性更强。

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