在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'？