frailtypack (version 2.2-12)

plot.nestedPenal: Plot Method for a Nested frailty model.

Description

Plots estimated baseline survival and hazard functions (output from an object of class `frailtyPenal' for nested frailty models). Confidence bands are allowed.

Usage

## S3 method for class 'nestedPenal':
plot(x, type.plot = "hazard", conf.bands=TRUE, ...)

Arguments

x

A nested model, more generaly an object of class frailtyPenal for Nested frailty models (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.

...

Other graphical parameters

Value

• Print a plot of the baseline survival or hazard functions with the confidence bands or not (conf.bands argument)

See Also

print.nestedPenal,summary.nestedPenal,frailtyPenal for Nested frailty models

Examples

data(dataNested)
  modNested<-frailtyPenal(Surv(t1,t2,event)~cluster(group)+subcluster(subgroup)+
  cov1+cov2,data=dataNested,n.knots=8,kappa1=50000)

# It takes around 90 minutes to converge (depends on the processor)  #

  plot(modNested,conf.bands=FALSE)