# NOT RUN {
library(ggplot2)
pcoxtheme()
data(heart, package="survival")
lam <- 0.02
alp <- 1
pfit <- pcoxtime(Surv(start, stop, event) ~ age + year + surgery + transplant
, data = heart
, lambda = lam
, alpha = alp
)
# Plot survival curves
psurv <- pcoxsurvfit(pfit)
plot(psurv)
# }
Run the code above in your browser using DataLab