R/plot_coh_pars.R

In ilc: Lee-Carter Mortality Models using Iterative Fitting Algorithms

plot_coh_pars <-
function(lca.obj){
  oldpar <- par(no.readonly=T)
  par(mfrow=c(1,2), mar=c(5, 5, 4, 1) + 0.1)
  plot(lca.obj$bx0, xlab='Age', ylab=expression(beta[x]^(0)), type='l')
  plot(lca.obj$itx, xlab='Year of birth', ylab=substitute(iota[t-x][' '] (a), list(a=lca.obj$adj)))
  title(paste(lca.obj$lab, names(lca.obj)[4], sep=' : '), outer=T, line=-1.5, cex.main=1.5)
  title(lca.obj$model, outer=T, line=-3, cex.sub=0.7, font.sub=3)
  invisible(par(oldpar))
}