Academic journal article
*Journal of Risk and Insurance*

# Portfolio Choice and Life Insurance: The CRRA Case

## Article excerpt

ABSTRACT

We solve a portfolio choice problem that includes life insurance and labor income under constant relative risk aversion (CRRA) preferences. We focus on the correlation between the dynamics of human capital and financial capital and model the utility of the family as opposed to separating consumption and bequest. We simplify the underlying Hamilton--Jacobi--Bellman equation using a similarity reduction technique that leads to an efficient numerical solution. Households for whom shocks to human capital are negatively correlated with shocks to financial capital should own more life insurance with greater equity/stock exposure. Life insurance hedges human capital and is insensitive to the family's risk aversion, consistent with practitioner guidance.

INTRODUCTION AND MOTIVATION

There is a glaring disconnect between the way most financial advisers sell life insurance and how they sell or promote investment products such as mutual funds, stocks, or bonds. Aside from the regulatory environment and the different licenses required--namely, securities as opposed to insurance--these two financial decisions are presented as if there were a "separation theorem" that justified their relative invariance. Moreover, although the concepts of risk tolerance, risk aversion, and utility are ubiquitous in the lingo and marketing material of the securities industry, the same terminology rarely enters the dialogue in the life insurance arena. This is quite odd since historically the economics of insurance was the breeding ground for much of the development in utility theory.

This observation is more than just anecdotal. The insurance literature's starting point for the optimal face amount of life insurance is the original work by Solomon Huebner in the early 1900s, based on the concept of a human life value (HLV). This idea is also at the heart of much forensic economics and litigation, where the courts must determine the value of a lost life, for the purposes of compensation. See Zietz (2003) or Todd (2004) for a comprehensive review of the academic literature on life insurance. Most financial planning and investment textbooks advocate something called a needs analysis, where the family determines how much would be required to meet long-term expenses and other financial goals if the primary source of income were lost. Other authors focus on the present value of lost wages. Regardless of the precise mechanics, rarely is this decision presented as a portfolio allocation problem akin to the investment in stocks and bonds. We believe it should. After all, life insurance is a hedge against the loss of human capital and portfolio hedging decisions should be made jointly, not independently.

Therefore, our objective in this article is to jointly analyze the decision of how much life insurance a family unit should have to protect against the loss of its breadwinner as well as how the family should allocate its financial resources between risk-free and risky assets vis-a-vis the dynamics of labor income and human capital in our (1) model.

We view this problem within the paradigm of portfolio choice and use the tools of financial economics pioneered by Merton (1969, 1971) to arrive at optimal controls for investment, consumption, and life insurance--but where the family unit is placed front and center.

A number of recent articles have extended the set of decisions included in the portfolio choice problem to highlight the interaction and the risks faced by the household, broadly defined. Thus, for example, Goetzmann (1993), Yao and Zhang (2005), and as Cocco (2005) focus on the role of the housing portfolio; Campbell and Cocco (2003) focus on optimal mortgage choices; Cairns, Blake, and Dowd (2006) examine portfolio choice in defined contribution pension plans; Sundaresan and Zapatero (1997) examine the role of defined benefit pensions; whereas Dybvig and Liu (2004) and Bodie et al. (2004) examine consumption and investment decisions in a life-cycle model with habit formation, stochastic opportunity set, stochastic wages, and labor supply flexibility. They vary retirement dates but do not endogenize it; Jagannathan and Kocherlakota (1996) and Viceira (2001) stress the impact of aging; Faig and Shum (2002) are motivated by the demand for illiquid assets; Koo (1998) and Davis and Willen (2000) model the role of labor income; Dammon, Spatt, and Zhang (2001) focus their attention on capital gains and income taxes; and Heaton and Lucas (2000) focus on the role of entrepreneurial risk. Others go back to basics and extend portfolio choice models to include more sophisticated (and realistic) processes for investment returns, like Chacko and Viceira (2005), or time-varying and mean-reverting risk premiums, like Kim and Omberg (1996) or Detemple, Garcia, and Rindisbacher (2003).

Although the above list is not exhaustive, the unifying theme of the burgeoning portfolio choice literature is that a singular personal financial feature is highlighted and carefully modeled to provide relevant financial economic insights related to that feature. Some of these articles test or calibrate their models against real-world data, whereas others end their contribution with a normative model. Most of them use Merton (1969, 1971) as a starting point.

Our objective is similar. We start with a traditional diffusion model of asset dynamics and constant relative risk aversion (CRRA) utility of consumption, but we focus on modeling labor income and life insurance purchases. A central feature of our model is the correlation between innovations to the labor income process and financial returns. We are specifically interested in the interaction between the utility-maximizing demand for life insurance and the optimal consumption pattern and asset allocation when there is a nonzero correlation between these two critical state variables. We deviate somewhat from previous models by focusing on the family unit that derives utility from consumption and the risk to this family of losing the wage/income source. The emphasis on the family as a unit, which survives independent of the life status of the breadwinner, alleviates the need for specifying a separate utility of bequest. See Zietz (2003) for a review of the literature and the references therein that analyze life insurance from the perspective of the family as opposed to the survivors.

Life insurance has been previously analyzed in a financial context. Yaari (1965), Hakansson (1969), Fischer (1973), Campbell (1980), Lewis (1989), Hurd (1989), Babbel and Ohtsuka (1989), and more recently, Purcal (2003) develop a number of insights into the demand for life insurance. These articles illustrate how the nature of mortality uncertainty, the utility and strength of bequest desire, risk aversion, increasing household wealth, and intertemporal rates of substitution impact the decision to purchase life insurance and the type of life insurance to buy. But these models do not give us much insight into the role of labor income dynamics and human capital over the life cycle. The same is true of the work by Richard (1975), which extended Merton (1969, 1971) to life insurance but with little numerical insight or analytic solutions, given the complexity of the problem. Along the same lines, Buser and Smith (1983) adopt a portfolio approach to life insurance but from a strictly one-period framework. Our contributions will be emphasized when we present the specifics of our model. Likewise, our work differs from the recent article by Chen et al. (2005), where the analysis is done in discrete time and via simulations. One of the main differences between the current work and most previous research is that we avoid specifying a utility of bequest yet are able to talk about the demand for life insurance. Our results confirm many of the earlier insights.

From a technical point of view, two state-variable problems--i.e., wage income and asset prices--normally lead to highly nonlinear Hamilton--Jacobi--Bellman (HJB) partial differential equations (PDEs) that are difficult to solve. Despite this nonlinearity, closed-form solutions to our model can be obtained in special cases when wage income and asset prices are perfectly (positively or negatively) correlated. Furthermore, when this correlation is between -1 and +1, we obtain a variety of efficient approximate solutions. One of the keys to our solution methodology, and indeed one of the contributions of this article, is to locate a similarity reduction, which reduces the dimension of the HJB equation.

