生存分析cox迴歸R語言與SAS以及STATA參數估計結果不同

進行生存分析的時候，單變量的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, ...)

在進行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)