[R] discrepancies between stata and r for a cox regression
mboutsen at ulb.ac.be
Mon May 18 17:50:55 CEST 2009
I would like to develop the use of R.
Trying R and more particulary the cox model, I am surprised by discrepancies between results with stata and R for a cox model
With the same data base, I get a hazard ratio (4.82) that is not the same obtained with stata (4.52)
You will find attached the file leukemia.dta I used (Stata)
Here are the codes for R
res <- coxph(Surv(TIME, STATUS)~TREAT, data=leukemia)
and here the codes for for stata
infix ID 2-3 TREAT 6 TIME 8-9 STATUS 11 LOGWBC 12-16 GENDER 18 using "g:rleukem.txt",clear
stset TIME, failure(STATUS==1)
SPSS and EPIinfo give the same HR than Stata
I tried with an other database without any problem
What would be the problem??? I changed of pc and versions of R (2.81 & 2.9.0) without any change. The means are the same for the two packages.
I saw a few posts with discrepancies but not with the same database
Thanks in advance
Department of Biostatistics
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
More information about the R-help