Next Article in Journal
Emotional Recognition and Its Relation to Cognition, Mood and Fatigue in Relapsing–Remitting and Secondary Progressive Multiple Sclerosis
Previous Article in Journal
Mapping the Scientific Research on Suicide and Physical Activity: A Bibliometric Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying Risk Factors for Hospitalization with Behavioral Health Disorders and Concurrent Temperature-Related Illness in New York State

by
Heather Aydin-Ghormoz
1,2,
Temilayo Adeyeye
1,2,*,
Neil Muscatiello
1,
Seema Nayak
1,
Sanghamitra Savadatti
1,2 and
Tabassum Z. Insaf
1,2
1
Center for Environmental Health, New York State Department of Health, Albany, NY 12208, USA
2
School of Public Health, University at Albany, Rensselaer, NY 12144, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(24), 16411; https://doi.org/10.3390/ijerph192416411
Submission received: 2 November 2022 / Revised: 1 December 2022 / Accepted: 1 December 2022 / Published: 7 December 2022

Abstract

:
Extreme temperature events are linked to increased emergency department visits, hospitalizations, and mortality for individuals with behavioral health disorders (BHD). This study aims to characterize risk factors for concurrent temperature-related illness among BHD hospitalizations in New York State. Using data from the NYS Statewide and Planning Research and Cooperative System between 2005–2019, multivariate log binomial regression models were used in a population of BHD hospitalizations to estimate risk ratios (RR) for a concurrent heat-related (HRI) or cold-related illness (CRI). Dementia (RR 1.65; 95% CI:1.49, 1.83) and schizophrenia (RR 1.38; 95% CI:1.19, 1.60) were associated with an increased risk for HRI among BHD hospitalizations, while alcohol dependence (RR 2.10; 95% CI:1.99, 2.22), dementia (RR 1.52; 95% CI:1.44, 1.60), schizophrenia (RR 1.41; 95% CI:1.31, 1.52), and non-dependent drug/alcohol use (RR 1.20; 95% CI:1.15, 1.26) were associated with an increased risk of CRI among BHD hospitalizations. Risk factors for concurrent HRI among BHD hospitalizations include increasing age, male gender, non-Hispanic Black race, and medium hospital size. Risk factors for concurrent CRI among BHD hospitalizations include increasing age, male gender, non-Hispanic Black race, insurance payor, the presence of respiratory disease, and rural hospital location. This study adds to the literature by identifying dementia, schizophrenia, substance-use disorders, including alcohol dependence and non-dependent substance-use, and other sociodemographic factors as risk factors for a concurrent CRI in BHD hospitalizations.

1. Introduction

Average annual temperatures in New York State (NYS) are expected to rise 2.2–5 °C by 2080 and the state can expect increased days of extreme heat in the warmer months over the next 60 years [1]. Most NYS regions are projected to have an average of 21–74 days over 32 °C by 2080 [1]. Despite this, NYS remains a region that experiences extreme cold events in the winter months, although these events are expected to decrease given the predicted rise in average temperatures [1]. By 2080, it is projected that NYS will have an average of 67–116 days at or below 0 °C [1]. NYS experiences large variations in regional temperatures and no single definition can be used to define extreme weather events across the state. In the United States from 2006 to 2010, 3340 people died from heat-related illness (HRI) and 6652 people died from cold-related illness (CRI) [2]. While previous studies also suggest that most temperature-associated mortality is a result of CRI, increasing ambient temperatures make mortality from HRI an emerging threat in NYS, especially for vulnerable populations [3].
Extreme heat and cold temperatures can impact health directly as well as exacerbate existing medical conditions and, in some instances, lead to death [3,4,5,6,7]. Prolonged exposure to high temperatures causes HRI, including heat stress, heat stroke and heat exhaustion, while prolonged exposure to low temperatures causes a variety of CRI, including hypothermia and frostbite. Both exposures can exacerbate medical conditions, such as respiratory and cardiovascular diseases. Several research studies from around the world have found associations between temperature increases or heatwaves and increased mental health emergency department (ED) visits and hospitalizations [8,9,10,11,12,13,14,15]. ED visits for behavioral health disorders, schizophrenia, mood disorders, neurotic disorders, and morbidity due to substance-use have also been shown to be associated with high temperature [16]. High ambient temperatures are also associated with increased mortality in populations with behavioral disorders, specifically substance-use disorders [17,18,19]. Behavioral health disorders (BHD) have also been found to be significant risk factors of HRI [11]. A study conducted in the U.S. in 2017 by Schmeltz and Gamble analyzed specific BHD diagnoses and the risk of HRI hospitalizations, as well as individual and environmental risk factors. Schmeltz and Gamble found that among BHD hospitalizations, the risk of a concurrent HRI was higher among patients with dementia, psychoactive substance abuse, and nondependent drug abuse [11]. Schmeltz and Gamble (2017) also analyzed risk factors by patient characteristics and showed an increased risk among African Americans, Native Americans, males, and adults over the age of 40 [11]. Additionally, those with lower income, the uninsured, rural patients, and those going to smaller hospitals had an increased risk of hospitalization for BHD with concurrent HRI [11].
The association between low ambient temperature and BHD-related ED visits, hospitalizations, and mortality is inconclusive in the literature. Research finds that mental health comorbid conditions are common in CRI mortality [20]. Additionally, antipsychotic medications, in the presence of multiple risk factors, have been observed to cause hypothermia in the first 7–10 days of initiation or when increasing the dosage [21]. Other research has suggested that cold temperatures are found to decrease adverse mental health outcomes [22,23].
There are a limited number of studies that have investigated, in a population of people hospitalized with BHD, which risk factors contributed to concurrent temperature-related illness (TRI). Due to climate, temperature, and demographic variation between US regions, there is also a need to conduct this study in NYS to identify state-specific risk factors. In this study, we aim to determine in a population of BHD hospitalizations which BHD have an increased risk of concurrent TRI in NYS. Furthermore, we seek to identify the demographic factors associated with an increased risk of a concurrent cold- or heat-related illnesses in NYS in the same population.

2. Materials and Methods

This study uses cross-sectional administrative data from 2005–2019 to identify possible risk factors in the population of BHD hospitalizations for TRIs. The primary data source for this study was hospitalization data from the NYS Statewide Planning and Research Cooperative System (SPARCS). NYS SPARCS provides several categories of statewide patient information data, including demographic data, such as age, sex, race, and ethnicity, insurance payor data, primary and secondary diagnoses, and hospital information for all hospital admissions in NYS excluding Federal and Veterans Health Administration hospitals, Indian Health Service hospitals, and institutionalized (e.g., prison) populations. The risk factors considered in this analysis were specific BHD diagnoses, including age, race, sex, insurance status, comorbidities, hospital size, and hospital location. Hospital size and location were included to assess risk in rural communities where mental health services are less accessible. Variables selected for analysis were based on prior U.S. national research on risk factors for BHD with concurrent HRI, specifically the research conducted by Schmeltz and Gamble [11]. This research was chosen as it provided national level estimates of risk factors for hospitalization with BHD and HRI that were useful to compare to findings in NYS.

2.1. Population

We analyzed administrative claims data of patients aged five and over who had been hospitalized in NYS between 2005 and 2019 with a primary or secondary diagnosis of BHD. We excluded BHD hospitalization under the age of five since diagnosis of BHD in very young children is difficult and potentially inaccurate. The classification of mental illness for this analysis was developed based on a review of the prior literature and expert guidance from the NYS Office of Mental Health.
This study includes a full analysis and two sub-analyses of risk factors for concurrent TRI in BHD hospitalizations. The full BHD analysis includes all mental disorder diagnoses in addition to dementia and the two sub-analyses are limited to mental disorder diagnoses (excluding dementia) and substance-use disorders diagnoses. The mental disorder categories for this study included depressive disorders, schizophrenia spectrum and other psychotic disorders, bipolar and related disorders, anxiety disorders, trauma and stressor related disorders, attention deficit hyperactivity disorder, disruptive disorders, impulse control disorders, conduct disorders, personality disorders, and obsessive compulsive and related disorders. The International Classification of Diseases, Ninth Revision (ICD-9) to Tenth Revision (ICD-10) crosswalk was developed using the Center for Medicare and Medicaid Services’ General Equivalence Mappings (see Supplementary Materials, Table S1 for the BHD Crosswalk) [24]. We included dementia in the full analysis, despite it often being considered a neurological disorder, as previous research demonstrated an association between hospitalization for dementia and concurrent HRI [11]. Additionally, cooler than average temperatures and high temperature variability have been linked to increased dementia hospitalizations [25]. The substance-use disorders categories include alcohol dependence, drug dependence, and non-dependent drug and alcohol use (see Table S1 for the BHD Crosswalk).