The method of similarity reductions is a technique for solving PDEs (see, e.g., Ockendon et al., 1999). The basic idea is to find a special combination of the variables so that a new PDE can be obtained, which is simpler than the original PDE. For example, in the classical one-dimensional diffusion equation, spatial (state) and time variables can be collapsed into a new variable by taking the ratio of the spatial variable and the square root of the time variable. This idea (or technique) has been used in the finance literature. For example, Merton (1969, 1971) solved the optimal-consumption and asset allocation problem by scaling the value function with a particular function of wealth. In the option pricing literature, the ratio of an asset and its average (maximum or minimum) value is used as a similarity variable for pricing lookback (and Asian) options (see, e.g., Wilmott, Dewynne, and Howison, 1993). Likewise, Kahl, Liu, and Longstaff (2003) used the ratio of asset and total wealth as a new variable, for solving an asset allocation problem with illiquid stock holdings. Thus, although locating similarity reductions is appealing, the practice is more an art than a science. We will elaborate on this in the "Solution Methodology" section.

Our main qualitative results are as follows. First, we find that the optimal amount of (face value) life insurance a family should have depends on volatility of the wage/income process as well as its correlation to investment returns. Higher volatility and higher correlations lead to lower, optimal, life insurance levels. Ceteris paribus, a tenured university professor needs more life insurance than a Wall Street investment banker, even if they both expect to earn the same wages over time. Furthermore, as one would expect, the risk profile of the wage/income process has a strong impact on the optimal mix between risky and risk-free assets for the family. This result should come as no surprise. The original work by Bodie, Merton, and Samuelson (1992), all the way to the recent work by Cocco, Gomez, and Maenhout (2005), have stressed the financial characteristics of human capital and how it serves as a substitute for financial capital. Our model easily reproduces this result and emphasizes the role of life insurance, which is a hedge for human capital. Finally, our model confirms the intuition shared by many financial advisors that the optimal amount of life insurance a family should hold is not sensitive to the risk aversion or risk tolerance of the family. With CRRA, the optimal demand for life insurance is insensitive to the degree of relative risk aversion, even though the family's financial portfolio is obviously quite sensitive to this parameter. This opens up the possibility of testing whether the actual holding of life insurance is, in fact, correlated with other proxies for a family's risk aversion. Our model provides justification for finding no such relationship.

The rest of the article is organized as follows. The setup of the problem is described in the "Model of the Family Life Cycle" section, followed by the setup of the HJB equation in the "The Hamilton--Jacobi--Bellman" section. Solution methodologies are discussed in the "Solution Methodology" section. In the "Special Case: [rho] = [+ or -] 1" section we focus on special cases when the wage and risky asset are perfectly correlated. The general case is discussed in the "General Case: -1 < 1" section. We further simplify the HJB equation using the similarity reduction technique, which is followed by the presentation of the numerical method for solving the reduced HJB equations, where the proper boundary conditions are also discussed. In the "Closed-Form Approximation Solution" section we present a closed-form approximate solution, which proves to be quite accurate in practice. In the "Numerical Results and Discussion" section, some results and the financial implications of these results are presented. We summarize the article in the "Conclusion." The derivations of the HJB equation and other nonessential derivations are placed in the appendices.

MODEL OF THE FAMILY LIFE CYCLE

The Input Variables

We will work with three dates of interest. The first date is the time horizon of the family [T.sub.H], which is exogenous and deterministic. The second date is the time of retirement [T.sub.R] < [T.sub.H]. The third date of interest is the death of the breadwinner, which is the end of the wage/income process taking place at a random stopping time, denoted by r. The wage/income process jumps to zero at the minimum of [T.sub.R] and [tau].

The family purchases short-term insurance on the life of the breadwinner, which is renegotiated and guaranteed renewable on an ongoing basis at a predetermined schedule that is driven by an instantaneous force of mortality (IFM) curve denoted by [[lambda].sub.y+t], where y is the age of the primary breadwinner at time zero. Let [I.sub.t] denote the insurance premium (in dollars) payable per unit time; a variable that is under the direct control of the family. One can think (2) of [I.sub.t] as a budget for insurance that induces a certain face value or death benefit [I.sub.t][[lambda].sub.y+t].

We let [M.sub.t] denote the market value of the family's assets including the value of all (risky) stocks and (risk-free) bonds on a mark-to-market basis. We assume that [M.sub.0] denotes the initial marketable wealth at time t = 0. The breadwinner works and converts labor and time into wages and income. A portion of this income is consumed and the remainder is saved in a diversified portfolio, consisting of risky stocks and safe bonds.

The variable [[alpha].sub.t] denotes the fractional allocation of the family's marketable wealth [M.sub.t] to the risky asset at time t. Thus, if [[alpha].sub.t] = 0, the family allocates all its marketable wealth--but no more--to risk-free bonds, and if [[alpha].sub.t] = 1, the family allocates all its marketable wealth--but no more--to risky stocks. The model also allows for [[alpha].sub.t] > 1 which would imply leverage. Thus, for example, [[alpha].sub.t] = 2 implies that 200 percent of net worth is invested in the risky stock. This is financed by borrowing 100 percent of wealth at some (constant) rate of interest, denoted by r. Recall that [[alpha].sub.t] is under direct control of the family and is one of the three choice variables in our model. We do not address liquidity constraints faced by younger families when borrowing. Many of our numerical results lead to large amount of leverage.

Let [c.sub.t] denote the instantaneousness consumption rate of the family per unit time. In general, our model is specified in real terms, and all parameters and choice variables will reflect this. The consumption rate is our third choice variable and [c.sub.t] is chosen to maximize the family's utility of consumption, which includes the breadwinner as well. We refer readers to the review article by Zietz (2003), and the various references therein, in which the goals of a family versus an individual is presented as a possible objective.

Even though a more realistic state-dependent instantaneous utility function can be used under our PDE-based methodology, we adopt the simpler CRRA utility function in this article. The precise functional form is

u(c) = 1/1 - [gamma] [c.sup.1 - [gamma]] (1)

for some positive coefficient of relative risk aversion [gamma] > 0.

Let [Z.sub.t] denote the market value of the risky asset (stock index, market portfolio) at time t, which follows

d [Z.sub.t] = [[mu].sub.m] [Z.sub.t]dt + [[sigma].sub.m] [Z.sub.t]d [B.sup.m.sub.t], (2)

where [[mu].sub.m] denotes the constant drift and am the volatility. The risk-free rate r, which is the rate at which the family can borrow money, must be lower than [[mu].sub.m] for economic equilibrium purposes.

Let [W.sub.t] denote the real wage/income rate of the family's breadwinner per unit time. This stochastic process is expected to increase over time and might be correlated with the risky asset. The wage process satisfies

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

where [tau] is the random time of death, [[mu].sub.w] denotes the constant drift, [[sigma].sub.w] denotes the constant volatility, and [B.sup.W.sub.t] denotes the Brownian motion driving the wages process [W.sub.t]. The Brownian motion [B.sup.m.sub.t] driving the risky asset is instantaneously correlated with the [B.sup.W.sub.t] driving the wage process via the relationship d <[B.sup.m.sub.t], [B.sup.W.sub.t]) = [rho][[sigma].sub.m][[sigma].sub.w]dt. This correlation variable [rho] is the primary focus of our numerical case study and results.

Note that if two individuals start working at time 0 and face the same ([[mu].sub.w], [[sigma].sub.w]) wage process, the expected present value of their income, i.e., the discounted value of their human capital, is identical. In other words, if both individuals start out earning $100,000 per year, and they expect it to grow by an inflation adjusted 5 percent each year, with a volatility of 3 percent, for example, then the discounted value of their human capital is identical. However, of course, the economic value of their human capital can be quite different. That will depend on the correlation between their wages and the market portfolio--and that is exactly our point. The amount of life insurance they will purchase should depend on the economic value of their human capital, not the discounted value of their human capital. Indeed, practitioners in the insurance industry tend to focus on the discounted value (i.e., seven times annual salary) in their insurance rules of thumb, whereas in truth, this might be an over/underestimate, depending on the correlation coefficient.

