进行生存分析的时候，单变量的cox回归三个软件参数估计都是一样的，在进行多变量生存分析的时候，R与SAS和STATA的结果不一致，SAS和STATA是一样的，刚开始还以为R的参数设置错了，后来查看帮助文档
Usage
coxph(formula, data=, weights, subset,
na.action, init, control,
ties=c("efron","breslow","exact"),
singular.ok=TRUE, robust=FALSE,
model=FALSE, x=FALSE, y=TRUE, tt, method, ...)
na.action, init, control,
ties=c("efron","breslow","exact"),
singular.ok=TRUE, robust=FALSE,
model=FALSE, x=FALSE, y=TRUE, tt, method, ...)
在进行cox回归的时候，有三种参数估计的方法， c("efron","breslow","exact")，我把method改成”breslow“，R的结果就和SAS一样了~
y <- Surv(data2$day,data2$censor==1)
coxmodel <- coxph(y~td+v0+tr1+tr2+tr3,data=data2,method=("breslow"))
#生存分析中，协变量较多的情况下要用breslow的方法，默认是efron，在单变量的情况下efron和breslow的方法没有区别
summary(coxmodel)
step(coxmodel)版权声明：本文为youruolinmc原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接和本声明。