The full text on this page is automatically extracted from the file linked above and may contain errors and inconsistencies.
Why Do Couples and Singles Save During Retirement? WP 21-09 Mariacristina De Nardi University of Minnesota, Federal Reserve Bank of Minneapolis, CEPR, and NBER Eric French University of Cambridge, CEPR, and IFS John Bailey Jones Federal Reserve Bank of Richmond Rory McGee University of Western Ontario and IFS Why Do Couples and Singles Save During Retirement? Mariacristina De Nardi, Eric French, John Bailey Jones and Rory McGee* May 27, 2021 Abstract While the savings of retired singles tend to fall with age, those of retired couples tend to rise. We estimate a rich model of retired singles and couples with bequest motives and uncertain longevity and medical expenses. Our estimates imply that while medical expenses are an important driver of the savings of middle-income singles, bequest motives matter for couples and highincome singles and generate transfers to nonspousal heirs whenever a household member dies. The interaction of medical expenses and bequest motives is a crucial determinant of savings for all retirees. Hence, to understand savings, it is important to model household structure, medical expenses, and bequest motives. * De Nardi: University of Minnesota, Federal Reserve Bank of Minneapolis, CEPR, and NBER. French: University of Cambridge, CEPR, and IFS. Jones: Federal Reserve Bank of Richmond. McGee: UWO and IFS. We thank Wilbert van der Klaauw, Monica Costa Dias, Tom Crossley, Ami Ko, Melanie Lührmann, Tatyana Koreshkova, Victor Rı́os-Rull and many seminar participants for helpful comments, and Karolos Arapakis, Justin Kirschner, Jeremy McCauley, Ali Merali and Rachel Rodgers for research assistance. We gratefully acknowledge support from NSF grant SES-2044748, the Social Security Administration through the Michigan Retirement Research Center (MRRC grant UM17-15), the Economic and Social Research Council (Centre for Microeconomic Analysis of Public Policy at the Institute for Fiscal Studies (RES-544-28-50001), and Grant Inequality and the insurance value of transfers across the life cycle (ES/P001831/1)). The views expressed in this paper are those of the authors and not necessarily those of the National Bureau of Economic Research, the Federal Reserve Bank of Minneapolis, the Federal Reserve Bank of Richmond, the IFS, or the CEPR. 1 1 Introduction The savings of couples and singles behave very differently over the retirement period. While the savings of singles tend to fall with age, those of retired couples tend to increase or stay constant until one of the spouses dies. Why is this the case? Do couples face different medical spending and longevity risks? How important are risks to the surviving spouse? Do couples have different attitudes toward bequests? Answering these questions helps us understand the savings of all retirees. Understanding why households save, in turn, is essential to predict how households would react to policy reforms and other changes. Previous work has mainly focused on singles. These studies find that medical expenses are an important reason why single retirees hold wealth, but that the desire to leave bequests might also be important. Much less is known about the reasons why retired couples save. Couples differ from singles in several notable ways. Being in a couple allows its members to pool their longevity and medical expense risks and potentially reduces medical spending due to spousal caregiving, but it also exposes each member to their spouse’s risks, including the income loss and high medical expenses that often accompany a spouse’s death. A couple also cares about leaving resources to the surviving spouse and potentially about leaving bequests to other heirs. Disentangling these motivations is important because couples account for about 50% of people age 70 and older and tend to be richer. We start our analysis by documenting couples’ and singles’ savings, medical expenditures, and longevity. We then construct a rich model of retiree saving that incorporates heterogeneity in life expectancy and medical expenses and includes the large jump in the medical expenses that is typically observed before a death. We model means-tested social insurance and bequest motives, both upon death of the last survivor and upon the death of one of the people in a couple, in which case net worth is optimally split between the surviving spouse and other heirs. We estimate our model using the Method of Simulated Moments (MSM) and Assets and Health Dynamics of the Oldest Old (AHEAD) data. Finally, we use our estimated model to evaluate the extent to which retirees’ savings are driven by medical expenses, bequest motives, and their interaction. Our paper provides multiple contributions. We establish new facts on the savings of singles and couples during retirement: While the savings of retired singles stay roughly constant or fall during retirement, those of couples increase as long as both spouses are alive. Net worth drops sharply (by $160,000 on average) when the first spouse dies. A large share of this drop is explained by high medical expenses before 2 death, but not the majority, which is accounted for by transfers to nonspousal heirs. By the time the second spouse dies, much of the couples’ wealth has vanished. We use our estimated model to quantify the effects of various saving motives during retirement for couples and singles. Our results for singles reinforce earlier findings: Most singles save mainly to self-insure against future medical expenses, but at the top of the permanent income (PI) distribution, bequests are also important, especially when interacting with medical expense risk. Our results for couples reveal several new findings. First, bequest motives have large effects on the savings of older couples in the top two PI terciles. That is, at age 84, couples in the highest PI tercile would hold 26% less wealth (that is, $325,000 instead of $437,000) if they did not have bequest motives. Couples in the middle PI tercile would hold about 14% less ($172,000 instead of almost $200,000). In contrast, couples at the bottom, just like singles at the bottom and median PI level, would barely change their wealth. Thus, more couples than singles save to leave bequests, and bequest motives affect their savings by a much larger extent. This is for at least two reasons. First, couples are wealthier, and we estimate bequests to be a luxury good. Second, couples not only receive direct utility from leaving bequests at the time of one spouse’s death, but also indirect utility from the expectation that the last surviving spouse will leave bequests at his or her death. Second, in the presence of bequest motives, medical expenses have a smaller effect on the savings of couples at all PI levels. For instance, at age 84, the median wealth of couples in the top PI tercile would fall by less than 10% (from $437,000 to $417,000) if they did not face any medical expenses. This suggests that medical expenses are only one of the reasons why they continue to save at very advanced ages. Third, while the savings of couples are driven more by bequests than by medical expenses, the interaction of medical expenses (including for the surviving spouse) and bequest motives has large effects. At age 84, the median wealth of couples in the top PI tercile would fall by more than 65% (from $437,000 to $150,000) if they did not face any medical expenses and had no bequest motives. Because wealth saved for bequests also helps insure against medical spending, the incremental effects of bequest motives or medical spending, even when summed, can be much smaller than their joint effect. This is due to the way bequest and precautionary saving motives interact. When medical expenses are uncertain, households that reduce their consumption to finance future medical expenses may die with unspent wealth. By giving these unspent funds value, bequest motives reduce the opportunity cost of self-insuring through saving. Moreover, medical expenses introduce uncertainty over whether a household will actually leave a bequest. This might generate precautionary saving aimed at increasing the likelihood that a bequest is 3 made. A third way in which these motives interact is that medical expenses have a negative wealth effect, which tends to decrease both savings and bequests. Fourth, couples save to provide resources to the surviving spouse after the death of the first spouse. This is important because there are considerable differences in the life expectancy of each spouse. At age 70, married women and men expect to live 15.8 years and 11.5 years, respectively. In addition, while the average life expectancy of the couple is 13.7 years, the oldest survivor expects to live for 17.9 years. During their remaining years alone, survivors face higher medical spending (including from an elevated risk of entering a nursing home) and receive lower income. As a result, couples wish to leave resources to the surviving spouse upon a partner’s death. This has a large effect on savings at all PI levels. If couples instead placed no weight on the survivor, at age 84 their median wealth in the middle PI tercile would fall by more than 17% (from $199,000 to $165,000). Studies omitting this risk will understate precautionary saving motives. Fifth, when one member of a high-PI couple dies, even though most of the estate is usually left to the surviving spouse, large bequests to nonspousal heirs are often distributed. Changes to saving motives thus affect both the saving of couples and singles and the division of wealth between the surviving spouse and nonspousal heirs. For example, when bequest motives are removed, couples save less, but when one spouse dies all of the estate is left to the survivor. In fact, we find that removing bequest motives often increases the amount of wealth held by widows and widowers. Taken together, our results show that bequest motives play a larger role for the savings of couples than those of singles and for the savings of high-PI households than those of low-PI households, and that the interaction of precautionary savings and bequests are most important for couples and high-PI households. In contrast, providing for the surviving spouse is important for all couples. We reach similar conclusions when we look at all households, couples and singles, together. For median wealth, and at lower percentiles, the effects of removing medical spending exceed the effects of eliminating bequest motives. At higher asset percentiles, and for the mean, bequest motives have larger effects. The interaction of bequest and precautionary saving motives is large across the entire asset distribution. The rest of the paper is organized as follows. Section 2 summarizes the findings of key related papers and motivates our main modeling choices. Section 3 briefly introduces our data and displays the key facts that we seek to understand. Section 4 introduces our model, and section 5 details our estimation procedure. Section 6 displays important features of the data affecting households’ saving decisions and evaluates our model fit. To build additional trust in the predictions of our estimated model, Section 7 evaluates the performance of our model in terms of important 4 additional aspects of the data that we do match at the estimation stage. Section 8 evaluates the impact of medical expenses and bequest motives for singles and couples, and Section 9 concludes. 2 Related literature and modeling choices The goal of this paper is to document and explain the differences in the saving behavior of couples and singles and their contribution to aggregate retirement savings. Although much of the structural work on retiree saving has focused on singles,1 a few papers, including Braun, Kopecky and Koreshkova [9] and Nakajima and Telyukova [43], include couples. Our paper builds on these contributions by modeling what happens when a spouse dies. As in Poterba, Venti, and Wise [49], we find that net worth falls when households lose a spouse and that this fall is too large to be attributable to medical and other end-of-life expenses. A novel aspect of our framework is that we allow households to leave bequests to nonspousal heirs when only one of the spouses dies. This feature proves important in matching the fall in wealth observed at the time of the first spouse’s death. Bequests to heirs also occur when the final member of a household dies. The desire to leave such bequests has received considerable attention as a potential explanation why households retain high levels of wealth at very old ages, as in Ameriks, Caplin, Laufer, and Van Nieuwerburg [4], Ameriks, Briggs, Caplin, Shapiro, and Tonetti [3], Dynan, Skinner, and Zeldes [23], and Lockwood [41]. In addition, De Nardi, French, and Jones [19] show that this kind of bequest motive helps fit the asset and Medicaid recipiency profiles of singles. Besides modeling household demographic transitions and bequest motives toward the surviving spouse and other heirs, we take into account health and medical expense risk, social insurance for couples and singles, and heterogeneity in life expectancy by PI and marital status. These features are important because households face potentially large out-of-pocket medical and nursing home expenses (Braun, Kopecky, and Koreshkova [9], French and Jones [27, 28], Palumbo [47], Feenberg and Skinner [25], and Marshall, McGarry, and Skinner [42]), which generate precautionary savings (as in Kopecky and Koreshkova [38], and Laitner, Silverman, and Stolyarov [39]). But these risks are partially insured by means-tested programs such as Medicaid and Supplemental Security Income, which provide strong saving disincentives (Hubbard, Skinner, and Zeldes [32], and De Nardi, French, and Jones [18]). Finally, previous 1 De Nardi, French, and Jones [20] provide a detailed literature review. 5 work has shown that high-income individuals live longer than low-income individuals (Attanasio and Emmerson [6], and Rı́os-Rull and Pijoan Mas) and that married people live longer than single people. 3 3.1 Savings, medical spending, and Medicaid The AHEAD dataset We use the AHEAD dataset, which begins with households of non-institutionalized individuals age 70 or older in 1993/94 and subsequently surveys their survivors every two years. Appendix A describes the details of our sample selection and data work. Three important decisions are that (i) to abstract from labor supply decisions, we only consider retired households; (ii) to be consistent with the transitions in the model we drop household who either get married or divorced during the sample period; and (iii) due to the well-known under-reporting of assets and medical spending in the initial wave, we drop the 1993/94 data. Our 1996 data include 4,634 households, of whom 1,388 are couples and 3,246 are singles. This represents 24,274 household-year observations for which at least one household member was alive. We distinguish three types of household: single man, single woman, and couple. Single people include those who are divorced, never married, or widowed when first observed, and those who become widowed over the sample period. 3.2 Permanent Income At retirement time, permanent income (PI) captures substantial household exante heterogeneity: Households with different PI ranks receive different flows of retirement income and face different processes for health, mortality and medical expenses. A contribution of this paper is distinguishing heterogeneity from risks. To estimate a household’s PI, we first sum all of its annuitized income sources (Social Security benefits, defined benefit pension benefits, veterans benefits and annuities).2 Because there is a roughly monotonic relationship between lifetime earnings and our annuity income measure, this measure is also a good indicator of its income during the working period. We then construct a PI measure comparable across households of different ages and sizes. To do so, we regress annuity income on a household fixed effect, dummies for household structure, a polynomial in age, 2 Since we model means-tested social insurance from SSI and Medicaid explicitly in our model (through a consumption floor), we do not include SSI transfers. 6 and interactions between these variables. The rank order of each household’s estimated fixed effect provides our measure of its PI. This is a time-invariant measure that follows the household even after one of its members dies. See Appendix B for details. 3.3 Savings, medical expenses, and bequest Our measure of wealth (or net worth, or assets) is the sum of all assets less mortgages and other debts.3 Figure 1 displays median net worth, conditional on birth cohort and PI tercile, for different configurations of couples and singles. We break the data into 4 cohorts containing people who in 1996 were ages 72-77, 7883, 84-89, and 90-102. For clarity we display data for two cohorts in this picture (Appendix I displays all of our data, as well as model fit). To construct these profiles, we calculate medians for each cohort for those alive in each calendar year (an unbalanced panel). The line for each cohort starts at the cohort’s average age in 1996. Panel (a) displays net worth for current singles, that is all people who are single at a given age and in a given cohort. It highlights that the savings of elderly singles depend on their income (which is predetermined at retirement). Individuals with the lowest lifetime incomes reach retirement with little wealth and run them down as they age. After that point, they rely on annuitized income and government insurance to fund their consumption and medical expenses. Elderly singles in the middle and top of the income distribution also run down their wealth as they age, but very slowly, and carry some wealth into very old ages. Panel (b) shows the net worth profiles that would arise if we were to limit the sample to only people who were initially single, without including those who become single upon being widowed. Because the newly single tend to have more wealth, the trajectories in Panel (a) slope down less than those in Panel (b); the two sets of trajectories are nonetheless quite similar. Panel (c) of Figure 1 plots the net worth of current couples, including couples as long as they remain intact but dropping them when they become single. Panel (c) reveals several interesting results. First, relative to singles in the same PI tercile, couples in the lowest PI tercile reach retirement with more wealth. Second, the asset profiles of the lower- and middle-income surviving couples display no decumulation. Third, couples in the top PI tercile increase their net worth until almost age 90. 3 The AHEAD has information on the value of housing and real estate, autos, liquid assets (which include money market accounts, savings accounts, Treasury bills, etc.), individual retirement accounts, Keoghs, stocks, the value of a farm or business, mutual funds, bonds, and other assets. 7 Median Net Worth by Cohort and Income: Data, Initial Singles 270 240 240 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) Median Net Worth by Cohort and Income: Data, Singles 270 210 180 150 120 90 60 30 0 72 210 180 150 120 90 60 30 76 80 84 88 92 96 0 72 100 76 80 84 Age 96 (a) (b) Median Net Worth by Cohort and Income: Data, Initial Couples 100 400 550 350 500 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 92 Median Net Worth by Cohort and Income: Data, Couples 600 450 400 350 300 250 200 150 100 300 250 200 150 100 50 50 0 72 88 Age 76 80 84 88 92 96 0 72 100 Age 76 80 84 88 92 96 100 Age (c) (d) Figure 1: Median wealth, AHEAD data. Panel (a): current singles. Panel (b): initial singles. Panel (c): current couples. Panel (d): initial couples. Each line shows medians for a cohort-PI cell over the period 1996-2014, plotted against average cohort ages. Thicker lines denote higher PI groups. 8 Panel (d) of Figure 1 reports median wealth for the population of those who are initially in a couple but might have subsequently lost a spouse, in which case we report the wealth of the surviving spouse. Relative to intact couples, initial couples in the higher PI terciles decumulate their wealth more quickly. Some of this decumulation may reflect the higher decumulation rates of singles, but some is also due to the loss of wealth that occurs when one of the spouses dies. To formalize and quantify this observation, we study the difference in net worth between two sets of couples: those in which one spouse dies and those that are similar in multiple dimensions 6 years prior to death, except for the fact that they do not experience a spousal death during a certain time period.4 This matching analysis accounts for observable initial heterogeneity across couples and thus identifies the impact of death and the declining health that precedes it. We use these data to estimate the following event study specification: ai,t = fi + 4 X gj + dj Di × 1{t − Ti = j} + ei,t (1) j=−4 where ai,t denotes the wealth of household i in calendar year t, and Di indicates whether i belongs to the group that does not lose a spouse (Di = 0) or to the group that does (Di = 1). We normalize the date of death, Ti , to occur at year 0 and follow the households for three waves before and two waves after the death. The coefficients of interest are the parameters {dj }, which show the extent to which wealth rise or fall for households who experience the death of a spouse, relative to households in the matched sample who experience only a placebo death at date Ti .5 Panel (a) of Figure 2 reports our estimated death-related asset decline (the values of dj ). On average, the wealth of households experiencing a death declines by an additional $160,000 compared to those not experiencing it. While the decline begins up to 6 years in advance of death, a substantial share of the decline occurs in the final two years. Part of this decline reflects a gradual worsening of health and a related increase in medical expenses. The rest reflects bequests made at the time of death of the first spouse and changes in consumption. The gap between the treatment and control groups continues to widen during the post-death years. This may reflect additional death-related effects, for instance the loss of a spouse’s income, or other differences between intact couples and singles. 4 We match on PI, initial wealth, and age. See Appendix C for details and Jones, De Nardi, French, McGee and Rodgers [35] for additional outcomes. 5 The coefficients gj and dj are defined relative to their value 6 years before death (the omitted indicator category). In interpreting the summation, recall that AHEAD interviews occur every other year. 9 -5 -200 1000s of 2014 Dollars -100 0 1000s of 2014 Dollars 0 5 10 15 100 Some of these asset declines are explained by high medical expenses prior to death. The AHEAD contains high-quality data on what the household spends outof-pocket on private insurance premia, drug costs, hospital stays, nursing home care, doctor visits, dental visits, and outpatient care, including those incurred during the last year of life. To construct our measure of out-of-pocket medical spending, we sum all of these components, along with funeral, other end-of-life expenses, and imputed Medicare insurance contributions. Panel (b) of Figure 2 shows that when a member of the household dies, there is a sharp increase in medical expenses: they are $8,000 larger during each of the two years preceding the death of a spouse, compared with similar couples experiencing no death. Over the 6 years preceding death, the total difference in medical expenses between these two households surpasses $27,000. Thus high end-of-life medical ex= 17% of the average fall in wealth. The effect for singles penses can explain 27,000 160,00 is even larger, totalling almost $40,000. -6 -4 -2 0 Years from Death Data 2 4 -6 95% CI -4 -2 0 Years from Death Data (a) 2 4 95% CI (b) Figure 2: Average changes around the death of a spouse, initial couples, AHEAD/MCBS data. Panel (a): Mean asset decline around the death of a spouse. Panel (b): Mean increase in the sum of out-of-pocket medical and death related expenses. Solid lines: estimates, shaded region: 95% confidence interval. Death dates are centered at year 0. If not medical and other end-of-life expenses, what explains the drop in wealth at the time of death of a spouse? When one spouse dies, a family member of the deceased, usually the spouse, is asked how the estate is split between the surviving spouse and other heirs, such as children. Comparing the wealth of the newly single with nonspousal transfers shows that most of the wealth is left to the surviving 10 spouse, but that, on average $79,000 are transferred to other heirs around the time of death of the first spouse. Transfers to other heirs at the time of death of a spouse are large and support the view that bequest motives are important. Although we do not require our model to match the decline in wealth around the death of a spouse or transfers to other heirs, our estimated structural model matches these aspects of the data very well, and thus helps validate its predictions. 3.4 Medicaid recipiency Medicaid Recipiency by Cohort and Income: Data, All 0.6 0.55 0.5 Medicaid Recipiency 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 72 76 80 84 88 92 96 100 Age Figure 3: Medicaid recipiency by age, income, and cohort. Because previous work has shown that means-tested social insurance programs affects saving, it is important to match the extent to which retirees use such programs. Figure 3 plots Medicaid recipiency rates conditional on age, PI tercile, and cohort for all households in our sample (the AHEAD has reliable measures of Medicaid recipiency, but not spending). It shows that the lowest-income retirees are most likely to end up on Medicaid and that, and consistently with the asset-tested nature of these programs (high-asset households are not eligible) and the asset profiles shown in Figure 1, Medicaid recipiency increases with age. 4 The model A retired household maximizes utility by choosing savings, current and future consumption, and, upon death of his or her spouse, the split of wealth between 11 himself or herself and other heirs. These choices also determine final bequests upon the death of the last household member. Households make these choices at household head age t, t = tr , tr + 1..., T + 1, where tr is the initial period that we consider and T is the maximum potential lifespan. For tractability, we assume that wives are always two years younger than their husbands, so that one age is sufficient to characterize a household. Consistent with the AHEAD data frequency, our time period is two years long. The household begins the period in one of the following possible marital statuses: a couple, a newly widowed man or woman, or a single man or single woman. Single men, single women, and couples use the cash on hand they have at the beginning of a period to consume and save. After that, mortality, health and medical expense shocks occur, income net of taxes is received, and government transfers take place. At this point, the amount of cash-on-hand that is carried into the next period is known, and the household enters next period. People who are newly widowed have an additional decision, which is how to divide their estate between themselves and bequests to others. Once that decision is made, they become single men or women, with wealth equal to the remainder of their estates. When the final member of a household dies, all remaining net worth goes to its heirs. Consistent with reality, our timing implies that medical costs associated with death are collected before any bequests can be made and that Medicaid pays bills that are incurred even when a patient dies. 4.1 Preferences The per-period utility functions for singles and couples are given by uS (c) = c1−ν , 1−ν uC (c) = 2 (c/η)1−ν , 1−ν 1 < η ≤ 2, ν ≥ 0, (2) respectively, where c is total consumption and the parameter η determines the extent to which couples enjoy economies of scale in the transformation of consumption goods to consumption services. The household weighs future utility with the factor β. The household derives utility θj (b) from leaving bequest b, where θ0 (b) is the utility from bequests when there are no surviving members in the household, while θ1 (b) is the utility from bequests when there is a surviving spouse. It takes the form θj (b) = φj (b + κj )1−ν , 1−ν (3) where κj determines the curvature of the bequest function, and φj determines its intensity. De Nardi [16] used this functional form to help explain savings over the 12 life cycle. It supports several interpretations of the “bequest motive:” dynastic or “warm glow” altruism (as in Becker and Tomes [7], Abel and Warshawsky [1], or Andreoni [5]); strategic motives (as in Bernheim, Schleifer and Summers [8] or Brown [10]); or some form of utility from wealth itself, as in (Carroll [12] and Hurd [33]). 4.2 4.2.1 Sources of uncertainty and budget constraints Income Because there is much less income uncertainty during retirement than during the working period, we simplify the model by assuming that the household’s non-asset income at time t, yt (·), is a deterministic function of the household’s permanent income, I, age, and family structure ft (single man, single woman, or couple). yt (·) = y(I, t, ft ). (4) We do not include received bequests as a source of income, because very few households aged 70 or older receive them. 4.2.2 Health and survival uncertainty: Health and survival are individual, rather than household-level, variables and we use gender, g ∈ {h, w}, to differentiate between men and women. A person’s health status, hsg , indicates whether he or she is in a nursing home, in bad health, or in good health. The transition probabilities for a person’s future health status depend on that person’s current health status, permanent income, age, gender and marital status πt (·) = Pr(hst+1 |hst , I, t, g, ft ). (5) Survival depends on the same variables. Let st (I, g, hst , ft ) denote the probability that an individual alive at age t survives to age t + 1. 4.2.3 Medical expense uncertainty We use mt+1 to denote the sum of medical spending that is either paid out-ofpocket by the household or covered by Medicaid between periods t and t + 1. While we treat this total as exogenous, the division of expenses between the household and Medicaid depends on the household’s financial resources, the total expenses that must be covered, and the level of the consumption floor. In other words, mt+1 gives 13 the household’s maximum possible medical spending obligation, but because of social insurance the amount faced by poorer households may be much smaller. We allow mt+1 to depend on the health status of each family member at both the beginning and end of the period, permanent income, age, household’s family structure (also differentiating single men and single women) at the beginning and end of each period, and an idiosyncratic component, ψt+1 . h w ln mt+1 = m(hsht , hsw t , hst+1 , hst+1 , I, t + 1, ft , ft+1 ) (6) h w + σ(hsht , hsw t , hst+1 , hst+1 , I, t + 1, ft , ft+1 ) × ψt+1 . We normalize the variance of ψt+1 to 1. Allowing medical expenses to depend on the household’s composition and health status at both the beginning of a period (which was realized at the very end of the previous period) and the period’s end (which will be carried over into the subsequent period) allows us to capture the jump in medical spending that occurs when a family member dies – that is ft changes – and to incorporate the impact of two subsequent health realizations on medical spending. The latter better helps us account for the cost of prolonged periods of bad health or nursing home stays. Our timing implies that the value of mt+1 is not known to the household when it decides how much to consume between periods t and t + 1. Following Feenberg and Skinner [25] and French and Jones [27], we assume that ψt+1 can be decomposed as ψt+1 = ζt+1 + ξt+1 , ξt+1 ∼ N (0, σξ2 ), (7) ζt+1 = ρm ζt + t+1 , N (0, σ2 ), (8) t+1 ∼ where ξt+1 and t+1 are serially and mutually independent. We discretize ξ and ζ, using the methods described in Tauchen [52]. 4.2.4 Budget constraints Let at denote net worth at the beginning of period t and r denote their constant pre-tax rate of return. Total post-tax income is given by Υ(r at + yt (·), τft ), with the vector τft summarizing the tax code, which depends on family structure. Define the resources available before government transfers as x et = at + Υ r at + yt (·), τft − mt , (9) To capture Medicaid and SSI, we assume that government transfers bridge the gap between a minimum consumption floor and the household’s financial resources, n o trt (e xt , ft ) = max 0, cmin (ft ) − x et , (10) 14 where we allow the guaranteed consumption level cmin to vary with family structure. To save on state variables we follow Deaton [17] and sum the household’s financial resources after government transfers into cash-on-hand: xt = at + Υ(r at + yt (·), τft ) − mt + trt (e xt , ft ). (11) Households divide their cash-on-hand between consumption and savings: at+1 = xt − ct , ct ∈ [cmin (ft ), xt ], ∀t. (12) (13) Equation (13) ensures that consumption is at least as high as the consumption floor and that savings are non-negative. Next period’s resources before transfers and next period’s cash-on-hand can then be expressed in terms of this period’s cash-on-hand and are, respectively, x et+1 = (xt − ct ) + Υ r (xt − ct ) + yt+1 (·), τft+1 − mt+1 , (14) xt+1 = xt − ct + Υ r (xt − ct ) + yt+1 (·), τft+1 − mt+1 + trt+1 (e xt+1 , ft+1 ). (15) 4.3 Recursive formulation Let ft = S indicate a single-person household. The value function for a single person of age t and gender g is g Vt (xt , hst , I, ζt ) = max uS (ct ) + βst (I, g, hst , S) ct g × Et Vt+1 (xt+1 , hst+1 , I, ζt+1 ) + β[1 − st (I, g, hst , S)]Et θ0 (xt+1 ) , (16) subject to Equations (4)-(8), and (13)-(15). A newly-single person – one who was part of a couple in the previous period and single now – distributes bequests towards other heirs before making savings and consumption decisions as a single person: n o ng g Vt (xt , hst , I, ζt ) = max θ1 (bt ) + Vt (xt − bt , hst , I, ζt ) , (17) bt subject to equation bt ∈ [0, xt − cmin (ft )], 15 (18) which prohibits the surviving spouse from using bequests to become eligible for government transfers. The value function for couples (ft = C) can be written as w h C Vt (xt ,hst , hst , I, ζt ) = max uC (ct ) ct h C h w + βst (I, w, hsw , C)s (I, h, hs , C)E V (x , hs , hs , I, ζ ) t t t+1 t t t+1 t+1 t+1 t+1 w h nw w w + βst (I, w, hst , C)[1 − st (I, h, hst , C)]Et Vt (xt+1 , hst+1 , I, ζt+1 ) h nh h h + β[1 − st (I, w, hsw , C)]s (I, h, hs , C)E V (x , hs , I, ζ ) t t t+1 t t t+1 t+1 t+1 h + β[1 − st (I, w, hsw (19) t , C)][1 − st (I, h, hst , C)]Et θ0 (xt+1 ) , subject to Equations (4)-(8) and (13)-(14). The dating of the continuation value for new widows, Vtnw (·), reflects that wives are one model period (two years) younger than their husbands. We solve our model numerically: see Appendix E for more details. 5 Estimation We adopt a two-step strategy to estimate the model. In the first step, we estimate or calibrate those parameters that, given our assumptions, can be cleanly identified outside our model. These include health transitions, out-of-pocket medical expenses, and mortality rates from raw demographic data. In addition, we fix the discount factor β at an annual value of 0.97, and we set the consumption floor for couples to be a multiple of the consumption floor for singles according to statutory rules for Medicaid and Supplemental Social Insurance. In the second step, we estimate the rest of the model’s parameters, which include risk aversion, the consumption equivalence scale, bequest parameters, and the consumption floor for singles, ∆ = ν, η, φ0 , φ1 , κ0 , κ1 , cmin (ft = S) , with the method of simulated moments (MSM), taking as given the parameters that were estimated in the first step. In particular, we find the parameter values that allow the simulated life-cycle decision profiles to “best match” (as measured by a GMM criterion function) those from the data. Appendix E provides more details on the mechanics of our MSM procedure. 16 Because our goal is to explain why retirees save so much and do so at rates that differ by income, we match moments of the distribution of net worth by cohort, age, and PI tercile. Because we wish to study differences in savings patterns of couples and singles, we match asset profiles for the singles and couples separately. Finally, because Medicaid is an important program insuring the medical expenses and consumption of the poor we also match Medicaid recipiency. More specifically, the moment conditions that comprise our estimator are given by 1. The 25th percentile, median, and 75th percentile of asset holdings by PI tercilecohort-year for current singles (including bequests). 2. The 25th percentile, median, and 75th percentile of asset holdings by PI tercilecohort-year for those who are currently couples with both members currently alive.6 3. Medicaid recipiency by PI tercile-cohort-year for all households currently alive. Appendix F contains a detailed description of these moment conditions, the weighting matrix in our GMM criterion function, and the asymptotic distribution of our parameter estimates. When there is a death in a couple, we drop the household from the current couples’ moments and add the surviving spouse to the appropriate moment condition for current singles. This is consistent with our assumption that singles differ only in their state variables, regardless of marital history. Such an assumption also holds well in the data, meaning that, conditional on our state variables the savings of recent singles are very similar to those of the people we were single in the past. (See Figure 1.) This is perhaps not surprising, given that the vast majority of single people in our sample were married at some point in the past. Our estimation targets do not include the falls in wealth that follow the first spouse’s death because we use them to validate our estimated model by applying the same estimation procedure to the AHEAD data and our model-simulated data. When estimating our model, we face two well-known problems. First, in a crosssection, older households were born in earlier years than younger households and, due to secular income growth, have lower lifetime incomes. Because of this, the asset levels of households in older cohorts will likely be lower as well. As a result, comparing older households born in earlier years to younger households born in later years leads to understated asset growth. We address this problem by starting 6 We also include bequests for the small number of couples where both spouses die between adjacent waves. 17 our simulations with initial conditions that come from the data (including wealth, income, age, and so on). This allow us to capture much of the heterogeneity across cohorts. Second, lower-income households and singles tend to die at younger ages than higher-income households and couples. The average survivor in a cohort thus has higher lifetime income, and thus more wealth, than the average deceased member of the same cohort. This “mortality bias” is more severe at older ages, when a greater share of the cohort members are dead. As a result, not accounting for mortality bias leads to overstated asset growth. We address this problem by modeling mortality as a function of key observables that include gender, permanent income and marital status. Our simulated profiles thus incorporate the same mortality bias as the data. 6 Estimation results In section 3 we report some of the key facts about savings and Medicaid recipiency that we require our estimated model to match. In this section we report some of the most relevant features of our first-step model estimates, and then discuss our second-step estimates. 6.1 First-step estimation results As they are key elements affecting saving behavior, the most important features of our first-step estimates pertain to life expectancy and nursing home risk, income and its drop when one of the spouses dies, and medical spending. Appendix G describes our first-step estimation procedures in detail. 6.1.1 Life expectancy and nursing home risk We estimate health transitions and mortality rates simultaneously by fitting the transitions observed in the AHEAD to a multinomial logit model. We allow the transition probabilities to depend on age, sex, marital status, current health status, PI, and interactions of these variables. Table 1 shows life expectancies at age 70 for singles and couples, respectively, that we obtain when we use our estimated transition probabilities to simulate demographic histories, beginning at age 70, for different gender-PI-health-family structure combinations. It shows that rich people, women, married people, and healthy people live much longer than their poor, male, single, and sick counterparts. For instance, a single man at the 10th permanent income percentile and in a nursing home expects to live only 3.0 more years, while a single woman at the 90th percentile and in good 18 Income Percentile Men Nursing Bad Home Health Good Nursing Health Home Women Bad Health Good Health All Singles 10 50 90 3.0 3.0 2.9 6.9 7.8 8.1 8.7 10.3 10.9 4.1 4.1 3.8 11.3 12.3 12.5 13.2 14.9 15.4 10.2 11.5 12.0 Couples 10 50 90 2.7 2.8 2.7 7.8 9.4 10.4 9.8 12.2 13.5 4.0 4.0 3.9 12.1 13.7 14.6 14.1 16.3 17.3 11.3 13.4 14.5 Single Men Married Men Single Women Married Women 9.0 11.5 13.9 15.8 Oldest Survivor∗ Probability that Oldest Survivor is Woman Notes: ∗ 17.9 63.7% Oldest survivor among households that were a couple at age 70. Table 1: Life expectancy in years, conditional on reaching age 70 health expects to live 15.4 more years. The far right column of the top two panels shows average life expectancy conditional on permanent income, averaging over both genders and health states. Singles at the 10th percentile of the permanent income distribution live on average 10.2 years, while singles at the 90th percentile live on average 12.0 years. People in couples live about 2 years longer than singles: single women live on average 13.9 years versus 15.8 for married women but, conditional on PI and health, the differences in longevity are much smaller. Thus, married people live longer than singles, but a significant part of the difference is explained by the fact that married people tend to have higher PI and to be in better health. The bottom part of Table 1 shows the expected years of remaining life for the oldest survivor in a household when both the man and the woman are 70. On average the last survivor lives 17.9 more 19 years. The woman is the oldest survivor 63.7% of the time. Income Percentile Men Bad Good Health Health Women Bad Good Health Health All Singles 10 50 90 23.6 22.8 20.3 25.3 24.8 22.8 35.8 35.5 32.2 37.9 38.2 35.8 32.8 32.5 30.1 Couples 10 50 90 17.3 16.6 14.6 19.2 18.8 16.8 34.4 34.1 31.4 37.0 37.3 34.5 28.7 28.7 26.3 Single Men Married Men Single Women Married Women 26.4 19.5 37.2 36.3 Table 2: Probability of ever entering a nursing home, conditional on being alive at age 70 Table 2 shows that single men and women face on average a 26% and 37% chance of being in a nursing home for an extended stay (at least 60 days in a year), respectively, while married men and women face on average a 20% and 36% chance of being in a nursing home for an extended stay. Married people are much less likely to transition into a nursing home at any age, but married people, especially women, often become single as their partner dies. Furthermore, married people tend to live longer than singles, so they have more years of life to potentially enter a nursing home. Compared to gender or marital status, permanent income and age 70 health have a smaller effect on ever being in a nursing home. This is because those with high permanent income, or in good health, are less likely to be in a nursing home at each age but they tend to live longer. 20 6.1.2 Income As noted above, we model non-asset income as a function of PI, age, and family structure. Figure 4 presents predicted income profiles for those at the 20th and 80th percentiles of the PI distribution. The average income of couples ranges from about $14,000 at the 20th percentile to over $30,000 at the 80th . As a point of comparison, median wealth holdings for these two groups at age 74 are $70,000 and $330,000, respectively. To evaluate how household income drops when one of the spouses dies, Figure 4 displays three scenarios for each PI level, all commencing from the income of a couple. Under the first scenario, the household remains a couple until age 100. Under the second one, the man dies at age 80. Under the third one, the woman dies at age 80. Our estimates imply that couples in which the husband dies at age 80 suffer a 40% decline in income, while couples in which the wife dies at 80 suffer a 30% decline in income. 20000 10000 0 2005 dollars 30000 Income, by Permanent Income Percentile, Initially Couples 70 80 90 100 age 20th percentile couple 20th percentile to female 20th percentile to male 80th percentile couple 80th percentile to female 80th percentile to male Figure 4: Income, conditional on permanent income and family structure. All households begin as couples, then either stay in a couple or switch to being to a single man or single woman at age 80 The income losses that occur at the death of a spouse reflect the fact that even though Social Security and private defined benefit pensions have survivors’ benefits, these benefits replace only a fraction of the deceased spouse’s income. By way of example, a wife can base her Social Security benefits either on her own Social Security contributions, or on those of her husband (assuming the husband is retired). Should the husband die, his benefits would cease. If the wife’s benefits were based on her husband’s contributions, after his death her individual benefit would double, from 21 50% to 100% of the husband’s benefit. In spite of this, the spouse would receive only (100%/[100% + 50%]) = 67% of the household’s original Social Security benefits.7 If the wife’s benefits were based on her own contributions, her individual benefit would increase only if her husband’s contributions were higher than hers. The relative decrease in household benefits would then depend on each spouse’s contributions. Our fixed effects regression estimates, which capture the average drop in income at the time of a spouse’s death, average over all these cases. 6.1.3 Medical spending Because our model explicitly accounts for Medicaid payments, which depend on a household’s savings, the measure of medical expenses we need is the sum of out-ofpocket spending by the household and Medicaid payments.8 Although the AHEAD contains detailed data on out-of-pocket medical spending, it does not contain Medicaid payments. A key empirical contribution of this paper is to construct the sum of these two components by combining out-of-pocket medical spending information in the AHEAD with Medicaid payment data from the Medicare Current Beneficiary Survey (MCBS). The MCBS contains extremely high quality administrative and survey information on both Medicaid payments and out-of-pocket medical spending (De Nardi, French, Jones, and McCauley [21]). Its main limitation for our purposes is that the data are collected at the individual, rather than the household level: although the MCBS contains marital status, it lacks information on the medical spending or health of the spouse. To exploit the strengths of each dataset, we use the conditional mean matching procedure described in Appendix D to impute a Medicaid payment for each individual in each period in the AHEAD. The procedure preserves both the mean and the distribution of combined medical spending, conditional on Medicaid recipiency, age, income, out-of-pocket spending and other health and medical utilization variables. The regression of Medicaid payments on these variables has an R2 statistic of 0.67, suggesting that our predictions are accurate. We model medical spending as a function of a polynomial in age, marital status, health of each spouse (at the beginning and the end of the period), interactions of these variables, and persistent and transitory spending shocks (see Equation (6)). 7 These calculations also operate in reverse, with “husband” and “wife” exchanged. To perform the calculations, we make several assumptions, including that both spouses begin receiving benefits at the normal retirement age. See https://socialsecurity.gov/planners/retire/yourspouse. html for more details of calculation of spousal benefits. 8 A large share of the elderly’s medical spending is covered by Medicare co-pays, which are not means-tested. Our measure is net of those co-pays. 22 To estimate these profiles, we use fixed effects rather than ordinary least squares for two reasons. First, differential mortality causes the composition of our sample to vary with age, whereas we are interested in how medical expenses vary for the same individuals as they grow older. Second, cohort effects are ubiquitous and using fixed effects we construct the profiles for our model using the mean fixed effect for the cohort with an average age of 74 in 1996. We estimate the medical spending persistence parameter ρm and the variance of the transitory and persistent medical spending shocks using a standard error components model. Our estimates imply that approximately 40% of the cross-sectional variation in log medical spending is explained by observables. Of the remaining cross-sectional variation, 40% comes from the persistent shock and 60% from the transitory shock. Our estimated value of ρm is 0.85. See Appendix G.4 for details. Mean OOP + Medicaid by Income Tercile: Model, Singles Mean OOP + Medicaid by Income Tercile: Model, Couples 28 55 50 top 24 OOP + Medicaid (000s of 2014 dollars) OOP + Medicaid (000s of 2014 dollars) 26 second third 22 20 18 16 14 12 10 8 top second third 40 35 30 25 20 15 10 6 4 72 45 76 80 84 88 92 96 100 Age 5 72 76 80 84 88 92 96 100 Age (a) (b) Figure 5: Mean medical spending. Panel (a): current singles. Panel (b): current couples. Figure 5 shows model-predicted mean medical spending for the AHEAD cohort aged 72-77 in 1996, where we simulate the households through age 100.9 Three general trends are apparent. First, medical expenses rise rapidly with age, in part because older individuals are more likely to reside nursing homes or die.10 Second, spending rises only modestly with PI, if at all. Although our underlying coefficients show that medical expenses rise with PI when health is held constant, lower-income households are often in worse health. Finally, the spending of singles (left panel) is roughly half that of couples (right panel). 9 Because of the structure of our medical spending model, which requires two periods of health realizations, the simulation results start at age 76. 10 Our simulations include end-of-life spending. 23 Figure 6 shows out-of-pocket medical spending. We find out-of-pocket expenses by simulating our estimated structural model and calculating Medicaid payments. Subtracting these payments from total medical spending (shown in Figure 5) yields out-of-pocket spending. Because Medicaid covers a larger share of medical expenses in poorer households, out-of pocket spending has a strong income gradient. Likewise, as medical expenses rise with age, the share covered by Medicaid rises as well. This leads out-of-pocket medical spending to rise more slowly with age than total spending, especially among the low-income. For instance, among singles in the bottom PI tercile, out-of-pocket medical spending stays between $3,700 and $5,400; among those at the top, spending rises by a factor of 3, from $5,600 to $18,000. Mean OOP Medex by Income Tercile: Model, Singles Mean OOP Medex by Income Tercile: Model, Couples 19 30 top second OOP Medex (000s of 2014 dollars) OOP Medex (000s of 2014 dollars) 17 32.5 third 15 13 11 9 7 5 3 72 27.5 top second third 25 22.5 20 17.5 15 12.5 10 7.5 76 80 84 88 92 96 100 5 72 76 Age 80 84 88 92 96 100 Age (a) (b) Figure 6: Mean out-of-pocket medical spending. Panel (b): current couples. Panel (a): current singles. To find the implications of our estimated spending process for lifetime medical spending, we update the calculations in Jones, De Nardi, French, McGee and Kirschner [34] to use our latest estimates. In particular, we simulate a large number of demographic and medical spending histories and calculate discounted sums at each age. Figure 7 shows that, as of age 70, single households will on average incur $63,000 of medical expenses over the remainders of their lives, and that nearly 10% of them will incur medical expenses in excess of $126,000. The corresponding values for initial couples are $129,000 for the mean and $234,000 for the 90th percentile. Some of this total may be picked up by Medicaid, but it is nonetheless significant. The spending statistics for singles rise with age until age 90. This is not merely an artifact of mortality bias. Although older people have less time to live, Figure 5 shows that medical expenses rise with age. In other words, at any point in time, most survivors have yet to incur the bulk of their lifetime expenses. As De Nardi, 24 French, and Jones [18] and others have pointed out, this “backloading” of medical expenses gives it the potential to be an important driver of saving. In contrast, as couples become singles their medical expenses fall. Remaining Lifetime Medical Spending of Surviving Households Medical Spending (000s of 2014 dollars) 250 225 200 175 90th %tile, Initial Couples 90th %tile, Initial Singles Means, Initial Couples Means, Initial Singles 150 125 100 75 50 25 0 70 74 78 82 86 Age 90 94 98 102 Figure 7: Mean and 90th percentile of remaining lifetime medical spending for surviving households, initial singles and initial couples 6.2 Second-step estimation results Table 3 reports our estimated preference parameters. Our estimate of ν, the coefficient of relative risk aversion, is 3.7, a value similar to that estimated for retired singles in De Nardi, French, and Jones [18], and to those typically used in the life cycle literature. Our estimate of consumption equivalence scale η is almost identical to the “modified” OECD scale. It also lies within the confidence interval estimated by Hong and Rı́os-Rull [31], who estimate it using data on life insurance holdings and a structural life-cycle model of consumption and saving decisions. After reviewing a variety of estimates, Fernandez-Villaverde and Krueger [26] argue in favor of similar economies of scale. In combination with our estimated value of ν, our estimated η implies that the consumption of a couple must be 1.63 times that of a single to equate the marginal utility of consumption across the two groups.11 Intuitively, η is identified 11 To see this note that the marginal utilities of singles and couples are 25 ∂uS (cS ) ∂c = (cS )−ν and ν: coefficient of RRA 3.698 (0.09073) η: consumption equivalence scale 1.514 (0.2036) φ0 : bequest intensity, single (in 000,000s) 133.3 (25.5) κ0 : bequest curvature, single (in 000s) 9,175 (1,020) φ1 : bequest intensity, surviving spouse 7,581 (1,411) κ1 : bequest curvature, surviving spouse (in 000s) 244.7 (25.77) cmin (f = 1): annual consumption floor, singles 4,108 (113.8) Notes: Standard errors in parentheses. We set cmin (f = 2) = 1.5 · cmin (f = 1). Table 3: Estimated second-step parameters from the extent to which couples save, relative to singles, to self-insure future risks. Our estimates for θ0 and κ0 govern the strength of terminal bequest motives. Panel (a) of Figure 8 plots their implications in terms of the share of resources that would be left as bequests by single households who know that they will die next period and by married couples who know that both members of their couple will die next period for sure. In that case, single households with less than $29,600 in resources will leave no bequests,12 while those having $100,000 will leave over 75% of their resources in bequests. In contrast, couples require almost $50,000 in resources before they choose to leave bequests to other heirs. These differences in behavior for singles and couples do not come from bequest parameters, which we constrain to be the same when there are no survivors. Rather, they are due to the difference in household size between couples and singles and to the presence of equivalence scales. Together, these parameters imply that bequests are luxury goods for singles (i.e., the threshold level of wealth above which households leave bequests is higher) and that they are even more of a luxury good, but are weaker for couples. Panel ∂uC (cC ) ∂c 12 C −ν 1/ν C 2 2 = η1−ν (c ) . Equating these marginal utilities yields ccS = η1−ν . This is an annual value. In our two-year framework, the threshold will (approximately) double. See Appendix H. 26 Share Allocated to Bequests Share Allocated to Bequests 1 .8 .6 .4 .2 0 0 50 100 150 200 1 .8 .6 .4 .2 0 0 50 Cash on Hand (1000s USD) Singles 100 150 200 Cash on Hand (1000s USD) Couples Singles (a) DFJ Lockwood ABCST L&T (b) Figure 8: Estimated Bequest Motives. Panel (a): expenditure share allocated to bequests for couples (solid) and singles (dashed) facing certain death in the next period. Panel (b): comparing our results for singles (dashed red), with those in De Nardi, French and Jones[18] (DFJ), Lockwood [41], Ameriks, Briggs, Caplin, Shapiro and Tonetti [3] (ABCST), and Lee and Tan [40] (L&T). (b) of Figure 8 compares our estimated bequest motives for singles with those from several previous papers and displays that our estimated bequest motives for singles are of similar strength and slightly more of a luxury good than those estimated by Lockwood [41] and Lee and Tan [40]. Our estimated consumption floor implies that the consumption of a single household is bounded below at $4,108 per year. The floor for couples is set to 150% of this value. Our estimated floor is best interpreted as an “effective” consumption floor that accounts for transactions costs, stigma, and other determinants of Medicaid usage. 6.3 Model fit We require our model to match the observed heterogeneity in savings and Medicaid recipiency by age, cohort, income, and wealth levels. More specifically, we target the 25th percentile, the median, and the 75th percentile of wealth, conditional on PI tercile, by cohort and age, for both couples and singles. We also target average Medicaid recipiency by cohort, age, and PI tercile. These targets help separately identify household precautionary savings, bequest motives, and the degree to which bequest motives are luxury goods. Figure 9 plots median net worth by age, PI, and birth cohort for current couples 27 (Panel (a)) and current singles (Panel (b)), in the data (solid line) and from our model (dashed line). For clarity, we only show targeted moments for two of our birth cohorts in each panel of these set of graphs, but the model fit is similarly good across cohorts (See Appendix I). Figure 9 highlights several important features of saving behavior during retirement. First, higher-PI households dissave more slowly than lower-PI households throughout retirement and, second, that this force is stronger for couples, as long both members of the couple stay alive. More specifically, singles in the lowest PI tercile have almost no wealth, while couples in the same PI tercile hold onto their retirement savings until age 90. Singles in the middle PI tercile start decumulating their wealth from the time that we start tracking them (age 75), while couples of the same age and PI tercile display flat asset trajectories over much of the period that we observe them. Finally, households in the highest PI tercile, and especially current couples, have highest savings during their retirement years. Our estimated model matches all of these aspects of the data. Median Net Worth: Data (Solid) vs. Model (Dashed), Singles Median Net Worth: Data (Solid) vs. Model (Dashed), Couples 270 600 550 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 240 210 180 150 120 90 60 30 0 72 500 450 400 350 300 250 200 150 100 50 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) (b) Figure 9: Median wealth. Panel (a): current singles. Panel (b): current couples. Solid lines: cohorts aged 72-77 and 84-89 in 1996. Dashed lines: model simulations. Figure 10 compares the Medicaid recipiency profiles generated by our model (dashed line) to those in the data (solid line) and shows that our model also matches the patterns of Medicaid usage across age and household PI. Matching these features of Medicaid recipiency means that our model is able to capture household reliance on social insurance, in addition to matching the risk of the catastrophic medical spending that leads higher-income households to utilize Medicaid. Matching Medicaid recipiency rates is important for identifying the consumption floor, which in 28 Medicaid Recipiency: Data (Solid) vs. Model (Dashed), All 0.8 0.7 Medicaid Recipiency 0.6 0.5 0.4 0.3 0.2 0.1 0 72 76 80 84 88 92 96 100 Age Figure 10: Medicaid Recipiency. Solid lines: cohorts aged 72-77 and 84-89 in 1996. Dashed lines: model simulations. turn is important to evaluate the role of precautionary savings because a low consumption floor implies high consumption risk, and thus a high role of precautionary savings. The role of precautionary savings is in turn important for understanding bequest motives, since if precautionary savings are high, the role for bequest motives is weaker. Appendix I shows all the moments we match and our model’s fit. Matching multiple asset percentiles across multiple PI terciles requires the model to reproduce saving behavior across a wide variety of income flows, household structures and initial asset holdings. High-PI households and couples save at higher rates than their lowPI and single counterparts, and this is especially visible at the 25th and 75th asset percentiles as shown in Figures 18 and 19. Among couples at the top of the PI distribution, wealth rises with age at the 75th percentile of wealth, but are constant at the 25th percentile. Among singles at the bottom of the PI distribution, wealth is constant at the 75th percentile, while at the 25th percentile it falls rapidly, or is 0 throughout. This variation disentangles the saving motives operating within our model. 7 Model validation To build additional trust in the predictions of our estimated model, we turn to comparing its implications with some aspects of the data that are important given our question but are not targeted by our estimation procedure. These aspects of the data are medical spending and net worth dynamics around the time of death of one 29 -5 1000s of 2014 Dollars 0 5 10 15 of the spouses in a couple. -6 -4 -2 0 Years from Death Data 95% CI 2 4 Simulations Figure 11: Average change in out-of-pocket medical expenses around death of a spouse, initial couples, relative to similar households with no death. Estimates (solid blue line) and their 95% confidence interval (shaded region) from the AHEAD data and differences-in-differences estimates (dashed red line) from model-simulated data. Death dates are centered at year 0. Our first validation exercise focuses on the rise in out-of-pocket medical spending around the death of a spouse and is shown in Figure 11. In this graph, we use the same matched-differences-in-differences estimation as described in Figure 2 for both actual and model-generated data. Unlike Figure 2, this figure excludes Medicaid payments and focuses on the out-of-pocket component because it is this part that captures the extent to which asset declines at death are driven by medical expenses. Figure 11 shows that our model generates a trajectory of out-of-pocket medical expenses around the time of death that matches almost exactly the one observed in the data. This is notable for two reasons. First, because we did not require the model to match it by construction. Second, because although the sum of out-of-pocket and Medicaid payments is exogenous in our model, a household’s out-of-pocket spending component is partly determined by their wealth, which determines Medicaid eligibility and thus the share of expenses paid by Medicaid. Consequently, the fact that our model closely replicates the dynamics of out-of-pocket medical spending through the dynamics of saving and means-tested transfers, demonstrates that we endogenously generate the correct exposure of wealth to medical expense risk. Our second validation exercise returns to the asset dynamics of couples when one of their spouses dies. Here, we display results not only for our full sample of matched initial couples (Panel (a)) but also for those in the top half of the 1996 30 100 50 1000s of 2014 Dollars -200 -100 0 0 -300 1000s of 2014 Dollars -150 -100 -50 -400 -200 -6 -4 -2 0 Years from Death Data 95% CI (a) Initial Couples 2 4 -6 Simulations -4 -2 0 Years from Death Data 95% CI 2 4 Simulations (b) Initial Couples, Above Median Wealth Figure 12: Average change in wealth around death of a spouse for initial couples (Panel (a)) and initial couples in the top half of the 1996 wealth distribution (Panel (b)). Differences-in-differences estimates (solid blue line) and their 95% confidence interval (shaded region) from the AHEAD data and differences-in-differences estimates (dashed red line) from model-simulated data. Death dates are centered at year 0. initial wealth distribution (Panel (b)). We add this panel because matching well the asset dynamics of wealthier households around death implies that our luxury good specification of bequest motives matches households’ saving behavior of wealthier households. Figure 12 shows that the net worth dynamics around a spousal death generated by our model lie well within the 95% confidence interval of their AHEAD data counterpart, even when we condition by household wealth. For the full sample of initial couples both the data and the simulations predict a decline of $160,000 while for the sample of wealthier households we replicate the decline of $255,000 around death. Our model thus captures well both medical spending and net worth dynamics around the time of death of a spouse across the wealth distribution, which implies that it also matches the extent to the decline in wealth at death which is due to changes in out-of-pocket medical expenses as opposed to bequests to other heirs than the spouse. These features of the data, in turn, have important implications about the relative importance of precautionary savings and bequest motives, forces that help identify risk aversion and bequest motive parameters (given the consumption floor, which is pinned down by Medicaid recipiency). 31 8 What drives savings? Median Net Worth: Experiment (Solid) vs. Baseline (Dashed), Initial Singles 250 500 Median Net Worth: Experiment (Solid) vs. Baseline (Dashed), Couples 225 450 200 400 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) To quantify the savings determinants of retired households, we perform several decomposition exercises, calculating the changes in total retiree wealth that occur as we switch off various saving motives. In doing so, we assume that household wealth accumulation prior to retirement is unchanged. Because of this, our decompositions provide a lower bound on the savings changes, compared to an environment in which people can re-optimize their savings at younger ages. Our approach allows us to focus on the drivers of post-retirement savings and to better understand policy reforms that are not announced far in advance. More specifically, to perform our decomposition, we fix our model parameters at their estimated values and change one feature of the model at a time, recompute the optimal saving decisions, and compare the resulting asset accumulation profiles to those generated by our baseline model. To make the differences clearer, we display the simulated median asset profiles for our youngest cohort, which has an average age of 74 in 1996. When couples lose a spouse and become singles, we drop them from the analysis. We do the same when singles die. Therefore, the profiles embed the same mortality bias that is in the data. 175 150 125 100 75 50 25 0 72 350 300 250 200 150 100 50 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) Singles (b) Couples Figure 13: Median wealth by permanent income tercile: baseline model (dashed lines) and model with no medical expenses (solid lines), initial singles (Panel (a)) and current couples (Panel (b)). We start by asking whether medical expenses are important drivers of saving behavior. To answer this question, we zero out both expected and realized medical expenses and report the corresponding asset trajectories. 32 Panel (a) of Figure 13 shows that the poorest singles, in the bottom PI tercile, save almost nothing to cover medical expenses: with little wealth, they choose to rely on Medicaid to pay for their medical needs. It is singles in the middle PI tercile who are most sensitive to medical spending: at age 84, removing medical expenses would cause their wealth to fall by nearly 40%, from $87,000 to $54,000. These individuals have enough wealth for medical expenses to pose a concern, but not enough for bequest motives to dominate. At the top PI tercile, eliminating medical spending has a smaller effect on savings. At age 84, their wealth would fall by 7%, from $238,000 to $221,000. The wealth they hold to leave bequests allows them to self-insure against medical expenses (as noted by Dynan, Skinner and Zeldes [23]). Interestingly, after age 90 the savings of richer singles are higher without medical expenses. This is because eliminating medical spending leaves households with more resources to consume or save. Panel (b) shows that the effect of medical spending on couples’ wealth is relatively small at every PI tercile, suggesting that medical expenses are only one of the reasons why they continue to hold wealth at very advanced ages. The largest effect occurs for couples in the top PI tercile, for whom net worth at age 84 declines from $437,000 to $417,000, by 5%. Comparing Panel (a) and Panel (b) shows that medical expenses have a larger effect on the savings of singles than of couples. Median Assets: Experiment (Solid) vs. Baseline (Dashed), Couples 500 225 450 200 400 Assets (000s of 2014 dollars) Assets (000s of 2014 dollars) Median Assets: Experiment (Solid) vs. Baseline (Dashed), Initial Singles 250 175 150 125 100 75 50 25 0 72 350 300 250 200 150 100 50 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) Singles (b) Couples Figure 14: Median wealth by permanent income tercile: baseline model (dashed lines) and model with no bequest motives (solid lines). Singles (Panel (a)) and couples (Panel (b)). We next consider the role of our estimated bequest motives by shutting down the utility from leaving any bequests (that is, we set φ0 = φ1 = 0) while holding all other 33 Median Net Worth: Experiment (Solid) vs. Baseline (Dashed), Initial Singles 250 500 Median Net Worth: Experiment (Solid) vs. Baseline (Dashed), Couples 225 450 200 400 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) preference parameters constant. Figure 14 reports the resulting effects on savings. For singles, eliminating bequest motives has little effect on the median savings of the bottom two PI terciles, although they do reduce median savings at the top tercile from $238,000 to $197,000 (17%) at age 84. For couples, bequest motives are more important and have a large effect on savings for the top two PI terciles. Panel (b) of Figure 14 shows while savings barely change for couples in the bottom third of the PI distribution, eliminating bequest motives reduces savings at age 84 from $199,000 to $172,000 (14%) in the middle tercile and from $437,000 to $325,000 (26%) in the top tercile. Thus, more couples than singles save to leave bequests, and bequest motives affect their savings by a much larger extent. This is for at least two reasons. First, couples are wealthier, and bequests are a luxury good. Second, couples not only receive direct utility from leaving bequests at the time of one spouse’s death, but also indirect utility from the expectation that the last surviving spouse will leave bequests at her or his death. 175 150 125 100 75 50 25 0 72 350 300 250 200 150 100 50 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) Singles (b) Couples Figure 15: Median wealth by permanent income tercile: baseline model (dashed lines) and model with no bequest motives and no medical expenses (solid lines). Singles (Panel (a)) and couples (Panel (b)). While medical expenses and bequest motives taken each in isolation have relatively modest effects on the saving behavior of most households, their interaction could have a much larger impact. To quantify the importance of these interactions, we eliminate both medical expenses and bequest motives. The resulting savings profiles are displayed in Figure 15. Because singles in the bottom and middle PI tercile are too poor for bequest motives to actively influ34 ence their saving behavior (as seen in Figure 14), the effects on their savings from eliminating both bequest motives and medical expenses are very similar to those in the case in which we only eliminate medical expenses (Figure 13). In contrast, for those in the top PI tercile, for whom the bequest motive is active, savings falls from $238,000 to $82,000 (by 66%). This drop is almost three times larger than the sum of the drops in the previous two experiments (24%). This is due to the way bequest and precautionary saving motives interact. When medical expenses are uncertain, households that reduce their consumption to finance future medical expenses may end up dying with unspent wealth. In the absence of bequest motives, these accidental bequests have no value. Bequest motives reduce the opportunity cost of self-insuring through saving.13 At the same time, medical expenses introduce uncertainty over the size of the household’s bequest. This might generate precautionary saving aimed at increasing the size of the bequest. Panel (b) of Figure 15 displays the corresponding results for couples. The savings of couples in the bottom PI tercile fall from $62,000 to $41,000 (by 34%) at age 84. For couples in the middle PI tercile, savings falls from $199,000 to $99,000 (by 50%) at age 84. Because bequest motives are active for this group, unlike for the corresponding singles, the interaction of bequest and precautionary saving motives is much larger for middle PI couples than singles. This interaction is also very strong for couples in the top PI tercile. Rather than accumulating wealth in retirement, couples at the top tercile would run down their wealth holdings from $437,000 to $150,000 (66%) by age 84. Finally, couples save to provide resources to the surviving spouse after the death of the first spouse. This is important because while at age 70 married men live on average 11.5 years and women live 15.8 years (or on average 13.7 years), the oldest survivor lives on average 17.9 years. During these remaining years the survivor faces lower income and high medical spending, including medical spending from an elevated risk of being in a nursing home when single. Studies omitting this risk will understate this precautionary savings motive. To quantify the effect of this force, Figure 16 compares our baseline savings with those implied by our model when we shut down the weight on the utility of the surviving spouse. It shows that the desire to leave resources to the surviving spouse upon a partner’s death has a large effect on savings at all PI levels. For instance, among couples in the lowest PI tercile, median wealth at age 84 would fall from $62,000 in the benchmark to $46,000 when there is no concern for the welfare of the surviving spouse, a 26% drop. Median-asset couples in the second PI tercile at age 13 A similar argument applies to saving against longevity risk and for the decision to self-insure instead of purchasing formal Long Term Care insurance (See Lockwood [41] for an example). 35 Median Net Worth: Experiment (Solid) vs. Baseline (Dashed), Couples 500 Net Worth (000s of 2014 dollars) 450 400 350 300 250 200 150 100 50 0 72 76 80 84 88 92 96 100 Age Figure 16: Median wealth for couples by permanent income tercile: baseline model (dashed lines) and model with no weight on the surviving spouse (solid lines). 84 hold $199,000 in the benchmark economy and $165,000 in this counterfactual, a 17% drop. At the highest PI tercile, median wealth falls from $437,000 to $368,000 at the same age, a 16% drop. Taken together, Figures 13-16 show that medical spending is a more important determinant of savings for singles than for couples. In contrast, bequest motives affect the savings of couples and high-PI households more than those of singles and lowPI households. The figures also show that the interaction of precautionary savings and bequests are most important for couples and high-PI households. Moreover, providing for the surviving spouse is important for all couples. For clarity, Figures 13-16 only display profiles for households whose composition does not change. Thus, couples do not experience a decline in net worth upon death of the first spouse and singles do not include the newly widowed, who tend to be richer. To assess the saving motives across all retired households simultaneously, Table 8 aggregates the preceding simulations.14 The first column of Table 8 shows summary statistics for our baseline specification. The distribution of wealth is right skewed, with the mean ($372,300) more than twice as large as the median ($150,100). The second column shows the effects of eliminating medical expenses. Consistent with Figure 13, medical expenses have little effect on the saving of the wealthiest: the 75th percentile falls by 2.7% and the 14 Formally, we take the same joint distribution of household initial conditions as in the data and simulate a panel accounting for the realisations of individual mortality, health and medical expenditure shocks. We use all retired households alive at each age in the simulated panel, except for their wealth when first observed which we take as given – and calculates summary statistics for wealth. 36 Percentage Change from Baseline 25th Percentile Median 75th Percentile Mean Baseline Wealth No Medical Expenses No Bequest Motives No Bequests or Medical Expenses No Weight on Surviving Spouse 50.0 150.1 391.9 372.3 -61.4% -20.9% -2.7% -3.0% 8.3% 6.9% -8.4% -17.4% -61.9% -45.6% -42.8% -44.2% -83.3% -36.33% -25.3% -27.4% Notes: Wealth in baseline simulations are expressed in thousands of 2014 dollars. Table 4: Aggregate Implications of Retiree Saving Motives mean, which is dominated by the richest households, falls only 3.0%. The effects are much larger in the lower half of the wealth distribution: the median falls by 20.9%. When medical expenses are eliminated, middle-wealth singles reduce their wealth and middle-wealth couples transfer a larger fraction of their estates from surviving spouses to nonspousal heirs. The effects reverse in the third column, where we eliminate bequest motives. Bequests are luxury goods and their elimination causes mean wealth to fall by 17.4%. In contrast, eliminating bequest motives leads median wealth to increase by 6.9%. Couples save less when bequest motives are removed, but when one spouse dies all of the estate is left to the surviving spouse. For middle-wealth households, the latter effect dominates. Taken together, the second and third column reinforce our earlier findings that medical expenses are more important for poorer households and singles, while bequest motives are more important for the rich. The fourth column also reinforces Figures 13-16, in that the effect of eliminating medical expenses and bequest motives together is large and exceeds the sum of the individual changes. This is a key finding because the previous literature has focused on the role of precautionary saving motives due to medical spending without a bequest motive. Given that our results present robust evidence of bequest motives, and that their interaction with medical expenses has a very large impact on savings, our results show that it is essential to model both bequest motives and medical spending. The final column of Table 8 shows that setting to zero the weight that couples 37 place on the surviving spouse also has large effects. This is not surprising because this experiment leads couples to bequeath all their wealth to nonspousal heirs upon the death of the first spouse, but highlights that it is important to model couples and spousal deaths to fully capture life cycle saving motives. 9 Conclusions This paper documents new facts about the retirement savings and risks of couples and singles. It then develops and estimates a rich structural model that contains these risks and explains their savings. In terms of new facts, we compare the saving behavior of couples and singles. Households with high PI, and especially couples, tend to accumulate wealth, whereas those with low PI, and especially singles, tend to decumulate it. Furthermore, upon the death of a spouse, there are large drops in wealth, and bequests are often distributed. To understand the relative roles of risks versus bequest motives as drivers of these savings decisions, it is crucial to properly measure the risks faced by couples and singles. While medical expenses for couples and singles are not very different on a per capita basis until someone dies, medical spending jumps at the time of a spouse’s death. At the same time, household income falls, exposing households to a substantial loss in financial resources upon a spouse’s death. We embed these estimated risks within a rich structural model of savings that incorporates heterogeneity in financial resources, life expectancy, and medical expense risks, including the large jump in medical expenses that is typically observed before a death. Our model also accounts for means-tested social insurance and bequest motives. It allows households to leave bequests upon death of the last survivor. In addition, and very importantly, households can choose to distribute bequests when the first member of the couple dies (net worth is split between the surviving spouse and other heirs). We estimate our model using MSM, targeting the savings of both couples and singles and Medicaid recipiency rates by cohort, PI, and wealth levels. Our model fits all of these important features of the data well. In addition, although we do not require our model to match the dynamics of out-of-pocket medical spending and wealth around the time of a spouse’s death, our model also matches those well. This validates key predictions of our model. Finally, we use our estimated model to measure the extent to which the savings of retirees are driven by medical expenses, bequest motives, and their interaction. We uncover several new findings. Medical expenses are more important for singles, 38 whereas bequest motives are more important for couples. This is due to two main reasons. First, couples not only receive direct utility from leaving bequests at the time of one spouse’s death, but also indirect utility from the expectation that the last surviving spouse will leave bequests upon their own death. Second, bequests are a luxury good, and couples are wealthier. We also find that introducing bequest motives reduces the effects of medical expenses on savings. This is because savings intended for medical spending that go unspent due to early death are still valued when bequest motives are present. This finding shows that we need to model both medical spending and bequests. Lastly, although medical spending is important, the desire to transfer resources to one’s heirs other than their spouse is key to understanding the wealth drops at the death of a first spouse. Hence, to understand the drivers of savings, it is essential to model the bequest motives and medical spending of both couples and singles. Why people save has important policy implications. Our results suggest that couples and high-PI singles can easily self-insure against medical spending risk because they save to leave bequests. Low-PI singles are well insured through Medicaid and do not need to save for medical expenses. In contrast, middle-PI singles set aside significant amounts of wealth for precautionary purposes because they are not rich enough to want to leave bequests but are too rich to qualify for Medicaid. It is thus the middle-PI singles who would respond most to changes in public health insurance. Medicaid also includes rules determining its relative generosity for couples and singles. Evaluating whether these rules are well designed requires carefully modeling and quantifying the saving motives of all retirees. 39 References [1] Andrew B Abel and Mark Warshawsky. Specification of the joy of giving: Insights from altruism. The Review of Economics and Statistics, 70(1):145, 1988. [2] Joseph G. Altonji and Lewis M. Segal. Small sample bias in gmm estimation of covariance structures. Journal of Business and Economic Statistics, 14(3):353– 366, 1996. [3] John Ameriks, Joseph Briggs, Andrew Caplin, Matthew D. Shapiro, and Chris Tonetti. Long-term care utility and late in life saving. Journal of Political Economy, 128(6):2375–2451, 2020. [4] John Ameriks, Andrew Caplin, Steven Laufer, and Stijn Van Nieuwerburgh. The joy of giving or assisted living? using strategic surveys to separate bequest and precautionary motives. Journal of Finance, 66(2):519–561, 2011. [5] James Andreoni. Giving with impure altruism: Applications to charity and ricardian equivalence. Journal of Political Economy, 97:1447–1458, 1989. [6] Orazio P. Attanasio and Carl Emmerson. Mortality, health status, and wealth. Journal of the European Economic Association, 1(4):821–850, 2003. [7] Gary S. Becker and Nigel Tomes. Human capital and the rise and fall of families. Journal of Labor Economics, 4(3):1–39, 1986. [8] B. Douglas Bernheim, Andrei Shleifer, and Lawrence H. Summers. The strategic bequest motive. Journal of Political Economy, 93(5):1045–1076, 1985. [9] R. Anton Braun, Karen A. Kopecky, and Tatyana Koreshkova. Old, sick, alone and poor: A welfare analysis of old-age social insurance programs. The Review of Economic Studies, 84(2):580–612, 2017. [10] Meta Brown. End-of-life transfers and the decision to care for a parent. mimeo, 2003. [11] Moshe Buchinsky. Recent advances in quantile regression models: A practical guideline for empirical research. Journal of Human Resources, 33:88–126, 1998. [12] Christopher D. Carroll. Why do the rich save so much? NBER Working Paper 6549, 1998. 40 [13] Gary Chamberlain. Comment: Sequential moment restrictions in panel data. Journal of Business & Economic Statistics,, 10(1):20–26, 1992. [14] Martin David, Roderick J. A. Little, Michael E. Samuhel, and Robert K. Triest. Alternative methods for cps income imputation. Journal of the American Statistical Association, 81(393):29–41, 1986. [15] James B. Davies and Antony F. Shorrocks. The distribution of wealth. In Anthony B. Atkinson and François Bourguignon, editors, Handbook of Income Distribution, pages 605–675. Handbooks in Economics, vol. 16. Amsterdam; New York and Oxford: Elsevier Science, North-Holland, 2000. [16] Mariacristina De Nardi. Wealth inequality and intergenerational links. Review of Economic Studies, 71(3):743–768, 2004. [17] Angus Deaton. Saving and liquidity constraints. Econometrica, 59(5):1221– 1248, 1991. [18] Mariacristina De Nardi, Eric French, and John Bailey Jones. Why do the elderly save? the role of medical expenses. Journal of Political Economy, 118(1):39–75, 2010. [19] Mariacristina De Nardi, Eric French, and John Bailey Jones. Medicaid insurance in old age. American Economic Review, 106(11):3480–3520, 2016. [20] Mariacristina De Nardi, Eric French, and John Bailey Jones. Savings after retirement: A survey. Annual Review of Economics, 8(1):177–204, 2016. [21] Mariacristina De Nardi, Eric French, John Bailey Jones, and Jeremy McCauley. Medical spending of the u.s. elderly. Fiscal Studies, 37(3-4):717–747, 2016. [22] Darrell Duffie and Kenneth J. Singleton. Simulated moments estimation of markov models of asset prices. Econometrica, 61(4):929–952, 1993. [23] Karen E. Dynan, Jonathan Skinner, and Stephen P. Zeldes. The importance of bequests and life-cycle saving in capital accumulation: A new answer. American Economic Review, 92(2):274–278, 2002. [24] Itzik Fadlon and Torben Heien Nielsen. Family health behaviors. American Economic Review, 109(9):3162–91, 2019. [25] Daniel Feenberg and Jonathan Skinner. The risk and duration of catastrophic health care expenditures. Review of Economics and Statistics, 76:633–647, 1994. 41 [26] Jesús Fernández-Villaverde and Dirk Krueger. Consumption over the Life Cycle: Facts from Consumer Expenditure Survey Data. The Review of Economics and Statistics, 89(3):552–565, 2007. [27] Eric French and John Bailey Jones. On the distribution and dynamics of health care costs. Journal of Applied Econometrics, 19(4):705–721, 2004. [28] Eric French and John Bailey Jones. The effects of health insurance and selfinsurance on retirement behavior. Econometrica, 79(3):693–732, 2011. [29] Miguel Gouveia and Robert P Strauss. Effective federal individual income tax functions: An exploratory empirical analysis. National Tax Journal, pages 317– 339, 1994. [30] Nezih Guner, Remzi Kaygusuz, and Gustavo Ventura. Income taxation of us households: Facts and parametric estimates. Review of Economic Dynamics, 17(4):559–581, 2014. [31] Jay H. Hong and José-Vı́ctor Rı́os-Rull. Life Insurance and Household Consumption. The American Economic Review, 102(7):3701–3730, 2012. [32] R. Glenn Hubbard, Jonathan Skinner, and Stephen P. Zeldes. Precautionary saving and social insurance. Journal of Political Economy, 103(2):360–399, 1995. [33] Michael D. Hurd. Mortality risk and bequests. Econometrica, 57(4):779–813, 1989. [34] John Bailey Jones, Mariacristina De Nardi, Eric French, Rory McGee, and Justin Kirschner. The lifetime medical spending of retirees. Economic QuarterlyFederal Reserve Bank of Richmond, 104(3):103, 2018. [35] John Bailey Jones, Mariacristina De Nardi, Eric French, Rory McGee, and Rachel Rodgers. Medical spending, bequests, and asset dynamics around the time of death. Economic Quarterly, Fourth Quarter:135–157, 2020. [36] Thomas F. Juster, James P. Smith, and Frank Stafford. The measurement and structure of household wealth. Mimeo, 1998. [37] Ruud H. Koning. http://www.xs4all.nl/ Kernel: A gauss library for kernel estimation. 42 [38] Karen Kopecky and Tatyana Koreshkova. The impact of medical and nursing home expenses and social insurance policies on savings and inequality. American Economic Journal: Macroeconomics, 6:29–72, 2014. [39] John Laitner, Dan Silverman, and Dmitriy Stolyarov. The role of annuitized wealth in post-retirement behavior. American Economic Journal: Macroeconomics, 10(3):71–117, 2018. [40] Siha Lee and Kegon T. K. Tan. Bequest Motives and the Social Security Notch. October 2019. [41] Lee M. Lockwood. Incidental bequests and the choice to self-insure late-life risks. American Economic Review, 108(9):2513–50, 2018. [42] Samuel Marshall, Kathleen M. McGarry, and Jonathan S. Skinner. The risk of out-of-pocket health care expenditure at the end of life. Working Paper 16170, National Bureau of Economic Research, 2010. [43] Makoto Nakajima and Irina A Telyukova. Home equity in retirement. International Economic Review, 61(2):573–616, 2020. [44] Whitney K. Newey. Generalized method of moments specification testing. Journal of Econometrics, 29(3):229–256, 1985. [45] Whitney K. Newey and Daniel L. McFadden. Large sample estimation and hypothesis testing. In Robert Engle and Daniel L. McFadden, editors, Handbook of Econometrics, Volume 4. Elsevier, Amsterdam, 1994. [46] Ariel Pakes and David Pollard. Simulation and the aysmptotics of optimization estimators. Econometrica, 57(5):1027–1057, 1989. [47] Michael G. Palumbo. Uncertain medical expenses and precautionary saving near the end of the life cycle. Review of Economic Studies, 66:395–421, 1999. [48] Jorn-Steffen Pischke. Measurement error and earnings dynamics: Some estimates from the PSID validation study. Journal of Business & Economics Statistics, 13(3):305–314, 1995. [49] James M. Poterba, Steven F. Venti, and David A. Wise. The composition and drawdown of wealth in retirement. Journal of Economic Perspectives, 25:95–117, 2011. 43 [50] James Powell. Estimation of semiparametric models. In Robert Engle and Daniel L McFadden, editors, Handbook of Econometrics, Volume 4. Elsevier, Amsterdam, 1994. [51] Susann Rohwedder, Steven J. Haider, and Michael Hurd. Increases in wealth among the elderly in the early 1990s: How much is due to survey design? Review of Income and Wealth, 52(4):509–524, 2006. [52] George Tauchen. Finite state markov-chain approximations to univariate and vector autoregressions. Economics Letters, 20(2):177–181, 1986. 44 Appendix A: Sample selection and data handling Our main dataset is the AHEAD data, which is a cohort within the Health and Retirement Study (HRS). To keep the dynamic programming problem manageable, we assume a fixed difference in age between spouses, and we take the average age difference from our data. In our sample, husbands are on average 3 years older than their wives; in the model we assume the difference is one period, or 2 years. To keep the data consistent with this assumption, we drop all households where the wife is more than 4 years older or 10 years younger than her husband. We do not use 1993/94 assets, nor medical expenses, due to underreporting issues (Rohwedder, Haider, and Hurd [51]) and thus begin with 6,047 households. Because we only allow for household composition changes through death, we drop households where an individual enters a household or an individual leaves the household for reasons other than death. Fortunately, attrition for reasons other than death is a minor concern in our data. We drop 401 households who get married, divorced, were same sex couples, or who report making other transitions not consistent with the model, 753 households who report earning at least $3,000 in any period, 171 households with a large difference in the age of husband and wife, and 87 households with no information on the spouse in a household. As a result, we are left with 4,634 households, of whom 1,388 are couples and 3,246 are singles. This represents 24,274 household-year observations where at least one household member was alive. An advantage of the AHEAD relative to other datasets is that it provides panel data on health status and nursing home stays. We assign individuals a health status of “good” if self-reported health is excellent, very good or good, and a health status of “bad” if self-reported health is fair or poor. We assign individuals to the nursing home state if they were in a nursing home at least 120 days since the last interview (on average 60 days per year) or if they spent at least 60 days in a nursing home before the next scheduled interview and died before that scheduled interview. The AHEAD has information on the value of housing and real estate, autos, liquid assets (which include money market accounts, savings accounts, T-bills, etc.), IRAs, Keoghs, stocks, the value of a farm or business, mutual funds, bonds, “other” assets and investment trusts less mortgages and other debts. One problem with asset data is that the wealthy tend to underreport their wealth in virtually all household surveys (Davies and Shorrocks [15]). This could lead us to understate asset levels at all ages. However, Juster, Smith, and Stafford [36] show that the wealth distribution of the AHEAD matches up well with aggregate values for all but the richest 1 % of households. Given that we match the quantiles (conditional on permanent income) 45 rather than means, underreporting at the very top of the wealth distribution should not affect our results. Social Security and pension benefits are included in our measure of annual income, hence differences in Social Security and pension wealth appear in our model as differences in the permanent income measure we use to predict annual income (and other state variables). Appendix B: Computing permanent income We assume that log income evolves as following ln yit = κ(t, fit ) + h(Ii ) + ωit , (20) where κ(t, fit ) is a flexible function of age t and family structure fit (i.e., couple, single man or single woman) and ωit represents measurement error. The variable Ii is the household’s percentile rank in the permanent income (PI) distribution. Since it is a summary measure of lifetime income at retirement, it should not change during retirement and is thus a fixed effect over our sample period. However, income could change as households age and potentially lose a family member. To estimate Equation (20) we first estimate the fixed effects model ln yit = κ(t, fit ) + αi + ωit , (21) which allows us to obtain a consistent estimate of the function κ(t, fit ). Next, note that as the number of time periods over which individual i is observed (denoted Ti ) becomes large, Ti Ti 1 X 1 X ln yit − κ(t, fit ) − ωit = ln yit − κ(t, fit ) = αi = h(Ii ). (22) plim Ti t=1 Ti →∞ Ti t=1 We thus calculate the PIPranking Ii for every household in our sample by taking the i percentile ranking of T1i Tt=1 [ln yit − κ b(t, fit )], where κ b(t, fit ) is the estimated value of κ(t, fit ) from Equation (21). Put differently, we take the mean residual per person from the fixed effects regression (where the residual includes the estimated fixed effect), then take the percentile rank of the mean residual per person to construct Ii . However, we also need to estimate the function h(Ii ), which converts the estimated index Ii back to a predicted level of income that can be used in the dynamic programming model. To do this we estimate the function [ln yit − κ(t, fit )] = h(Ii ) + ωit 46 (23) where the function h(Ii ) is a flexible functional form. In practice we model κ(t, fit ) as a third order polynomial in age, dummies for family structure, and family structure interacted with an age trend. When estimating Equation (23), we replace the function κ(t, fit ) with its estimated value. We model h(Ii ) as a fifth order polynomial in our measure of permanent income percentile. Given that we have, for every member of our sample, t, fit , and estimates of Ii and the functions κ(., .), h(.), we can calculate the predicted value ln ybit = κ b(t, fit )+ b h(Ibi ). It is ln ybit that we use when simulating the model for each household. A regression of ln yit on ln ybit yields a R2 statistic of .74, suggesting that our predictions are accurate. Appendix C: Matching methodology We match each household that suffers a death to a household that did not experience a death at the same date or within the window shown in our graphs, but did experience a death within 6-10 years. (Fadlon and Nielson [24] impose a similar restriction in their analysis of health shocks.) More specifically, for each household we take a matched control household from the set of households who have the same initial (1996) household composition, have a PI percentile within 15 percentage points and where their 1996 wealth differs by no more than 2.5% (or $5,000 in levels). We then use a random number generator to select a single control household from this set giving each household an equal probability of being selected. Our event study specification then uses a sample of 476 households who experience the death of a spouse (2,383 household-year observations) matched to 476 households who do not experience a death in the sample window (2,624 household-year observations). Jones, De Nardi, French, McGee and Rodgers [35] present more results on this. Appendix D: Imputing Medicaid payments and outof-pocket medical expenses Our goal is to estimate the data generating process for the sum of Medicaid payments and out-of-pocket expenses: this is the variable ln(m) in Equation (6) of the main text. The AHEAD data contain information on out-of-pocket medical spending, but not on Medicaid payments. Fortunately, the Medicare Current Beneficiary Survey (MCBS) has extremely high quality information on both Medicaid payments and out-of-pocket medical spending. One drawback of the MCBS, however, is that although it has information on marital status and household income, it does not include information on the medical spending or health of the spouse. 47 We use a two-step imputation procedure to exploit the strenghts of both datasets. First step of imputation procedure We use the MCBS to infer Medicaid payments for recipients, conditional on observable variables that exist in both the MCBS and the AHEAD datasets. Let oopit denote out-of-pocket medical expenses for AHEAD person i at time t, and let M cdit denote the dollar value of Medicaid payments. To impute M cdit , we follow David, Little, Samuhel and Triest [14] and French and Jones [28] and use the following predictive mean matching regression approach. First, using every member of the MCBS sample with a positive Medicaid indicator (i.e., a Medicaid recipient), we regress the variable of interest M cd on the vector of observable variables z, yielding M cd = zβ + ε. Second, for each individual j in the [ MCBS we calculate the predicted value M cdjt = zjt β̂, and for each member of the [ sample we calculate the residual ε̂jt = M cdjt − M cdjt . Third, we sort the predicted [ value M cdjt into deciles and keep track of all values of ε̂jt within each decile. The variables zjt include nursing home status, number of nights spent in a nursing home, an age polynomial, total household income, marital status, self-reported health, race, visiting a medical practitioner (doctor, hospital or dentist), out-ofpocket medical spending, education and death of an individual. Because the measure of medical spending in the AHEAD is medical spending over two years, we take two-year averages of the MCBS data to be consistent with the structure of the AHEAD. The regression of M cdjt on zjt yields a R2 statistic of .67, suggesting that our predictions are accurate. Second step of imputation procedure: Next, for every individual i in the AHEAD sample with a positive Medicaid [ indicator, we impute M cdit = zit β̂ using the value of β̂ estimated using the MCBS. Then we impute εic for each member of the AHEAD sample by finding a random [ [ individual j in the MCBS with a value of M cdjt in the same decile as M cdit in the ] [ AHEAD, and set εit = ε̂jt . The imputed value of M cdit is M cdit = M cdit + εit , and ] the imputed value of the sum mit is m e it = oopit + M cdit . As David, Little, Samuhel and Triest [14] point out, our imputation approach is equivalent to hot-decking when the “z” variables are discretized and include a full set of interactions. The advantages of our approach over hot-decking are two-fold. First, many of the “z” variables are continuous, and it seems unwise to discretize them. Second, we use a large number of observable variables “z” because we find that adding extra variables greatly improving goodness of fit when imputing Medicaid payments. Even a small number of variables generates a large number of hot-decking cells because hot-decking uses a full set of interactions. Thus, in this context, hot 48 decking is too data intensive. We predict Medicaid payments for 3,756 household-wave observations with Medicaid recipiency in the AHEAD and report the results of this imputation exercise in Table 5. Households where the last surviving member died between the previous and current waves of the sample have the largest imputed Medicaid payments. In contrast, couples have the smallest Medicaid payments per individual, but Medicaid payments at the household level are larger than for either single men or women. The imputed Medicaid payments for each spouse in a couple are approximately equal and the results for couples are not driven by the expenditures of only the husband or wife. In Table 5, new widows and widowers (those whose spouse has died between the current and preceding waves of the AHEAD sample) are included in the rows for single men and single women – on average the Medicaid payments of dead spouses are less than 20% of the single households’ Medicaid payments. The distribution of imputed Medicaid expenditures in the AHEAD is close to that in the MCBS. The mean imputed Medicaid expenditure, M cdit , for Medicaid recipients in the AHEAD is $14,050. This is lower than the corresponding value of $16,000 in the MCBS. This difference is due to the distribution of observable variables, z, in both the MCBS and AHEAD. In particular, the number of nights in a nursing home is 20% lower in the AHEAD than in the MCBS, which is in part due to the fact that AHEAD respondents were not in a nursing home when they answered the survey for the first time. Family Structure Dead Couples Single Men Single Women Number of Medicaid Eligible Households Mean Standard Deviation 1,040 287 351 2,078 21,800 15,300 11,600 12,300 31,500 25,700 21,500 21,800 Table 5: Imputed Medicaid payments for Medicaid beneficiaries in AHEAD 49 Appendix E: Computing optimal decision rules and MSM estimates We compute the value functions by backward induction. We start with singles. At time-T , we find the value function and decision rules by maximizing Equation (16), subject to the relevant constraints, using VTg+1 = θ0 (xt − ct ), g = h, w. This yields the value function VTg . We then find the value function and decision rules at time T − 1 by solving Equation (16) with VTg . Continuing this backward induction yields decision rules for periods T − 2, T − 3, ..., tr . We find the value function for new widows and widowers, Vtng , and the split of the estate between spousal and nonspousal heirs, by solving Equation (17), subject to the relevant constraints. These calculations utilize the value function for singles, Vtg , described immediately above. We find the time-T value function for couples, VTc by maximizing Equation (19), subject to the relevant constraints and the value function for the singles, and setting the time-T + 1 continuation value to θ0 (xt − ct ). This yields the value function VTc and the decision rules for time T. We then find the decision rules at time T − 1 by solving Equation (19) using VTc and VTng , g = h, w. Continuing this backward induction yields decision rules for time T − 2, T − 3, ..., 1. More specifically, to solve our value functions, we discretize the persistent component and the transitory components of the health shock, ξ and ζ, using the methods described in Tauchen [52]. We assume a finite number of permanent income categories and we break the space of cash-on-hand into a finite number of grid points. We use linear interpolation within the grid and linear extrapolation outside of the grid to evaluate the value function at cash-on-hand values that we do not directly compute. The end result is a set of decision rules for each possible combination of cash-on-hand, PI, health status, and persistent health shock (ζ). The mechanics of our MSM approach are as follows. We compute life-cycle histories for a large number of artificial households. Each of these households is endowed with an initial value of the state vector (t, ft , xt , I, hsht , hsw t ) that is drawn from the data for 1996. In addition, each household is assigned the entire health and mortality history that is recorded for this same household in the AHEAD data. This way we generate attrition in our simulations that exactly follows the attrition in the data (including by initial wealth and mortality). The simulated medical expenditure shocks ζ and ξ are Monte Carlo draws from discretized versions of our estimated shock processes. Our decision rules, combined with the initial conditions and simulated shocks, allows us to simulate a household’s wealth, medical expenses, health, and mortality. 50 We construct life-cycle profiles from the artificial histories in the same way that we compute them from the real data. We use these profiles to construct moment conditions and evaluate the match using our GMM criterion. As done when constructing the figures from the AHEAD data, we drop cells with fewer than 10 observations from the moment conditions. We search over the parameter space for the values that minimize the criterion. Appendix F details our moment conditions, the weighting matrix in our GMM criterion function, and the asymptotic distribution of our parameter estimates. Appendix F: Moment conditions and the asymptotic distribution of parameter estimates We estimate the parameters of our model in two steps. In the first step, we estimate the vector χ, the set of parameters than can be estimated without explicitly using our model. In the second step, we use the method of simulated moments (MSM) to estimate the remaining parameters, which are contained in the M × 1 vector ∆ = ν, η, φ0 , φ1 , κ0 , κ1 , cmin (ft = S) . ˆ of the “true” preference vector ∆0 is the value of ∆ that minOur estimate, ∆, imizes the (weighted) distance between the estimated life cycle profiles for wealth and Medicaid recipiency in the data and those generated by the model. For each calendar year t ∈ {t0 , ..., tT } = {1996, 1998, ..., 2014},15 we match the 25th, 50th and 75th percentile of wealth for 3 permanent income terciles in 4 birth year cohorts, both for singles and couples. Because the 1996 (period-t0 ) distribution of simulated wealth is bootstrapped from the 1996 data distribution, we match wealth for the period 1998 to 2014, hence 9 waves. In addition, we require each cohortincome-age cell have at least 15 observations to be included in the GMM criterion. As a result, we end up with 414 asset targets. We construct the moment conditions for wealth by building on French and Jones [28] (useful references include Buchinsky [11] and Powell [50]). Suppose that household i of family type f belongs to birth cohort c and permanent income tercile p. Let aqf cqt (∆, χ) denote the model-predicted qth asset quantile for households in households i’s group at time t.16 Assuming that observed wealth have a continuous conditional density, af cqt will satisfy q Pr ait ≤ af cpt (∆0 , χ0 ) | f, c, p, t, household i observed at t = πq , 15 Because we do not allow for macro shocks, in any given cohort, t is used only to identify the individual’s age. 16 The elements of χ include the permanent income boundaries. 51 where πq is the probability value associated with quantile q – when q denotes a median, πq would equal 1/2. The preceding equation can be rewritten as a moment condition. In particular, applying the indicator function produces E 1{ait ≤ aqf cpt (∆0 , χ0 )} − πq | f, c, p, t, household i observed at t = 0. Letting Ip denote the values contained in the pth permanent income tercile, we can convert this conditional moment equation into an unconditional one (e.g., Chamberlain [13]): E [1{ait ≤ aqf cpt (∆0 , χ0 )} − πq ] × 1{fit = f } × 1{ci = c} × 1{Ii ∈ Ip } × 1{household i observed at t} t = 0, (24) for q ∈ {.25, .5, .75}, f ∈ {single, couple}, c ∈ {1, 2, 3, 4}, p ∈ {1, 2, 3},and t ∈ {t1 , t2 ..., tT }. We also match Medicaid recipiency rates. We divide individuals into 4 cohorts, match data from 9 waves, and stratify the data by permanent income, but combine couples and singles. Let ucpt (∆, χ) denote the model-predicted recipiency rate for households in cohort c and permanent income tercile p at time t. Let uit be the {0, 1} indicator that equals 1 when household i receives Medicaid. The associated moment condition is E uit − ucpt (∆0 , χ0 ) × 1{ci = c} × 1{Ii ∈ Ip } × 1{individual i observed at t} t = 0 (25) for c ∈ {1, 2, 3, 4}, p ∈ {1, 2, 3},and t ∈ {t1 , t2 ..., tT }. To summarize, the moment conditions used to estimate our model consist of: the moments for asset quantiles described by Equation (24) and the moments for the Medicaid recipiency rates described by Equation (25). We have a total of J = 505 moment conditions. Suppose we have a dataset of I independent individuals that are each observed at up to T separate calendar years. Let ϕ(∆; χ0 ) denote the J-element vector of moment conditions described immediately above, and let ϕ̂I (.) denote its sample c I denote a J × J weighting matrix, the MSM estimator ∆ ˆ is given analog. Letting W by I c I ϕ̂I (∆; χ0 ), argmin ϕ̂I (∆; χ0 )0 W 1+τ ∆ 52 where τ is the ratio of the number of observations to the number of simulated observations. In practice, we estimate χ0 using the approach described in the main text. Computational concerns, however, compel us to treat χ0 as known in the analysis that follows. Under the regularity conditions stated in Pakes and Pollard [46] and Duffie and ˆ is both consistent and asymptotically normally Singleton [22], the MSM estimator ∆ distributed: √ ˆ − ∆0 I ∆ N (0, V), with the variance-covariance matrix V given by V = (1 + τ )(D0 WD)−1 D0 WSWD(D0 WD)−1 , where S is the variance-covariance matrix of the data; D= ∂ϕ(∆; χ0 ) ∂∆0 (26) ∆=∆0 c I }. is the J×M gradient matrix of the population moment vector; and W = plimI→∞ {W Moreover, Newey [44] shows that if the model is properly specified, I ˆ χ0 )0 R−1 ϕ̂I (∆; ˆ χ0 ) ϕ̂I (∆; 1+τ χ2J−M , where R−1 is the generalized inverse of R = PSP, P = I − D(D0 WD)−1 D0 W. c I converges to S−1 , The asymptotically efficient weighting matrix arises when W the inverse of the variance-covariance matrix of the data. When W = S−1 , V simplifies to (1 + τ )(D0 S−1 D)−1 , and R is replaced with S. Even though the optimal weighting matrix is asymptotically efficient, it can be biased in small samples. (See, for example, Altonji and Segal [2].) Hence, we use a “diagonal” weighting matrix, as suggested by Pischke [48]. This diagonal weighting scheme uses the inverse of the matrix that is the same as S along the diagonal and has zeros off the diagonal of the matrix. We estimate D, S, and W with their sample analogs. For example, our estimate of S is the J × J estimated variance-covariance matrix of the sample data. When 53 estimating this matrix, we use sample statistics, so that apqt (∆, χ) is replaced with the sample median for group pqt. One complication in estimating the gradient matrix D is that the functions inside the moment condition ϕ(∆; χ) are non-differentiable at certain data points; see equation (24). This means that we cannot consistently estimate D as the numerical derivative of ϕ̂I (.). Our asymptotic results therefore do not follow from the standard GMM approach, but rather the approach for non-smooth functions described in Pakes and Pollard [46], Newey and McFadden [45] (section 7), and Powell [50]. To find D, it is helpful to rewrite equation (24), for the case of medians, as: Pr fit = f & ci = c & Ii ∈ Ip & individual i observed at t × # "Z median af cpt (∆0 ,χ0 ) 1 = 0. (27) f ait f, c, Ii ∈ Ip , t dait − 2 −∞ It follows that the rows of D for the median asset moments are given by Pr fit = f & ci = c & Ii ∈ Iq & individual i observed at t × ∂amedian (∆0 ; χ0 ) f cpt median f af cpt f, c, Ii ∈ Ip , t × . (28) ∂∆0 In practice, we find f amedian f cpt |f, c, p, t , the conditional p.d.f. of wealth evaluated at the median amedian f cpt , with a kernel density estimator written by Koning [37]. The derivatives for the 25th and 75th percentiles are found in analogous fashion. Appendix G: First stage estimates In the first-step estimation procedure, we estimate or calibrate those parameters that, given our assumptions, can be cleanly identified outside our model. In particular, we estimate health transitions and mortality rates from raw demographic data, estimate income profiles from the AHEAD, and calibrate the returns to saving and income taxes using external sources. We use our measure of total medical expenses (which combines data from the AHEAD and the MCBS) to estimate the process for medical expenses as a function of exogenous household-level state variables which are observed in the AHEAD. We outline the estimation and calibration of these objects next. 54 G.1 Health and mortality Each person’s health status, hsg , has four possible values: dead; in a nursing home; in bad health; or in good health. We allow transition probabilities for an individual’s health depend on his or her current health, age, household composition f , permanent income I, and gender g.17 It follows that the elements of the health transition matrix are given by πi,j,k (t, fit , Ii , g = Pr hsgi,t+2 = k hsgi,t = j; t, fi,t , Ii , g , (29) with the transitions covering a two-year interval, as the AHEAD interviews every other year.18 We estimate health/mortality transition probabilities by fitting the transitions observed in the AHEAD to a multinomial logit model.19 We enumerate each possible health status as follows 0 Dead 1 Nursing Home (30) hs = 2 Bad Health 3 Good Health. Our multinomial logit assumption gives the following expression for health and mortality transitions exit βk πi,j,k (t, fit , Ii , g = P3 , (31) xit βh h=0 e with βk denoting the coefficient vector for future outcome k. The coefficient for death, β0 , is normalized to zero. We jointly estimate health transitions and survival 17 We do not allow health transitions to depend on medical spending. The empirical evidence on whether medical spending improves health, especially at older ages, is surprisingly mixed; see, for example, the discussion in De Nardi, French and Jones [19]. Likely culprits include reverse causality – sick people have higher expenditures – and a lack of insurance variation – almost every retiree gets Medicare. 18 As discussed in De Nardi, French and Jones [19], one can fit annual models of health and medical spending to the AHEAD data. The process becomes significantly more involved, however, especially when accounting for the dynamics of two-person households. 19 We do not control for cohort effects. Instead, our estimates are a combination of period (crosssectional) and cohort probabilities. While our AHEAD sample covers 18 years, it is still too short to track a single cohort over its entire post-retirement lifespan. This may lead us to underestimate the lifespans expected by younger cohorts as they age. Nevertheless, lifespans have increased only modestly over the sample period. Accounting for cohort effects would have at most a modest effect on our estimates. 55 probabilities at the individual level using a maximum likelihood estimator. We allow the transition probabilities to depend on xit which includes age, sex, current health status, marital status, and permanent income. In particular we include a third order age polynomial, indicators for gender and marital status (interacted with an age quadratic), single woman (interacted with a first order age polynomial), contemporaneous indicators for health (interacted with age), and a quadratic in permanent income (interacted with a first order polynomial in age and marital status). Using the estimated transition probabilities, we simulate demographic histories, beginning at age 70, for different gender-PI-health combinations. We report moments of these histories in Section 6.1.1. G.2 Income Our procedure for estimating the income profiles described in section 6.1.2 is part of the procedure for inferring permanent income described in Appendix B. We set the annual rate of return on saving to 4%. G.3 The tax system We parameterize the tax system using the functional form developed by Gouveia and Strauss [29]. Average tax rates for total pre-tax income, y, are given by −1 t(y, τft ) = bft [1 − (sft y pft + 1) pft ], (32) where bft , sft and pft are parameters of the tax system that can differ between couples and singles, and τft = bft , sft , pft is their union. We use the estimates of τ for married and singles households with no children provided by Guner et al. [30]. After-tax income Υ(y, τft ) is then Υ(y, ft ) = (1 − t(y, τft )) × y. (33) To reflect differences in the generosity of Medicaid, and other means tested social insurance programs available to the elderly, between couples and singles we impose that the consumption floor for couples is 150% of the value we estimate for singles. G.4 Medical spending Our medical spending measure is the sum of expenditures paid out-of-pocket plus those paid by Medicaid (see Appendix D for its construction). The AHEAD’s 56 measure of medical spending is backward-looking: in each wave, the household is asked about the expenses it incurred since the previous interview. Let mi,t denote the expenses incurred by household i between ages t − 1 and t. We observe household’s health at the beginning and the end of this interval, that is, at the time of the interview conducted at age t − 1 and at the time of the interview conducted at age t. We assume that medical expenses depend upon a household’s PI, its family structure at both t − 1 and t, the health of its members at both dates, and an idiosyncratic component ψi,t : h w ln mi,t = m(hshi,t−1 , hsw i,t−1 , hsi,t , hsi,t , fi,t−1 , fi,t , Ii , t) + υi,t , (34) h w υi,t = σ(hshi,t−1 , hsw i,t−1 , hsi,t , hsi,t , fi,t−1 , fi,t , Ii , t) × ψi,t . (35) We normalize the variance of ψ to 1, hence σ 2 (·) gives the conditional variance of υ. Including current and lagged family structure indicators allows us to account for the jump in medical spending that occurs in the period a family member dies. Likewise, including health indicators for both periods allows us to distinguish persistent health episodes from transitory ones. To show how we estimate m(·) and σ(·), write Equation (6) as ln mit = x1i β1 + x2it β2 + ϑi + ςit , (36) where x1i denotes a vector of time-invariant variables, x2it denotes a vector of timevarying variables, ϑi is an unobserved person-specific term, and ςit captures any remaining variation. We assume that E(ςit | ϑi ) = 0. We estimate Equation (36) in three steps. First, we regress log medical spending on the time-varying factors in Equation (34), namely age, household structure, and health, and interaction terms (such as gender and PI interacted with the time varying variables) using a fixed effects estimator. In particular we regress log medical spending on a fourth order age polynomial, indicators for single man (interacted with an age quadratic), single woman (interacted with an age polynomial), the contemporaneous and lagged values of indicators for {man in bad health, married man in a nursing home, single man in a nursing home, woman in bad health, married woman in a nursing home, single woman in a nursing home}, whether the man died (interacted with age and permanent income), whether the woman died (interacted with age, and permanent income). Because fixed effects regression cannot identify the effects of time-invariant factors, which are combined into the estimated fixed effects, in the second step we collect the residuals from the first regression, inclusive of the estimated fixed effects, and regress them on the time-invariant factors, namely a quadratic in permanent 57 income and a set of cohort dummies. The level of m(·) is set to be consistent with the outcomes of the cohort aged 72-76 in 1996. A key feature of our spending model is that the conditional variance and the conditional mean of medical spending depends on demographic and socioeconomic factors, through the function σ(·) shown in Equation (35). In the third step of our estimation procedure, we use the coefficients for m(·) in hand, found in the previous 2 (·), we square two steps, to back out the residual υ from Equation (6). To find σ[ 2 the residuals and regress υ on the demographic and socioeconomic variables in Equation (35). We assume that ψi,t can be decomposed as ψi,t = ζi,t + ξi,t , ξi,t ∼ N (0, σξ2 ), ζi,t = ρm ζi,t−2 + i,t , i,t ∼ N (0, σ2 ), (37) (38) where ξi,t and i,t are serially and mutually independent. With the variance of ψi,t normalized to 1, σξ2 can be interpreted as the fraction of idiosyncratic variance due to transitory shocks. We estimate the parameters of Equations (7) and (8) using a standard error components method. Although the estimation procedure makes no assumptions on the distribution of the error terms ψi,t , we assume normality in the simulations. [27] shows that if the data are carefully constructed, normality captures well the far right tail of the medical spending distribution.20 Approximately 40% of the cross-sectional variation in log medical spending is explained by observables, which are quite persistent. Of the remaining cross-sectional variation, 40% comes from the persistent shock ζ and 60% from the transitory shock ξ. In keeping with the results in Feenberg and Skinner [25], French and Jones [27] and De Nardi, French and Jones [18], we estimate substantial persistence in the persistent component, with ρm = 0.85. 20 To help us match the distribution of medical spending, we bottom code medical spending at 10% of average medical spending. French and Jones [27] also bottom code the data to match the far right tail of medical spending. Because we include Medicare B payments in our medical spending measure, which most elderly households pay, for the vast majority of households these bottom coding decisions are not important. 58 Appendix H: The marginal propensity to consume Consider an individual with no income and no medical expenses. In the last period of life, this individual solves 1−ν 1 1 1−ν c + βφ0 (1 + r)(a − c) + κ0 , max c 1−ν 1−ν where c denotes consumption and a denotes beginning-of-period-wealth. The firstorder condition for an interior solution is −ν c−ν = βφ0 (1 + r) (1 + r)(a − c) + κ0 , Define ϕ = [βφ0 (1 + r)]1/ν . (39) We can then solve for bequests, b = (1 + r)(a − c) : 1+r b (ϕa − κ0 ) . = b + κ0 ⇒ b = ϕ a− 1+r 1+r+ϕ (40) It immediately follows that for intentional bequests, the marginal propensity to b ∂ bequeath out of wealth, ∂a , is 1+r MP B = ϕ . 1+r+ϕ (41) With non-negative bequests, however, it follows from Equation (40) that the propensity to bequeath is identically zero unless κ0 a>a= . (42) ϕ It is widely accepted that M P B (or M P C = 1 − M P B) and a provide a better characterization of the bequest motive than φ0 and κ0 (see, e.g., De Nardi et al., 2010 or Lockwood, 2018). When one expects to die at a future date rather than next period, the effective threshold where the bequest motive becomes operative will be higher than a. There is a one-to-one mapping between the vectors (M P B, a) and (φ0 , κ0 ), and thus if one of these objects can be identified, the other one can be as well. Using equations (42), we know κ0 = ϕa. It then follows from equations (39) and (41) that: ϕ = (1 + r) MP B 1−M PB ϕν MP B ⇒ φ0 = = (1 + r) β(1 + r) 1 − MP B 59 ν 1 . β(1 + r) We also need to map parameters between models with one-year and two-year period lengths. An important intermediate step in this process is to consider the sequence of 1-year problems where the individual lives two periods and then dies. We will use hats to designate parameters for this one-year model, and unmarked b + parameters for the two-year framework. To simplify the algebra, suppose that β(1 r̂) = 1, so that b c is the same in the two periods of life. The present value budget constraint is b̂ 2 + r̂ + , a=b c 1 + r̂ (1 + r̂)2 which can be rewritten as b c= 1 + r̂ 2 + r̂ b̂ a− . (1 + r̂)2 (43) ∂ b̂ \ If the marginal propensity to bequeath out of starting wealth is M P B = ∂a , (1+r̂)2 the marginal propensity to consume on a annual basis, found by differentiating equation (43), is: i 1 + r̂ h ∂ĉ \ \ = 1−M PB . (44) M PC = ∂a 2 + r̂ Next, consider the asset threshold b a. Recall that in the final period of life, wealth below b a are consumed. Let e a denote the asset threshold two periods prior to death. By definition, when wealth two periods prior to death equal e a, b̂ = 0. Equation (43) then implies that the asset threshold one period prior to death is 1 + r̂ 1 + r̂ e a = e a. b a = (1 + r̂) e a− (45) 2 + r̂ 2 + r̂ \ To complete the mapping, suppose that 1 + r = (1 + r̂)2 , M P B = M P B, and a=e a. Equations (44) and (45) then imply that 1 + r̂ \ M PC = M P C, 2 + r̂ 1 + r̂ b a= a. 2 + r̂ (46) (47) In other words, when we operate at a one-year frequency, the marginal propensity to consume and the asset threshold are (more or less) half their two-year counterparts. 60 Appendix I: Additional model fits Figure 17 compares the Median asset profiles generated by the model (dashed line) to those in the data (solid line) for all of the birth cohorts and permanent income terciles that we match in estimation. We reproduce profiles for the cohorts reported in the main text, but also show profiles for the cohorts that were omitted. In all cohorts, the model matches differences in the level and trajectory of wealth by age and household permanent income. Median Net Worth: Data (Solid) vs. Model (Dashed), Couples 550 550 500 500 450 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) Median Net Worth: Data (Solid) vs. Model (Dashed), Couples 600 450 400 350 300 250 200 150 100 350 300 250 200 150 100 50 50 0 72 400 76 80 84 88 92 96 0 72 100 76 80 84 Age 92 96 100 Median Net Worth: Data (Solid) vs. Model (Dashed), Singles 200 240 180 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) Median Net Worth: Data (Solid) vs. Model (Dashed), Singles 270 210 180 150 120 90 60 30 0 72 88 Age 160 140 120 100 80 60 40 20 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age Figure 17: Median wealth. Top panel: current couples. Bottom panel current singles: Solid lines: cohorts aged 72-77, 84-89 in 1996 in the left-hand panel and ages 78-83 and 90+ in 1996 in the right-hand side panel. Dashed lines: model simulations. Figure 18 plots the corresponding data for the 75th percentile of the asset distribution. Comparing the 75th wealth percentiles to the median ones towards the beginning of the retirement period, when we start tracking our sample, shows that 61 75th %tile Net Worth: Data (Solid) vs. Model (Dashed), Couples 75th %tile Net Worth: Data (Solid) vs. Model (Dashed), Couples 1400 1100 1000 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 1200 1000 800 600 400 200 900 800 700 600 500 400 300 200 100 0 72 76 80 84 88 92 96 0 72 100 76 80 84 Age 88 92 96 100 Age 75th %tile Net Worth: Data (Solid) vs. Model (Dashed), Singles 75th %tile Net Worth: Data (Solid) vs. Model (Dashed), Singles 700 600 550 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 600 500 400 300 200 100 500 450 400 350 300 250 200 150 100 50 0 72 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age Figure 18: 75th wealth percentile. Top panel: current couples. Bottom panel current singles: Solid lines: cohorts aged 72-77, 84-89 in 1996 in the left-hand panel and ages 78-83 and 90+ in 1996 in the right-hand side panel. Dashed lines: model simulations. the 75th asset percentiles are about twice the medians for the two highest PI groups and much larger than twice the median for the lowest PI tercile. In addition, the 75th asset percentiles show even more evidence of wealth accumulation during retirement than do median wealth, both for couples and singles. Finally, the 75th percentiles suggest that current couples tend to save more during their retirement than singles, even given similar initial wealth. As with the medians, the model does a good job of matching the 75th percentile of asset holdings across PI groups and ages, both for singles and couples. For example, it reproduces the flatter profiles of lower PI groups, while still generating eventual wealth holdings of more than $1,000,000 for the top terciles. Figure 19 plots the corresponding data for the 25th percentile of the asset distribution. Comparing the 25th wealth percentiles to the medians towards the beginning 62 25th %tile Net Worth: Data (Solid) vs. Model (Dashed), Couples 25th %tile Net Worth: Data (Solid) vs. Model (Dashed), Couples 300 270 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 240 250 200 150 100 50 210 180 150 120 90 60 30 0 72 76 80 84 88 92 96 0 72 100 76 80 84 Age 25th %tile Net Worth: Data (Solid) vs. Model (Dashed), Singles 96 100 90 100 80 90 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) 92 25th %tile Net Worth: Data (Solid) vs. Model (Dashed), Singles 110 80 70 60 50 40 30 20 70 60 50 40 30 20 10 10 0 72 88 Age 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age Figure 19: 25th wealth percentile. Top panel: current couples. Bottom panel current singles: Solid lines: cohorts aged 72-77, 84-89 in 1996 in the left-hand panel and ages 78-83 and 90+ in 1996 in the right-hand side panel. Dashed lines: model simulations. of the retirement period, when we start tracking our sample, shows that the 25th asset percentiles are about half the median ones for the two highest PI groups, with the exception of the lowest PI tercile for singles, where it is almost always 0. In addition, the 25th asset percentiles show even more evidence that couples and singles’ wealth evolves differently as the wealth holdings of singles decline. As with the 75th percentile, it seems that current couples tend to save more during their retirement than singles, even given similar initial wealth. Our model, in addition to matching the median and the 75th percentile of asset holdings by age, PI, and cohort, also does a good job of generating the differential deccumulation we see for the 25th percentile. Finally, Figure 20 compares the Medicaid recipiency profiles generated by the model (dashed line) to those in the data (solid line) for the same birth cohorts and permanent income terciles and shows that our model matches important patterns of 63 Medicaid Recipiency: Data (Solid) vs. Model (Dashed), All Medicaid Recipiency: Data (Solid) vs. Model (Dashed), All 0.8 0.65 0.6 0.7 0.55 0.5 Medicaid Recipiency Medicaid Recipiency 0.6 0.5 0.4 0.3 0.2 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.1 0.05 0 72 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) (b) Figure 20: Medicaid Recipiency. Solid lines: cohorts aged 72-77, 84-89 in 1996 in panel (a) and ages 78-83 and 90+ in 1996 in panel (b). Dashed lines: model simulations. Medicaid usage. We reproduce the cohorts reported in the main text and add the cohorts not shown there. The model matches difference in the level of usage by age and household permanent income for all cohorts. It does especially well in capturing the run-up in Medicaid use by wealthier households at very old ages. Because individuals who are single, poor or sick die at younger ages, in an unbalanced panel the sample composition changes over time, introducing mortality bias. Figure 21 compares the mortality bias found in the AHEAD data (left panel) with that predicted by the model (right panel). The solid lines show median wealth for all households (both singles and couples) observed at a given point in time, even if all members died in a subsequent wave, that is, the unbalanced panel. The dashed lines show median wealth for the subsample of households with at least one member still alive in the final wave, that is, the balanced panel. They show that the asset profiles for households who survived to the final wave (the balanced panel) have much more of a downward slope. The difference between the two sets of profiles confirms that people who died during our sample period tended to have lower wealth than the survivors. Our model matches the extent of the mortality bias found in the data, without being required to do so by construction. As in the data, restricting the profiles to long-term survivors reveals much more asset decumulation. 64 Median Net Worth: Everyone in Sims (Solid) vs. Survivors (Dashed), All 400 400 350 Net Worth (000s of 2014 dollars) Net Worth (000s of 2014 dollars) Median Net Worth: Everyone in Data (Solid) vs. Survivors (Dashed), All 450 350 300 250 200 150 100 250 200 150 100 50 50 0 72 300 76 80 84 88 92 96 100 Age 0 72 76 80 84 88 92 96 100 Age (a) (b) Figure 21: Mortality Bias. Left panel: data. Right panel: model simulations. Solid lines: all household in the data at a given point in time. Dashed lines: households with at least one member surviving to 2014. 65