r - 设置参考组，变化是使用 R 的 cox 回归中的值

Question

我是R中的新手，我想进行cox回归分析，如图所示，如何设置参考组并将HR更改为1，其他是参考组的具体值？

非常感谢，圣诞快乐！

Answer 1

无需担心协变量的哪个值是“参考”值，而是直接估计感兴趣的数量。例如，您可以估计并获得针对 x=7 的所有风险比的置信限。这是使用 Rrms包时的简单方法，它有一个函数，它是 R包函数cph的前端。survivalcoxph

require(rms)
dd <- datadist(mydata); options(datadist='dd')
# Example model with 3 predictors, interaction with the variable
# x1 that you are interested in.  So hazard ratios will depend on x2
# which we set to 10.  Omit x2 from list() if want to use the default
# value (median/mode of x2).  The following obtains hazard ratios for
# x1 = 6 7 8 9 against x1 = 7.
f <- cph(Surv(time, event) ~ x1 * x2 + x3, data=mydata)
contrast(f, list(x1=6 : 9, x2=10), list(x1=7, x2=10))
# Type ?contrast.rms to see how to plot these with confidence limits

要获得生存概率，请使用该survest函数，但请注意，您的列标题具有误导性，因为没有一个生存概率。