#  Figure 19.11. 


library(statmod)
x.values = seq(0, 3, by = 0.01)
h1<-dinvgauss(x.values, 1, 1)/(1-pinvgauss(x.values, 1, 1))
h2<-dinvgauss(x.values, 0.49, 1)/(1-pinvgauss(x.values, 0.49, 1))
h3<-dinvgauss(x.values, 1.69, 1)/(1-pinvgauss(x.values, 1.69, 1))

#  Set up graphics device and graphical parameters.
dev.new(width = 7, height = 7)
par(oma = c(1, 2, 1, 1))
axis.size <- 1.5
axis.values <- seq(from = 0, to = 3, by = 0.5)

plot(x.values, h1, type = "l", xlim = c(0, 3), ylim = c(0, 3.05),
     xlab = "Time", ylab = "Hazard",
     xaxt = "n", yaxt = "n", lwd = 2, cex.lab = 1.5)


lines(x.values, h2, lwd = 2)
lines(x.values, h3, lwd = 2)
axis(1, at = axis.values, cex.axis = axis.size)
axis(2, at = axis.values, cex.axis = axis.size)

dev.print(device = postscript, "19.10.eps", horizontal = TRUE)