在 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”数据框中添加“危险比”(hr)的第 5 列?另外,有没有办法预测 EVENT2 概率而不是“hr”?
最佳答案
predict.coxph
函数允许您生成几种不同“类型”的输出。其中之一是“预期”,这可能就是您所说的“概率”。这不是一个真正的概率,因为当相对风险、“基线危险”和观察次数很高时,数字有时会超过 1.0。
“类型”的“风险”选项返回风险比。
有一个 survfit.coxph 可以让人们计算预测的存活率。它返回的对象同时具有 surv
和一个 cumhaz
列出组件。
你可能想试试这个:
sme$cumhaz <- survfit(fit, newdata=sme)$cumhaz
关于r - 在 R 中的 Cox 回归后将预测的危险比列添加到数据帧,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/30516275/