Description Usage Arguments Value Author(s) References See Also Examples

Plots flexible hazard ratio curves allowing non-linear relationships between continuous predictors and survival. To better understand the effects that each continuous covariate has on the outcome, results are expressed in terms of hazard ratio curves, taking a specific covariate value as reference. Confidence bands for these curves are also derived.

1 2 3 |

`x` |
An object of class HR |

`predictor` |
Variable named in the formula or included as a predictor in the |

`prob` |
Value between 0 and 1. If |

`pred.value` |
Value from the variable |

`conf.level` |
Level of confidence. Defaults to 0.95 (corresponding to 95%). |

`round.x` |
Rounding of numbers in the plot. |

`ref.label` |
Label for the reference covariate. By default is the name of the covariate. |

`col` |
Vector of dimension 3 for the colors to plot. |

`main` |
These arguments to title have useful defaults here. |

`xlab` |
The range of x and y values with sensible defaults. |

`ylab` |
The range of x and y values with sensible defaults. |

`lty` |
Vector of dimension 2 for the line type. |

`xlim` |
The range of x and y values with sensible defaults. |

`ylim` |
The range of x and y values with sensible defaults. |

`xx` |
Vector of values (from the variable |

`...` |
Other arguments. |

No value is returned.

Artur Araújo and Luís Meira-Machado

Carmen Cadarso-Suarez, Luis Meira-Machado, Thomas Kneib and Francisco Gude. Flexible hazard ratio curves for continuous predictors in multi-state models: a P-spline approach. Statistical Modelling, 2010, 10:291-314

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ```
# Example 1
library(survival)
data(whas500)
fit <- coxph(Surv(lenfol, fstat)~age+hr+gender+diasbp+pspline(bmi)+pspline(los),
data=whas500, x=TRUE)
hr1 <- smoothHR(data=whas500, coxfit=fit)
plot(hr1, predictor="bmi", prob=0, conf.level=0.95)
# Example 2
hr2 <- smoothHR( data=whas500, time="lenfol", status="fstat", formula=~age+hr+gender+diasbp+
pspline(bmi)+pspline(los) )
plot(hr2, predictor="los", pred.value=7, conf.level=0.95, xlim=c(0,30), round.x=1,
ref.label="Ref.", xaxt="n")
xx <- c(0, 5, 10, 15, 20, 25, 30)
axis(1, xx)
``` |

```
Loading required package: survival
Loading required package: splines
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.