A Bayesian joint model for zero‐inflated integers and left‐truncated event times with a time‐varying association: Applications to senior health care

• Published in: Statistics in medicine Vol. 40; no. 1; pp. 147 - 166
• Main Authors: Piulachs, Xavier, Andrinopoulou, Eleni‐Rosalina, Guillén, Montserrat, Rizopoulos, Dimitris
• Format: Journal Article
• Language: English
• Published: England Wiley Subscription Services, Inc 15.01.2021
• Subjects: age‐related mortality hazard; Emergency medical care; left truncation; morbidity; Older people; Parameter estimation; private health care; shared‐parameter joint model; time‐varying association; zero inflation; age-related mortality hazard; private health care; morbidity; left truncation; time-varying association; shared-parameter joint model; zero inflation
• ISSN: 0277-6715; 1097-0258; 1097-0258
• DOI: 10.1002/sim.8767

Population aging in most industrialized societies has led to a dramatic increase in emergency medical demand among the elderly. In the context of private health care, an optimal allocation of the medical resources for seniors is commonly done by forecasting their life spans. Accounting for each subject's particularities is therefore indispensable, so the available data must be processed at an individual level. We use a large and unique dataset of insured parties aged 65 and older to appropriately relate the emergency care usage with mortality risk. Longitudinal and time‐to‐event processes are jointly modeled, and their underlying relationship can therefore be assessed. Such an application, however, requires some special features to also be considered. First, longitudinal demand for emergency services exhibits a nonnegative integer response with an excess of zeros due to the very nature of the data. These subject‐specific responses are handled by a zero‐inflated version of the hierarchical negative binomial model. Second, event times must account for the left truncation derived from the fact that policyholders must reach the age of 65 before they may begin to be observed. Consequently, a delayed entry bias arises for those individuals entering the study after this age threshold. Third, and as the main challenge of our analysis, the association parameter between both processes is expected to be age‐dependent, with an unspecified association structure. This is well‐approximated through a flexible functional specification provided by penalized B‐splines. The parameter estimation of the joint model is derived under a Bayesian scheme.