[R] Relative Risk/Hazard Ratio plots for continuous variables
Frank E Harrell Jr
f.harrell at Vanderbilt.Edu
Tue May 25 14:41:23 CEST 2010
On 05/25/2010 07:28 AM, Laura Bonnett wrote:
> Dear all,
> I am using Windows and R 2.9.2 for my analyses. I have a large dataset and
> I am particularly interested in looking at time to an event for a continuous
> variable. I would like to produce a plot of log(relative risk) or relative
> risk (also known as hazard ratio) against the continuous variable.
Please use correct terminology. A risk is a probability (except perhaps
in finance) whereas a hazard is a rate (instantaneous conditional risk).
What you want is relative hazard.
> I have spent a long time looking for advice on how to do this but my search
> has proved fruitless - sorry if I've missed something obvious. It seems
> that there are options such as muhaz, survfit, coxph and cph that may enable
> some plots to be produced but none that specifically look at the relative
> risk one.
> In addition to the survival analysis, I have incorporated the mfp function
> (from package mfp).
> I currently use code such as,
> coxfit1<- coxph(Surv(rtime,rcens)~cts,data=data1)
> plot(coxfit1) nor plot(coxfit2) produce the relevant relative risk vs.
> continuous variable that I need.
Replace Design with the new rms package and see
The plot you need is a basic one provided by rms (or its ancestor Design).
> Can anyone help?
> Thank you,
> [[alternative HTML version deleted]]
> R-help at r-project.org mailing list
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
Frank E Harrell Jr Professor and Chairman School of Medicine
Department of Biostatistics Vanderbilt University
More information about the R-help