Plots output from a phreg regression

Source: R/plot.phreg.R

Plot(s) of the hazard, density, cumulative hazards, and/or the survivor function(s) for each stratum.

# S3 method for phreg
plot(
  x,
  fn = c("haz", "cum", "den", "sur"),
  main = NULL,
  xlim = NULL,
  ylim = NULL,
  xlab = "Duration",
  ylab = "",
  col,
  lty,
  printLegend = TRUE,
  score = 1,
  fig = TRUE,
  ...
)

Arguments

x

A phreg object
fn

Which function should be plotted? Default is the hazard function(s).
main

Header for the plot
xlim

x limits
ylim

y limits
xlab

x label
ylab

y label
col

Color(s) for the curves. Defaults to black.
lty

Line type for the curve(s). Defaults to 1:(No. of strata).
printLegend

Logical, or character ("topleft", "bottomleft", "topright" or "bottomright"); if TRUE or character, a legend is added to the plot if the number of strata is two or more.
score

Multiplication factor for the hazard function.
fig

logical, should the graph be drawn? If FALSE, data is returned.
...

Extra parameters passed to 'plot' and 'lines'.

Value

No return value if fig = TRUE, otherwise the cumulative hazards function (coordinates), given fn = "cum".

Note

Reference hazard is given by the fit; zero for all covariates, and the reference category for factors.

See also

phreg

Author

Göran Broström

Examples


y <- rllogis(40, shape = 1, scale = 1)
x <- rep(c(1,1,2,2), 10)
fit <- phreg(Surv(y, rep(1, 40)) ~ x, dist = "loglogistic")
plot(fit)