- 1Department of Human Development and Family Science, Aging Studies Institute, Center for Aging and Policy Studies, and the Lerner Center for Public Health Promotion and Population Health, Syracuse University, Syracuse, NY, United States
- 2Department of Sociology, Aging Studies Institute, Center for Aging and Policy Studies, and the Lerner Center for Public Health Promotion and Population Health, Syracuse University, Syracuse, NY, United States
- 3Center for Demography and Ecology, and Center for Demography of Health and Aging, University of Wisconsin-Madison, Madison, WI, United States
- 4University of Alabama at Birmingham, Birmingham, AL, United States
Background: Recent efforts have been made to collect data on neighborhood-level attributes and link them to longitudinal population-based surveys. These linked data have allowed researchers to assess the influence of neighborhood characteristics on the health of older adults in the US. However, these data exclude Puerto Rico. Because of significantly differing historical and political contexts, and widely ranging structural factors between the island and the mainland, it may not be appropriate to apply current knowledge on neighborhood health effects based on studies conducted in the US to Puerto Rico. Thus, we aim to (1) examine the types of neighborhood environments older Puerto Rican adults reside in and (2) explore the association between neighborhood environments and all-cause mortality.
Methods: We linked data from the 2000 US Census to the longitudinal Puerto Rican Elderly Health Conditions Project (PREHCO) with mortality follow-up through 2021 to examine the effects of the baseline neighborhood environment on all-cause mortality among 3,469 participants. Latent profile analysis, a model-based clustering technique, classified Puerto Rican neighborhoods based on 19 census block group indicators related to the neighborhood constructs of socioeconomic status, household composition, minority status, and housing and transportation. The associations between the latent classes and all-cause mortality were assessed using multilevel mixed-effects parametric survival models with a Weibull distribution.
Results: A five-class model was fit on 2,477 census block groups in Puerto Rico with varying patterns of social (dis)advantage. Our results show that older adults residing in neighborhoods classified as Urban High Deprivation and Urban High-Moderate Deprivation in Puerto Rico were at higher risk of death over the 19-year study period relative to the Urban Low Deprivation cluster, controlling for individual-level covariates.
Conclusions: Considering Puerto Rico's socio-structural reality, we recommend that policymakers, healthcare providers, and leaders across industries to (1) understand how individual health and mortality is embedded within larger social, cultural, structural, and historical contexts, and (2) make concerted efforts to reach out to residents living in disadvantaged community contexts to understand better what they need to successfully age in place in Puerto Rico.
1. Introduction
The twenty-first century in the Commonwealth of Puerto Rico (hereafter, Puerto Rico)—an unincorporated United States (US) territory—is an era characterized by rapid population aging, reductions in social and economic resources, rampant disparities in access to adequate healthcare, and the ongoing reconstruction of the built environment post-Hurricanes Irma and María (1–6). The constellation of these factors infers that many older Puerto Rican adults may lack access to resources, services, and contexts considered necessary for promoting healthy aging.1 In order to understand contemporary conditions in Puerto Rico, it is important to consider how historical contexts contribute to health inequities over time, particularly for older adults at increased risk for poor health and mortality.
Researchers have argued that social, political, and economic inequalities in Puerto Rico derive from the impacts of US colonialism—a structural and social determinant of health (7, 8). One significant impact of US colonialism was the transition of Puerto Rico from a rural agricultural society to an urban industrial society in the early twentieth century (9). This transition brought public health benefits, including improved sanitation practices and housing conditions, the creation of local health boards and hospitals, and increased access to primary education. However, urbanization in Puerto Rico also led to widening economic and racial disparities that resulted in unfavorable neighborhood and living conditions among socially marginalized individuals (e.g., poor and Black Puerto Ricans) (10).
For example, San Juan, the capital of Puerto Rico, has a long history of continuous urban growth and economic development. Under US control, San Juan experienced substantial modernization, including changes in land use efficiency and aggregation of local areas that connected land use with global-scale factors. Notably, new and growing opportunities in the San Juan wage labor market were a major driver for rural-dwelling Puerto Ricans to relocate in the early twentieth century; this rural-to-urban migration affected the subsequent development and preservation of several neighborhoods in the metropolitan area (11). Due to their lower socioeconomic position, Puerto Ricans from rural areas were forced to reside in poor and disadvantaged, communities in San Juan, such as La Perla (12). In addition to rural-urban migration patterns, rapid population growth and efforts to mirror the US model of suburbanization were additional factors that influenced variations in the investment of resources across neighborhoods in San Juan throughout the twentieth century that contributed to contemporary residential segregation patterns (13). For example, a study examining residential segregation in the San Juan-Bayamón metropolitan area, the most racially diverse metropolitan area in Puerto Rico, found that neighborhoods with a higher percentage of Black residents were associated with lower socioeconomic status (14).
In addition, a study focusing on the socioeconomic features of neighborhoods to assess health disparities in Puerto Rico found that municipalities (considered county-equivalents by the US Census) with a low socioeconomic position (SEP) were linked to higher cancer-related mortality rates (15). Importantly, the study showed that more deprived municipalities of Puerto Rico were in the island's central region.2 In contrast, less deprived municipalities were concentrated in the San Juan metropolitan area. This suggests that residents living in municipalities with lower SEPs may lack access to healthcare services and health-promoting resources due to economic, environmental, and physical barriers that impact health and increase the risk of mortality. However, these findings are conditional based on the assumptions made regarding area-based socioeconomic status. Better inference of neighborhood effects would require a more nuanced approach on how specific constructs of the neighborhood environment are measured (e.g., census tract vs. census block group) and how they influence health and the risk of mortality (16).
Recent efforts have been made to collect data on neighborhood-level attributes and link them to longitudinal population-based surveys [e.g., the Health and Retirement Study Contextual Data Resource (HRS-CDR)] (17). These linked data have allowed researchers to assess the influence of neighborhood characteristics on the health of older adults in the US. However, these data only include the contiguous US and exclude Alaska, Hawai'i, and the five permanently inhabited US territories, including Puerto Rico. Because of significantly differing historical and political contexts, and widely ranging structural factors between Puerto Rico and the US mainland, it is not appropriate to apply current knowledge on neighborhood health effects based on studies conducted in the contiguous US to Puerto Rico. In addition, despite Puerto Rico's status as an unincorporated US territory, its social and economic contexts are more like Latin American and Hispanic-Caribbean countries than the US, which may lead to substantially different risk factors for poor health and mortality.
In this study, we aim to highlight multilevel perspectives and analyses of social determinants of health among older adults residing in Puerto Rico. We address a gap in the literature by using longitudinal data from the Puerto Rican Elderly Health Conditions Project (PREHCO) linked with 2000 US Census data to (1) examine the types of neighborhood environments older Puerto Rican adults reside in and, (2) explore the association between neighborhood environment and all-cause mortality.
2. Background
It is widely recognized that physical and social environments influence health behaviors, health outcomes, and mortality in the US. Although the neighborhood environment affects the health of people of all ages, the effects of the neighborhood environment may be accentuated among older adults as they are more likely than younger adults to have spent decades in the same community, have decreased physical mobility and cognitive functioning, and rely more on community resources for social integration and support (18). The combination of these factors may result in an early onset of age-related diseases (19), reduced life expectancy (20), and an increased risk of all-cause mortality (21–23). Notably, a vast array of research has shown that individuals residing in neighborhoods with greater deprivation have poorer health behaviors (24), lack access to preventive health services (25), are exposed to chronic stress and pollutants (26), experience greater biological weathering (27), have worse health outcomes (28), and experience higher mortality rates (29). In many of these studies, neighborhood deprivation is based on socioeconomic contextual variables or indices related to income, education, employment, and housing, typically at the census tract level. Although these socioeconomic indicators have different meanings for older adults, it is noteworthy that the influence of socioeconomic deprivation persists in the oldest ages (30). Indeed, several studies suggest a cumulative effect of disadvantage across the lifespan that results in poor health and an increased risk of mortality (31). However, there is limited knowledge of how these multilevel processes influence population health and mortality in Puerto Rico due to the lack of data infrastructure to support these inquiries.
2.1. Neighborhood socioeconomic context
Various theoretical perspectives and conceptual frameworks have been put forth to explain why the neighborhood socioeconomic context (NSEC) plays a vital role in poor health outcomes and mortality risk. For instance, the ecological framework with a life course perspective would suggest that individuals living in disadvantaged NSECs are more likely to have a low socioeconomic position themselves due to constrained opportunity structures (22, 32, 33). Individuals who spend their early life in lower-income neighborhoods have less access to quality education than their peers residing in higher-income communities. This limits opportunities to obtain higher levels of education and marketable job skills and reduces lifetime earnings (34, 35). Thus, the importance of neighborhood context as a fundamental cause of mortality cannot be overlooked (36), particularly given the vast literature documenting how education shapes access to resources that promote better health and an individual's exposure to multiple health risks (37).
Another theoretical consideration is the systemic perspective, which infers that the NSEC affects the social, service, and physical environments of communities shared by residents. Namely, neighborhoods characterized by low socioeconomic levels are linked to underinvestment in health-promoting resources, such as lack of green and recreational spaces, adequate public transportation, affordable and high-quality grocery stores, and access to medical and social services (23). For example, individuals residing in high-poverty neighborhoods are less likely to have access to recreational opportunities to walk and exercise and are more likely to live in food swamps3 (38, 39). Not being able to engage in healthy behaviors due to these structural challenges can increase the likelihood of early disease onset, reduce active life expectancy (e.g., physical mobility), and increase the risk of mortality. Overall, the emphasis of the NSEC on health is important from a public health perspective since resource-poor environments can be potentially addressed through community-level interventions, including investments in public education, transportation, expansion of door-to-door services (e.g., Meals on Wheels), and affordable and quality housing to name a few.
Although research has overwhelmingly demonstrated that the NSEC is a crucial determinant of health, other neighborhood-related factors interplay with the NSEC, such as a neighborhood's age structure, racial composition, residential stability, and family structure that shape opportunities and health-enhancing resources made available for residents across communities. We provide a summary of how each of these neighborhood-level determinants potentially influences health outcomes and the risk of mortality.
2.2. Neighborhood age structure
The age structure of a neighborhood may be particularly important to older adults who age in place as it may influence the provision of health services and facilities (including Medicaid reimbursements), perceptions of neighborhood safety, and opportunities for social engagement (40). Previous research has shown that neighborhoods with a high concentration of older adults are associated with better health among older adults, including those who are socioeconomically disadvantaged (41, 42). Evidence suggests that the presence of older adults in the community facilitates social integration and cultivates social ties, mutual support, social cohesion, and perceived safety (43), which is independently associated with various population-level health outcomes, including mortality (44). Several pathways have been hypothesized on how aspects of the social environment may influence health and mortality, including the impact of health behaviors and physiology (e.g., allostatic load) (44, 45). Specifically, individuals with positive social ties are less likely to engage in smoking and drinking and are more likely to receive preventive health screenings (e.g., cancer screenings). In contrast, socially isolated individuals are more likely to have weakened immune function, cardiovascular disease, and cognitive impairment. Older adults with chronic health conditions, disabilities, who live alone, and have reduced social networks are at an increased risk of social isolation, which has been shown to negatively impact health and mortality.
With the population of Puerto Rico is rapidly aging—due to a combination of outmigration among younger cohorts of adults, declining fertility, and increased longevity—these demographic changes will challenge the ability of Puerto Rico and local communities to meet the growing demands of older adults, including care and quality of life, that may further strain the collective (and scarce) resources available (1, 3, 5, 46). Specifically, increases in poverty and declining economic conditions across the archipelago, changes in the family structure, and the limited availability (and proximity) of individuals and/or services to provide long-term care for older adults in Puerto Rico (due in part to out migration of family and professionals) may result in poor health and an increased risk of mortality. Older Puerto Ricans, cognizant of these social realities have expressed concerns with loss of family cohesion and intergenerational support due to their children's search of economic opportunities outside of Puerto Rico (47). This suggests that places in Puerto Rico that have a larger concentration of older adults, particularly in rural areas, may not have the resources necessary for older adults to successfully age in place.
2.3. Neighborhood racial composition
Neighborhood racial composition has been shown to be associated with poor health and an increased risk of mortality among older adults due in part to exposure of institutionalized and systemic anti-Black racism across the life course (48–50). A large body of research shows that Black (including African American and Afro-Latino) individuals in the US overwhelmingly reside in residentially segregated neighborhoods that are characterized by concentrated economic disadvantage, which is often associated with disinvestment of municipal resources (e.g., high-quality medical care), poorly maintained infrastructures (e.g., sidewalks and green spaces), and densely populated and subpar housing quality (51–53). These conditions stem from racial capitalism and environmental racism that intentionally create the underdevelopment of non-White spaces (54). The purposeful underdevelopment of these communities results in unequal exposure to contextual health-related risks that over time exact wear and tear on the body, which contributes to a process of “weathering,” leading to physiological dysregulation, the early onset of disease and disability, and ultimately mortality (55).
Although Puerto Rico appears to have a more flexible attitude toward race (i.e., the concept of “racial democracy”) than the US, there is ample evidence documenting that racial minorities, immigrants (e.g., Dominican immigrants), and phenotypically dark-skinned individuals in Puerto Rico are stigmatized, discriminated against, and experience more socioeconomic disadvantage than their more socially advantaged counterparts (14, 56–59). Notably, Black communities in Puerto Rico4 are largely located along the coastal regions of the Puerto Rican archipelago—a legacy of plantation slavery—and are regions that exhibit lower levels of education, lower median household income, lower median housing values, and higher rates of poverty and unemployment relative to predominantly White communities in Puerto Rico (60). Indeed, for Black Puerto Ricans, systemic and institutional racism across generations and across the lifespan have led to the inequitable access of social, educational, and material resources that have direct (e.g., access to health care) and indirect (e.g., stress and psychosocial resources) effects on health and mortality.
A community-based study of Puerto Rican adults aged 25–55 years in Guayama, Puerto Rico (a southeastern coastal town) found that respondents that are culturally defined as negro (Black) have higher systolic blood pressure (SBP) and diastolic blood pressure (DBP) than those who are classified as blanco (White) or trigueño (racially mixed)5 (61). Additionally, Black Puerto Ricans who occupy higher socioeconomic status (SES) positions exhibit higher SBP and DBP relative to their Black counterparts in low SES contexts (61). The authors posit that Black Puerto Ricans' chronic exposure to institutional and interpersonal discrimination may be linked to their adverse cardiovascular responses (i.e., high blood pressure). Thus, deeply embedded, and multiple dimensions of racism in Puerto Rico are associated with the pronounced residential segregation of Black Puerto Ricans that results in constrained access to resources and opportunities which affect health and mortality.
2.4. Neighborhood residential stability
Living in residentially stable neighborhoods is theorized to promote the health and wellbeing of its residents as it facilitates the development of interpersonal bonds and ties (i.e., social cohesion) that individuals can draw on in times of need (i.e., social support) and may encourage healthy behaviors, and extend longevity. However, a study by Ross et al. (62) found that residential stability was only associated with enhanced psychological wellbeing among residents in affluent neighborhoods. In contrast, residential stability did not benefit the mental health of residents in impoverished communities. Ross et al. posit that living in a poor, stable neighborhood does not confer mental health advantages since residents of these environments do not have the instrumental and material resources needed to mitigate the high levels of disorder in their communities. For example, the chronic stress associated with living long-term in a neighborhood where the streets are dirty, noisy, and dangerous repeatedly activates the stress response, which can contribute to blood pressure and brain changes associated with mental and physical health outcomes (63). Thus, the effects of residential stability need to be considered in the context of a neighborhood's economic resources available.
Data from the U.S. Census and Puerto Rican Community Surveys show that Puerto Rican have high residential stability (64); however, no study, to our knowledge, has examined whether neighborhood-level variation in residential stability is beneficial or detrimental to the health of older adults in Puerto Rico. The scant research that does exist on island-born Puerto Ricans residing in the mainland U.S. has shown that living in ethnically dense, low NSECs reported worse physical health than island-born Puerto Ricans living in other types of NSECs (65). Individuals residing in ethnic enclaves tend to share common sociocultural characteristics (e.g., language and cultural background) and have strong social ties with community members, which have been found to be beneficial for health and mortality. However, enclaves that are formed involuntarily due to housing discrimination may not offer opportunities necessary for economic development at the individual and community levels. Given the high rates of poverty across the archipelago, we can infer that residential stability may not confer health benefits for Puerto Ricans who are living in disadvantaged NSECs.
2.5. Neighborhood family structure
Research on the association between neighborhood family structure and mortality is scarce; however, neighborhood family structure is related to the formation of social ties, which has been shown to have a robust association with extended longevity (66). For example, residents in neighborhoods with high family dissolution (e.g., single-parent households) have lower participation rates in formal voluntary organizations and local affairs. These forms of participation provide opportunities for individuals to integrate within the larger community—additionally, neighborhoods with a high percentage of individuals living alone present opportunities for crime. Individuals who live alone are more likely to go outside alone, which increases the likelihood of a targeted crime (e.g., robbery). These incidents are more likely to instill perceptions of neighborhood disorder that may contribute to the dissolution of social ties and an increased risk of mortality.
Traditionally, Puerto Ricans are very family oriented, embody familism,6 and their families encompass extended and non-blood relatives (e.g., godparents and informally adopted children). The traditional structure of family dynamics in Puerto Rico has historically benefited older family members who often rely on family-based care. Recent research shows that intergenerational co-residence (e.g., children living with their older parents) is associated with increased functional and health support among older adults in Puerto Rico (67). However, the outmigration of younger Puerto Ricans to the US mainland, has led to a significant reduction in the number of family members available to provide care for older adult family members. Moreover, with increasing numbers of Puerto Ricans migrating in search of economic and educational opportunities, we can expect a higher risk of social isolation and lower social participation among older adults, which may be detrimental to mental and physical health (3). Thus, we can expect that communities in Puerto Rico with a high proportion of older adults that live alone and have a high proportion of single-parent households may be associated with worse health and an increased risk of mortality.
2.5.1. The present study
There is compelling theoretical and empirical evidence illustrating how various dimensions of the neighborhood environment co-occur and/or interact to influence the risk of mortality. Given the limited knowledge on the types of residential environments that older Puerto Ricans reside in, it is important to characterize the places where they live based on the factors discussed above. Previous research has shown that using latent class (or profile) models offers an efficient and statistically robust means of summarizing many indicators that constitute neighborhood risks and resources that are not captured by continuous scales or indices (68, 69). We intend to employ this method to classify how various neighborhood characteristics cluster together to create distinct neighborhood typologies that capture risk for all-cause mortality.
3. Materials and methods
3.1. Data
3.1.1. Individual-level data
This study used data from the Puerto Rico Elderly Health Conditions Project (PREHCO), a representative longitudinal cohort study of community-dwelling Puerto Ricans aged 60 and older residing on the archipelago's main island that began in May 2002, with follow-up interviews completed in 2006–2007 and 2021–2022 (the data and documentation are not yet publicly available) (70). Response rates for the first two waves of PREHCO are high (>90.0%). The 4,291 respondents included in the PREHCO baseline sample were derived from a multistage, stratified sample of older adults, including oversampling in regions heavily populated by Afro-descendant individuals (e.g., residents in Loíza) and individuals over 80 years of age. Face-to-face interviews were conducted with each respondent in Spanish or with a proxy if a respondent had cognitive limitations. Additional information on the study and its design is provided elsewhere (71–73).
PREHCO obtained mortality information on respondents using a combination of the National Death Index (NDI) mortality data and PREHCO-identified deaths using reports by family members or the Puerto Rican death registry. Respondents were matched to the National Death Index (NDI) from their first PREHCO interview in 2002–2003 to December 2020, using the available matching variables in the PREHCO study, including social security number (SSN), name (first, middle, father's last name and/or mother's last name), birth date (month and year), and sex (female or male). We would like to note that many Puerto Ricans use two surnames, which adds to the difficulty in NDI matching. Thus, the investigators examined different combinations of respondents' last names to increase the likelihood of a positive match for those with two last names. Additional deaths were identified through November 2021 using family reports or the Puerto Rican death registry. The data file comprising the currently restricted PREHCO mortality database contains the PREHCO respondent's case identification number, the mortality status of the respondent (presumed dead or alive), year of death, month of death, day of death (for some), and cause of death (for most respondents). Two thousand eight hundred and thirty-two all-cause presumed deaths were identified from the cohort of 4,291 PREHCO respondents.
3.1.2. Neighborhood-level data
Data on baseline neighborhood characteristics were constructed from the 2000 Decennial US Census at the block-group level downloaded from Social Explorer and were linked with the PREHCO data (74). Census block groups typically include 600 to 3,000 people and is the smallest geographical unit for which the US Census Bureau publishes sample data. PREHCO respondents were linked to their affiliated census block group by linking their records in the public-use PREHCO to the restricted-use PREHCO geographic data file. These data were then merged with the 2000 US Census data using Federal Information Processing Standard (FIPS) codes to link the files. Out of the 2,477 unique census block group identifiers for Puerto Rico in 2000,7 we identified 233 unique census block groups in which PRECHO respondents resided at the time of the baseline interview, with 1–47 observations in each block group.
3.1.3. Sample selection
The baseline PREHCO cohort sample consisted of 4,291 unique respondents aged 60 and older. Given the design of the present study, we focused on individuals who were able to complete the full interview at baseline (n = 3,713). Respondents that needed proxies to do the interview were not asked health-related questions relevant to the present study (n = 578). Furthermore, we excluded respondents (n = 24) in neighborhoods with <5 individuals in any given block group to minimize statistical bias (75). Lastly, we excluded ~6% of participants (n = 220) due to missingness on baseline covariates. The variables with the highest prevalence of missing values were body mass index (BMI; 5%) and receipt of government-related income and services (1%). The final analytical sample included 3,469 participants.
Participants excluded from the analytical sample were more likely to be older (76.8 vs. 70.3 years), less likely to be married or partnered (38.8 vs. 53.2%), reported lower levels of education (6.1 vs. 8.3 years), and were more likely to receive government-related income and services (see Supplementary Table 1). Additionally, excluded participants were less likely to be obese, current smokers, and physically active. Excluded participants were also more likely to report chronic conditions and disability. We caution readers that the health profiles of excluded participants may be underestimated since proxy interviews were not asked all of the health-related questions. Thus, our analytical sample of community-dwelling older Puerto Ricans is relatively healthier than the general population of older adults in Puerto Rico. A detailed scheme showing the exclusion criteria and the analytic sample inclusion is provided in Supplementary Figure 1.
Additionally, given that measures included in our analysis are time varying, we briefly describe changes in sample characteristics for Wave 2 of PREHCO. From our analytical sample of 3,469 participants identified in Wave 1, 941 respondents (27%) did not have information reported in Wave 2 relevant to the analysis, including 226 proxy interviews, 27 respondents that became institutionalized, 319 that were lost at follow-up, and 369 respondents that were reported dead. To keep all respondents in our analysis, we conducted multiple imputation using chained equations (MICE) for missing data at Wave 2 using the mi suite of commands in Stata (76). We used the distribution of the observed data from Waves 1 and 2 to estimate a set of plausible values for the missing data in Wave 2. We then used Bodner's approach of generating the number of imputed data sets equivalent to the total percent missing and Rubin's rule for combining across data sets (77–82).
3.2. Measures
3.2.1. Mortality
The outcome of interest is all-cause mortality from May 2002 to November 2021. We calculated the time to censoring or death from the year of the interview to the year of death or censoring. For those who did not die in the interval, the censoring date was November 2021. We used years as the time metric.
3.2.2. Individual-level characteristics
Three groups of individual-level variables were considered as potential confounders in examining the role of neighborhood context and all-cause mortality—demographic, socioeconomic, and health characteristics.
3.2.2.1. Demographic variables
Age is measured in continuous years. We also included an age squared term, so we can model more accurately the effect of age rather than assuming the effect is linear for all ages. Sex was dichotomized as male or female. Marital status was dichotomized as married or partnered vs. never married, widowed, separated, or divorced. A dichotomous indicator for whether the respondent had moved from their main residence reported at baseline was also included.
3.2.2.2. Socioeconomic variables
Educational attainment is measured as continuous years of education completed. Given that PREHCO has limited variables for assessing individual income (e.g., not having combined household annual income or poverty thresholds) and wealth (e.g., not having a standardized measure of all assets and debt), we used indirect measures of income, including whether a respondent reports having difficulty paying for their daily necessities (categorized as never, sometimes, and often), whether they receive income from social welfare or the department of the family8 (yes/no), whether they receive income from the nutritional assistance program (yes/no), and whether they have government-sponsored health insurance (excluding Medicare; yes/no). Given the strong association between individual-level socioeconomic position and mortality, it is crucial to adjust for individual socioeconomic measures to ensure the validity of neighborhood-level factors (83).
3.2.2.3. Health characteristics
We included indicators related to the morbidity process such as health behaviors, health conditions, and disability (84). Health behaviors included obesity, current smoking status, and physical activity. Dichotomous indicators were used to classify respondents as obese (i.e., body mass index of ≥30 kg/m2), for whether the respondent reported being a current smoker at the time of the interview (yes/no), and whether they engaged in either moderate or vigorous physical activity at least three times per week (yes/no).
Health conditions included cardiometabolic diseases, other chronic conditions, and severe depression. Cardiometabolic diseases were a sum of whether the respondent self-reported heart problems (e.g., coronary heart disease, congestive heart failure, and heart attack), stroke, hypertension (including medication use) and diabetes (including medication use), ranging from 0 to 4. Other chronic conditions were a sum of self-reported cancer, lung disease, and arthritis, ranging from 0 to 3. We used the geriatric depression scale in its 15-item version (GSD-15) to classify respondents as having severe depression (85). Possible scores range from 0 (no depression) to 15 points (severe depression). Respondents were classified as having severe depression if they scored ≥10 points.
Disability was based on whether a respondent reported limitations in activities of daily living (ADLs) and instrumental activities of daily living (IADLs). ADLs are a continuous measure ranging from 0 to 6 and included difficulty with bathing, eating, dressing, walking across a room, getting in and out of bed, and using the toilet (86). IADLs are a continuous measure ranging from 0 to 7 and included difficulty with using the telephone, managing transportation, buying food or clothing, preparing meals, doing household tasks, taking medications, and managing finances.
3.2.3. Neighborhood-level characteristics
We included variables at the block group level that are theoretically related to and have been identified in previous studies as being associated with all-cause mortality. Neighborhood characteristics included 19 indicators related to the neighborhood constructs of socioeconomic status, household composition, minority status, and housing and transportation. These indicators included the proportion of the population living in a rural area,9 Black residents, residents aged ≥65 years, older adults living alone, residents that lived in the same house past 5 years (residential stability), residents with <9 years of education, residents aged ≥16 years unemployed, residents aged ≥16 years employed in management, professional, and related occupations, households with income ≥$40,000, households with interest, dividend, or rental income, households with public assistance income, residents below 150% of the poverty threshold, single-parent households with children <18 years of age, renter-occupied housing units, residents living in crowded housing units, occupied housing units without complete plumbing, occupied housing units without a telephone, occupied housing units without a motor vehicle, and homes valued ≥150 k.
3.3. Statistical analysis
A latent profile analysis (LPA) was conducted using the gsem feature on Stata to characterize the types of neighborhood environments that older Puerto Ricans resided in at baseline (88). LPA is a semi-parametric finite mixture model that identifies homogenous subgroups based on common characteristics, creating mutually exclusive and exhaustive latent classes. LPA sorts data using posterior probabilities that calculate the probability of membership in each latent class given. Unlike other agglomerative approaches, such as cluster analysis and factor analysis, LPA is a non-parametric statistical technique that relaxes assumptions about normality and linearity in the variables used in the analyses, making LPA ideal for analyzing neighborhood-level characteristics since they do not have normal distributions. We selected the class solution that best represented the data using a combination of model fit statistics, the interpretability of the classes that emerged, and sample size per class once combined with the PREHCO data. When comparing class solutions based on model fit statistics, generally, lower values of the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) are preferred (89); and entropy with values approaching 1, indicating a clear delineation of classes, are preferred (90).
Next, we described the characteristics of the PREHCO analytic sample by each neighborhood cluster that emerged. Means and percentages were calculated using the xtsum and xttab features in Stata to account for the multilevel design and repeat observations.
Lastly, we estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for all-cause mortality by applying a multilevel mixed-effects parametric survival model with a Weibull distribution and Berndt–Hall–Hall–Hausman (BHHH) optimization algorithm using the mestreg feature in Stata. We modeled our data with a three-level hierarchical structure: respondents (level 1) nested within each wave (level 2) and census block groups (level 3). Time-to-event was defined as the elapsed time, in years, from the baseline interview to the date of death or the end of the study follow-up, whichever came first. When we fitted a model, we included the neighborhood clusters and controlled for individual-level demographic variables: sex, age, age squared, and marital status (Model 1). Next, we proceeded to add individual-level socioeconomic indicators: education, income from social welfare, income from the nutritional assistance program, and government-sponsored health insurance (Model 2). Lastly, we added individual-level health characteristics: obesity, smoking, physical activity, cardiometabolic conditions, other chronic conditions, severe depression, and disability (i.e., ADLs and IADLs; Model 3).
All data wrangling, visualization, and analyses were conducted in Stata/MP version 17.0 (91). The data were weighted using PREHCO-provided sampling weights to ensure the representativeness of the PREHCO survey and to account for the sampling design to get reliable statistical estimates. The study protocol was deemed exempt by the Institutional Review Board at Syracuse University.
4. Results
4.1. Neighborhood clusters derived from the LPA
Latent profile models were fit based on 19 block-group level indicators using the 2,477 observations (i.e., unique block groups) available in the 2000 US Census for Puerto Rico, ranging from two to seven classes. Based on the model fit statistics, sample size, and accounting for interpretability, we chose the five-class model as having the best fit for further analysis (see Supplementary Tables 2–7). A summary of the five-class solution of neighborhood clusters is presented in Table 1.
We labeled the first cluster Urban Low Deprivation (Class 1), representing 7.9% of census block groups in Puerto Rico (N = 236). This cluster was characterized by block groups that were almost all urban, had the lowest proportion of Black individuals present, the highest proportion of older adults present, very favorable socioeconomic conditions, stable family structure, and favorable housing features relative to the other classes.
The second cluster was labeled Urban Moderate-Low Deprivation (Class 2) and represented 32.3% census block groups in Puerto Rico (N = 783). This cluster was characterized by block groups that were like the previous neighborhood cluster but notably had lower socioeconomic conditions, family structures that were somewhat less stable, and less favorable housing conditions compared to the first neighborhood cluster.
We labeled the third cluster Rural Moderate Deprivation (Class 3), representing 5.9% of census block groups in Puerto Rico (N = 143). This cluster was characterized by block groups that were predominantly rural, had a low proportion of Black individuals present, the lowest proportion of older adults living alone, unfavorable socioeconomic conditions, stable family structure, and unfavorable housing conditions relative to previous classes.
We labeled the fourth cluster Urban Moderate-High Deprivation (Class 4), representing 47.0% of census block groups in Puerto Rico (N = 1,149). This cluster was characterized by block groups that were predominantly urban, a higher proportion of older adults living alone, less favorable socioeconomic conditions, family structures that were less stable, and less favorable housing conditions relative to previous classes.
The final cluster represented 6.8% of census block groups in Puerto Rico (N = 166) and was labeled Urban High Deprivation (Class 5). This cluster was characterized by block groups that were urban, had the highest proportion of Black individuals present, the lowest proportion of older adults present yet the highest proportion of older adults living alone, very unfavorable socioeconomic conditions, unstable family structure, and unfavorable housing conditions relative to the other classes.
To better contextualize where these neighborhood clusters are geographically located in Puerto Rico, we provide a map of the neighborhood clusters identified in Puerto Rico by census block group (Figure 1). Neighborhoods that were classified as Urban Low Deprivation and Urban High Deprivation were mainly found in the municipalities of San Juan (the largest municipality), Ponce (the largest municipality outside the San Juan area), and Mayagüez (the largest municipality on the west side of the island). Neighborhoods characterized as Urban Low-Moderate Deprivation tended to be clustered outside larger municipalities (e.g., outside of San Juan). Neighborhoods characterized as Urban High-Moderate Deprivation and Rural Moderate Deprivation were distributed across the archipelago. Notably, neighborhoods in the Rural Moderate Deprivation cluster tended to be in the mountainous regions of the archipelago (i.e., the central part of Puerto Rico).
Figure 1. The distribution of neighborhood clusters by year 2000 census block groups in Puerto Rico. Data source: 2000 U.S. Decennial Census.
4.2. Characteristics of older Puerto Ricans by neighborhood cluster
The summary statistics of the PREHCO study sample by neighborhood cluster are presented in Table 2. We find that PREHCO respondents who resided in neighborhoods classified as Urban Low Deprivation (Class 1; n = 224), Urban Low-Moderate Deprivation (Class 2; n = 1,153), and Rural Moderate Deprivation (Class 3; n = 153) had a lower proportion of deaths over the study period relative to those residing in more disadvantaged neighborhood contexts. Older Puerto Ricans residing in the most advantaged neighborhood contexts included a higher proportion of female respondents, were older, less likely to move residences between waves, more educated, did not suffer from economic deprivation, and were less likely to report cardiometabolic conditions and disability. In contrast, respondents residing in the Urban High Deprivation (Class 5; n = 281) cluster had a higher proportion of individuals who died over the study period. The composition of this neighborhood cluster included a lower proportion of females, were younger, were the least likely to be married or partnered, more likely to move residences between waves, were less educated, suffered from economic deprivation, and were more likely to be classified with severe depression. Respondents in the Rural Moderate Deprivation (Class 3; n = 153) and Urban High-Moderate Deprivation (Class 4; n = 1,658) neighborhood clusters had similar demographic compositional profiles. However, respondents in the Rural Moderate Deprivation cluster had the lowest years of education attained relative to the other neighborhood clusters and had worse socioeconomic profiles relative to the Urban High-Moderate Deprivation cluster. Moreover, respondents in the Rural Moderate Deprivation cluster had relatively healthier behavioral profiles (e.g., more physically active, and lower proportion of obese individuals and current smokers) compared to the Urban High-Moderate Deprivation cluster.
Table 2. Characteristics of observations included in multilevel analysis by neighborhood cluster, PREHCO 2002–2007.
4.3. Association of neighborhood clusters with all-cause mortality
The results of the fitted multilevel survival models are summarized in Table 3. Hazard ratios (HR) are presented with 95% confidence intervals (CI). Hazard ratios >1 indicate that the mortality hazard is increasing, whereas hazard ratios <1 indicate that the mortality hazard is decreasing. The results of Model 1 (our base model) show that neighborhood clusters are associated with an increased hazard in all-cause mortality among older Puerto Ricans. Older adults that resided in the Urban Low-Moderate Deprivation [HR: 2.94; 95% CI (1.33, 6.49)], Rural Moderate Deprivation [HR: 2.60; 95% CI (1.10, 6.13)], Urban High-Moderate Deprivation [HR: 3.55; 95% CI (1.58, 7.94)], and Urban High Deprivation [HR: 5.59; 95% CI (2.24, 13.96)] clusters at baseline had higher mortality rates over the study period relative to the Urban Low Deprivation cluster. We also observed that female and married or partnered respondents had lower mortality rates over the study period, and that increasing age was associated with higher mortality rates, which is consistent with results from studies in high- and middle-income countries.
Table 3. All-cause mortality estimated from multilevel survival models of older Puerto Rican adults (n = 3,469).
Controlling for individual-level socioeconomic characteristics (Model 2) reduced the HR gradient of all the neighborhood clusters associated with all-cause mortality observed in Model 1. For example, adjusting for individual-level socioeconomic characteristics decreased the HR by ~32–37% for the Urban Low-Moderate Deprivation, Urban High-Moderate Deprivation, and Urban High Deprivation clusters but they remained significantly associated with all-cause mortality. Conversely, adjusting for individual-level socioeconomic characteristics reduced the Rural Moderate Deprivation cluster to non-significance [HR: 1.83; 95% CI (0.78, 4.30)]. Furthermore, our results indicate that higher levels of education and receiving nutritional assistance was associated with lower mortality over the study period, whereas reporting government-sponsored health insurance was associated with higher mortality over the study period.
Additionally controlling for individual-level health characteristics (Model 3) further reduced (changed) the HR for all the neighborhood clusters. The Urban High-Moderate Deprivation and Urban High Deprivation clusters exhibited an ~20–25% decrease (change) in the HR and were still significantly associated with all-cause mortality. For the Urban Low-Moderate Deprivation cluster, adjusting for individual-level health characteristics reduced the association to non-significance [HR: 1.86; 95% CI (0.92, 3.78)]. We also found that current smoker status, and reporting affirmative to individual items for cardiometabolic disease, other chronic conditions, and IADL limitations increased the hazard by 53, 19, 33, and 12%, respectively. In contrast, respondents that reported engaging in physical activity decreased the hazard by 28%.
Post-estimation tests of coefficients from the final model indicated that older Puerto Ricans residing in the Urban High Deprivation cluster were at the highest risk of death over the study period compared to all the other neighborhood clusters in Puerto Rico. Smoothed hazard estimates of the risk of mortality by neighborhood cluster demonstrating this are shown in Figure 2.
5. Discussion
Using a population-based sample of community-residing individuals aged 60 and older in Puerto Rico, this study builds on prior literature documenting the effect of neighborhood environments on all-cause mortality among older adults. Using latent profile analysis to classify neighborhoods based on indicators related to the constructs of socioeconomic status, household composition, minority status, and housing and transportation resulted in five neighborhood clusters with varying patterns of social (dis)advantage: Urban Low Deprivation, Urban Low-Moderate Deprivation, Rural Moderate Deprivation, Urban High-Moderate Deprivation, and Urban High Deprivation. Our results show that older Puerto Ricans residing in neighborhoods classified as Urban High Deprivation and Urban High-Moderate Deprivation in Puerto Rico (over half of our analytical sample) exhibited an increased risk of mortality over the 19-year study period after adjustment for individual-level covariates. This suggests that a high concentration of unsupportive contexts for healthy aging increases the risk of premature death. This finding is consistent with other studies in the US and Latin America that have found exposure to disadvantaged neighborhood contexts to be a robust predictor of poor health outcomes and increased risk of mortality (27, 28, 92).
In contrast, residing in neighborhoods classified as Rural Moderate Deprivation and Urban Low-Moderate Deprivation was associated with all-cause mortality among older adults, however the association was attenuated once individual-level socioeconomic factors and health-related characteristics were accounted for, respectively. Previous research has shown that individuals residing in rural communities in the US tend to be less educated, have higher rates of poverty, and are less likely to have health insurance than those residing in urban communities (93). In the case of older Puerto Ricans that reside in Rural Moderate Deprivation contexts, our results indicate that the socioeconomic composition of individuals residing within these communities is a more important risk factor for all-cause mortality than the deprivation that exists at the community level. Furthermore, we can infer that older adults with socioeconomic or material advantages living in these communities can alleviate some of the adverse effects and exposures associated with these environments, which may be a family-level social selection mechanism that is unaccounted for in this study (94). It is possible that individuals with economic advantages residing in rural areas in Puerto Rico have been there for generations and chose to stay for reasons related to social, cultural, human, and financial capital (95). For older adults in Urban Low-Moderate Deprivation neighborhood contexts, we can infer that these individuals may self-select into neighborhoods with access to a wealth of social and material resources, such as having access to preventive health care services, and access to medical care that allows for the management of age-related diseases, which can increase longevity.
With the combination of rapid aging and compounding disasters in Puerto Rico, it is imperative to document and account for multilevel determinants of mortality for older adults across later stages of the life course. From a risk environment perspective, there is a need to delineate the environmental factors associated with the risk of mortality, such as the types of environments (e.g., physical, social, economic, and policy) and level of environmental influence (micro and macro), because understanding the places in which harm is produced and reduced offers a broader vision for intervention (96). For instance, a recent review found that the long-term impacts of air pollution, heavy metals, chemicals, ambient temperature, noise, radiation, and urban residential surroundings are associated with increased mortality (97). Since aging is an active response to “weathering,” we must consider how these environmental exposures are related to increases in inflammation, metabolic dysregulation, and genetic damage across the life span, increasing mortality risk. Specific to older adults, as their biological capacity declines with normal aging, the effects of deleterious environmental exposures may be exacerbated among individuals who enter the later stage of the life course with pre-existing health conditions and disabilities (98). Indeed, the biophysiological mechanisms underlying the neighborhood-mortality association are just beginning to be elucidated. Nonetheless, evidence does show that there are links between social factors, physiological dysregulation, and adult mortality (99). Future data collection efforts of older adults in Puerto Rico should include measures that represent multiple regulatory physiological systems (e.g., cardiovascular, metabolic, and immune) to comprehensively capture neighborhood influences on biology, and their contribution to health and mortality risks.
Considering Puerto Rico's socio-structural reality—including high levels of poverty, a deficient infrastructure, a fragile healthcare system, the dismantling of the public education system, and hazardous environmental exposures—a health disparities framework was established to reflect historical and sociocultural influences of the Puerto Rican population (100). We can draw on this framework to highlight how present disparities are rooted in historical, cultural, political, and economic factors that influence biology and behaviors and to illustrate the complex relationship between the neighborhood environment and mortality. For example, a recent study found that Puerto Rican adults residing in San Juan had multiple lifestyle risk factors and cardiometabolic conditions and recommended targeted efforts to improve the health care system and material resources among socially disadvantaged populations (101). While increasing material resources among older residents in the most disadvantaged neighborhood contexts may ease some of the challenges of aging in place, it does not get at the systemic causes of these challenges. For instance, the ports of Puerto Rico are controlled by mainland US agencies, leading to the high costs of (healthy) food on the archipelago (100). As a result, some older adults may forgo eating foods that may improve or better manage their health and decrease their mortality risk since they must make constrained choices on what to spend their limited incomes on. Thus, we recommend that policymakers, health care providers, and leaders across industries to (1) understand how individual health and mortality is embedded within larger social, cultural, structural, and historical contexts, and (2) make concerted efforts to reach out to residents living in disadvantaged community contexts to understand better what they need to successfully age in place in Puerto Rico. A study of residents in La Perla (an informal shantytown in San Juan with a high proportion of older adults) found that despite living in socially and economically disadvantaged residential environments, the residents reported high residential satisfaction because they built their neighborhood environment according to their community needs and have a network of support (102). This suggests that community engagement is essential to identify the health and social needs of Puerto Rican older adults and improve health in neighborhoods directly affected by inequities (103).
5.1. Limitations
Several limitations of this research should be acknowledged. First, we must recognize the physical resilience and robustness of Puerto Ricans who survived to older ages (i.e., aged ≥ 60 years) who were able to participate in the PREHCO study. Previous research has found that survival bias (or, selective survival) can attenuate associations between harmful exposures and age-related diseases, suggesting that the effects of harmful neighborhood environments may not be as pronounced among older adults and are likely underestimated (104).
Second, there are limitations associated with the operationalization of neighborhoods. We selected the smallest census unit for which we could obtain data—census block groups—to conceptualize neighborhoods in this study, an improvement from previous studies that have used census tracts as a neighborhood unit. However, recent research has emerged on the importance of activity spaces—defined as the places individuals encounter due to their day-to-day activities, which may not necessarily include their residential areas (105). Older adults may have activity spaces in more favorable or less advantageous environments relative to their residential settings that affect resources, exposures, benefits, and risks that have multifaceted effects on health and mortality. Future data collection efforts should consider capturing mobility and location information on older adults in Puerto Rico.
Third, using LPA to classify neighborhood clusters depends on the measures included to identify class types. Our findings may be biased by the exclusion of neighborhood characteristics important for distinguishing underlying neighborhood clusters, such as the built environment (e.g., availability of green spaces), availability of health care (e.g., number of physicians and number of facilities), neighborhood crime (e.g., violent offenses), and air pollution (e.g., PM2.5), which we lacked data on, to determine whether the identification of neighborhood clusters is improved. Nonetheless, we included multiple neighborhood variables across multiple neighborhood constructs that have been used in previous studies of all-cause mortality.
Fourth, as with any observational study, this study has unmeasured potential confounders that limit causal inference. For example, due to the limited measures related to income and wealth available in PREHCO, we could not examine if the influence of the neighborhood context differed by individual-level socioeconomic status (SES; e.g., low vs. moderate vs. high SES). Previous research has shown that death rates were higher among low SES individuals residing in high SES neighborhoods (92, 106). This suggests that there are potentially other subpopulations not captured in this study who are at higher risk for death.
Finally, we did not examine residential trajectories over time, which is especially relevant for Puerto Rico given the budget crisis, the great recession, the debt crisis, and Hurricanes Irma and María that may have resulted in increases in spatial inequality. PREHCO has publicly available data for two waves (2002–2003 and 2006–2007). The third wave of surviving respondents of PREHCO will be publicly available soon, and the fourth wave of data collection will begin later this year. These data will allow the creation of a longitudinal database to examine residential trajectories over time and their association with mortality.
Despite these limitations, our study makes several contributions on the role of neighborhoods on older adult mortality. First, we focus on older adults in Puerto Rico—a segment of the US and Latino population that is overlooked in US-based neighborhoods research and aging research more broadly. Second, we used latent profile analysis to summarize multiple indicators that constitute neighborhood risks and resources that are not captured by continuous scales or indices. Using this approach to identify neighborhood clusters associated with an increased risk of death in late life may help inform “upstream” points for structural interventions that can extend healthy years of life among older adults who have had adverse experiences throughout their life course. Third, we used longitudinal data to help establish causal inference. Using multilevel methods and longitudinal data, we assessed the temporal relationship of the association between the neighborhood context at baseline and 19-year all-cause mortality, controlling for possible confounders, allowing for more robust causal inferences. This investigation serves as a foundation to highlight a multilevel perspective of social determinants of health in Puerto Rico. Collectively, we must reframe the narrative on the aging process in Puerto Rico to understand the interplay that historical, environmental, social, behavioral, and biological factors have on health and mortality in this rapidly aging population. Through these efforts, we can identify opportunities to assess and improve the health and wellbeing of older Puerto Rican adults.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://doi.org/10.3886/ICPSR34596.v1 for PREHCO and https://www.socialexplorer.com/ for the 2000 US Decennial Census.
Ethics statement
The original two waves of the PREHCO study complied with all the IRB requirements at the University of Wisconsin-Madison and the University of Puerto Rico. The use of the NDI mortality data with the PREHCO study complied with IRB requirements at the University of Alabama-Birmingham. The use of the geographic data with the PREHCO study complied with IRB requirements at Syracuse University.
Author contributions
CG conceptualized and designed the study, organized and conducted the statistical analysis, interpreted the results and findings, prepared all data visualization, and wrote the manuscript. MG assisted with the interpretation and validation of results and revised the manuscript critically for important intellectual content. MM created the mortality data for PREHCO for use in this study and revised the manuscript. MC assisted with the conceptualization of the study, revised the manuscript, and acquired financial support for PREHCO. All authors approve the submitted version of this manuscript and agree to be accountable for the content of the work.
Funding
The authors gratefully acknowledge the funding of this research by the National Institute on Aging of the National Institutes of Health (R01AG064769, R01AG064769-03S1, and L30AG074407) and the services and support of the Center for Aging and Policy Studies, awarded to Syracuse University, in consortium with Cornell University and the University at Albany, by the National Institute on Aging of the National Institutes of Health (P30AG066583). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.995529/full#supplementary-material
Footnotes
1. ^Healthy aging here is defined as good health, high mental and physical functioning, and active involvement in life.
2. ^The central region of Puerto Rico has many mountainous places and a higher proportion of individuals living in rural areas.
3. ^Food swamps describe the combination of food deserts and a high-density of stores and restaurants that offer high-calorie fast food and junk food.
4. ^Predominantly Black communities in Puerto Rico are located in the municipalities of Loíza, Yabucoa, and Patillas to name a few.
5. ^Trigueño in English literally means “wheat colored” and is used to (1) describe a light-skinned person with a slightly tanned complexion; or (2) describe someone with European, African, and Native American heritage. Being classified as trigueño may be socially advantageous since it does not carry the stigma of blackness.
6. ^Familism emphasizes obligation and duties of families to one another.
7. ^Every census block group has a unique 12-digit FIPS code.
8. ^The Puerto Rico department of the family specializes in individual and family social services.
9. ^The US Census Bureau classifies urban/rural areas based on population thresholds, density, and land-use characteristics. Typically, rurality is designated as all population, housing, and territory not included within an urbanized area or urban cluster (≥2,500 people) (87).
References
1. Pérez C, Ailshire JA. Aging in Puerto Rico: A comparison of health status among Island Puerto Rican and Mainland US older adults. J Aging Health. (2017) 29:1056–78. doi: 10.1177/0898264317714144
2. García C, Rivera FI, Garcia MA, Burgos G, Aranda MP. Contextualizing the COVID-19 era in Puerto Rico: Compounding disasters and parallel pandemics. J Gerontol Ser B. (2021) 76:e263–7. doi: 10.1093/geronb/gbaa186
3. Matos-Moreno A, Verdery AM, Mendes de Leon CF, De Jesús-Monge VM, Santos-Lozada AR. Aging and the left behind: Puerto Rico and its unconventional rapid aging. Gerontologist. (2022) 2022:gnac082. doi: 10.1093/geront/gnac082
4. Keenan JM, Hauer ME. Resilience for whom? Demographic change and the redevelopment of the built environment in Puerto Rico. Environ Res Lett. (2020) 15:e074028. doi: 10.1088/1748-9326/ab92c2
5. Lozada ARS, Estrada ALV. The Population Decline of Puerto Rico: An Application of Prospective Trends in Cohort-Component Projections. (2015). Available online at: http://rgdoi.net/10.13140/RG.2.1.2895.5927 (accessed May 16, 2022).
6. Ramos JGP, Garriga-López A, Rodríguez-Díaz CE. How is colonialism a sociostructural determinant of health in Puerto Rico? AMA J Ethics. (2022) 24:305–12. doi: 10.1001/amajethics.2022.305
7. Valle AJ. Race and the empire-state: Puerto Ricans' unequal U.S. Citizenship Sociol Race Ethn. (2019) 5:26–40. doi: 10.1177/2332649218776031
8. Czyzewski K. Colonialism as a broader social determinant of health. Int Indig Policy J. (2011) 2:7337. doi: 10.18584/iipj.2011.2.1.5
9. Colón-López A, García C. 20th century Puerto Rico and later-life health: The association between multigenerational education and chronic conditions in Island-Dwelling older adults. J Aging Health. (2022) 2022:89826432210975. doi: 10.1177/08982643221097532
10. Zulawski A. Environment, urbanization, and public health: The bubonic plague epidemic of 1912 in San Juan, Puerto Rico. Lat Am Res Rev. (2022) 53:500–16. doi: 10.25222/larr.424
11. Grau HR, Aide TM, Zimmerman JK, Thomlinson JR, Helmer E, Zou X. The ecological consequences of socioeconomic and land-use changes in postagriculture Puerto Rico. Bioscience. (2003) 53:1159. doi: 10.1641/0006-3568(2003)053[1159:TECOSA]2.0.CO;2
12. Urban F. La Perla – 100 years of informal architecture in San Juan, Puerto Rico. Plan Perspect. (2015) 30:495–536. doi: 10.1080/02665433.2014.1003247
13. Azar D, Rain D. Identifying population vulnerable to hydrological hazards in San Juan, Puerto Rico. GeoJournal. (2007) 69:23–43. doi: 10.1007/s10708-007-9106-8
14. Denton NA, Villarrubia J. Residential segregation on the Island: The role of race and class in Puerto Rican neighborhoods. Sociol Forum. (2007) 22:51–76. doi: 10.1111/j.1573-7861.2006.00004.x
15. Torres-Cintrón M, Ortiz AP, Ortiz-Ortiz KJ, Figueroa-Vallés NR, Pérez-Irizarry J, Díaz-Medina G, et al. Using a socioeconomic position index to assess disparities in cancer incidence and mortality, Puerto Rico, 1995-2004. Prev Chronic Dis. (2012) 9:E15. doi: 10.5888/pcd9.100271
16. Diez Roux AV. Estimating neighborhood health effects: The challenges of causal inference in a complex world. Soc Sci Med. (2004) 58:1953–60. doi: 10.1016/S0277-9536(03)00414-3
17. Dick C. The health and retirement study: Contextual data augmentation. For Health Econ Pol. (2022) 2021:68. doi: 10.1515/fhep-2021-0068
18. Glass TA, Balfour JL. Neighborhoods, aging, and functional limitations. In: Kawachi I, Berkman LF, editors, Neighborhoods and Health. Oxford: Oxford University Press (2003). p. 303–34. Available online at: http://www.oxfordscholarship.com/view/10.1093/acprof:oso/9780195138382.001.0001/acprof-9780195138382-chapter-14 (accessed April 17, 2020).
19. Clarke PJ, Weuve J, Barnes L, Evans DA, Mendes de Leon CF. Cognitive decline and the neighborhood environment. Ann Epidemiol. (2015) 25:849–54. doi: 10.1016/j.annepidem.2015.07.001
20. Gill TM, Zang EX, Murphy TE, Leo-Summers L, Gahbauer EA, Festa N, et al. Association between neighborhood disadvantage and functional well-being in community-living older persons. J Am Med Assoc Intern Med. (2021) 181:1297. doi: 10.1001/jamainternmed.2021.4260
21. Geronimus AT. The weathering hypothesis and the health of African-American women and infants: Evidence and speculations. Ethn Dis. (1992) 2:207–21.
22. Lawton MP, Nahemow L. Ecology and the aging process. In: Eisdorfer C, Lawton MP, editors, The Psychology of Adult Development and Aging. Washington, DC: American Psychological Association (1973). p. 619–74. Available online at: http://content.apa.org/books/10044-020 (accessed June 2, 2020).
23. Robert SA. Socioeconomic position and health: The independent contribution of community socioeconomic context. Annu Rev Sociol. (1999) 25:489–516. doi: 10.1146/annurev.soc.25.1.489
24. Stimpson JP, Ju H, Raji MA, Eschbach K. Neighborhood deprivation and health risk behaviors in NHANES III. Am J Health Behav. (2007) 31:215–22. doi: 10.5993/AJHB.31.2.10
25. Kirby JB, Kaneda T. Neighborhood socioeconomic disadvantage and access to health care. J Health Soc Behav. (2005) 46:15–31. doi: 10.1177/002214650504600103
26. Ailshire JA, García C. Unequal places: The impacts of socioeconomic and race/ethnic differences in neighborhoods. Generations. (2018) 42:20–7.
27. Lei MK, Beach SRH, Simons RL. Biological embedding of neighborhood disadvantage and collective efficacy: Influences on chronic illness via accelerated cardiometabolic age. Dev Psychopathol. (2018) 30:1797–815. doi: 10.1017/S0954579418000937
28. Powell WR, Buckingham WR, Larson JL, Vilen L, Yu M, Salamat MS, et al. Association of neighborhood-level disadvantage with Alzheimer disease neuropathology. J Am Med Assoc Netw Open. (2020) 3:e207559. doi: 10.1001/jamanetworkopen.2020.7559
29. Ward-Caviness CK, Pu S, Martin CL, Galea S, Uddin M, Wildman DE, et al. Epigenetic predictors of all-cause mortality are associated with objective measures of neighborhood disadvantage in an urban population. Clin Epigenetics. (2020) 12:44. doi: 10.1186/s13148-020-00830-8
30. Yen IH, Michael YL, Perdue L. Neighborhood environment in studies of health of older adults. Am J Prev Med. (2009) 37:455–63. doi: 10.1016/j.amepre.2009.06.022
31. Dannefer D. Cumulative advantage/disadvantage and the life course: Cross-fertilizing age and social science theory. J Gerontol B Psychol Sci Soc Sci. (2003) 58:S327–37. doi: 10.1093/geronb/58.6.S327
32. Moore KD. An ecological framework of place: Situating environmental gerontology within a life course perspective. Int J Aging Hum Dev. (2014) 79:183–209. doi: 10.2190/AG.79.3.a
33. Elder GH, Rockwell RC. The life-course and human development: An ecological perspective. Int J Behav Dev. (1979) 2:1–21. doi: 10.1177/016502547900200101
34. Cutler DM, Lleras-Muney A. Understanding differences in health behaviors by education. J Health Econ. (2010) 29:1–28. doi: 10.1016/j.jhealeco.2009.10.003
35. Frisvold D, Golberstein E. School quality and the education–health relationship: Evidence from Blacks in segregated schools. J Health Econ. (2011) 30:1232–45. doi: 10.1016/j.jhealeco.2011.08.003
36. Link BG, Phelan J. Social conditions as fundamental causes of disease. J Health Soc Behav. (1995) 35:80. doi: 10.2307/2626958
37. Mirowsky J, Ross CE. Education, Social Status, and Health. 1st ed. London: Routledge (2017). Available online at: https://www.taylorfrancis.com/books/9781351328074 (accessed July 15, 2022).
38. Cooksey-Stowers K, Schwartz M, Brownell K. Food swamps predict obesity rates better than food deserts in the United States. Int J Environ Res Public Health. (2017) 14:1366. doi: 10.3390/ijerph14111366
39. Perry A, Patlak M, Ramirez AG, Gallion KJ. Making Healthy Food and Beverages the Affordable, Available, Desired Choices Among Latino Families. San Antonio, TX: Salud America! (2015).
40. Cagney KA. Neighborhood age structure and its implications for health. J Urban Health. (2006) 83:827–34. doi: 10.1007/s11524-006-9092-z
41. Kubzansky LD, Subramanian SV, Kawachi I, Fay ME, Soobader MJ, Berkman LF. Neighborhood contextual influences on depressive symptoms in the elderly. Am J Epidemiol. (2005) 162:253–60. doi: 10.1093/aje/kwi185
42. Subramanian SV, Kubzansky L, Berkman L, Fay M, Kawachi I. Neighborhood effects on the self-rated health of elders: Uncovering the relative importance of structural and service-related neighborhood environments. J Gerontol B Psychol Sci Soc Sci. (2006) 61:S153–60. doi: 10.1093/geronb/61.3.S153
43. Cagney KA, Cornwell EY. Place, aging, and health. In: Future Directions for the Demography of Aging: Proceedings of a Workshop. National Academies Press (US) (2018). Available online at: https://www.ncbi.nlm.nih.gov/books/NBK513076/ (accessed June 17, 2022).
44. Seeman TE, Crimmins E. Social environment effects on health and aging: Integrating epidemiologic and demographic approaches and perspectives. Ann N Y Acad Sci. (2006) 954:88–117. doi: 10.1111/j.1749-6632.2001.tb02749.x
45. Berkman LF. Social integration, social networks, and health. In: Encyclopedia of Health and Behavior. Thousand Oaks, CA: SAGE Publications, Inc. (2004). Available online at: http://sk.sagepub.com/reference/behavior/n192.xml (accessed April 16, 2020).
46. Matos-Moreno A, Santos-Lozada AR, Mehta N, Mendes de Leon CF, Lê-Scherban F, De Lima Friche AA. Migration is the driving force of rapid aging in Puerto Rico: A research brief. Popul Res Policy Rev. (2022) 41:801–10. doi: 10.1007/s11113-021-09683-2
47. Zsembik BA, Bonilla Z. Eldercare and the changing family in Puerto Rico. J Fam Issues. (2000) 21:652–74. doi: 10.1177/019251300021005007
48. Charles CZ. The dynamics of racial residential segregation. Annu Rev Sociol. (2003) 29:167–207. doi: 10.1146/annurev.soc.29.010202.100002
49. Massey DS. American apartheid: Segregation and the making of the underclass. Am J Sociol. (1990) 96:329–57. doi: 10.1086/229532
50. Sewell AA. The racism-race reification process: A mesolevel political economic framework for understanding racial health disparities. Sociol Race Ethn. (2016) 2:402–32. doi: 10.1177/2332649215626936
51. Williams DR, Rucker TD. Understanding and addressing racial disparities in health care. Health Care Financ Rev. (2000) 21:75–90.
52. Williams DR, Collins C. Racial residential segregation: A fundamental cause of racial disparities in health. Public Health Rep. (2001) 116:404–16. doi: 10.1016/S0033-3549(04)50068-7
53. Williams DR, Lawrence JA, Davis BA. Racism and health: Evidence and needed research. Annu Rev Public Health. (2019) 40:105–25. doi: 10.1146/annurev-publhealth-040218-043750
54. Seamster L, Purifoy D. What is environmental racism for? Place-based harm and relational development. Environ Sociol. (2021) 7:110–21. doi: 10.1080/23251042.2020.1790331
55. Crimmins EM, Seeman TE. Integrating biology into the study of health disparities. Popul Dev Rev. (2004) 30:89–107.
56. Caraballo-Cueto J, Godreau IP. Colorism and health disparities in home countries: The case of Puerto Rico. J Immigr Minor Health. (2021) 23:926–35. doi: 10.1007/s10903-021-01222-7
57. Liu MM, Crowe M, Telles EE, Jiménez-Velázquez IZ, Dow WH. Color disparities in cognitive aging among Puerto Ricans on the archipelago. SSM Popul Health. (2022) 17:100998. doi: 10.1016/j.ssmph.2021.100998
58. Rodriguez-Díaz CE, Lewellen-Williams C. Race and racism as structural determinants for emergency and recovery response in the aftermath of Hurricanes Irma and Maria in Puerto Rico. Health Equity. (2020) 4:232–8. doi: 10.1089/heq.2019.0103
59. Lloréns H, García-Quijano C, Godreau IP. Racismo En Puerto Rico: Surveying perceptions of racism. Cent J. (2017) 29:154–83.
60. U.S. Census Bureau. Quick Facts: Puerto Rico. Washington, DC (2021). Available online at: https://www.census.gov/quickfacts/fact/table/PR/IPE120221 (accessed July 2022).
61. Gravlee CC, Dressler WW, Bernard HR. Skin color, social classification, and blood pressure in Southeastern Puerto Rico. Am J Public Health. (2005) 95:2191–7. doi: 10.2105/AJPH.2005.065615
62. Ross CE, Reynolds JR, Geis KJ. The contingent meaning of neighborhood stability for residents' psychological well-being. Am Sociol Rev. (2000) 65:581. doi: 10.2307/2657384
63. Yang L, Zhao Y, Wang Y, Liu L, Zhang X, Li B, et al. The effects of psychological stress on depression. Curr Neuropharmacol. (2015) 13:494–504. doi: 10.2174/1570159X1304150831150507
64. U.S. Census Bureau. QuickFacts: Puerto Rico. Washington, DC. (2020). Available online at: https://www.census.gov/quickfacts/fact/table/PR/AGE775218 (accessed July 2022).
65. Roy AL, Hughes D, Yoshikawa H. Intersections between nativity, ethnic density, and neighborhood SES: Using an ethnic enclave framework to explore variation in Puerto Ricans' Physical Health. Am J Community Psychol. (2013) 51:468–79. doi: 10.1007/s10464-012-9564-0
66. Sampson RJ. Neighborhood family structure and the risk of personal victimization. In: Byrne JM, Sampson RJ, editors, The Social Ecology of Crime. New York, NY: Springer New York (1986). p. 25–46. Available online at: http://link.springer.com/10.1007/978-1-4613-8606-3_2 (accessed July 15, 2022).
67. Quashie NT, Andrade FCD, Meltzer G, García C. Living arrangements and intergenerational support in Puerto Rico: Are fathers disadvantaged? J Gerontol Ser B. (2022) 2022:gbac044. doi: 10.1093/geronb/gbac044
68. Weden MM, Escarce JJ, Lurie N. Technical Detail and Appendices for a Study of Neighborhood Archetypes for Population Health Research. (RAND Labor and Population working paper series). Report No.: WR-805. Santa Monica, CA: RAND Corporation (2010). Available online at: https://www.rand.org/pubs/working_papers/WR805.html (accessed July 2022).
69. Humphrey J, Lindstrom M, Barton K, Shrestha P, Carlton E, Adgate J, et al. Social and environmental neighborhood typologies and lung function in a low-income, urban population. Int J Environ Res Public Health. (2019) 16:1133. doi: 10.3390/ijerph16071133
70. Palloni A, Luisa Davila A, Sanchez-Ayendez M. Puerto Rican elderly: Health Conditions (PREHCO) Project, 2002–2003, 2006–2007. Inter-university Consort Polit Soc Res. (2013) 2013:ICPSR34596. doi: 10.3886/ICPSR34596
71. McEniry M, Palloni A. Early life exposures and the occurrence and timing of heart disease among the older adult Puerto Rican population. Demography. (2010) 47:23–43. doi: 10.1353/dem.0.0093
72. Palloni A, McEniry M, Dávila AL, Gurucharri AG. The influence of early conditions on health status among elderly Puerto Ricans. Soc Biol. (2005) 52:132–63. doi: 10.1080/19485565.2005.9989106
73. McEniry M, Palloni A, Davila AL, Gurucharri AG. Early life exposure to poor nutrition and infectious diseases and its effects on the health of older Puerto Rican adults. J Gerontol B Psychol Sci Soc Sci. (2008) 63:S337–48. doi: 10.1093/geronb/63.6.S337
74. U.S. Census Bureau. SF3 Prefixed Tables Summary File 3 (SF 3) - Sample Data. Social Explorer (2000). Available online at: https://www.socialexplorer.com/reports/socialexplorer/en/report/183c339a-e43d-11ec-9428-d3638f2d61a2 (accessed July 2022).
75. Clarke P. Theory and methods: When can group level clustering be ignored? Multilevel models versus single-level models with sparse data. J Epidemiol Community Health. (2008) 62:752–8. doi: 10.1136/jech.2007.060798
76. White IR, Royston P, Wood AM. Multiple imputation using chained equations: Issues and guidance for practice. Stat Med. (2011) 30:377–99. doi: 10.1002/sim.4067
77. Bodner TE. What improves with increased missing data imputations? Struct Eq Model. (2008) 15:651–75. doi: 10.1080/10705510802339072
78. Bennett DA. How can I deal with missing data in my study? Aust N Z J Public Health. (2001) 25:464–9. doi: 10.1111/j.1467-842X.2001.tb00659.x
79. Graham JW, Cumsille PE, Shevock AE. Methods for handling missing data. In: Handbook of Psychology: Research Methods in Psychology. Vol 2 2nd, ed. Hoboken, NJ: John Wiley & Sons Inc (2013). p. 109–41. doi: 10.1002/9781118133880.hop202004
80. Schafer JL. Multiple imputation: A primer. Stat Methods Med Res. (1999) 8:3–15. doi: 10.1191/096228099671525676
81. Rubin DB. Multiple Imputation for Nonresponse in Surveys. New York, NY: John Wiley & Sons (2004).
82. Rubin DB. Multiple imputation after 18+ years. J Am Stat Assoc. (1996) 91:473–89. doi: 10.1080/01621459.1996.10476908
83. Meijer M, Röhl J, Bloomfield K, Grittner U. Do neighborhoods affect individual mortality? A systematic review and meta-analysis of multilevel studies. Soc Sci Med. (2012) 74:1204–12. doi: 10.1016/j.socscimed.2011.11.034
84. Crimmins EM, Kim JK, Vasunilashorn S. Biodemography: New approaches to understanding trends and differences in population health and mortality. Demography. (2010) 47:S41–64. doi: 10.1353/dem.2010.0005
85. Yesavage JA, Sheikh JI. 9/Geriatric Depression Scale (GDS): Recent evidence and development of a shorter version. Clin Gerontol. (1986) 5:165–73. doi: 10.1300/J018v05n01_09
86. Katz S. Assessing self-maintenance: Activities of daily living, mobility, and instrumental activities of daily living. J Am Geriatr Soc. (1983) 31:721–7. doi: 10.1111/j.1532-5415.1983.tb03391.x
87. Ratcliffe M, Burd C, Holder K, Fields A. Defining Rural at the U.S. Census Bureau. American Community Survey and Geography Brief. Washington, DC: U.S. Census Bureau (2016). Available online at: https://www.census.gov/content/dam/Census/library/publications/2016/acs/acsgeo-1.pdf (accessed July 2022).
88. Collins LM, Lanza ST. Latent Class and Latent Transition Analysis. Wiley Series in Probability Statistics. Hoboken, NJ: John Wiley & Sons, Inc. (2009). Available online at: http://doi.wiley.com/10.1002/9780470567333 (accessed March 28, 2020).
90. Celeux G, Soromenho G. An entropy criterion for assessing the number of clusters in a mixture model. J Classif. (1996) 13:195–212. doi: 10.1007/BF01246098
92. Yen IH, Kaplan GA. Neighborhood social environment and risk of death: Multilevel evidence from the Alameda county study. Am J Epidemiol. (1999) 149:898–907. doi: 10.1093/oxfordjournals.aje.a009733
93. Centers for Disease Control and Prevention. About Rural Health. Washington, DC (2017). Available online at: https://www.cdc.gov/ruralhealth/about.html (accessed June 2, 2020).
94. Jokela M. Are neighborhood health associations causal? A 10-year prospective cohort study with repeated measurements. Am J Epidemiol. (2014) 180:776–84. doi: 10.1093/aje/kwu233
95. Flora CB, Flora JL, Gasteyer SP. Rural Communities: Legacy and Change. 4th ed. London: Routledge (2018). Available online at: https://www.taylorfrancis.com/books/9780429963247 (accessed November 26, 2022).
96. Rhodes T. The “risk environment”: A framework for understanding and reducing drug-related harm. Int J Drug Policy. (2002) 13:85–94. doi: 10.1016/S0955-3959(02)00007-5
97. Rojas-Rueda D, Morales-Zamora E, Alsufyani WA, Herbst CH, AlBalawi SM, Alsukait R, et al. Environmental risk factors and health: An umbrella review of meta-analyses. Int J Environ Res Public Health. (2021) 18:704. doi: 10.3390/ijerph18020704
98. Geller AM, Zenick H. Aging and the environment: A research framework. Environ Health Perspect. (2005) 113:1257–62. doi: 10.1289/ehp.7569
99. Yang YC, McClintock MK, Kozloski M, Li T. Social isolation and adult mortality: The role of chronic inflammation and sex differences. J Health Soc Behav. (2013) 54:183–203. doi: 10.1177/0022146513485244
100. Lafarga Previdi I, Vélez Vega CM. Health disparities research framework adaptation to reflect Puerto Rico's socio-cultural context. Int J Environ Res Public Health. (2020) 17:8544. doi: 10.3390/ijerph17228544
101. Mattei J, Tamez M, Ríos-Bedoya CF, Xiao RS, Tucker KL, Rodríguez-Orengo JF. Health conditions and lifestyle risk factors of adults living in Puerto Rico: A cross-sectional study. BMC Public Health. (2018) 18:491. doi: 10.1186/s12889-018-5359-z
102. Caldieron J. Residential satisfaction in La Perla Informal Neighborhood, San Juan, Puerto Rico. OIDA Int J Sustain Dev. (2011) 2:77–84.
103. García AA, West Ohueri C, Garay R, Guzmán M, Hanson K, Vasquez M, et al. Community engagement as a foundation for improving neighborhood health. Public Health Nurs. (2021) 38:223–31. doi: 10.1111/phn.12870
104. McGuinness MB, Karahalios A, Kasza J, Guymer RH, Finger RP, Simpson JA. Survival bias when assessing risk factors for age-related macular degeneration: A tutorial with application to the exposure of smoking. Ophthalmic Epidemiol. (2017) 24:229–38. doi: 10.1080/09286586.2016.1276934
105. Cagney KA, York Cornwell E, Goldman AW, Cai L. Urban mobility and activity space. Annu Rev Sociol. (2020) 46:623–48. doi: 10.1146/annurev-soc-121919-054848
Keywords: neighborhood characteristics, mortality, Puerto Rican adults, latent variable analysis, multilevel survival analysis, PREHCO, older adults, social determinants of health (SDOH)
Citation: García C, Garcia MA, McEniry M and Crowe M (2023) The neighborhood context and all-cause mortality among older adults in Puerto Rico. Front. Public Health 11:995529. doi: 10.3389/fpubh.2023.995529
Received: 16 July 2022; Accepted: 14 February 2023;
Published: 09 March 2023.
Edited by:
Michel Guillot, Pennsylvania Hospital, United StatesReviewed by:
Alexis Santos, The Pennsylvania State University (PSU), United StatesIrma Elo, University of Pennsylvania, United States
Copyright © 2023 García, Garcia, McEniry and Crowe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Catherine García, Y2dhcmNpMjQmI3gwMDA0MDtzeXIuZWR1