2.2. Outcomes

The two main outcomes of interest in this study were hospitalization with BHD and a concurrent HRI diagnosis during the warmer months (May-September), and hospitalization with BHD and a concurrent CRI during the colder months (October–April). TRI data were obtained using hospital administrative coding for HRI and CRI from the NYS SPARCS dataset. Among hospitalizations with BHD and a concurrent HRI diagnosis, HRI was classified as heat stroke, sunstroke, heat syncope, heat exhaustion and other heat effects, accidents caused by excessive heat due to weather conditions, and exposure to excessive natural heat and sunlight (ICD-9: 992, E900.0, E900.9; ICD-10: T67, X30, X32). Approximately 95% of HRI diagnoses were made between May and September. Patients hospitalized with BHD and a concurrent HRI were compared to patients hospitalized with only BHD and no HRI during the warmer months. Among hospitalizations with BHD and a concurrent CRI diagnosis, CRI was classified as hypothermia, frostbite with tissue necrosis, trench foot, chilblains, other effects of reduced temperature, accident due to excessive cold due to weather conditions, and injury by extremes of cold (ICD-9: 991, E901.0, E901.9, E988.3; ICD-10: T68, T69, X31, T33, T34). Approximately 82% of CRI diagnoses were made between October and April. Patients hospitalized with BHD and a concurrent CRI were compared to patients hospitalized with only BHD and no CRI during the colder months.

2.3. Covariates

2.3.1. Individual Level Covariates

These included demographics (age, sex, and race/ethnicity), health insurance and concurrent comorbidities. We categorized insurance payor into Medicare, Medicaid, private/HMO, uninsured, or other (Worker’s Compensation, Civilian Health, and Medical Program of the Uniformed Service (CHAMPUS), Title V, other government programs, and non-government programs). Comorbidities for TRIs include diabetes, cardiovascular disease, respiratory disease, and obesity (See Table S2 for comorbidities crosswalk) [20,26,27].

2.3.2. Hospital and Area Level Covariates

The hospital size variable describes urbanicity and health resources. It was developed for this study using a classification from the Healthcare Cost and Utilization Project (HCUP) Statistical Brief using number of hospital beds, location, and teaching status to determine size [28]. The hospital location variable describes urbanicity of the census tract it lies within and was developed using the Rural-Urban Commuting Area (RUCA) codes—Classification B. RUCA codes provide a classification of U.S. census tracts using measures of population density, urbanization, and daily commuting [29]. The most recent RUCA codes in use are based on the 2010 decennial census and the 2006–2010 American Community Survey. Schmeltz and Gamble included an income variable in their analysis using Census tract income quartiles based on residential ZIP codes. We did not create an income variable from Census tract income quartiles since patients with address information missing in SPARCS would have been excluded from complete case analysis. Given the high prevalence of BHD and exposure to extreme temperatures in homeless populations, the exclusion of patients with no addresses would not be useful to this analysis. SPARCS does include a homelessness variable; however, it was considered unreliable since there was a significant difference between the ICD-9 (5.9%) and ICD-10 (0.4%) classification years (p < 0.0001), and the prevalence of homelessness in the BHD population is expected to be higher.

2.4. Statistical Analysis

A univariate analysis was performed on all individual and hospital level variables; frequencies and distributions of all relevant variables were assessed for all BHD hospitalizations, for BHD with concurrent CRI, and for BHD with concurrent HRI. Subsequently, a bivariate analysis using the Chi-squared test was conducted to assess demographic differences between hospitalizations with and without the outcome of interest. Frequency counts and percentages of BHD with concurrent CRI and HRI were also conducted. A bivariate log-binomial regression was conducted to estimate relative risks (RR) of specific mental health diagnoses for BHD hospitalization with a concurrent CRI or HRI [30,31,32,33]. Cross-sectional research with a binary outcome such as this study often uses logistic regression; however, using log-binomial regression avoids overestimation of the prevalence odds ratio by directly estimating the prevalence ratio [30]. Where BHD counts were less than 15, we ran exact logistic regressions to confirm the findings. The results were very similar to the log-binomial RRs, with incrementally tighter precision. Subsequently, multivariable log-binomial regression was conducted using predictive models and complete case analysis to assess the individual and regional risk factors for hospitalization with BHD and a concurrent TRI.
log(πTRI) = β0 + β1 Age + β2 Sex + β3 Race + β4 Payor + β5 Comorbid Conditions + β6 Hospital Size + β7 Hospital Location
Statistical analyses for this study were performed using SAS software (SAS 9.4, SAS Institute, Cary, NC, USA).

3. Results

NYS is located in the Northeastern part of the United States and experiences wide variations in regional temperatures. The temperature distribution during the warm (May through September) and cold months (October through April) of the study period (2005–2019) varied widely (Table S3). In summary, the average maximum temperature for the warm months in NYS was 24.72 °C (76.50 °F) while the cold months was 7.99 °C (46.38 °F). The range of maximum temperature during the study period was from −0.67 °C to 39.21 °C (30.79 °F to 102.58 °F) for warm months and −25.06 °C (−13.11 °F) through 31.63 °C (88.93 °F) for the cold months. The average maximum heat index during the warm months was 25.53 °C (77.95 °F) with a range of −0.67 °C to 50.28 °C (30.79 °F to 122.50 °F), while during the cold months, the maximum heat index ranged from −25.06 °C to 34.67 °C (−13.11 °F to 94.41 °F), with an average maximum heat index of 7.46 °C (45.43 °F).
The study population consists of 12,272,393 BHD hospitalizations from 2005 to 2019 among NYS residents, including those of New York City. There was a total of 18,872,852 diagnoses of the selected BHD in our analysis. About 62% of hospitalizations had one BHD diagnosis, 22% had two BHD diagnoses, and 11% had three diagnoses (Table S4). Table 1 provides a summary of the demographic distribution of the study population. The mean age in the BHD hospitalizations study population was 56 years, with about 41% of hospitalizations in the 40–64 age category. About 51% of the study population were female. The race and ethnicities of the study population were non-Hispanic white (57%), Hispanic/Latino (12%), non-Hispanic Black (19%), Asian/Pacific Islander (2%), and other (11%); other races included Native American, Non-Hispanic multiracial, and other races. The insurance payor status in the study population was Medicare (35%), Medicaid (37%), private insurance (19%), and uninsured (8%). The comorbidity frequency was diabetes (23%), cardiovascular disease (22%), respiratory disease (35%), and obesity (10%). The majority of hospitalizations were in large hospitals (54%) or in urban locations (92%) (Table 1).
Among BHD hospitalizations, significantly decreased risks of concurrent TRI were observed for several diagnoses including depressive disorders, anxiety disorders, and bipolar disorders compared to BHD hospitalizations without these diagnoses (Table 2). Multivariate log-binomial models characterizing the individual and regional risk factors for a concurrent TRI among BHD hospitalizations are detailed in Table 3 and Table 4.

3.1. BHD Hospitalizations with Concurrent HRIs

