我在做癌症survival analysis的cure model,用R function: smcure(),
具体code如下:
smcure(Surv(RFS,Recurrence)~SNP+age65+Var1+Var2+Var3,cureform=~SNP+age65+Var1+Var2+Var3,data=data_R,model="ph",Var=TRUE,nboot=2000)
RFS是时间,Recurrence 是binary,即0,1。5个变量都是binary的,即0,1。这个运行时间超级长,几个小时,因为nboot比较大。结果由两个部分组成:Cure probability model 和 Failure time distribution model。
我想问问有没有人懂这个结果,如何解释这两个models, 这两个结果跟SAS的结果INCEPART和SURVPART是一样的吗?
还有,我怎么能得到变量的odds ratio跟hazard ratio estimation 和 confidence intervals?