用coxphf预测

时间:2016-03-12 03:58:18

标签: r statistics survival-analysis

我想知道是否有人实施过"预测" R中使用Firth校正(包coxphf)进行cox回归的函数。

使用coxphf进行的预测与使用coxph的预测相同吗?

谢谢!

以下是一个例子:

if

1 个答案:

答案 0 :(得分:2)

在搜索包含名为coxphf的函数的软件包并在包中找到相同名称的软件包之后,我安装了软件包并运行了帮助页面?coxphf中的第一个示例。然后我看到函数的结果是一个继承自coxph的对象,所以predict.coxph函数应该提供结果,它确实:

> library(coxphf); time<-c(1,2,3)
> cens<-c(1,1,1)
> x<-c(1,1,0)
> sim<-cbind(time,cens,x)
> sim<-data.frame(sim)
> myfit <- coxphf(sim, formula=Surv(time,cens)~x)
> class(myfit)
[1] "coxphf" "coxph" 
> predict(myfit)
[1]  0.7673587  0.7673587 -1.5347175

据推测,pkg:coxphf的作者认为这是合适的。否则添加类&#34; coxph&#34;本来是不负责任的。