# NOT RUN {
#Long computing times
# }
# NOT RUN {
CBDfit <- fit(cbd(), data = central2initial(EWMaleData),
ages.fit = 55:89)
CBDResBoot <- bootstrap(CBDfit, nBoot = 500)
plot(CBDResBoot)
plot(CBDResBoot, parametricbx = FALSE, probs = seq(2.5, 97.5, 2.5))
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab