Greetings, I am new to R and have been comparing CPH survival analysis hazard ratios between R and SAS PhReg. The binary covariates' HRs are the same, however the continuous variables, for example age, have quite different HRs although in the same direction. SAS PhReg produces HRs which are the change in risk for every one increment change in the independent variable. How do I interpret the HRs produced by R? Thanks much, C Colleen Ross, MS Clinical Research Unit Kaiser Permanente Colorado (303) 614-1244 NOTICE TO RECIPIENT: If you are not the intended recipient ...{{dropped}}
Colleen.Ross at kp.org wrote:> Greetings, > > I am new to R and have been comparing CPH survival analysis hazard ratios > between R and SAS PhReg. The binary covariates' HRs are the same, however > the continuous variables, for example age, have quite different HRs > although in the same direction. SAS PhReg produces HRs which are the > change in risk for every one increment change in the independent variable. > How do I interpret the HRs produced by R? Thanks much, C > >I'm not aware of peculiarities. You're not giving us much to go on though. In fact, not even the function used to fit the model with.> ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >Exactly...
On Wed, 10 Jan 2007, Colleen.Ross at kp.org wrote:> I am new to R and have been comparing CPH survival analysis hazard ratios > between R and SAS PhReg. The binary covariates' HRs are the same, however > the continuous variables, for example age, have quite different HRs > although in the same direction. SAS PhReg produces HRs which are the > change in risk for every one increment change in the independent variable. > How do I interpret the HRs produced by R? Thanks much, C >What function did you use to fit the model in R? If you used coxph(), in the survival package then you should get the same answers as SAS. If you used cph() in the design package then the output says what the increment is that correponds to the quoted hazard ratio, and the default is the interquartile range. -thomas