我正在使用 cph 包装器 {rms} 和 rcs 运行带有样条的 coxph 模型，如下所示：
y <- cph(Surv(time, event) ~ rcs(data$z, df=4) + Age_recruitment + factor(eth) + bmi + factor(dm) , na.action=na.omit, data=data)
恐怕我什至无法发布数据样本，因为它是敏感的健康信息。data$z 是一个连续变量。我收到如下错误：
Error in nact$nmiss : $ operator is invalid for atomic vectors
我理解这意味着 nact 是一个向量，因此不能与 $ 一起使用 - 但不确定为什么会发生这种情况，因为它是 cph 调用的“内部” - 即不是我创建的东西。如果出现问题，我尝试删除样条元素，但会发生同样的错误。任何人都可以为我阐明这一点吗？
我的目标是每次使用不同的 RCS 运行一系列 coxph 模型并有意义地输出结果。
谢谢。