There were 1914 BHD hospitalizations with concurrent HRIs from 2005–2019 during the months of May through September; and this varied between 50–200 cases per year and did not exhibit a discernible trend over time during the study period (Figure S1a). The characteristics of the BHD hospitalizations with concurrent HRI are summarized in Table 1 and were significantly different from those with BHD alone for several variables including age, sex, race, insurance payer, and hospital size (p < 0.0001–0.01), (Table 1). Less than 2% of BHD hospitalizations with HRI and a comorbidity occurred in the same individual.
In Table 2, we examine the frequency and distribution of the selected BHD with concurrent TRI diagnoses. Of the 1914 BHD hospitalizations with a concurrent HRI, there were 2661 BHD diagnoses. Mental illness accounted for 48% of BHD with concurrent HRI diagnoses, including depressive disorders (18%), schizophrenia spectrum and other psychotic disorders (13%), and anxiety disorders (8%). Dementia accounted for 18% of BHD with concurrent HRI diagnoses. The other 34% of the diagnoses were substance-use disorders; of which the majority was accounted for by non-dependent drug use (70%). Log-binomial regression analysis showed that dementia (RR 1.65; 95% Confidence Interval (CI): 1.49, 1.83) and schizophrenia spectrum and other psychotic disorders (RR 1.24; 95% CI: 1.10, 1.39), including schizophrenia (RR 1.38; 95% CI:1.19, 1.60), were associated with an increased risk for HRI among BHD hospitalizations. In a sub-analysis of multiple diagnoses among BHD hospitalizations, we did not observe an increased risk of HRI in those with dementia and schizophrenia, dementia and substance-use disorders, and schizophrenia and substance-use disorders compared to those without these multiple diagnoses (Table S5).
Risk factors for hospitalization for BHD and concurrent HRI include age, specifically 40–64 years (RR 1.46; 95% CI: 1.27, 1.69), 65–74 years (RR 2.26; 95% CI: 1.87, 2.73), and 75 years and over (RR 2.62; 95% CI: 2.19, 3.13), male gender (RR 1.87; 95% CI: 1.69, 2.06), non-Hispanic Black race (RR 1.34; 95% CI: 1.18, 1.51) and medium hospital size (RR 1.18; 95% CI: 1.07, 1.32). Significantly decreased risks of hospitalization for MBD and concurrent HRI were associated with cardiovascular disease (RR 0.83; 95% CI: 0.74, 0.93) and respiratory disease (RR 0.86 95% CI: 0.78, 0.95). In the sub-analysis of hospitalization with a mental disorder and concurrent HRI, older age and non-Hispanic Black persisted as risk factors. In a sub-analysis of hospitalization with substance-use disorders and concurrent HRI, older age and male gender persisted as risk factors. Additionally, being uninsured, having other insurance status, and the presence of respiratory disease were found to be risk factors. A significantly decreased risk of hospitalization for substance-use disorders and concurrent HRI was associated with cardiovascular disease (RR 0.80; 95% CI: 0.66, 0.97) (Table 3).

3.2. BHD Hospitalizations with Concurrent CRIs

There were 7738 BHD hospitalizations with concurrent CRIs from 2005–2019 during the months of October to April. BHD hospitalizations with concurrent CRI diagnoses trended upward over the study period, increasing from an average of 400 to over 1000 per year, or from 5% to 11% (Figure S1b). The characteristics of those hospitalized with concurrent CRI are also summarized in Table 1 and were significantly different from those with BHD alone for all variables (p < 0.0001–0.05). In BHD hospitalizations, 12% with CRI and diabetes, 8% with CRI and respiratory disease, 7% with CRI and CVD, and 2% with CRI and obesity occurred in the same individual with two to three hospitalizations over the study period.
Of the 7738 hospitalizations of BHD with a concurrent CRI, there were 11,449 BHD diagnoses. Mental illness accounted for 44% of BHD with concurrent CRI diagnoses, including depressive disorders (15%), schizophrenia spectrum and other psychotic disorders (12%), and anxiety disorders (7%). Dementia accounted for 16% of BHD with concurrent CRI diagnoses. The other 40% of the diagnoses were substance-use disorders, primarily non-dependent drug use (56%). A log-binomial regression analysis showed that dementia (RR 1.52; 95% CI:1.44, 1.60), schizophrenia spectrum and other psychotic disorders (RR 1.22; 95% CI: 1.15, 1.30), which included schizophrenia (RR 1.41; 95% CI:1.31, 1.52) were associated with an increased risk of CRI among BHD hospitalizations. Substance-use disorders (RR 1.62; 95% CI:1.55, 1.70), which included alcohol dependence (RR 2.10; 95% CI:1.99, 2.22), and non-dependent drug/alcohol use (RR 1.20; 95% CI:1.15, 1.26) were also associated with an increased risk of CRI among BHD hospitalizations. In a sub-analysis of multiple diagnoses among BHD hospitalizations, we observed an increased risk of concurrent CRI in those with dementia and schizophrenia, dementia and substance-use disorders, and schizophrenia and substance-use disorders compared to those without these multiple diagnoses (Table S5).
Risk factors for hospitalization for BHD and concurrent CRI include age, specifically 40–64 years (RR 1.53; 95% CI: 1.43, 1.64), 65–74 years (RR 1.55; 95% CI: 1.40, 1.71), and 75 years and over (RR 1.78; 95% CI: 1.63, 1.95), male gender (RR 2.06; 95% CI: 1.96, 2.17), non-Hispanic Black race (RR 1.23; 95% CI: 1.15, 1.30), being on Medicare (RR 1.51; 95% CI: 1.38, 1.64), Medicaid (RR 1.71; 95% CI:1.59, 1.85), or uninsured status (RR 1.76; 95% CI: 1.59, 1.94), the presence of respiratory disease (RR 1.08; 95% CI: 1.03, 1.14), and hospitalization in a rural location, including a large rural city or town (RR 1.24; 95% CI: 1.13,1.36), or a small rural town (RR 1.54; 95% CI:1.35, 1.75). Significantly decreased risks of hospitalization for BHD and concurrent CRI were associated with diabetes (RR 0.89; 95% CI 0.84, 0.95), cardiovascular disease (RR 0.68; 95% CI: 0.64, 0.72) and obesity (RR 0.49; 95% CI: 0.43, 0.54). Sub-analysis of the risk factors for hospitalization with a mental disorder and concurrent CRI showed similar risk associations to the full model. Sub-analysis of risk factors for hospitalization with substance-use disorders and concurrent CRI included the following risk factors: being aged 5–17 years and increasing age, being on Medicare, Medicaid, or having uninsured status, and the presence of respiratory disease. Significantly decreased risks of hospitalization for substance-use disorders and concurrent CRI were associated with diabetes, cardiovascular disease, and obesity (Table 4).

4. Discussion

4.1. Mental Health Disorders

In the population of BHD hospitalizations in NYS, several mental health diagnoses and dementia were found to be risk factors for concurrent TRI. BHD hospitalized patients with dementia diagnoses were observed to be more likely to have a concurrent TRI compared to those without dementia diagnoses. BHD hospitalized patients with schizophrenia spectrum diagnoses and other psychotic disorders, particularly schizophrenia, were observed to be more likely to have a concurrent TRI compared to those without these diagnoses. Similar to findings in this study, Schmeltz and Gamble found dementia to be a risk factor for HRI [11]. People who have dementia or schizophrenia can have impaired decision making, communication, and comprehension. These factors may contribute to failure to stay adequately hydrated in extreme heat events, to dress appropriately for the weather, to seek shelter, to recognize the signs of TRI, and be vulnerable to prolonged exposure to extreme weather conditions. Dementia most commonly presents in older populations who are at increased risk of hyperthermia, hypothermia, and other TRI due to age related processes such as decreased metabolic rate, diminished sweating capacity, and inability to thermoregulate [34,35]. Additionally, hyperthermia and hypothermia are both known side-effects of antipsychotic medication used in the treatment of schizophrenia and dementia due to thermoregulatory disturbance [21,36,37].

4.2. Substance Abuse Disorders

Among BHD hospitalizations, patients with substance-use disorder diagnoses, alcohol dependence and non-dependent drug or alcohol use were observed to be more likely to have a concurrent CRI compared to patients without these diagnoses. Alcohol and drugs are known to impair mental state [35], behavioral responses, and can increase the risk of cold exposure [35]. Several substances, including alcohol, cocaine, and other stimulants have been found to affect thermoregulation [38,39].

4.3. Demographics

Among those hospitalized for BHD, older age groups, males, and non-Hispanic Blacks were associated with increased risk for a concurrent HRI or CRI. The HRI findings support those observed by Schmeltz and Gamble in their national research on risk factors for BHD with concurrent HRI [11]. Older age is a known risk factor for TRI, due to thermoregulatory deterioration and existing morbidities [34,35,40]. Older age remained a risk factor for concurrent TRI in the mental illness sub-analysis and for concurrent CRI in the substance-use disorders sub-analyses. Males with BHD hospitalization were observed to have an increased risk of concurrent TRI compared to females with BHD hospitalization; and the risk remained in the sub-analyses of cases diagnosed with mental illness and substance-use disorders. Males have previously been found to be at higher risk of hyperthermia and hypothermia compared to females [41,42]. The disparity in TRI risk may be attributed to a difference in risk-taking behavior, substance-use, occupational exposure, and utilization of healthcare services [11,43,44].
In the study population, having an additional concurrent diagnosis of respiratory disease was another risk factor for concurrent CRI compared to those BHD hospitalized without respiratory disease, and the risk remained in sub-analyses of both mental and substance-use disorders. Patients hospitalized with substance-use disorders had an increased risk of concurrent HRI if they had a concurrent respiratory infection. Respiratory diseases are found to be more prevalent in homeless populations and highly transmissible in homeless shelters [45,46]. Previous research also supports an association between exposure to extreme cold and increased risk of respiratory infection [20,47]. Unfortunately, in our analysis, data on the housing status of patients was not available, but future analyses on this risk factor relative to housing status is warranted. Protective effects for hospitalization with BHD and concurrent TRI associated with comorbidities such as diabetes, cardiovascular disease, and obesity need to be explored further, but may be due to the fact that those who are sick with multiple comorbidities have a lower risk of exposure, as they may not be outdoors for long periods of time. The observed decreased risk could also be due to the underreporting of secondary diagnoses, where comorbidity diagnoses take precedence in treatment and are recorded and reported instead of TRI diagnoses.
In the population of people hospitalized with BHD, results indicate that Non-Hispanic Blacks had a higher risk of concurrent TRI compared with Whites hospitalized with BHD, and these findings persisted in the mental illness sub-analysis. Research on Medicare claims administrative data found that Black Medicare beneficiaries had significantly higher risk of inpatient and outpatient visits for hyperthermia and hypothermia [42]. Non-Hispanic Blacks may be at increased risk of BHD with concurrent TRI due to demographic disparities in the social determinants of health, such as differential environmental conditions including housing and neighborhood characteristics and access to heat and air-conditioning [48,49]. Non-Hispanic Blacks experience disparities in access and utilization of health care services and mental health medications, resulting in higher rates of inpatient care for mental health services [50]. Non-Hispanic Blacks with mental illness are more likely to be prescribed anti-psychotics in general, as well as at higher doses [51]. Also, Non-Hispanic Blacks are disproportionately affected by homelessness and therefore exposure to extremes of temperature [52].