As mentioned earlier, the IFM (or hazard rate) is denoted by [[lambda].sub.y+t], where y is the initial age. The conditional probability of survival, from age y to age y + t, under the law of mortality defined by [[lambda].sub.y+t] can be computed via

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

where the notation on the left-hand side of the equation is standard in the actuarial literature. If the future-lifetime random variable is exponentially distributed, we have that Pr[[tau] [greater than or equal to] s] = [e.sup.-[lambda]s], and the hazard rate is constant at all ages [[lambda].sub.y+t] = [lambda]. In the general case, the function [[lambda].sub.y+t] is expected to increases with time (age). In our model, the family unit knows exactly how much they will have to pay for insurance--regardless of how much and when they want to purchase it--from the current time t, to the time of retirement [T.sub.R]. Thus, in addition to precluding whole life and other more complicated forms of insurance, we do not allow for stochastic mortality rates or antiselection effects that might complicate the insurance purchase problem. For now, we ignore insurance loadings, commissions, and profits, which can easily be handled by working with a loaded hazard rate [[??].sub.y+t] instead of a biological mortality rate [[lambda].sub.y+t].

The Financial Wealth Dynamics

Based on the wage/income process Wt and the evolution of the risky asset price [Z.sub.t], the family budget constraint for the wealth process [M.sub.t] follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

for t < min [[tau], [T.sub.R]], and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

for t > min [[tau], [T.sub.R]].

At the instant of death t = [tau], we add the death benefit by distinguishing between the value of marketable wealth one instant prior to death and one instant after. More precisely,

[M.sub.[tau]+] = [M.sub.[tau]-] + [I.sub.[tau]]/[[lambda].sub.y+[tau]], (7)

where [I.sub.[tau]]/[[lambda].sub.y+[tau]] is the death benefit or the face value of the insurance policy at time [tau] and [[lambda].sub.y+[tau]] is the IFM, or hazard rate.

Finally, the family-unit or household objective function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (8)

where [delta] is the subjective discount rate and u(c) denotes the instantaneous utility of consumption at time t. The family is searching for an asset allocation strategy [[alpha].sub.s], an insurance buying strategy [I.sub.s], and a consumption strategy [c.sub.s] that maximize the discounted value of utility of consumption between time t and the terminal horizon of the family unit.

As explained and justified earlier, we have "eliminated" the need for a utility of bequest by treating the family as one unit that is dependent on the breadwinner for their source of income. Thus, in contrast to the insurance models of Yaari (1965), Fischer (1973), or Campbell (1980), the primary breadwinner is not forced to decide how much "they love their family," relative to how much they "love themselves." Rather, the family allocates consumption across the entire horizon of the family, while protecting the wage/income flow using insurance. This assumption alleviates the need to measure the strength of bequest but assumes that the family derives the exact same level of utility, whether the primary breadwinner is alive or not.

THE HJB EQUATION

With the model's foundations and formulation in place, let

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

be the indirect utility function. Assume we can find the optimal control that satisfies Equation (9), then J satisfies the following HJB equation (see Appendix A for the detailed derivation):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (10)

where the [PHI] satisfies a simpler (fundamental) HJB equation but without wage income and life insurance:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

For CRRA utility, the analytic form of [PHI](M, t; T) can be obtained as

[PHI](M, t; T) = h(t; T)[M.sup.1-[gamma]]/1 - [gamma]' (12)

where h(t; T) satisfies the following ordinary differential equation

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

Hereafter, the prime symbol is used to denote the derivative with respect to time. The closed-form solution of Equation (13), with zero utility at time T, can be obtained as

