5

在 R 中运行 Cox PH 回归后，我需要在数据框中添加预测风险比列。数据框是面板数据，其中 numgvkey 如果公司标识符和年龄是时间标识符。您可以从此链接下载一小部分日期： https ://drive.google.com/file/d/0B8usDJAPeV85VFRWd01pb0h1MDA/view?usp=sharing

我做了以下事情：

library(survival)
library(readstata13)
sme <- read.dta13("sme.dta")
reg<-coxph(Surv(age,EVENT2)~L1FETA+frailty(numgvkey), ties=c("efron"),  data=sme)
summary(reg)
hr <- predict(reg, type="risk")

如何在我的 'sme' 数据框中添加第 5 列“危险比”（hr）？另外，有什么方法可以预测 EVENT2 的概率而不是“hr”？

4

1 回答 1

6

该predict.coxph功能允许您生成几种不同“类型”的输出。其中之一是“预期”，这可能是您所说的“概率”。这并不是一个真正的概率，因为当相对风险、“基线风险”和观察时间很高时，数字有时会超过 1.0。

“类型”的“风险”选项返回风险比。

有一个 survfit.coxph 允许计算预测的生存率。它返回的对象同时 surv具有cumhaz列表组件。

你可能想试试这个：

sme$cumhaz <- survfit(fit, newdata=sme)$cumhaz
于 2015-05-28T22:32:19.543 回答