4.4. Insurance/Access to Care

We observed that among BHD hospitalizations, being on Medicare, Medicaid, and uninsured status were risk factors for concurrent CRI compared to BHD hospitalized patients with private insurance, and this risk persisted in sub-analyses of both mental and substance-use disorders. Patients hospitalized with substance-use disorders had an increased risk of concurrent HRI if they were uninsured or had an “other“ insurance status compared to those with private insurance. Medicare recipients who are 65 years and over are therefore a vulnerable population for CRI. Patients with low socioeconomic status are more likely to be uninsured or on Medicaid and homeless people have a low prevalence of insurance coverage due to barriers to healthcare access [53,54,55]. These vulnerable populations have decreased access to healthcare, resources, and potentially shelter, heating, and air-conditioning. Although we were not able to evaluate homelessness in our study, BHD are prevalent in homeless populations where barriers to healthcare and resource access, particularly shelter and air conditioning during extreme temperature events, are common [56,57]. These factors render the homeless population extremely vulnerable to TRIs, underlining the importance of public health responses when extreme temperatures are expected, especially in extreme cold events where mortality is even greater [56]. For example, a study conducted in New York City observed that approximately one-fourth of those hospitalized for CRI were homeless, and over one-third of those who died from CRI were homeless [20]. In our study population, several diagnoses had decreased risks of hospitalization with a concurrent TRI, such as depressive disorder, bipolar disorders, and anxiety disorders. This decreased risk may possibly be attributed to individuals with these BHD being more likely to stay indoors and not be exposed to extreme temperatures from factors such as impaired decision making, communication, and comprehension. Further research is needed to explore these associations.

4.5. Hospital Size & Location

In the study population, concurrent HRI hospitalizations were most likely to be in medium size hospitals compared to large hospitals, while concurrent CRI hospitalization was more likely in rural hospital locations, especially small rural towns, compared to urban BHD hospitalizations. Schmeltz and Gamble similarly observed that BHD patients in medium sized hospitals in the south and midwestern United States had higher risk of hospitalization with concurrent HRI [11]. Schemltz and Gamble also found a higher risk of concurrent HRI in rurally located hospitals; however, in our study, hospital location did not play a significant role in the risk of concurrent HRI among BHD cases. Further analysis could help determine the factors that may contribute to an increased risk of HRI associated with medium sized hospitals. Potential contributary factors to increased CRI in rural areas may be a decreased level of access to health resources and shelter, in addition to a greater likelihood of engagement in outdoor activities. Additionally, in NYS, most rural areas are upstate, which is much colder than downstate.

4.6. Strengths & Limitations

Many previous studies focus on the association between temperature and BHD morbidity and mortality, but to our knowledge this is the first study to assess risk factors in BHD hospitalizations for both concurrent heat- and cold-related illnesses in NYS. An important strength of this study is the use of administrative data specific to NYS to perform an in-depth exploration of several BHD. This allowed us to make NYS-specific individual-level analyses and inferences using a large sample size. Cross sectional research such as this is also useful for establishing preliminary findings and creating a direction for future research.
This study has some limitations. Firstly, we are unable to attribute any adverse temperature-related health outcome to specific heat/cold events. We used administrative billing data to define temperature-related hospitalizations. As such, we may only have captured cases severe enough to need hospitalization. Due to differential data entry from provider to provider, there is also the potential for misclassification. In addition, the underreporting of secondary diagnoses can lead to misclassification, which may result in the observed paradoxical decreased risk from comorbidities that are known to increase the risk of TRI [58,59]. We did not have a reliable indicator of homelessness in the SPARCS dataset. Without a clear indication of homelessness in the study population, we had to make assumptions based on the known prevalence of mental health disorders and substance-use issues in this population. An additional limitation, not specific to our study, is the change in ICD-9 to ICD-10 classification systems that led to differential diagnostic coding, particularly with substance-use disorders, and schizophrenia and other psychotic disorders. From the results of a log-binomial regression sensitivity analysis for substance-use disorders and schizophrenia and other psychotic disorders, we observed that the direction of associations for the risk factors were not affected before and after 2015, when the shift from ICD-9 to ICD-10 coding in SPARCS occurred (results not shown) indicating that the difference in classification did not affect our risk factor analyses. Finally, 18% of CRI cases occurred outside of the selected cold months of October–May. We limited the analysis to season periods to avoid overlap but recommend that future studies use data for the whole year to capture all cases.

5. Conclusions

Our results agree with prior risk factor findings for BHD hospitalization with a concurrent HRI. This study adds to the literature by identifying dementia, schizophrenia, substance-use disorders, including alcohol dependence and non-dependent substance-use, increasing age, male gender, non-Hispanic Black and other race, Medicare, Medicaid, uninsured status, and respiratory disease as risk factors for BHD hospitalization with a concurrent CRI. Study results can be used to aid public health initiatives from state and local health departments and the community to develop TRI adaptation strategies tailored to the needs of the most vulnerable populations. TRIs are preventable health events that can benefit from targeted programs aimed at preparing the elderly for extreme temperature events, providing shelter for homeless populations, and assisting vulnerable groups with information and resources during extreme temperature events.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph192416411/s1, Figure S1a: Trends in Behavioral Health Disorders with Concurrent Heat-Related Illness; Figure S1b: Trends in Behavioral Health Disorders with Concurrent Cold-Related Illness; Table S1: Behavioral Health Disorders Diagnoses ICD-9/ICD-10 Crosswalk; Table S2: Temperature Related Illness Comorbidity ICD9/ICD10 Crosswalk; Table S3: Distribution of select climate parameters (2005–2019) in New York State; Table S4: Number of BHD Diagnoses Per Patient; Table S5: Hospitalization counts and risk ratios (RRs) of hospitalizations due to multiple primary or secondary diagnosis of behavioral health disorders (BHD) with a concurrent heat-related illness (HRI) (May-September) or cold-related illness (CRI) (October-April), 2005–2019; Table S6: Total Inpatient Hospitalizations in New York State (2005–2019).

Author Contributions

