Modeling Multiple Time-Varying Related Groups: A Dynamic Hierarchical Bayesian Approach With an Application to the Health and Retirement Study

Abstract As the population of the older individuals continues to grow, it is important to study the relationship among the variables measuring financial health and physical health of the older individuals to better understand the demand for healthcare, and health insurance. We propose a semiparametric approach to jointly model these variables. We use data from the Health and Retirement Study which includes a set of correlated longitudinal variables measuring financial and physical health. In particular, we propose a dynamic hierarchical matrix stick-breaking process prior for some of the model parameters to account for the time dependent aspects of our data. This prior introduces dependence among the parameters across different groups which varies over time. A Lasso type shrinkage prior is specified for the covariates with time-invariant effects for selecting the set of covariates with significant effects on the outcomes. Through joint modeling, we are able to study the physical health of the older individuals conditional on their financial health, and vice-versa. Based on our analysis, we find that the health insurance (medicare) provided by the government (of the United States) to the older individuals is very effective, and it covers most of the medical expenditures. However, none of the health insurances conveniently cover the additional medical expenses due to chronic diseases like cancer and heart problem. Simulation studies are performed to assess the operating characteristics of our proposed modeling approach. Supplementary materials for this article, including a standardized description of the materials available for reproducing the work, are available as an online supplement.

[1]  Kiranmoy Das,et al.  A Bayesian approach of analysing semi-continuous longitudinal data with monotone missingness , 2020, Statistical Modelling.

[2]  Scott Lindroth,et al.  Dynamic Nonparametric Bayesian Models for Analysis of Music , 2010 .

[3]  John D. Lafferty,et al.  Dynamic topic models , 2006, ICML.

[4]  A. L. Pretorius,et al.  Bayesian estimation in animal breeding using the Dirichlet process prior for correlated random effects , 2003, Genetics Selection Evolution.

[5]  Partha Deb,et al.  Are primary care physicians, public and private sector specialists substitutes or complements? Evidence from a simultaneous equations model for count data. , 2008, Journal of health economics.

[6]  Arnaud Doucet,et al.  Generalized Polya Urn for Time-varying Dirichlet Process Mixtures , 2007, UAI.

[7]  Jayaran Sethuramant A CONSTRUCTIVE DEFINITION OF DIRICHLET PRIORS , 1991 .

[8]  Sarah Brown,et al.  Debt and distress: Evaluating the psychological cost of credit , 2005 .

[9]  G. Casella,et al.  The Bayesian Lasso , 2008 .

[10]  A. Kapteyn,et al.  Health, Wealth, and the Role of Institutions , 2003, The Journal of Human Resources.

[11]  Peter E. Rossi,et al.  A Bayesian analysis of the multinomial probit model with fully identified parameters , 2000 .

[12]  James Waterman Glover United States Life Tables , 2013 .

[13]  P. Deb,et al.  Demand for Medical Care by the Elderly: A Finite Mixture Approach , 1997 .

[14]  T. Ferguson A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[15]  D. Dunson,et al.  Bayesian Selection and Clustering of Polymorphisms in Functionally Related Genes , 2008 .

[16]  Kiranmoy Das,et al.  A semiparametric approach to simultaneous covariance estimation for bivariate sparse longitudinal data , 2014, Biometrics.

[17]  D. Dunson Bayesian dynamic modeling of latent trait distributions. , 2006, Biostatistics.

[18]  C. Antoniak Mixtures of Dirichlet Processes with Applications to Bayesian Nonparametric Problems , 1974 .

[19]  A. Soest,et al.  Health and Wealth of Elderly Couples: Causality Tests Using Dynamic Panel Data Models , 2004, Journal of health economics.

[20]  L. Carin,et al.  The Matrix Stick-Breaking Process , 2008 .

[21]  Michael I. Jordan,et al.  Hierarchical Dirichlet Processes , 2006 .

[22]  M. Daniels,et al.  A Nonparametric Prior for Simultaneous Covariance Estimation. , 2013, Biometrika.

[23]  Husˇkova´ Marie,et al.  On sequentially adaptive asymptotically efficient rank statistics , 1985 .

[24]  Michael J Daniels,et al.  A Flexible Bayesian Approach to Monotone Missing Data in Longitudinal Studies With Nonignorable Missingness With Application to an Acute Schizophrenia Clinical Trial , 2015, Journal of the American Statistical Association.

[25]  Shane T. Jensen,et al.  Semiparametric Bayesian Modeling of Income Volatility Heterogeneity , 2011 .

[26]  Jim E. Griffin,et al.  Stick-breaking autoregressive processes , 2011 .

[27]  Doug Brown,et al.  Primary Care Physicians , 2020, Encyclopedia of Behavioral Medicine.

[28]  J. Sethuraman A CONSTRUCTIVE DEFINITION OF DIRICHLET PRIORS , 1991 .

[29]  H. Schmitz,et al.  Broke, Ill, and Obese: Is There an Effect of Household Debt on Health? , 2014 .

[30]  Michael I. Jordan,et al.  Variational inference for Dirichlet process mixtures , 2006 .

[31]  Kiranmoy Das,et al.  A Bayesian two-stage regression approach of analysing longitudinal outcomes with endogeneity and incompleteness , 2019 .

[32]  Pulak Ghosh,et al.  Estimating Healthcare Demand for an Aging Population: A Flexible and Robust Bayesian Joint Model , 2011 .

[33]  Lancelot F. James,et al.  Gibbs Sampling Methods for Stick-Breaking Priors , 2001 .

[34]  Karin Meyer,et al.  Random regressions to model phenotypic variation in monthly weights of Australian beef cows , 2000 .

[35]  D. V. Dyk,et al.  A Bayesian analysis of the multinomial probit model using marginal data augmentation , 2005 .

[36]  N M Laird,et al.  Missing data in longitudinal studies. , 1988, Statistics in medicine.