frailtypack (version 2.4.1)

plot.frailtyPenal: Plot Method for a Shared frailty model.

Description

Plots estimated baseline survival and hazard functions from an object of class 'frailtyPenal'. Confidence bands are allowed.

Usage

## S3 method for class 'frailtyPenal':
plot(x, type.plot = "hazard", conf.bands=TRUE, pos.legend="topright",
 cex.legend=0.7, main, color, ...)

Arguments

x

A shared frailty model, more generally a frailtyPenal class object (output from calling frailtyPenal function).

type.plot

a character string specifying the type of curve. Possible value are "hazard", or "survival". The default is "hazard". Only the first words are required, e.g "haz", "su"

conf.bands

Logical value. Determines whether confidence bands will be plotted. The default is to do so.

pos.legend

The location of the legend can be specified by setting this argument to a single keyword from the list '"bottomright"', '"bottom"', '"bottomleft"', '"left"', '"topleft"', '"top"', '"topright"', '"right"' and '"center"'. The default is '"topright"'

cex.legend

character expansion factor *relative* to current 'par("cex")'. Default is 0.7

main

title of plot

color

color of the curve (integer or color name in quotation marks)

...

other unused arguments

Value

• Print a plot of a shared frailty model.

See Also

frailtyPenal

Examples

data(readmission)

# Shared frailty model

modSha <- frailtyPenal(Surv(time,event)~as.factor(dukes)+cluster(id),
n.knots=10,kappa1=10000,data=readmission,Frailty=TRUE,hazard="Splines")

plot(modSha,type="surv",conf=FALSE)

# Cox proportional hazard model using Penalized likelihood

modCox <- frailtyPenal(Surv(time,event)~as.factor(dukes),n.knots=10,
kappa1=10000,data=readmission,Frailty=FALSE,hazard="Splines")

plot(modCox)

# no confidence bands
plot(modSha,conf.bands=FALSE)
plot(modCox,conf.bands=FALSE)