## exact to the year of death like we have in CenSoc. num <- flexsurv::pgompertz(y + 1, shape = beta, rate = rate.vec) - flexsurv::pgompertz(y, shape = beta, rate ...