我想知道是否有人实施过"预测" R中使用Firth校正(包coxphf)进行cox回归的函数。
使用coxphf进行的预测与使用coxph的预测相同吗?
谢谢!
以下是一个例子:
if
答案 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;本来是不负责任的。