Conceptualization, H.A.-G.; methodology, H.A.-G., T.A. and T.Z.I.; validation, H.A.-G. and T.A.; formal analysis, H.A.-G.; data curation, H.A.-G.; writing—original draft preparation, H.A.-G.; writing—review and editing, T.A., N.M., S.N., S.S. and T.Z.I.; supervision, T.Z.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CDC grant funding provided to NYS as under the “Building Resilience Against Climate Effects: Implementing and Evaluating Adaptation Strategies that Protect and Promote Human Health” (project # NUE1EH001453) and “Enhancing Innovation and Capabilities of the Environmental Public Health Tracking Network” (project # NUE1EH001341) awards.

Institutional Review Board Statement

The New York State Department of Health Institutional Review Board (NYS DOH IRB) reviewed and approved this study with a waiver of participants’ informed consent on 22 October 2021 (Board Ref #: 21-016).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data for this study came from the NYS Statewide Planning and Research Cooperative System (SPARCS) database. This data is not publicly available as it would compromise individual privacy.

Acknowledgments

The authors would like to extend their thanks to Marc Manseau, Emily Leckman-Westin, Joshua Berezin, and Molly Finnerty from the New York State Office of Mental Health for their expert advice on this study. The authors would also like to thank Amandus Ankobil for his literature review on temperature-related illness and behavioral health disorders.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NYSERDA. Responding to Climate Change in New York State (ClimAID); New York State Research and Development Authority: Albany, NY, USA, 2014.
  2. Berko, J.; Ingram, D.D.; Saha, S.; Parker, J.D. Deaths Attributed to Heat, Cold, and Other Weather Events in the United States, 2006–2010. Natl. Health Stat. Rep. 2014, 30, 1–15. [Google Scholar]
  3. Gasparrini, A.; Guo, Y.; Hashizume, M.; Lavigne, E.; Zanobetti, A.; Schwartz, J.; Tobias, A.; Tong, S.; Rocklov, J.; Forsberg, B.; et al. Mortality Risk Attributable to High and Low Ambient Temperature: A Multicountry Observational Study. Lancet 2015, 386, 369–375. [Google Scholar] [CrossRef] [PubMed]
  4. Adeyeye, T.E.; Insaf, T.Z.; Al-Hamdan, M.Z.; Nayak, S.G.; Stuart, N.; DiRienzo, S.; Crosson, W.L. Estimating Policy-Relevant Health Effects of Ambient Heat Exposures Using Spatially Contiguous Reanalysis Data. Environ. Health 2019, 18, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. CDC Heat Related Illness. Available online: https://www.cdc.gov/niosh/topics/heatstress/heatrelillness.html (accessed on 8 July 2021).
  6. Herity, B.; Daly, L.; Bourke, G.J.; Horgan, J.M. Hypothermia and Mortality and Morbidity. An Epidemiological Analysis. J. Epidemiol. Community Health 1991, 45, 19–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. CDC Cold Related Illnesses. Available online: https://www.cdc.gov/niosh/topics/coldstress/coldrelatedillnesses.html (accessed on 5 August 2021).
  8. Chen, N.T.; Lin, P.H.; Guo, Y.L. Long-term Exposure to High Temperature Associated with the Incidence of Major Depressive Disorder. Sci. Total Environ. 2019, 659, 1016–1020. [Google Scholar] [CrossRef]
  9. Sherbakov, T.; Malig, B.; Guirguis, K.; Gershunov, A.; Basu, R. Ambient Temperature and Added Heat Wave Effects on Hospitalizations in California from 1999 to 2009. Environ. Res. 2018, 160, 83–90. [Google Scholar] [CrossRef]
  10. Trang, P.M.; Rocklov, J.; Giang, K.B.; Kullgren, G.; Nilsson, M. Heatwaves and Hospital Admissions for Mental Disorders in Northern Vietnam. PLoS ONE 2016, 11, e0155609. [Google Scholar] [CrossRef]
  11. Schmeltz, M.T.; Gamble, J.L. Risk Characterization of Hospitalizations for Mental Illness and/or Behavioral Disorders with Concurrent Heat-Related Illness. PLoS ONE 2017, 12, e0186509. [Google Scholar] [CrossRef]
  12. Peng, Z.; Wang, Q.; Kan, H.; Chen, R.; Wang, W. Effects of Ambient Temperature on Daily Hospital Admissions for Mental Disorders in Shanghai, China: A Time-Series Analysis. Sci. Total Environ. 2017, 590–591, 281–286. [Google Scholar] [CrossRef]
  13. Liu, X.; Liu, H.; Fan, H.; Liu, Y.; Ding, G. Influence of Heat Waves on Daily Hospital Visits for Mental Illness in Jinan, China-A Case-Crossover Study. Int. J. Environ. Res. Public Health 2018, 16, 87. [Google Scholar] [CrossRef] [Green Version]
  14. Lee, S.; Lee, H.; Myung, W.; Kim, E.J.; Kim, H. Mental Disease-Related Emergency Admissions Attributable to Hot Temperatures. Sci. Total Environ. 2018, 616–617, 688–694. [Google Scholar] [CrossRef]
  15. Basu, R.; Gavin, L.; Pearson, D.; Ebisu, K.; Malig, B. Examining the Association Between Apparent Temperature and Mental Health-Related Emergency Room Visits in California. Am. J. Epidemiol. 2018, 187, 726–735. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, X.; Lavigne, E.; Ouellette-kuntz, H.; Chen, B.E. Acute Impacts of Extreme Temperature Exposure on Emergency Room Admissions Related to Mental and Behavior Disorders in Toronto, Canada. J. Affect. Disord. 2014, 155, 154–161. [Google Scholar] [CrossRef]
  17. Marzuk, P.M.; Tardiff, K.; Leon, A.C.; Hirsch, C.S.; Portera, L.; Iqbal, M.I.; Nock, M.K.; Hartwell, N. Ambient Temperature and Mortality from Unintentional Cocaine Overdose. JAMA 1998, 279, 1795–1800. [Google Scholar] [CrossRef]
  18. Bohnert, A.S.; Prescott, M.R.; Vlahov, D.; Tardiff, K.J.; Galea, S. Ambient Temperature and Risk of Death from Accidental Drug Overdose in New York City, 1990–2006. Addiction 2010, 105, 1049–1054. [Google Scholar] [CrossRef] [Green Version]
  19. Auger, N.; Bilodeau-Bertrand, M.; Labesse, M.E.; Kosatsky, T. Association of Elevated Ambient Temperature with Death from Cocaine Overdose. Drug Alcohol Depend. 2017, 178, 101–105. [Google Scholar] [CrossRef]
  20. Lane, K.; Ito, K.; Johnson, S.; Gibson, E.A.; Tang, A.; Matte, T. Burden and Risk Factors for Cold-Related Illness and Death in New York City. Int. J. Environ. Res. Public Health 2018, 15, 632. [Google Scholar] [CrossRef] [Green Version]
  21. Zonnenberg, C.; Bueno-de-Mesquita, J.M.; Ramlal, D.; Blom, J.D. Hypothermia due to Antipsychotic Medication: A Systematic Review. Front. Psychiatry 2017, 8, 165. [Google Scholar] [CrossRef] [Green Version]
  22. Mullins, J.T.; White, C. Temperature and Mental Health: Evidence from the Spectrum of Mental Health Outcomes. J. Health Econ. 2019, 68, 102240. [Google Scholar] [CrossRef] [Green Version]
  23. Chan, E.Y.Y.; Lam, H.C.Y.; So, S.H.W.; Goggins, W.B.; Ho, J.Y.; Liu, S.; Chung, P.P.W. Association between Ambient Temperatures and Mental Disorder Hospitalizations in a Subtropical City: A Time-Series Study of Hong Kong Special Administrative Region. Int. J. Environ. Res. Public Health 2018, 15, 754. [Google Scholar] [CrossRef] [Green Version]
  24. CMS ICD-9-CM to and from ICD-10-CM and ICD-10-PCS Crosswalk or General Equivalence Mappings. Available online: https://www.nber.org/research/data/icd-9-cm-and-icd-10-cm-and-icd-10-pcs-crosswalk-or-general-equivalence-mappings (accessed on 15 November 2021).
  25. Wei, Y.; Wang, Y.; Lin, C.-K.; Yin, K.; Yang, J.; Shi, L.; Li, L.; Zanobetti, A.; Schwartz, J.D. Associations between Seasonal Temperature and Dementia-Associated Hospitalizations in New England. Environ. Int. 2019, 126, 228–233. [Google Scholar] [CrossRef] [PubMed]
  26. Kenny, G.P.; Yardley, J.; Brown, C.; Sigal, R.J.; Jay, O. Heat Stress in Older Individuals and Patients with Common Chronic Diseases. CMAJ Can. Med. Assoc. J. J. De L’association Med. Can. 2010, 182, 1053–1060. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Rawshani, A.; Rawshani, A.; Franzén, S.; Sattar, N.; Eliasson, B.; Svensson, A.M.; Zethelius, B.; Miftaraj, M.; McGuire, D.K.; Rosengren, A.; et al. Risk Factors, Mortality, and Cardiovascular Outcomes in Patients with Type 2 Diabetes. N. Engl. J. Med. 2018, 379, 633–644. [Google Scholar] [CrossRef] [PubMed]
  28. Tian, W. An All-Payer View of Hospital Discharge to Postacute Care, 2013: Statistical Brief #205. In Healthcare Cost and Utilization Project (HCUP) Statistical Briefs; Agency for Healthcare Research and Quality (US): Rockville, MD, USA, 2006. [Google Scholar]
  29. Morrill, R.; Cromartie, J.; Hart, G. Metropolitan, Urban, and Rural Commuting Areas: Toward a Better Depiction of the United States Settlement System. Urban Geogr. 1999, 20, 727–748. [Google Scholar] [CrossRef] [Green Version]
  30. Barros, A.J.; Hirakata, V.N. Alternatives for Logistic Regression in Cross-Sectional Studies: An Empirical Comparison of Models that Directly Estimate the Prevalence Ratio. BMC Med. Res. Methodol. 2003, 3, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Spiegelman, D.; Hertzmark, E. Easy SAS Calculations for Risk or Prevalence Ratios and Differences. Am. J. Epidemiol. 2005, 162, 199–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Richardson, D.B.; Kinlaw, A.C.; MacLehose, R.F.; Cole, S.R. Standardized Binomial Models for Risk or Prevalence Ratios and Differences. Int. J. Epidemiol. 2015, 44, 1660–1672. [Google Scholar] [CrossRef] [Green Version]
  33. Petersen, M.R.; Deddens, J.A. A Comparison of Two Methods for Estimating Prevalence Ratios. BMC Med. Res. Methodol. 2008, 8, 9. [Google Scholar] [CrossRef] [Green Version]
  34. Millyard, A.; Layden, J.D.; Pyne, D.B.; Edwards, A.M.; Bloxham, S.R. Impairments to Thermoregulation in the Elderly During Heat Exposure Events. Gerontol. Geriatr. Med. 2020, 6, 2333721420932432. [Google Scholar] [CrossRef]
  35. Hypothermia-Related Mortality--Montana, 1999-2004. MMWR. Morb. Mortal. Wkly. Rep. 2007, 56, 367–368.
  36. Hermesh, H.; Shiloh, R.; Epstein, Y.; Manaim, H.; Weizman, A.; Munitz, H. Heat Intolerance in Patients with Chronic Schizophrenia Maintained with Antipsychotic drugs. Am. J. Psychiatry 2000, 157, 1327–1329. [Google Scholar] [CrossRef]
  37. Best, K.; Tran, D.H.; Bulte, C.; Verceles, A.C.; Michael, M.B. Antipsychotic Medication-Induced Hyperthermia Leading to Cerebrovascular Accident: A Case Report. Cureus 2021, 13, e18651. [Google Scholar] [CrossRef]
  38. Sprague, J.E.; Riley, C.L.; Mills, E.M. Body Temperature Regulation and Drugs of Abuse. Handb. Clin. Neurol. 2018, 157, 623–633. [Google Scholar]
  39. Van Amsterdam, J.; Brunt, T.M.; Pierce, M.; van den Brink, W. Hard Boiled: Alcohol Use as a Risk Factor for MDMA-Induced Hyperthermia: A Systematic Review. Neurotox. Res. 2021, 39, 2120–2133. [Google Scholar] [CrossRef]
  40. Nayak, S.G.; Shrestha, S.; Kinney, P.L.; Ross, Z.; Sheridan, S.C.; Pantea, C.I.; Hsu, W.H.; Muscatiello, N.; Hwang, S.A. Development of a Heat Vulnerability Index for New York State. Public Health 2018, 161, 127–137. [Google Scholar] [CrossRef]
  41. MMWR Heat-Related Deaths—United States, 1999–2003. Available online: https://www.cdc.gov/mmwr/preview/mmwrhtml/mm5529a2.htm (accessed on 6 July 2022).
  42. Noe, R.S.; Jin, J.O.; Wolkin, A.F. Exposure to Natural Cold and Heat: Hypothermia and Hyperthermia Medicare Claims, United States, 2004–2005. Am. J. Public Health 2012, 102, e11–e18. [Google Scholar] [CrossRef]
  43. Centers for Disease Control and Prevention. Hypothermia-Related Deaths--United States, 1999-2002 and 2005. MMWR. Morb. Mortal. Wkly. Rep. 2006, 55, 282–284. [Google Scholar]
  44. McHugh, R.K.; Votaw, V.R.; Sugarman, D.E.; Greenfield, S.F. Sex and Gender Differences in Substance Use Disorders. Clin. Psychol. Rev. 2018, 66, 12–23. [Google Scholar] [CrossRef]
  45. Plevneshi, A.; Svoboda, T.; Armstrong, I.; Tyrrell, G.J.; Miranda, A.; Green, K.; Low, D.; McGeer, A. Population-Based Surveillance for Invasive Pneumococcal Disease in Homeless Adults in Toronto. PLoS ONE 2009, 4, e7255. [Google Scholar] [CrossRef] [Green Version]
  46. Burki, T.K. Homelessness and Respiratory Disease. Lancet. Respir. Med. 2013, 1, 767–768. [Google Scholar] [CrossRef]
  47. Mourtzoukou, E.G.; Falagas, M.E. Exposure to Cold and Respiratory Tract Infections. Int. J. Tuberc. Lung Dis. Off. J. Int. Union Against Tuberc. Lung Dis. 2007, 11, 938–943. [Google Scholar]
  48. Yardley, J.E.; Sigal, R.J.; Kenny, G.P. Heat Health Planning: The Importance of Social and Community Factors. Glob. Environ. Change-Hum. Policy Dimens. 2011, 21, 670–679. [Google Scholar] [CrossRef]
  49. Dialesandro, J.; Brazil, N.; Wheeler, S.; Abunnasr, Y. Dimensions of Thermal Inequity: Neighborhood Social Demographics and Urban Heat in the Southwestern U.S. Int. J. Environ. Res. Public Health 2021, 18, 941. [Google Scholar] [CrossRef] [PubMed]
  50. APA Mental Health Disparities: African Americans. Available online: https://www.psychiatry.org/psychiatrists/cultural-competency/education/mental-health-facts (accessed on 25 February 2022).
  51. Cook, T.B.; Reeves, G.M.; Teufel, J.; Postolache, T.T. Persistence of Racial Disparities in Prescription of First-Generation Antipsychotics in the USA. Pharmacoepidemiol. Drug Saf. 2015, 24, 1197–1206. [Google Scholar] [CrossRef]
  52. Henry, M. Azim Shivji, Tanya de Sousa, and Rebecca Cohen, Abt Associates Inc The 2015 Annual Homeless Assessment Report (AHAR) to Congress. Available online: https://www.huduser.gov/portal/sites/default/files/pdf/2015-AHAR-Part-1.pdf (accessed on 16 June 2022).
  53. Fryling, L.R.; Mazanec, P.; Rodriguez, R.M. Barriers to Homeless Persons Acquiring Health Insurance Through the Affordable Care Act. J. Emerg. Med. 2015, 49, 755–762.e2. [Google Scholar] [CrossRef]
  54. Smith, J.; Medalia, C. Health Insurance Coverage in the United States: 2014; US Census Bureau: Washington, DC, USA, 2015; pp. 60–253.
  55. Baggett, T.P.; O’Connell, J.J.; Singer, D.E.; Rigotti, N.A. The unmet health care needs of homeless adults: A national study. Am. J. Public Health 2010, 100, 1326–1333. [Google Scholar] [CrossRef]
  56. Zhang, P.; Wiens, K.; Wang, R.; Luong, L.; Ansara, D.; Gower, S.; Bassil, K.; Hwang, S.W. Cold Weather Conditions and Risk of Hypothermia Among People Experiencing Homelessness: Implications for Prevention Strategies. Int. J. Environ. Res. Public Health 2019, 16, 3259. [Google Scholar] [CrossRef] [Green Version]
  57. Gutwinski, S.; Schreiter, S.; Deutscher, K.; Fazel, S. The prevalence of mental disorders among homeless people in high-income countries: An updated systematic review and meta-regression analysis. PLoS Med. 2021, 18, e1003750. [Google Scholar] [CrossRef]
  58. Johnson, E.K.; Nelson, C.P. Values and Pitfalls of the Use of Administrative Databases for Outcomes Assessment. J. Urol. 2013, 190, 17–18. [Google Scholar] [CrossRef]
  59. Virnig, B. Strengths and Limitations of CMS Administrative Data in Research. Available online: https://resdac.org/articles/strengths-and-limitations-cms-administrative-data-research (accessed on 25 May 2022).
Table 1. Patient and hospital characteristics for all hospitalizations with behavioral health disorders (BHD), BHD and heat-related illnesses (HRI) (May-September), and BHD with cold-related illnesses (CRI) (October–April), 2005–2019.
Table 1. Patient and hospital characteristics for all hospitalizations with behavioral health disorders (BHD), BHD and heat-related illnesses (HRI) (May-September), and BHD with cold-related illnesses (CRI) (October–April), 2005–2019.
BHDBHD and HRI 1BHD and CRI 2
Total N12,272,39319147738
Age in years, Mean (SD)56 (20.9)61.3 (19.5)58.9 (20.0)
Age Categories in years, %
5–17309,307 (2.5)24 (1.3)113 (1.5)
18–392,481,211 (20.2)269 (14.1)1205 (15.6)
40–645,070,216 (41.3)758 (39.6)3517 (45.5)
65–741,579,414 (12.9)303 (15.8)988 (12.8)
75+2,832,245 (23.1)560 (29.3)1915 (24.8)
Sex, %
Male6,063,067 (49.4)1213 (63.4)5143 (66.5)
Female6,209,140 (50.6)701 (36.6)2595 (33.5)
Missing186 (.002)--
Race, %
Non-Hispanic White6,934,606 (56.6)1015 (53.2)4157 (53.8)
Hispanic/Latino1,454,637 (11.9)229 (12.0)743 (9.6)
Non-Hispanic Black2,303,695 (18.8)417 (21.9)1799 (23.3)
Asian/Pacific Islander200,237 (1.6)31 (1.6)107 (1.4)
Other 31,352,740 (11.1)216 (11.3)918 (11.9)
Missing26,478 (0.2)6 (0.3)14 (0.1)
Payer, %
Medicare4,040,182 (34.7)717 (39.3)2644 (36.0)
Medicaid4,341,858 (37.2)632 (34.6)3129 (42.6)
Private/HMO2,258,459 (19.4)298 (16.3)879 (11.9)
Uninsured894,255 (7.7)155 (8.5)645 (8.8)
Other 4124,859 (1.1)24 (1.3)48 (0.7)
Missing612,780 (5.0)88 (4.6)393 (5.1)
Comorbidities 5, %
Diabetes2,763,430 (22.5)435 (22.7)1562 (20.2)
Cardiovascular Disease2,706,675 (22.1)637 (22.8)1406 (18.2)
Respiratory Disease4,328,933 (35.3)604 (31.6)2886 (37.3)
Obesity1,239,795 (10.1)157 (8.2)343 (4.4)
Hospital Size 6, %
Small2,390,623 (19.5)347 (18.1)1390 (18.0)
Medium3,261,862 (26.6)586 (30.6)2100 (27.1)
Large6,619,908 (53.9)981 (51.3)4248 (54.9)
Hospital Location 7, %
Urban11,205,295 (91.8)1769 (92.4)6980 (90.2)
Large Rural City/Town722,465 (5.9)98 (5.1)504 (6.5)
Small Rural Town282,392 (2.3)47 (2.5)254 (3.2)
Missing62,241 (0.5)--
1 Heat-Related Illnesses (HRI) ICD-9 992, E900.0, E900.9; ICD-10, T67, X30, X32. 2 Cold-Related Illnesses (CRI) ICD-9 991, E901.0, E901.9, E988.3; ICD-10 T68, T69, X31, T33, T34. 3 Other Race includes Native American, Non-Hispanic multiracial, and other races. 4 Other includes Workers’ Compensation, CHAMPUS, Title V, other government programs, and non-government programs.5 Comorbidities do not add up to 100%; they are not mutually exclusive and not all patients have comorbidities. 6 Hospital Size Small: Rural 1–49 beds, Urban Non-Teaching 1–124 beds, Urban Teaching 1–249 beds; Medium: Rural 50–99 beds, Urban Non-Teaching 125–199 beds, Urban Teaching 250–424 beds; Large: Rural 100+ beds, Urban Non-Teaching 200+ beds, Urban Teaching 425+ beds. https://www.ncbi.nlm.nih.gov/books/NBK373736/table/sb205.t3/ (accessed on 14 December 2021). 7 Rural-urban commuting area (RUCA) codes were used to categorize rural and urban areas of NYS.
Table 2. Diagnosis counts and risk ratios (RRs) of hospitalizations due to a primary or secondary diagnosis of behavioral health disorders (BHD) with a concurrent heat-related illness (HRI) (May–September) or cold-related illness (CRI) (October–April), 2005–2019.
Table 2. Diagnosis counts and risk ratios (RRs) of hospitalizations due to a primary or secondary diagnosis of behavioral health disorders (BHD) with a concurrent heat-related illness (HRI) (May–September) or cold-related illness (CRI) (October–April), 2005–2019.
BHD 1HRI 2CRI 3
DescriptionCount%RR 495% CICount%RR 495% CI
Behavioral Health Disorder Diagnoses2661100 11,449100
Dementia47117.71.65(1.49, 1.83)188216.41.52(1.44, 1.60)
Depressive Disorders48918.40.63(0.57, 0.70)172015.00.53(0.50, 0.56)
Schizophrenia Spectrum and Other Psychotic Disorders35413.31.24(1.10, 1.39)140112.21.22(1.15, 1.30)
Schizophrenia19956.21.38(1.19, 1.60)81057.81.41(1.31, 1.52)
Bipolar and Related Disorders1214.50.77(0.64, 0.93)4674.10.76(0.69, 0.83)
Anxiety Disorders2158.10.53(0.46, 0.61)7456.50.45(0.41, 0.48)
Trauma and Stressor-Related Disorders511.90.66(0.50, 0.87)2252.00.73(0.64, 0.83)
Attention Deficit Hyperactivity Disorder250.90.81(0.54, 1.20)790.70.62(0.50, 0.78)
Disruptive, Impulse Control, and Conduct Disorders100.40.32(0.17, 0.60)1040.90.84(0.70, 1.02)
Personality Disorders110.40.19(0.10, 0.34)1471.30.63(0.54, 0.74)
Obsessive Compulsive and Related Disorders140.50.82(0.49, 1.39)760.71.16(0.93, 1.46)
Substance-Use Disorders90033.80.93(0.85, 1.02)460340.21.62(1.55, 1.70)
Alcohol Dependence13715.20.61(0.51, 0.72)155233.72.10(1.99, 2.22)
Drug Dependence10711.90.50(0.41, 0.61)4349.40.52(0.47, 0.57)
Non-Dependent Drug/Alcohol Use63270.21.10(1.00, 1.21)258256.11.20(1.15, 1.26)
1 ICD codes: See Crosswalk in Table S1. 2 Heat-Related Illnesses (HRI) ICD-9 992, E900.0, E900.9; ICD-10, T67, X30, X32. 3 Cold-Related Illnesses (CRI) ICD-9 991, E901.0, E901.9, E988.3; ICD-10 T68, T69, X31, T33, T34. 4 Risk ratios (relative risk) compared hospitalizations for BHD with concurrent TRI to hospitalizations for BHD without TRI, 2005–2019. The italics is category subsets.
Table 3. Multivariable models of risk factors for hospitalization for all behavioral health disorders, mental illness, and substance-use disorders with a concurrent heat-related illness (May–September), 2005–2019.
Table 3. Multivariable models of risk factors for hospitalization for all behavioral health disorders, mental illness, and substance-use disorders with a concurrent heat-related illness (May–September), 2005–2019.
CharacteristicsAll Behavioral Health DisordersMental Illness 1Substance-Use
Disorders
RR95% CIRR95% CIRR95% CI
Age Categories
5–17 0.85(0.56, 1.28)1.16(0.66, 2.03)1.45(0.68, 3.07)
18–391.006-1.00-1.00-
40–641.46(1.27, 1.69)2.27(1.80, 2.87)1.25(1.05, 1.48)
65–742.26(1.87, 2.73)3.49(2.63, 4.63)1.63(1.23, 2.16)
75+2.62(2.19, 3.13)3.38(2.56, 4.46)2.03(1.45, 2.83)
Sex
Male1.87(1.69, 2.06)1.71(1.50, 1.96)2.66(2.23, 3.17)
Female1.00-1.00-1.00-
Race
Non-Hispanic White1.00-1.00-1.00-
Hispanic/Latino1.14(0.98, 1.33)1.03(0.82, 1.29)1.18(0.96, 1.45)
Non-Hispanic Black1.34(1.18, 1.51)1.30(1.09, 1.56)1.17(0.98, 1.40)
Asian/Pacific Islander1.10(0.77, 1.59)1.16(0.68, 1.98)1.43(0.85, 2.41)
Other 21.08(0.92, 1.26)0.92(0.73, 1.17)1.06(0.85, 1.32)
Payer
Medicare0.97(0.83, 1.13)1.10(0.88, 1.37)0.94(0.74, 1.20)
Medicaid0.98(0.85, 1.13)1.10(0.90, 1.36)0.99(0.82, 1.19)
Private/HMO1.00-1.00-1.00-
Uninsured1.22(1.00 1.48)0.91(0.65, 1.27)1.45(1.13, 1.84)
Other 31.30(0.86, 1.97)0.61(0.25, 1.50)1.95(1.24, 3.08)
Comorbidities
Diabetes0.91(0.82, 1.03)1.05(0.90. 1.23)0.91(0.76, 1.09)
Cardiovascular Disease0.83(0.74, 0.93)0.89(0.75, 1.05)0.80(0.66, 0.97)
Respiratory Disease0.86(0.78. 0.95)0.90(0.78, 1.04)1.27(1.10, 1.47)
Obesity 1.00(0.84, 1.18)1.00(0.80, 1.24)1.03(0.80, 1.33)
Hospital Size 4
Small0.98(0.86, 1.12)0.90(0.75, 1.07)0.82(0.67, 1.01)
Medium1.18(1.07, 1.32)1.00(0.86, 1.17)1.16(1.00, 1.35)
Large1.00-1.00-1.00-
Hospital Location 5
Urban 1.00-1.00-1.00-
Large Rural City/Town0.95(0.77, 1.18)0.84(0.62, 1.13)1.03(0.76, 1.39)
Small Rural Town1.07(0.79, 1.45)1.06(0.69, 1.64)1.07(0.68, 1.68)
1 Includes depressive disorders, schizophrenia and other psychotic disorders, bipolar and related disorders, anxiety disorders, trauma and stressor related disorders, attention deficit hyperactivity disorders, disruptive, impulse control, and conduct disorders, personality disorders, and obsessive compulsive and related disorders. Excludes dementia since it is not a mental disorder. 2 Other Race includes Native American, Non-Hispanic multiracial and other races. 3 Other includes Workers’ Compensation, CHAMPUS, Title V, other government programs, and non-government programs. 4 Hospital Size Small: Rural 1–49 beds, Urban Non-Teaching 1–124 beds, Urban Teaching 1–249 beds; Medium: Rural 50–99 beds, Urban Non-Teaching 125–199 beds, Urban Teaching 250–424 beds; Large: Rural 100+ beds, Urban Non-Teaching 200+ beds, Urban Teaching 425+ beds. https://www.ncbi.nlm.nih.gov/books/NBK373736/table/sb205.t3/ (accessed on 14 December 2021). 5 Rural-urban commuting area (RUCA) codes were used to categorize rural and urban areas of NYS. 6 Risk Ratio (RR) of 1.00 indicates reference category.
Table 4. Multivariable model of risk factors for hospitalization for all behavioral health disorders, mental illness, and substance-use disorders with a concurrent cold-related illness (October–April), 2005–2019.
Table 4. Multivariable model of risk factors for hospitalization for all behavioral health disorders, mental illness, and substance-use disorders with a concurrent cold-related illness (October–April), 2005–2019.
CharacteristicsAll Behavioral Health DisordersMental Illness 1Substance-Use
Disorders
RR95% CIRR95% CIRR95% CI
Age Categories
5–17 0.76(0.62, 0.92)0.53(0.40, 0.72)2.34(1.79, 3.07)
18–391.006-1.00-1.00-
40–641.53(1.43, 1.64)1.48(1.34, 1.63)1.48(1.37, 1.60)
65–741.55(1.40, 1.71)1.24(1.08, 1.43)1.68(1.49, 1.90)
75+1.78(1.63, 1.95)1.34(1.17, 1.53)1.64(1.40, 1.92)
Sex
Male2.06(1.96, 2.17)1.89(1.76, 2.03)2.29(2.13, 2.47)
Female1.00-1.00-1.00-
Race
Non-Hispanic White1.00-1.00-1.00-
Hispanic/Latino0.80(0.74, 0.87)0.60(0.52, 0.70)0.74(0.67, 0.82)
Non-Hispanic Black1.23(1.15, 1.30)1.28(1.17, 1.41)0.91(0.84, 0.98)
Asian/Pacific Islander0.80(0.65, 0.98)0.72(0.51, 1.01)0.95(0.73, 1.23)
Other 21.07(0.99, 1.15)0.91(0.81, 1.03)1.03(0.94, 1.14)
Payer
Medicare1.51(1.38, 1.64)1.50(1.33, 1.70)1.43(1.27, 1.60)
Medicaid1.71(1.59, 1.85)1.57(1.40, 1.75)1.66(1.52, 1.82)
Private/HMO1.00-1.00-1.00-
Uninsured1.76(1.59, 1.94)1.45(1.24, 1.71)1.78(1.58, 2.00)
Other 30.84(0.63, 1.12)0.75(0.48, 1.16)0.78(0.54, 1.11)
Comorbidities
Diabetes0.89(0.84, 0.95)0.97(0.89, 1.06)0.91(0.83, 0.99)
Cardiovascular Disease0.68(0.64, 0.72)0.72(0.65, 0.79)0.60(0.54, 0.66)
Respiratory Disease1.08(1.03, 1.14)1.14(1.06, 1.22)1.18(1.11, 1.26)
Obesity 0.49(0.43, 0.54)0.50(0.43, 0.58)0.47(0.40, 0.55)
Hospital Size 4
Small0.92(0.86, 0.98)0.94(0.86, 1.03)0.78(0.72, 0.84)
Medium0.95(0.90, 1.01)0.93(0.85, 1.01)0.86(0.80, 0.92)
Large1.00-1.00-1.00-
Hospital Location 5
Urban 1.00-1.00-1.00-
Large Rural City/Town1.24(1.13, 1.36)1.25(1.10, 1.43)1.00(0.88, 1.14)
Small Rural Town1.54(1.35, 1.75)1.58(1.30, 1.91)1.05(0.86, 1.28)
1 Includes depressive disorders, schizophrenia and other psychotic disorders, bipolar and related disorders, anxiety disorders, trauma and stressor related disorders, attention deficit hyperactivity disorders, disruptive, impulse control, and conduct disorders, personality disorders, and obsessive compulsive and related disorders. Excludes dementia since it is not a mental disorder. 2 Other Race includes Native American, Non-Hispanic multiracial and other races. 3 Other includes Workers’ Compensation, CHAMPUS, Title V, other government programs, and non-government programs. 4 Hospital Size Small: Rural 1–49 beds, Urban Non-Teaching 1–124 beds, Urban Teaching 1–249 beds; Medium: Rural 50–99 beds, Urban Non-Teaching 125–199 beds, Urban Teaching 250–424 beds; Large: Rural 100+ beds, Urban Non-Teaching 200+ beds, Urban Teaching 425+ beds. https://www.ncbi.nlm.nih.gov/books/NBK373736/table/sb205.t3/ (accessed on 14 December 2021). 5 Rural-urban commuting area (RUCA) codes were used to categorize rural and urban areas of NYS. 6 Risk Ratio (RR) of 1.00 indicates reference category.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Aydin-Ghormoz, H.; Adeyeye, T.; Muscatiello, N.; Nayak, S.; Savadatti, S.; Insaf, T.Z. Identifying Risk Factors for Hospitalization with Behavioral Health Disorders and Concurrent Temperature-Related Illness in New York State. Int. J. Environ. Res. Public Health 2022, 19, 16411. https://doi.org/10.3390/ijerph192416411

AMA Style

Aydin-Ghormoz H, Adeyeye T, Muscatiello N, Nayak S, Savadatti S, Insaf TZ. Identifying Risk Factors for Hospitalization with Behavioral Health Disorders and Concurrent Temperature-Related Illness in New York State. International Journal of Environmental Research and Public Health. 2022; 19(24):16411. https://doi.org/10.3390/ijerph192416411

Chicago/Turabian Style

Aydin-Ghormoz, Heather, Temilayo Adeyeye, Neil Muscatiello, Seema Nayak, Sanghamitra Savadatti, and Tabassum Z. Insaf. 2022. "Identifying Risk Factors for Hospitalization with Behavioral Health Disorders and Concurrent Temperature-Related Illness in New York State" International Journal of Environmental Research and Public Health 19, no. 24: 16411. https://doi.org/10.3390/ijerph192416411

APA Style

Aydin-Ghormoz, H., Adeyeye, T., Muscatiello, N., Nayak, S., Savadatti, S., & Insaf, T. Z. (2022). Identifying Risk Factors for Hospitalization with Behavioral Health Disorders and Concurrent Temperature-Related Illness in New York State. International Journal of Environmental Research and Public Health, 19(24), 16411. https://doi.org/10.3390/ijerph192416411

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop