Life Insurance Contracts & Actuarial Present Value Valuations
Stochastic present value loss variable Z for life insurance contracts payable at the moment of death or at the end of the year of death; Continuous whole life insurance bar_A_x: Actuarial Present Value (APV / net single premium), integral formulation, and second moment ^2 bar_A_x; Discrete whole life insurance A_x: derivation using commutation functions and summation over life table cohorts; Term insurance contracts bar_A_x:n|^1 and A_x:n|^1: net single premium formulas and expiration without value; Pure endowment contracts A_x:n|^E = n_E_x = v^n n_p_x and endowment insurance contracts A_x:n| = A_x:n|^1 + n_E_x; Deferred insurances, varying insurance benefits (increasing whole life (IA)_x), and m-thly payable death benefits A_x^(m).
§4.1 The Stochastic Present Value Loss Variable for Life Insurance
A life insurance policy is a legally binding contract under which an insurer promises to pay a designated capital sum (the death benefit or sum insured $b_t$) upon the death of the insured life $(x)$, contingent upon death occurring within an agreed policy horizon.
Because human mortality is a stochastic biological phenomenon, both the timing of the claim payment and the number of years over which the insurer can invest premiums are random variables.
Let $T_x$ denote the future lifetime of $(x)$. In a continuous insurance contract paying a unit benefit of $\$1$ at the exact moment of death $T_x$, the present value of the benefit payment is a continuous random variable $Z$:
The expected value of $Z$, denoted $\mathbb{E}[Z]$, is the Actuarial Present Value (APV)—classically termed the Net Single Premium (NSP). The net single premium represents the exact capital lump sum the insurer must collect at policy inception to fund the promised death benefit, ignoring administrative expenses and profit loadings.
In discrete insurance contracts where claims are aggregated and paid at the end of the policy year of death, death occurs in year $K_x + 1$, where $K_x = \lfloor T_x \rfloor$ is the curtate future lifetime. The discrete present value random variable is:
§4.2 Continuous Whole Life Insurance: $\bar{A}_x$
A Continuous Whole Life Insurance pays a benefit of $1$ unit at the exact moment of death of $(x)$, whenever death occurs ($0 < T_x < \infty$).
Actuarial Present Value ($\bar{A}_x$)
The APV is the expected value of $Z = e^{-\delta T_x}$:
Second Moment and Variance of $Z$
To quantify the financial underwriting risk of the contract, the insurer must compute the variance of the present value random variable $\text{Var}(Z)$:
Notice that $\mathbb{E}[e^{-2\delta T_x}]$ has the exact same mathematical form as $\bar{A}_x$, but evaluated at a force of interest equal to twice the original force ($2\delta$). In actuarial notation, this second moment is denoted ${}^2\bar{A}_x$:
Therefore:
Under constant force of mortality $\mu$ and constant force of interest $\delta$:
§4.3 Discrete Whole Life Insurance: $A_x$ and Commutation Functions
In a Discrete Whole Life Insurance, a unit benefit of $1$ is paid at the end of the year of death: $Z = v^{K_x + 1}$.
Actuarial Present Value ($A_x$)
Summing over all possible discrete integer death years $k = 0, 1, 2, \dots$:
The second moment is:
Classical Commutation Functions: $D_x, C_x, M_x$
Before digital computers, evaluating summations across life table columns was labor-intensive. Actuaries streamlined calculations by pre-computing Commutation Functions combining interest discounting and cohort survival numbers:
1. Discounted Cohort Survivors ($D_x$):
2. Discounted Cohort Deaths ($C_x$):
3. Accumulated Discounted Deaths ($M_x$):
Multiplying $A_x$ by $v^x \ell_x / (v^x \ell_x) = 1$:
The entire infinite summation collapses into the ratio of two single numbers looked up from a commutation table!
§4.4 Term Insurance, Pure Endowment, and Endowment Insurance
Not all life insurance policies cover the entire lifespan. Commercial insurance packages combine term coverage with survival endowments:
1. Term Life Insurance ($A_{x:\overline{n}|}^1$ and $\bar{A}_{x:\overline{n}|}^1$)
A Term Life Insurance pays a benefit of $1$ unit if and only if $(x)$ dies within a fixed policy term of $n$ years. If the insured survives past $n$ years, the policy expires with zero payout.
- Continuous:
- Discrete:
2. Pure Endowment ($A_{x:\overline{n}|}^{\;\;1}$ or ${}_n E_x$)
A Pure Endowment contract pays a benefit of $1$ unit at time $n$ if and only if $(x)$ survives to attained age $x + n$. If $(x)$ dies prior to time $n$, zero benefit is paid.
The present value random variable is:
The Actuarial Present Value is denoted ${}_n E_x$:
The second moment is:
3. Endowment Insurance ($A_{x:\overline{n}|}$ and $\bar{A}_{x:\overline{n}|}$)
An Endowment Insurance combines term insurance and pure endowment: it pays $1$ unit at the time of death if death occurs before $n$, OR pays $1$ unit at maturity $n$ if $(x)$ survives to time $n$:
§4.5 Deferred and Varying Life Insurance Contracts
1. Deferred Life Insurance (${}_{m|} A_x$)
A deferred insurance contract provides whole life protection commencing only after an initial deferral period of $m$ years. If $(x)$ dies during the first $m$ years, no benefit is paid.
2. Arithmetically Increasing Whole Life Insurance ($(IA)_x$)
An increasing insurance contract pays a death benefit of $k + 1$ units if death occurs in policy year $k + 1$ (benefit is $\$1$ in year 1, $\$2$ in year 2, $\$3$ in year 3, etc.):
Defining the second-order commutation function $R_x = \sum_{y=x}^\infty M_y$:
§4.6 Approximating Claims Payable at the Moment of Death
In real life insurance practice, death benefits are paid immediately upon receipt of proof of death (typically within a few weeks of mortality), rather than delayed until the end of the policy anniversary year.
To approximate continuous claims $\bar{A}$ from discrete life table values $A$, actuaries assume the Uniform Distribution of Deaths (UDD) over each single year of age:
Since $\frac{i}{\delta} = \frac{e^\delta - 1}{\delta} \approx 1 + \frac{\delta}{2} \approx 1 + \frac{i}{2}$, paying benefits immediately at death rather than at year-end increases the net single premium by approximately half a year's interest ($\sim 2\%$ to $4\%$).
For an endowment insurance:
Notice that the factor $i/\delta$ applies only to the term insurance death component, never to the pure endowment maturity component ${}_n E_x$ (which is always paid exactly at integer time $n$).
§4.7 Continuous-Time Markov Multi-State Models and Kolmogorov Forward Equations
In modern life and health contingencies, individuals transition dynamically between multiple health states. Standard examples include:
1. The Disability Income Model (Three States): State 0 (Active / Healthy), State 1 (Temporarily or Permanently Disabled), State 2 (Dead / Absorbing).
2. Continuing Care Retirement Communities (CCRCs): Independent Living, Assisted Living, Skilled Nursing, and Death.
1. Transition Intensities and Chapman-Kolmogorov Equations
Let $S(t) \in \{0, 1, \dots, J\}$ denote the state occupied by a life aged $x$ at duration $t \ge 0$. The process is assumed to satisfy the continuous-time Markov property: the future transition probabilities depend solely on the current occupied state, independent of past state occupancy history.
The instantaneous transition intensity from state $i$ to state $j$ ($i \ne j$) at age $x + t$ is:
The total exit intensity leaving state $i$ is:
Let $p_{x}^{ij}(s, t) = P(S(t) = j \mid S(s) = i)$ denote the transition probability from state $i$ at duration $s$ to state $j$ at duration $t$ ($s \le t$). By the law of total probability, the transition matrix satisfies the Chapman-Kolmogorov Equation:
2. The Kolmogorov Forward Differential Equations
Differentiating $p_{x}^{ij}(s, t)$ with respect to terminal duration $t$ yields the system of Kolmogorov Forward Differential Equations:
In matrix notation, letting $\mathbf{P}(s, t) = [p_{x}^{ij}(s, t)]$ and defining the transition rate generator matrix $\mathbf{Q}(t)$ with off-diagonal entries $q_{ij}(t) = \mu_{x+t}^{ij}$ and diagonal entries $q_{ii}(t) = -\mu_{x+t}^{i \bullet}$:
Subject to the initial condition $\mathbf{P}(s, s) = \mathbf{I}$ (the identity matrix).
Solving this coupled system of linear ordinary differential equations numerically (using Runge-Kutta fourth-order integration) permits exact valuation of complex health insurance riders, disability waiver-of-premium benefits, and critical illness policies.
§4.8 Dependent Decrement Mechanics and Uniform Distribution of Decrements
In multiple decrement theory, an individual aged $x$ can exit an active status due to $m$ competing causes of decrement (such as death, disability, voluntary withdrawal, or retirement). Actuaries must distinguish between:
1. The Multiple Decrement Table (Dependent Rates): Where all causes operate simultaneously in the real insured population, yielding observed decrement probabilities $q_x^{(j)}$.
2. Associated Single Decrement Tables (Independent Rates): Theoretical probabilities $q_x^{\prime (j)}$ reflecting the rate at which cause $j$ would operate if all other $m - 1$ competing causes were completely eliminated.
1. The Fundamental Integral Relations
The dependent probability of exiting due to cause $j$ over the unit age interval is:
where the total survival probability under all decrements combined is:
assuming decrements act independently.
The independent single decrement probability is:
2. Conversions under Uniform Distribution of Decrements
When converting between dependent probabilities $q_x^{(j)}$ and independent probabilities $q_x^{\prime (j)}$, actuaries invoke fractional distribution assumptions.
Assumption A: Uniform Distribution of Decrements in the Multiple Decrement Table: Assume that in the combined table, exits from each cause $j$ occur uniformly across the year:
Under this condition, the independent rate for cause 1 in a two-decrement model is given exactly by:
if forces of decrement are proportional, or by the standard binomial approximation:
Actuarial Interpretation: An individual who exits due to cause 2 is exposed to cause 1 for an average of half a year ($1/2$). Therefore, the effective exposure base for cause 1 is reduced by $\frac{1}{2} q_x^{(2)}$, requiring the dependent rate $q_x^{(1)}$ to be inflated by $1 / (1 - \frac{1}{2} q_x^{(2)})$ to recover the true net rate $q_x^{\prime (1)}$.
Assumption B: Uniform Distribution of Decrements in Each Associated Single Decrement Table: If each independent decrement $j$ is uniformly distributed in its own single decrement table (${}_t q_x^{\prime (j)} = t \cdot q_x^{\prime (j)}$), the exact dependent rate in a two-decrement model is:
These conversion identities ensure mathematical consistency when building pension actuarial valuation software from separate mortality, disability, and employee turnover studies.
§4.9 Joint Life Modeling with Copulas and Common Shock Models
Classical joint life actuarial mathematics assumes that the future lifetime random variables $T_x$ and $T_y$ of two spouses are mutually independent:
In empirical demographic reality, spouses share common socioeconomic environments, dietary habits, infectious disease exposures, and psychological distress. Empirical studies document two pronounced phenomena:
1. Positive Lifetime Association: Long-term married couples exhibit significant positive correlation in life expectancy ($\rho(T_x, T_y) > 0$).
2. The Broken-Heart Phenomenon (Bereavement Effect): Upon the death of one spouse, the surviving spouse experiences an acute spike in mortality hazard rate $\mu$, particularly during the first six to twelve months of widowhood.
1. Common Shock Models
The simplest structural method to introduce positive mortality dependence is the Marshall-Olkin Common Shock Model.
Let each individual's death be driven by independent idiosyncratic fatal shocks plus a simultaneous common shock:
- Individual $(x)$ is subject to an idiosyncratic Poisson shock process with hazard intensity $\lambda_1$.
- Individual $(y)$ is subject to an idiosyncratic Poisson shock process with hazard intensity $\lambda_2$.
- Both individuals are jointly exposed to an environmental common shock process (e.g., motor vehicle crash, natural catastrophe, toxic exposure) with hazard intensity $\lambda_{12}$.
The individual total mortality forces are:
The joint survival function is:
For equal durations ($t = s$):
Under independence ($\lambda_{12} = 0$), $S(t, t) = e^{-(\mu_x + \mu_y)t}$. Because $\lambda_{12} > 0$, the probability that both spouses survive simultaneously is strictly higher than predicted by the independence model!
2. Archimedean Copula Framework
By Sklar's Theorem, any multivariate joint survival function $S(t, s)$ can be uniquely decomposed into its marginal survival functions $S_X(t) = {}_t p_x$ and $S_Y(s) = {}_s p_y$ via a survival copula $C(u, v)$:
A widely utilized parametric copula in actuarial science is the Clayton Copula:
where $\theta$ measures the degree of lower tail dependence. As $\theta \to 0$, $C_\theta(u, v) \to u \cdot v$ (independence).
Another primary family is the Frank Copula:
which exhibits symmetric dependence and captures married couple survival patterns with exceptional empirical fidelity. Incorporating copula dependence prevents life insurers from mispricing joint life and last survivor annuities.
§4.10 Multi-State Actuarial Cash Flow Valuation: Active-Disabled-Death Continuum
To price complex health contingencies (such as long-term care policies, critical illness insurance, and permanent disability income riders), actuaries combine the transition probabilities $p_{x}^{ij}(s, t)$ derived from Kolmogorov forward equations with continuous financial discounting.
1. Actuarial Present Value of Benefit Cash Flows
Consider an active life aged $x$ at time $t = 0$ in State 0 (Active). The policy provides:
- An annuity benefit paid continuously at rate $b(t)$ while the policyholder is in State 1 (Disabled).
- A lump sum death benefit $S(t)$ paid immediately upon transition from any live state to State 2 (Dead).
The Actuarial Present Value (APV) of the disability annuity benefit is:
The Actuarial Present Value of the death benefit is:
The term inside the brackets represents the total probability density of entering the absorbing death state at duration $t$, decomposing into deaths occurring directly from the active state ($p_{x}^{00} \mu^{02}$) and deaths occurring from the disabled state ($p_{x}^{01} \mu^{12}$).
2. Premium Equivalence in Multi-State Models
Let the policyholder pay a level annual premium continuously at rate $\bar{P}$ while active (State 0). Premium payments cease immediately upon disability or death (a built-in Waiver of Premium Benefit).
The Actuarial Present Value of premium inflows is:
By the Equivalence Principle, equating expected premium inflows to expected benefit outflows:
This general formulation provides the mathematical foundation for modern health insurance reserving across all international regulatory jurisdictions.
§4.11 Reversionary Annuities and Contingent Insurance Valuation
Actuaries frequently design contracts whose benefit payouts depend on the specific order in which two or more lives die.
1. The Reversionary Annuity ($\bar{a}_{x|y}$)
A Reversionary Annuity is issued to two lives $(x)$ and $(y)$. It pays a continuous annuity at rate 1 per year to $(y)$ following the death of $(x)$, provided $(y)$ is still alive. The payments cease upon the subsequent death of $(y)$.
Let $I_y(t)$ be the indicator that $(y)$ is alive at time $t$, and let $I_x(t)$ be the indicator that $(x)$ is alive at time $t$. Payment occurs at time $t$ if and only if $(y)$ is alive and $(x)$ is dead:
Taking expectations discounted by $e^{-\delta t}$:
Integrating over all future durations $t \in [0, \infty)$:
This identity proves that a reversionary annuity can be synthesized as a long position in a single life annuity on $(y)$ combined with a short position in a joint life annuity on $(xy)$.
2. Contingent Life Assurances ($\bar{A}_{xy}^1$ and $\bar{A}_{xy}^2$)
A Contingent Life Assurance pays a benefit upon the death of one life, conditional upon the survival status of the other life:
- $\bar{A}_{xy}^1$: Pays 1 at the moment of death of $(x)$, provided that $(x)$ dies before $(y)$.
- $\bar{A}_{xy}^2$: Pays 1 at the moment of death of $(x)$, provided that $(x)$ dies after $(y)$.
The net single premium for $\bar{A}_{xy}^1$ is:
For $\bar{A}_{xy}^2$:
Summing the two contingent assurances:
Because $(x)$ must either die before $(y)$ or die after $(y)$, the sum of the two mutually exclusive contingent assurances reconstitutes the unconstrained single life whole life assurance on $(x)$.
§4.12 Joint-and-Survivor Annuity Reductions and Cash Refund Guarantees
In defined benefit pension schemes and individual commercial retirement planning, retiring couples rarely select a joint life annuity that terminates on the first death, because the surviving spouse would be left without income. Nor do they typically select an unreduced 100% joint-and-survivor annuity, which carries high premium costs.
Instead, the vast majority of retirement plans provide Joint-and-Survivor Annuities with Benefit Reductions.
1. The Joint-and-Survivor Reduction Model
Consider a retirement annuity issued to $(x)$ (the primary retiree) and $(y)$ (the contingent spouse). The contract pays:
- Benefit rate 1 per year while both $(x)$ and $(y)$ are alive.
- Benefit rate $\theta_1$ per year to $(x)$ if $(y)$ dies first ($0 < \theta_1 \le 1$).
- Benefit rate $\theta_2$ per year to $(y)$ if $(x)$ dies first ($0 < \theta_2 \le 1$).
- Zero payments after both lives have died.
In standard pension design, commonly selected reduction factors are $\theta_1 = 1.0$ (no reduction to the primary worker) and $\theta_2 = 0.50$ ($50\%$ joint-and-survivor) or $\theta_2 = 0.667$ ($66.7\%$ joint-and-survivor).
Let $I_x(t)$ and $I_y(t)$ be the survival indicators. The payment rate at duration $t$ is:
Rearranging terms:
Taking expectations and discounting:
For the symmetric case where the benefit reduces to fraction $\theta$ upon either spouse's death ($\theta_1 = \theta_2 = \theta$):
If $\theta = 1/2$ (a $50\%$ reduction upon the first death):
Remarkably, the joint life status $\bar{a}_{xy}$ cancels out completely, and the annuity value is simply the arithmetic average of the two single life annuities!
2. Cash Refund and Installment Refund Annuities
Retirees are often reluctant to purchase life annuities because of the perceived risk that if they die shortly after retirement, the insurer retains the remaining capital. Insurers resolve this behavioral friction by offering Refund Annuities:
- Cash Refund Annuity: If the annuitant dies before cumulative annuity payments equal the initial purchase premium $P$, the remaining balance $P - \sum \text{Payments}$ is refunded immediately as a lump sum death benefit.
- Installment Refund Annuity: If total payments are less than $P$, the periodic annuity payments continue to the designated beneficiary until cumulative payments equal $P$.
Actuarial valuation splits the refund annuity into a life annuity plus a decreasing term insurance benefit on $(x)$, providing capital certainty to risk-averse pensioners.
Step-by-Step Actuarial Statistics Solutions
Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.
An extract from an actuarial valuation table compiled at an effective annual interest rate of $i = 5.00\%$ ($i = 0.05$, $v = 1/1.05 \approx 0.952381$) reports the following commutation values:
- $D_{40} = 15,000$
- $D_{50} = 8,800$
- $D_{60} = 4,800$
- $M_{40} = 4,200$
- $M_{50} = 3,300$
- $M_{60} = 2,400$
- $R_{40} = 84,000$
A life aged $40$ purchases various life insurance policies. Compute:
- The Net Single Premium for a Whole Life Insurance of $100,000\text{ BDT}$ payable at the end of the year of death ($100,000 A_{40}$).
- The Net Single Premium for a 20-Year Term Insurance of $100,000\text{ BDT}$ ($100,000 A_{40:\overline{20}|}^1$).
- The Net Single Premium for a 20-Year Pure Endowment of $100,000\text{ BDT}$ ($100,000 {}_{20} E_{40}$).
- The Net Single Premium for a 20-Year Endowment Insurance of $100,000\text{ BDT}$ ($100,000 A_{40:\overline{20}|}$).
- The Net Single Premium for an Increasing Whole Life Insurance paying $10,000\text{ BDT}$ in Year 1, $20,000\text{ BDT}$ in Year 2, increasing by $10,000\text{ BDT}$ each year thereafter ($10,000 (IA)_{40}$).
Using $A_x = M_x / D_x$:
For a sum insured of $100,000\text{ BDT}$:
For a 20-year term from age 40, terminal age is $40 + 20 = 60$:
For a sum insured of $100,000\text{ BDT}$:
For a maturity benefit of $100,000\text{ BDT}$:
An endowment insurance is the sum of term insurance and pure endowment:
Equivalently, via commutation functions:
For a sum insured of $100,000\text{ BDT}$:
Using $(IA)_x = R_x / D_x$:
For a scale benefit of $10,000\text{ BDT}$ per step:
An insurance company sells identical independent 10-year term insurance policies to a portfolio of $N = 100$ lives, each aged $x = 50$. Each policy pays a lump-sum death benefit of $B = 250,000\text{ BDT}$ at the end of the year of death if the insured dies within 10 years.
The actuarial life table and interest rate parameters are:
- Effective annual interest rate: $i = 6.00\%$ ($i = 0.06$)
- Single policy term actuarial values:
Let $S = \sum_{j=1}^{100} Z_j$ denote the aggregate present value of claims across the entire portfolio, where $Z_j = B \cdot v^{K_{50, j} + 1} \mathbb{I}(K_{50, j} < 10)$.
Perform the following calculations:
- Compute the expected value $\mathbb{E}[Z]$ and standard deviation $\sigma_Z$ for a single policy.
- Compute the expected aggregate claims $\mathbb{E}[S]$ and standard deviation $\sigma_S$ for the portfolio of 100 policies.
- Compute the coefficient of variation ($CV = \sigma / \mathbb{E}$) for a single policy and for the 100-policy portfolio.
- Using the Central Limit Theorem, calculate the total capital reserve $C$ the insurer must establish at policy inception so that the probability of aggregate claims exceeding reserves is no more than $1.0\%$ ($P(S > C) \le 0.01$, with $Z_{0.01} = 2.326$).
- Compute the required contingency risk loading as a percentage of expected claims.
For a single policy with benefit $B = 250,000\text{ BDT}$:
The variance of a single policy loss variable is:
For $N = 100$ independent, identically distributed policies:
- Single policy:
- 100-policy portfolio:
Pooling 100 independent risks reduces relative volatility by a factor of $\sqrt{100} = 10\times$, illustrating the foundational economic law of insurance risk pooling!
By the Central Limit Theorem, aggregate claims $S$ asymptotically follow a normal distribution $S \sim \mathcal{N}(\mathbb{E}[S], \sigma_S^2)$:
The insurer must maintain $2,453,000\text{ BDT}$ in reserves to ensure solvency at a $99\%$ confidence level.
The insurer must collect a $50.95\%$ safety buffer over expected claims to protect against extreme statistical mortality clusters in a small portfolio of 100 lives.
A continuous whole life insurance policy pays $100,000\text{ BDT}$ at the moment of death of a life aged $x = 60$. Mortality follows De Moivre's Law with terminal age $\omega = 100$. The effective annual interest rate is $i = 6.00\%$ ($i = 0.06$).
Perform the following calculations:
- Determine the force of interest $\delta$ and the discount factor $v$.
- Compute the exact Actuarial Present Value $\bar{A}_{60}$ using the continuous integral definition under De Moivre's Law.
- Compute the discrete Actuarial Present Value $A_{60}$ assuming payments at the end of the year of death.
- Compute the UDD approximation $\bar{A}_{60} \approx \frac{i}{\delta} A_{60}$ and determine the approximation error relative to the exact integral.
- Compute the second moment ${}^2\bar{A}_{60}$ and the standard deviation of the present value loss variable.
For $x = 60$ and $\omega = 100$:
Under De Moivre's Law, the future lifetime $T_{60}$ is uniformly distributed on $[0, 40]$:
The financial interest parameters are:
Calculate $e^{-40 \delta}$:
For a $100,000\text{ BDT}$ benefit:
Under De Moivre's Law, $q_{60+k} = \frac{1}{40 - k}$, and ${}_k| q_{60} = \frac{1}{40}$ for each integer year $k = 0, 1, \dots, 39$:
Calculate $a_{\overline{40}|0.06}$:
The UDD multiplier is:
Approximating $\bar{A}_{60}$:
Notice that because De Moivre’s Law has a constant death density $f(t) = 1/40$, deaths are exactly uniformly distributed over each single year of age! Consequently, the UDD approximation is not merely an approximation—it is mathematically exact to 6 decimal places!
Evaluate at $2\delta = 2 \times 0.0582689 = 0.1165378$:
The variance is:
The standard deviation of the present value of the death claim is $24,994\text{ BDT}$.