## Importing log-mortality data from Portugal:
data(PT)
Y <- PT
## Fitting the model
fit = blc(Y = Y, M = 100, bn = 20)
## Improvement:
improvement(fit)
improvement(fit, prob = 0.9) #90% credible intervals
Run the code above in your browser using DataLab