# NOT RUN {
data(ap)
(h0 = uh(.2, NULL, NULL, NULL, NULL, 15, 1)) # Uniform hazard
plot(h0, ylim=c(0,.3))
logLikuh(h0, ap)
r = Uhaz(ap, deg=2)
r$ll
logLikuh(r$h, ap)
plot(r$h, add=TRUE, col="red3")
# }
Run the code above in your browser using DataLab