考克斯按年分层

Question

我认为这应该是非常简单的事情，但我无法理解它。

我有以下代码：

library(survival)
cox <- coxph(Surv(SURV, DEAD)~YEAR, data)
summary(cox)

但我希望将结果分解为各个年份。

这是SPSS语法和解决方案的样子：

COXREG surv /STATUS=dead(1) /CONTRAST (year)=Indicator(1)
 /METHOD=ENTER year /PRINT=CI(95)
 /CRITERIA=PIN(.05) POUT(.10) ITERATE(20).
EXECUTE.

和STATA中的相同：

xi: stcox i.year

这是

的输出

str(data)

Answer 1

你没有向我们展示str（数据）或如何构建一个可重复的例子给出了“数据”。我怀疑“YEAR”将成为一个数字向量。如果它是一个因子变量，你会看到一个截距和n-1系数。因此，Interecpt系数将与“年”相同，而其他系数将与年（n）值匹配。你告诉SPSS引擎“年份”是一个“指标”，但是你没有给R引擎提供同样的礼貌。

试试这个：

data$year.ind <- factor(data$year)  # equivalent of SPSS INDICATOR
                                    # or SAS /CLASS
cox.mdl <- coxph(Surv(SURV, DEAD)~YEAR, data)
as.matrix(coef(coc.mdl)
summary(cox.mdl)

Answer 2

R经常拆分计算和显示结果以允许更多自由。我假设您需要coxph (?predict.coxph)的预测函数。

文档页面底部有一些示例，很可能是您想要的

predict(cox, type="terms")