hp.lifetab: Heligman-Pollard life table conversion

In HPbayes: Heligman Pollard mortality model parameter estimation using Bayesian Melding with Incremental Mixture Importance Sampling

Description

Generates a life table from the age-specific probabilities of death resulting from the estimation of the eight parameters of the Heligman-Pollard mortality model

Usage

hp.lifetab(hpp, nax, age = seq(0, 85, 1), l0 = 1e+05, 
with.CI=FALSE, CI=95)

Arguments

hpp	A vector of lenght 8 containing the Heligman-Pollard model parameter values
nax	A vector of equal length to age containing the appropriate nax values for each age interval in the life table
age	A vector containing the ages at which to estimate the life table columns
l0	The radix of the life table (Defaults to 100,000)
with.CI	Logical. If TRUE, produces two additional life tables each generated from the lower CI bound values of nqx and the upper bound values of nqx respectively
CI	Defines the width of the confidence interval. A value of CI=95 creates a 95 percent confidence interval with the upper bound defined as the 97.5th percentile nqx value for each age and the lower bound defined as the 2.5th percentile.

Value

A life table containing the columns "Age" (the starting age of the interval), "nax" (the set of nax values passed as an argument), "nqx" (the age-specific probabilites of dying in the interval x to x+n), "npx" (the probability of surviving the interval age x to x+n or 1-npx), "ndx" (the number of deaths in the interval), "lx" (the number of persons who have survived up to age x), "nLx" (the total number of person years lived within the interval x to x+n), "Tx" (the total number of person years left to live above age x+n), "ex" (life epectancy at age x).

References

Heligman, Larry and John H. Pollard. 1980 "The Age Pattern of Mortality." Journal of the Institute of Actuaries 107:49–80.

Examples

##requires a set of Heligman-Pollard parameters##
theta <- cbind(0.06008, 0.31087, 0.34431, 0.00698,
           1.98569, 26.71071, 0.00022, 1.08800)
age <- seq(0, 85, 1)
nax <- rep(.5, length(age))
HP.lt <- hp.lifetab(hpp=theta, nax=nax)

HPbayes documentation built on May 2, 2019, 5:53 a.m.