h(t; T) = [e.sup.-[delta]t] [([e.sup.[xi](T - t) - 1/[xi]).sup.[gamma]], (14)

where

[xi] = - [delta]/[gamma] + 1 - [gamma]/[gamma] (r + [([[mu].sub.m] - r).sup.2]/2[gamma][[sigma].sup.2.sub.m]). (15)

After applying the first-order conditions, the HJB equation for J can be rewritten as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

The terminal condition is

J(M, W, [T.sub.R]) = [PHI](M, [T.sub.R]; [T.sub.H]). (17)

Finally, the optimal control [c.sup.*.sub.t] (consumption), [[alpha].sup.*.sub.t] (allocation to risky equity), and [I.sup.*.sub.t] (amount spent on life insurance purchases) are given by

[c.sup.*.sub.t] = [([e.sup.[delta]t][J.sub.M]).sup.1/sub.[gamma]], (18)

[[alpha].sup.*.sub.t] ([[mu].sub.m] -r)[J.sub.M] + [rho][[sigma].sub.m][[sigma].sub.w] [WJ.sub.WM]/ [[sigma].sup.*.sub.m][MJ.sub.MM] (19)

[I.sup.*.sub.t] = [[([J.sub.M]/h(t; [T.sub.H]).sup.1/[gamma]] - M] [[lambda].sub.y+t]. (20)

SOLUTION METHODOLOGY

From here our article proceeds in the following manner. Remember that at this point, all we have is a two-dimensional PDE, representing the two uncertain variables--wages and investment returns--for the solution we have denoted by J; the indirect utility function. The optimal consumption, asset allocation, and insurance are all functions of J and its various derivatives.

To our knowledge, it is impossible to obtain a closed-form expression for J under the general conditions. However, there are some special and unique cases in which J can be obtained explicitly. For example, the well-known Merton (1969, 1971) analytic solution is obtained when the wage parameters are all set to zero. Another of those special cases is when the correlation between shocks to wages and investment returns are perfect. This effectively creates only one source of economic uncertainty in the model, and in this case too we derive explicit answers for the optimal consumption, asset allocation, and insurance.

In the general case, however, when the correlation between the random variables is not perfect, which means that we truly have two financial sources of uncertainty, the situation is more complicated. Nevertheless, we are able to employ a so-called similarity reduction technique, which is explained later on, that helps us to reduce the dimensionality of the problem, leaving us with only a one-dimensional PDE.

Of course, even though we are left with a one-dimensional PDE for J, we are still not "out of the woods" because even this (simpler) PDE does not have an analytic solution. At that point our narrative proceeds along two distinct paths.

First, we employ numerical methods to solve the new (albeit reduced dimension) PDE by brute force. More specifically, we use the method of lines, which converts the PDE into a system of ODEs that are then solved using numerical integration techniques. Although this provides us with reasonable solutions for any given set of parameters, it is time consuming and hence prone to the usual stability problems. It is not ideal for full-scale real-time implementation unless a more efficient numerical scheme can be found and implemented.

Then, as an alternative path to the above computationally intensive approach, we formulate an analytic approximation to the implied utility. This approximating function for J--which we can solve more easily--is quite accurate when compared with numerical results for nonperfect wage/returns correlations. It is confirmed against the case of perfect correlation. This analytic approximation functions at its best when the underlying economic parameters fall in a certain region, which we describe later.

Finally, we provide some case studies using our analytic approximation, since the parameter values we are interested in fall within the required range. We believe that one of the contributions of this article is our methodological approach to solving the main PDE, which is a problem that vexes many other articles in the literature. We now describe the highlights of the solution and refer the reader to the appendices for more details.

Special Case: [rho] = [+ or -] 1

As shown in Appendix A, the solution of the HJB Equation (16) takes the form

J(M, W, t) = h(t)[(M + k(t)W).sup.1-[gamma]]/1 - [gamma], (21)

where h is given by (14) and k satisfies

k'/k + 1/k + [eta] = 0, (2)

with

[eta] = [[mu.sub.w] - r [[lambda].sub.y+t] - [beta]([[mu.sub.m] -r) (23)

and

[beta] = [rho][[sigma].sub.w]/[[sigma].sub.m. (24)

Using the terminal condition k([T.sub.R]) = 0, k can be solved as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

When [[lambda].sub.y+t] = [lambda], i.e., a constant mortality rate, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

Remark: Since [beta] = [+ or -] [[sigma].sub.w]/[[sigma].sub.m] for [rho] = [+ or -] 1, [eta] is the only parameter affected by the value of [rho]. For given values of W/M and t, the solution is solely determined by the value of [beta] or [eta], while other parameter values are fixed. This relationship illustrates how risk aversion does not impact the optimal face value of life insurance.

General Case: -1 < [rho] < 1

For the general case of correlated wage and risky asset, a closed-form solution cannot be obtained. However, we can define a new variable x = W/M and use this as our "similarity variable" without the restriction of a simple analytical form. Even though our model itself is valid when M is negative, for technical reasons we will restrict our discussion (and examples) to the case of M > 0 for the rest of this article. We note that this is not a serious restriction since total household wealth is usually positive in reality. In the rare event that M is negative, we can either use the closed-form approximate solution (see the "Closed-Form Approximate Solution" section) or solve the two-dimensional Equation (16) using numerical methods described in Wang (2006).

Similarity Reduction. We now show that the two-dimensional HJB equation can be reduced by one dimension via the similarity reduction technique. We start by letting

J(M, W, t) = [M.sup.1 - [gamma]]/1 - [gamma] F(x, t).

This leads to a new PDE,

[F.sub.t] + [c.sub.0](1 - [gamma])F + [c.sub.1]x[F.sub.x] + [c.sub.2][x.sub.2] [F.sub.xx] = 0, (27)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (28)

The terminal condition is now

F(x, [T.sub.R]) = h([T.sub.R], [T.sub.H]). (29)

The coefficients [c.sub.0], [c.sub.1], [c.sub.2], and the terminal conditions are functions of x only, so that the HJB equation permits the similarity solution, which is a function of x and t.

Numerical Method. First, we truncate the semi-infinite domain (x > 0) into a finite one (0 < x < X) and set up a uniform mesh with size [delta]x = X/N and grid points [x.sub.j] = j[delta] x for j = 0, 1, ..., N. Then, we use the method of lines to solve (27); i.e., we approximate F by [F.sub.j] on the grid points and its spatial derivatives by finite differences, which leads to a system of ordinary differential equations with respect to time on grid point [x.sub.j]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

for j = 1, 2, ..., N - 1. Here [c.sub.0], [c.sub.1], and [c.sub.2] are evaluated on the grid points, based on the formula in (28). Note that the first upwind method is used for the first-order derivative term in x.

Boundary Conditions. We observe that [x.sub.0] = 0 corresponds to W = 0, i.e., no wage income. Therefore, the solution is of Merton type (with the addition of insurance),

with F not a function of x. In such cases, F satisfies the ordinary differential equation

[F.sub.t] + [c.sub.0](1 - [gamma])F = 0

or

[F'.sub.0] + [c.sub.0](1 - [gamma])[F.sub.0] = 0. (31)

At [x.sub.N] = X (the truncation of the domain), the situation is more complicated. Motivated by the earlier closed-form solution, we postulate

[F.sup.1/1 - [gamma]] ~ 1 + kx,

which is a linear function of x, therefore [([F.sup.1/1 - [gamma]]).sub.xx] ~ 0. We approximate this by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (32)

The system of ordinary differential Equation (30) with boundary conditions (31) and (32) is then solved by an appropriate numerical integration technique. For example, we can use a semi-implicit method, where the coefficients [c.sub.0], [c.sub.1], and [c.sub.2] are evaluated at the previous time level and the spatial derivatives are treated implicitly. In this article, we use a built-in Matlab solver ode15s.m. (3)

Closed-Form Approximate Solution

Motivated by the solution for the special case [rho] = [+ or -] 1, we now seek a closed-form approximate solution for the general case -1 < [rho] < 1. To proceed formally, we define a new parameter

[epsilon] = 1 - [[rho].sup.2] [([gamma][[sigma].sub.m][[sigma].sub.w]/[[mu].sub.m] - r).sup.2] (33)

and consider the case when [epsilon] [much less than] 1. Note that the special cases considered earlier correspond to [epsilon] = 0.

From an economic point of view, one can think of [member of] as a measure of "incompleteness" of the market. Clearly, when either [rho] = 1 or -1 or [[sigma].sub.m] = 0 or [[sigma].sub.w] = 0, the value of [epsilon] = 0, and we are in a world in which there is only one source of economic risk. Likewise, when [gamma] = 0, the individual is risk neutral, i.e., various risk-states of nature are treated the same and [member of] = 0 as well. This can be viewed (with a bit of creativity and imagination) as yet another form of market completeness.

Based on the procedure outlined in Appendix B, we can write our solution in the form J = [J.sub.0] + O([epsilon]). For [J.sub.0], we have the following explicit form:

[J.sub.0](M, W, t) = [M.sup.1 - [gamma]]/1 - [gamma] h[(1 + kx).sup.1-[gamma]], (34)

where h and k are given earlier. For a constant mortality rate, they can be solved explicitly as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (35)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (36)

with [xi] and [eta] defined earlier by Equations (15) and (23). Recall that a constant mortality rate [[lambda].sub.t] = [lambda] implies that future lifetime is exponentially distributed.

NUMERICAL RESULTS AND DISCUSSION

We now present a case study for a realistic set of economic parameters. The parameter values we have chosen for our results are as follows: [[mu].sub.m] = 7 percent; [[sigma].sub.m] = 20 percent; r = 2 percent; [[mu].sub.m] = 1 percent; [[sigma].sub.m] = 5 percent; [M.sub.0] = 1, 200; [W.sub.0] = 50; [T.sub.R] = 10, 20, 30; [T.sub.H] = 40, 50, 60; [delta] = 2 percent; y = 35, 45, 55.

The financial justification for these numbers are as follows. We assume that the family can invest their financial wealth in a broadly diversified index-fund that is expected to earn an inflation-adjusted arithmetic mean of [[mu].sub.m] = 7 percent per annum, with a volatility of [[sigma].sub.m] = 20 percent. These numbers are consistent with other asset allocation research in the literature (e.g., Campbell and Viciera, 2002; Chen et al., 2005). The risk-free rate of r = 2 percent after inflation is also consistent with current conditions in the inflation-adjusted bond market. Wages grow in real terms by [[mu].sub.m] = 1 percent per annum approximately, with a standard deviation (4) of [[sigma].sub.w] = 5 percent. This is consistent with wage inflation being higher than price inflation but also allows for shocks to wages, such as unemployment or business cycle effects. Finally, we assume the family's discount rate--or subjective rate of time preference--is [delta] = r = 2 percent.

In terms of other embedded parameters, we assume that insurance prices (and mortality rates) are driven by a Gompertz law of mortality, which, at age 45, starts at [[lambda].sub.45+t] = exp{(45 + t - 86.3)/9.5}/9.5. These numbers are consistent with survival rates implicit in pension-based mortality tables, and the analytic law of mortality has been used in a variety of actuarial and insurance articles (e.g., Frees, Carriere, and Valdez, 1996).

[FIGURE 1 OMITTED]

Our first order of business is to confirm that for these parameters values, our analytic approximation for [rho] [not equal to] 1 is valid and accurate. Figure 1 provides a visual confirmation of this fact by comparing the optimal insurance purchase [I.sub.y+t] as a function of the ratio between wages [W.sub.0]and marketable wealth [M.sub.0], when the correlation parameter is [rho] = 0.5.

Indeed, based on the above-mentioned parameter values, we compute the value of [epsilon] = 0.35(1 - [[rho].sup.2]), which is small compared with unity. The numerical solution was obtained using 100 grid points on a spatial interval x = [0, 0.25] and a time step size of [10.sup.-5] on the interval of t = [0, 20]. Our numerical experiments reveal that the computation becomes unstable at a larger time step due to oscillations in the optimal value of [[alpha].sup.*.sub.t]. This severe stability constraint may be caused by the nonlinear nature of the boundary conditions. Despite the fact that small time step sizes have to be used, our numerical method is still competitive, since we have reduced the dimension of the HJB, compared with the method in Purcal (2003), where a two-dimensional HJB was solved directly.

Given the accuracy of the closed-form approximate solution for the particular set of parameter values used, we turn our attention to the financial implications of the solution. In particular, we compute a number of cases that demonstrate the impact of risk aversion and wage correlation on the optimal asset allocation and demand for term life insurance.

Table 1 displays the optimal allocation to stocks and the optimal face amount (death benefit) for life insurance, assuming the primary breadwinner has [T.sub.R] = 20 years to retirement and the family horizon is [T.sub.H] = 50 years. This individual is y = 45 years old, earns [W.sub.0]= 50,000 per year in real terms, and currently has [M.sub.0] = 200,000 saved.

To be consistent with how insurance is discussed in practice, we display the optimal face value of life insurance [I.sub.t.sup.*]/[[lambda].sub.y+t] instead of the amount spent on life insurance [I.sub.t.sup.*] or the ratio of wealth insurance [I.sub.t.sup.*]/[M.sub.t].

The main qualitative insights from our model are as follows. First, the optimal amount of life insurance and allocation to stocks depend on the correlation between wage shocks and the investment return on stocks. As we mentioned in the earlier part of the article, by varying only the correlation between wage shocks and market shocks, we are essentially comparing the demand for insurance and asset allocation for two individuals that have identical discounted value of wages (since the wage process parameters are the same). The economic value of their human capital is obviously different, as evidenced by the fact that the amount of life insurance is lower when the correlation is greater than zero. In other words, the lower amount of life insurance is a proxy for a reduced economic value of human capital. Indeed, if the economic value of human capital was the same, the amount of insurance would be the same as well.

For example, a 45-year-old who earns $50,000 per annum and who already has $200,000 in accumulated savings will optimally purchase between $762,000 and $965,000 of life insurance, depending on the correlation between wages (i.e., human capital) and markets (financial capital). If this individual is working in a job/career whose wages profile is countercyclical to the financial market ([rho] = -1), then the optimal amount of insurance is $964,000. On the other hand, if the wage profile is perfectly correlated with financial markets, the optimal amount of life insurance is a lower $762,000. The economic intuition for the impact of wage-market correlations on the optimal demand, for insurance, as in Chen et al. (2005), is that when human capital is highly correlated with financial capital, the utility-adjusted value (i.e., economic value) of human capital is much lower, and hence the family requires less life insurance. Also, when the wage process is highly correlated with the returns from the risky asset, the allocation to financial risky asset is reduced to counteract the existing (and implicit) allocation to risky assets within the wage process. A number of models developed in the cited articles (e.g., Bodie, Merton, and Samuelson, 1992) confirm this. Notice, however, that the optimal face value life insurance, [I.sub.t.sup.*]/[[lambda].sub.y+t], does not depend on the risk aversion of the family. For example, at zero wage-to-market correlation ([rho] = 0), the optimal face value is approximately $855,000, which is more than 16 times the breadinner's annual wage. This number is somewhat higher than the often-heard rule of thumb that people should have 5-8 times their annual income in life insurance.

Some might find it puzzling that our utility-based model does not predict even higher levels of life insurance for families that are more risk averse. On a theoretical level, it is important to remember that we did utilize fair actuarial pricing. In practice though, this invariance is fairly consistent with practitioner intuition, when it comes to life insurance. Namely, a family should be protected with a minimal level of insurance, regardless of how risk tolerant the family (or breadwinner) consider themselves. Of course, risk tolerance and risk aversion have a very strong impact on the demand for the risky asset [alpha].sub.t], as one would expect. Families that are more risk tolerant will allocate more financial and marketable wealth to risky equities, in highly leveraged amounts, as evidenced by our numerical values of [[alpha].sub.t] [much greater than] 1. It is not clear how more general HARA utility, highly loaded insurance premiums as well as liquidity and borrowing constraints would impact these results. Preliminary and ongoing research reported in Huang and Milevsky (in press) indicates that risk aversion does impact the demand for insurance in these more general cases, yet the practical numerical impact is relatively small.

Table 2 provides similar results but with the key difference that the individual is now y = 55 years old, and therefore [T.sub.R] = 10 years from retirement and [T.sub.H] = 40 years from the terminal horizon of the family unit. All other parameters are exactly the same so that we can focus on the impact of "aging" on the optimal face value (i.e., death benefit) of life insurance and the optimal asset allocation.

Table 2 reveals substantially lower values for the optimal asset allocation proportions [[alpha].sub.t], as well as the face amount of life insurance [I.sub.t.sup.*]/[[lambda].sub.y+t]. Intuitively, the individual is closer to retirement, with only 10 more remaining years of work and wages. In this case, the optimal amount of life insurance lies between $430,000 and $480,000, depending on the wage-market correlation [rho]. This is only eight times the annual wage of [W.sub.0]= 50,000. The optimal amount of life insurance coverage closely tracks the expected discounted value of future wages. The breadwinner has 10 more years of labor income and buys eight times annual income of life insurance. In Table 1 the breadwinner had 20 more years of labor income and he/she purchased 16 times annual income. In addition, once the family unit is 10 years away from retirement, the amount of financial wealth invested in the risky asset is greatly reduced. In fact, when the wage-market correlation [rho] > 0, the family eliminates leverage and stops borrowing, as evidenced by at < 1 values.

Finally, Table 3 takes us back to age 35, which is [T.sub.R] = 30 years prior to retirement and [T.sub.H] = 60 years prior to the end of the family's planning horizon. One other change in Table 3 is that we have reduced the initial (current) wealth [M.sub.0] = 1 in contrast to the [M.sub.0] = 4,[W.sub.0] = 100 of the previous two tables. This particular parameter assumption is motivated by the low levels of wealth one would expect to observe among (young) 35-year-old individuals.

In this third example, the amount of life insurance purchased by the family falls in the millions of dollars, ranging from $1.03 million to $1.45 million, depending on the correlation variable [rho]. With 30 years to retirement, the family purchases almost 20 times annual wages in life insurance. These numbers are consistent with a human capital perspective that protects the family against the loss of the income source.

Another important insight from Table 3 is the (abnormally) high amounts allocated to the risky investment asset [[alpha].sub.t]. Regardless of the level of risk aversion, whether it is a high [gamma] or a low [gamma], the family places thousands of a percent in risky stocks. This is financed by borrowing thousands of a percent, a.k.a. shorting the risky free bond in our model. Clearly, these very high numbers are unrealistic (in practice) and at the very least are unfeasible, given the various institutional restrictions on borrowing with very little collateral. However, we do not believe this (odd) result represents a flaw or problem with our model, since this tends to plague most portfolio choice and asset allocation models in which risk aversion is low, borrowing is unlimited and time horizons are long.

CONCLUSION

This article focused on a subset of portfolio choice problems, where the emphasis was placed on the demand for life insurance, as a function of labor income. Our primary research question was to investigate the interaction between the demand for life insurance, the demand for risky assets, and the optimal level of consumption, as a function of one's human capital, all within a CRRA framework.

From a technical point of view, we show that the underlying HJB equation can be simplified using a similarity reduction technique, which then allows for the implementation of an efficient numerical method. Furthermore, for realistic financial economic parameter values, a closed-form approximation is derived, which greatly simplifies the calculations. Various tests confirm the robustness and accuracy of the analytic approximation compared with the full numerical solutions. Ongoing research by Huang and Milevsky (in press) is extending the application of this similarity reduction technique to more general HARA utility functions and a broader universe of mortality-contingent claims, such as retirement pensions and longevity insurance.

A variety of financial case studies illustrating our algorithm are also provided within the article. Our main qualitative and practical result is that households, whose primary breadwinner's wages are negatively correlated with financial market returns should optimally purchase more life insurance and can afford to take more risky positions with their financial portfolio. However, we also find that the optimal face value of life insurance is insensitive to the family's risk aversion.

From a practical perspective, our model validates a number of rules of thumb, used by financial planners when making portfolio and investment recommendations for their wealth management clients. First, we find that younger investors, i.e., people who are farther away from the date at which their employment wage process will hit zero, are more able to tolerate financial risk and will therefore invest more in the risky asset compared with the risk-free bond. In fact, we find a very high (optimal) leverage ratio that can range in the hundreds of a percent, for families whose risk aversion is low. Second, we find that households should insure against the loss of their primary breadwinner's wages relatively independent of how risk tolerant or risk averse the family defines itself. This is consistent with the HLV concept, originally introduced by Solomon Huebner in the early 1900s. According to this approach to financial planning and insurance, there is a universal multiple of wages that a family should insure against losing, regardless of whether this particular family wants to gamble (i.e., is risk averse.).

Finally, we find that the demand for life insurance should depend on the riskiness of the human capital (wage process) as well as its correlation with financial returns. High anticipated levels of wage volatility and/or high levels of correlation induce a reduction in the demand for life insurance. In other words, the discounted HLV alone is not enough to determine the optimal amount of life insurance. The HLV alone might over/underestimate the demand for life insurance when the remainder of the lifecycle portfolio is taken into account. Indeed, it is the economic value of human capital that is the quantity that must be insured.

APPENDIX A: DERIVATION OF THE HJB EQUATION

We assume in the time period h, the survival probability of the individual is h [p.sub.y+t], the death probability is h [q.sub.y+t] = 1 - h [p.sub.y+t], the current age of the individual is y, and the current time is t.

Let h be an arbitrarily small time increment; we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A1)

where [PHI](*, *) is the optimal objective function when there is no wage income and no need for paying insurance premium. By Ito's formula, we have

d J([M.sub.t], [W.sub.t], t) = ([J.sub.t] + [A.sup.J.sub.t])dt + [[sigma].sub.m] [M.sub.t] J [M.sup.d] [B.sup.m.sub.t] + [[sigma].sub.W] [M.sub.t] J [W.sup.d] [B.sub.t.sup.W], (A2)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A3)

Integrating Equation (A2), we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A4)

Similarly, for [PHI], we have

d[PHI]([M.sub.t], t) = ([[PHI].sub.t] + [[B.sup.[PHI].sub.t]]) dt + [[sigma].sub.m][[alpha].sub.t][M.sub.t][PHI][M.sub.d][B.sup.M.sub.t], (A5)

where

[[B.sup.[PHI].sub.t] (([micro].sub.m] - r)[M.sub.t][[alpha].sub.t] + r [M.sub.t] - [c.sub.t]) [J.sub.M] + 1/2 [[sigma].sup.2.sub.m] [M.sup.2.sub.t] J MM.

Integrating Equation (A4) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A6)

Combining (A1), (A5), and (A3), we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A7)

Dividing by h and rearranging Equation (A6), as h [right arrow] 0, we obtain

J([M.sub.t], [W.sub.t], t) [[lambda].sub.y+t] [greater than or equal to] [J.sub.t] + [A.sup.J.sub.t] + [PHI] ([M.sub.t] + [I.sub.t]/[[lambda].sub.y+t], t) [[lambda].sub.y+t] + [e.sup.-[[delta]t u([c.sub.t]). (A8)

The equality holds when we take the optimal control:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A9)

Using Equations (A3) and (A7), after rearrangement, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Here [PHI] satisfies the following HJB (without wage and insurance):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A11)

In order to simplify the presentation, we now drop the subscript t from Mt and Wt.

APPENDIX B: DERIVATION OF CLOSED-FORM SOLUTION

Motivated by the closed-form solution for constant wage Merton (1971), we seek the solution of the HJB Equation (16) in the form

J(M, W, t) = h(t) = h(t)[(M + k(t))W).sup.1-[gamma]]/1 -[gamma]. (B1)

Note that when W = 0 (no wage income), this solution reduces to the Merton formula (12). Simple calculation and algebraic manipulation lead to the following expressions

[I.sup.*.sub.t]/[lambda]M = kx, (B2)

[[alpha].sup.*.sub.t] = [[mu].sub.m] - r/[gamma][[sigma].sup.2.sub.m] (1 + kx)--[[beta]kx, (B3)

[[c.sup.*.sub.t]/M = [e.sup.[delta]/[gamma]t][h.sup. -1/[gamma] (1 + kx), (B4)

where [beta] = [rho][sigma] w/[[sigma].sub.m] and x = W/M. From (16), we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (B5)

When [rho] = [+ or -] 1, the last term in the equation above drops out, and when 1 + kx [not equal to] 0, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (B6)

This can be viewed as a linear algebraic equation for x and for the solution to exist, we must have the coefficients vanish for all x, which leads to the following two first-order ordinary differential equations for h and k:

[1/1 - [gamma]] [h'/h] + [gamma]/1-[gamma] [e.sup.[delta]/[gamma]t] [h.sup.-1/[gamma]] + r + [([[mu].sub.m] - r).sup.2]/2[gamma][[sigma].sup.2.sub.m] = 0 (B7)

k'/k + 1/k + [[mu].sub.w] - r - [[lambda].sub.y + t] - [beta] ([[mu].sub.m] - r) = 0. (B8)

Note that (B7) is the same as (13), and h is simply the Merton solution given by (14). Furthermore, our solution reduces to (12) at t = [T.sub.R], which yields the terminal conditions for k,

k([T.sub.R]) = 0. (B9)

Remark: When (B5) is viewed as a quadratic equation of x, the solution (B1) exists when this quadratic equation is of a special form (with a common factor 1 + kx) since we only have two degrees of freedom (h and k). Another special case when the quadratic equation has a common factor 1 + kx is when [[sigma].sub.w] = 0, in which case we can also eliminate a common factor 1 + kx. Merton's solution of constant wage (Merton, 1971) belongs to this case when [[mu].sub.w] is also zero.

The Case of -1 < [rho] < 1

We assume that there exists an asymptotic expansion of the value function in the form J = [J.sub.0] + [epsilon] [J.sub.1], where [epsilon] = (1 - [[rho].sup.2] [([[gamma].sup.[sigma] [m.sup.[sigma]w/[[mu].sup.m] - r]).sup.2]

[J.sub.0] (M, W, t) = [M.sup.1-[gamma]]/1 - [gamma] h(t)[(1 + k(t)x).sup.1-[gamma]], [J.sub.1] (M, W, t) = [M.sup.1-[gamma]]/1 - [gamma]) [F.sub.1](x,t).

In the following, we show that h and k are the same functions obtained in the "Special Case: [rho] = [+ or -] 1" section valid for -1 [less than or equal to] [rho] [less than or equal to] 1 as long as is small.

We proceed as before, by substituting J into the HJB Equation (16) and collect the terms of the same order of [epsilon]. At the zero order, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (B10)

which is (B6) multiplied by 1 + kx. At the first order of [epsilon], we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B11)

where [[??].sub.0], [[??].sub.1] and [[??].sub.2] are similarly defined as [c.sub.0], [c.sub.1], and [c.sub.0]in the "General Case: -1 < 1" section, with F replaced by h[(1 + kx).sup.1-[gamma]] G is a function of [F.sub.1], as well as h and k.

For practically relevant parameter values, as shown by the numerical examples, the zero-order approximation gives sufficiently accurate results. In principle, if we are interested in a more accurate approximation, we could always use the numerical method outlined earlier.

REFERENCES

Babbel, D. F., and E. Ohtsuka, 1989, Aspects of Optimal Multi Period Life Insurance, Journal of Risk and Insurance, 56: 460-481.

Bodie, Z., J. B. Detemple, S. Otruba, and S. Walter, 2004, Optimal Consumption-Portfolio Choices and Retirement Planning, Journal of Economic Dynamics and Control, 28: 1013-1226.

Bodie, Z., R. C. Merton, and W. Samuelson, 1992, Labor Supply Flexibility and Portfolio Choice in a Life Cycle Model, Journal of Economic Dynamics and Control, 16: 427-49.

Buser, S. A., and M. L. Smith, 1983, Life Insurance in a Portfolio Context, Insurance: Mathematics and Economics, 2: 147-157.

Cairns, A., D. Blake, and K. Dowd, 2006, Stochastic Lifestyling: Optimal Dynamic Asset Allocation for Defined Contribution Pension Plans, Journal of Economic Dynamics and Control., 30: 843-877.

Campbell, J. Y., and J. F. Cocco, 2003, Household Risk Management and Optimal Mortgage Choice, Quarterly Journal of Economics, 118: 1149-1194.

Campbell, J. Y., and L. Viciera, 2002, Strategic Asset Allocation: Portfolio Choice for Long Term Investors (Oxford, UK: Oxford University Press).

Campbell, R. A., 1980, The Demand for Life Insurance: An Application of the Economics of Uncertainty, The Journal of Finance, 35: 1155-1172.

Chen, P., R. Ibbotson, M. A. Milevsky, and K. Zhu, 2005, Human Capital, Asset Allocation and Life Insurance, Financial Analysts Journal, 62(1): 97-109.

Cocco, J. F., 2005, Portfolio Choice in the Presence of Housing, Review of Financial Studies, 18: 535-567.

Cocco, J. F., F. J. Gomes, and P. J. Maenhout, 2005, Consumption and Portfolio Choice Over the Life cycle, Review of Financial Studies, 18: 491-533.

Dammon, R. M., C. S. Spatt, and H. H. Zhang, 2001, Optimal Consumption and Investment With Capital Gains Taxes, Review of Financial Studies, 14: 583-616.

Davis, S., and P. Willen, 2000a, Using Financial Assets to Hedge Labor Income Risks: Estimating the Benefits, Working Paper, University of Chicago.

Davis, S., and P. Willen, 2000b, Occupational-level Income Shocks and Asset Returns: Their Covariance and Implications for Portfolio Choice, Working Paper, University of Chicago.

Detemple, J., R. Garcia, and M. Rindisbacher, 2003, A Monte Carlo Method for Optimal Portfolios, Journal of Finance, 58: 401-446.

Dybvig, P. H., and H. Liu, 2004, Lifetime Consumption and Investment: Retirement and Constrained Borrowing, Working Paper, Washington University.

Faig, M., and P. Shum, 2002, Portfolio Choice in the Presence of Personal Illiquid Projects, Journal of Finance, 57: 303-328.

Fischer, S., 1973, A Life Cycle Model of Life Insurance Purchases, International Economic Review, 14: 132-152.

Frees, E. W., J. Carriere, and E. Valdez, 1996, Annuity Valuation With Dependent Mortality, Journal of Risk and Insurance, 63: 229-261.

Goetzmann, W. N., 1993, The Single Family Home in the Investment Portfolio, Journal of Real Estate Finance and Economics, 6: 201-222.

Hakansson, N. H., 1969, Optimal Investment and Consumption Strategies Under Risk, an Uncertain Lifetime, and Insurance, International Economic Review, 10: 443-466.

Heaton, J., and D. J. Lucas, 2000, Portfolio Choice and Asset Prices: The Importance of Entrepreneurial Risk, Journal of Finance, 55: 1163-1198.

Hurd, M., 1989, Mortality Risks and Bequests, Econometrica, 57: 779-813.

Huang, H., and M. A. Milevsky, In press, Portfolio Choice and Mortality-Contingent Claims: The General HARA Case, Journal of Banking and Finance.

Jagannathan, R., and N. R. Kocherlakota, 1996, Why Should Older People Invest Less in Stocks Than Younger People? Federal Reserve Bank of Minneapolis Quarterly Review, 20: 11-23.

Kahl, M., J. Liu, and F. Longstaff, 2003, Paper Millionaires: How Valuable Is Stock to a Stockholder Who Is Restricted From Selling It? Journal of Financial Economics, 67: 385-410.

Kim, T. S., and E. Omberg, 1996, Dynamic Non-myopic Portfolio behavior, Review of Financial Studies, 9: 141-161.

Koo, H. K., 1998, Consumption and Portfolio Selection With Labor Income: A Continuous Time Approach, Mathematical Finance, 8: 49-65.

Lewis, F. S., 1989, Dependents and the Demand for Life Insurance, American Economic Review, 79: 452-467.

MaCurdy, T. E., 1982, The Use of Time-Series Processes to Model the Error Structure of Earnings in a Longitudinal Data Analysis, Journal of Econometrics, 18: 83-114.

Merton, R. C., 1969, Lifetime Portfolio Selection Under Uncertainty: The Continuous Time Case, Review of Economic Studies, 51: 247-257.

Merton, R. C., 1971, Optimum Consumption and Portfolio Rules in a Continuous Time Model, Journal of Economic Theory, 3: 373-413.

Ockendon, J. R., S. D. Howison, A. Lacey, and A. Movchan, 1999, Applied Differential Equations. (Oxford, UK: Oxford University Press).

Purcal, S., 2003, A Stochastic Control Model for Individual Asset Liability Management, Working Paper, University of New South Wales.

Richard, S. F., 1975, Optimal Consumption, Portfolio and Life Insurance Rules for an Uncertain Lived Individual in a Continuous Time Model, Journal of Financial Economics, 2: 187-203.

Sundaresan, S., and F. Zapatero, 1997, Valuation, Optimal Asset Allocation and Retirement Incentives of Pension Plans, Review of Financial Studies, 10: 631-660.

Viceira, L., 2001, Optimal Portfolio Choice for Long-Horizon Investors With Nontradable Labor Income, Journal of Finance, 55: 1163-1198.

Wang, J., 2006, Numerical PDE Techniques for Personal Finance and Insurance Problems, Ph.D. Dissertation, York University, Toronto.

Wilmott, P., J. N. Dewynne, and S. D. Howison, 1993, Option Pricing: Mathematical Models and Computation (Oxford, UK: Oxford Financial Press).

Yaari, M. E., 1965, Uncertain Lifetime, Life Insurance and the Theory of the Consumer, Review of Economic Studies, 32: 137-150.

Yao, R., and H. H. Zhang, 2005, Optimal Consumption and Portfolio Choices With Risky Housing and Borrowing Constraints, Review of Financial Studies, 18: 197-239.

Zietz, E. N., 2003, An Examination of the Demand for Life Insurance, Risk Management and Insurance Review, 6: 159-191.

(1) And although we refer to this type of protection as life insurance, our framework could also apply to (permanent) disability and other income replacement insurance. We would like to thank one of the JRI reviewers for pointing out the applicability of our model beyond traditional life insurance.

(2) For example, if the IFM curve at age y = 35 is [[lambda].sub.35] = 0,001, and the family spends I = $50 on life insurance, this entitles the family to a death benefit of 50/0.001 = $50,000.

(3) We have experimented with different solvers and the solutions are almost identical.

(4) We use 5 percent for the standard deviation of wages even though there is a lack of consensus in the literature regarding the precise magnitude of this variable. For example, MaCurdy (1982) estimates the error structure for wages and earnings of prime-age males, using data from the Michigan Panel of Income Dynamics, and estimates numbers much higher than 5 percent.

Huaxiong Huang is a professor in the Department of Mathematics and Statistics. Moshe A. Milevsky is the executive director of the IFID Center and Associate Professor in the Schulich School of Business. Jin Wang was a research associate at the IFID Centre, York University, 4700 Keele St., Toronto, Canada, when this article was written. M. Milevsky can be reached at milevsky@yorku.ca. The authors would like to acknowledge helpful comments from the JRI editor and two anonymous referees.

TABLE 1 Life Insurance (in Thousands) and Allocation to Stocks 20 Years Away From Retirement, at Age 45 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = -1.0 -0.5 [gamma] = 1.5 [[alpha].sup.*.sub.t] 6.06 5.18 Low [I.sup.*.sub.t]/[[lambda].sub.y+t] $964.7 $907.7 [gamma] = 3 [[alpha].sup.*.sub.t] 3.63 2.88 Medium [I.sup.*.sub.t]/[[lambda].sub.y+t] $964.7 $907.7 [gamma] = 4.5 [[alpha].sup.*.sub.t] 2.82 2.11 High [I.sup.*.sub.t]/[[lambda].sub.y+t] $964.7 $907.7 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = [rho] = 0.0 +0.5 +1.0 [gamma] = 1.5 4.40 3.69 3.06 Low $855.2 $806.8 $762.1 [gamma] = 3 2.20 1.59 1.05 Medium $855.2 $806.8 $762.1 [gamma] = 4.5 1.47 0.894 0.384 High $855.2 $806.8 $762.1 Note: Parameter assumptions: [[mu].sub.m] = 7%, [[sigma].sub.m] = 20%, r = 2%, [[micro].sub.w] = 1%, [[sigma].sub.w] = 5%, [M.sub.O] = 200, [W.sub.O] = 50, [T.sub.R] = 20, [T.sub.H] = 50, [delta] = 2%, y = 45, m = 86.3, b = 9.5. TABLE 2 Life Insurance (in Thousands) and Allocation to Stocks 10 Years Away From Retirement, at Age 55 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = -1.0 -0.5 [gamma] = 1.5 [[alpha].sup.*.sub.t] 3.46 3.09 Low [I.sup.*.sub.t]/[[lambda].sub.y+t] $484.8 $470.2 [gamma] = 3 [[alpha].sup.*.sub.t] 2.03 1.69 Medium [I.sup.*.sub.t]/[[lambda].sub.y+t] $484.8 $470.2 [gamma] = 4.5 [[alpha].sup.*.sub.t] 1.56 1.22 High [I.sup.*.sub.t]/[[lambda].sub.y+t] $484.8 $470.2 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = [rho] = 0.0 +0.5 +1.0 [gamma] = 1.5 2.73 2.40 2.09 Low $456.1 $442.6 $429.6 [gamma] = 3 1.37 1.06 0.775 Medium $456.1 $442.6 $429.6 [gamma] = 4.5 0.911 0.616 0.337 High $456.1 $442.6 $429.6 Note: Parameter assumptions: [[mu].sub.m] = 7%, [[sigma].sub.m] = 20%, r = 2%, [[mu].sub.w] = 1%, [[sigma].sub.w] = 5%, [M.sub.O] = 200, [W.sub.O] = 50, [T.sub.R] = 10, [T.sub.H] = 40, [delta] = 2%, y = 55, m = 86.3, b = 9.5. TABLE 3 Life Insurance (in Thousands) and Allocation to Stocks 30 Years Away From Retirement, at Age 35 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = -1.0 -0.5 [gamma] = 1.5 [[alpha].sup.*.sub.t] 1,570 1,270 Low [I.sup.*.sub.t]/[[lambda].sub.y+t] $1,450 $1,320 [gamma] = 3 [[alpha].sup.*.sub.t] 966 717 Medium [I.sup.*.sub.t]/[[lambda].sub.y+t] $1,450 $1,320 [gamma] = 4.5 [[alpha].sup.*.sub.t] 765 533 High [I.sup.*.sub.t]/[[lambda].sub.y+t] $1,450 $1,320 Correlation Between Wages and Stock Market Risk Aversion [rho] = [rho] = [rho] = 0.0 +0.5 +1.0 [gamma] = 1.5 1,010 789 599 Low $1,210 $1,110 $1,030 [gamma] = 3 505 325 171 Medium $1,210 $1,110 $1,030 [gamma] = 4.5 337 170 28.8 High $1,210 $1,110 $1,030 Note: Parameter assumptions: [[mu].sub.m] = 7%, [[sigma].sub.m] = 20%, r = 2%, [[mu].sub.w] = 1%, [[sigma].sub.w] = 5%, [M.sub.O] = 200, [W.sub.O] = 50, [T.sub.R] = 30, [T.sub.H] = 60, [delta] = 2%, y = 35, m = 86.3, b = 9.5.…