Introduction

Attention-deficit/hyperactivity disorder (ADHD), the most common childhood neurodevelopmental disorder, is a growing health crisis in the United States. ADHD is characterized by long-lasting negative consequences, including increased risks for substance abuse and addiction, anti-social and criminal behavior, lower academic attainment, and significantly greater personal and overall economic adversity; an overall estimated incremental median cost of $949.24 per year is associated with ADHD cases, with the average direct expenditure increasing in younger age groups up to preschool age1,2,3,4,5,6. Moreover, childhood ADHD prevalence has been increasing since the late 1990s and early 2000s7,8,9, resulting in an estimated 8.7% (5.3 million) of U.S. children aged 3–17 years diagnosed with ADHD in 2018/2019 [National Survey of Children’s Health (NSCH) 2018–19]. Symptoms in diagnosed adolescents and young adults have worsened during the COVID-19 pandemic, but it is not known if this trend will continue post-COVID-1910,11,12,13. ADHD outcome data should be promptly reported and analyzed to recognize any potential need for population-specific interventions. Notably, ADHD diagnosis is difficult, and the symptoms are nonspecific14,15,16. Therefore, while current trends indicate a growing number of ADHD cases during the pandemic10,11,12,13, possible false positive ADHD diagnosis due to increased severity of ADHD-associated symptoms stipulates effective diagnosis and treatment of ADHD.

Based on representative samples of U.S. physician office visits, rates of ADHD diagnosis and medication treatment have been significantly increasing17,18,19,20,21. The combined effect of several newly approved ADHD drugs, the 5th edition revision of the Diagnostic and Statistical Manual of Mental Disorders manual (DSM-5), and updated ADHD guidelines for children, teens, and adults may have contributed to the growing number of ADHD cases22,23,24. The DSM-5 and American Academy of Pediatrics (AAP) guidelines revised the ADHD diagnosis age criteria by raising the age of symptom onset from seven to twelve years old and outlining pediatric guidelines for children aged 4–18 years22,24. Furthermore, the updated DSM-5 guidelines require adolescents to display only five, instead of six, ADHD symptoms22.

National 2016 prevalence estimates indicate that 62.0% of children 2–17-year-old diagnosed with ADHD take medication, accounting for 5.1% of all U.S. children within that age range25. The AAP recommends preschool-age children (4–5-year-old) receive behavioral therapy administered by a parent or teacher as the first treatment option, with medication prescribed only if behavioral therapy does not reduce impairment, and that elementary-school-age children (6–11-year-old) and adolescents (12–18-year-old) receive a combination of medication and behavioral therapy24. Consistent with the AAP guidelines, the Healthy People 2020 developmental objectives include increasing the proportion of ADHD-diagnosed preschool-age children who receive behavioral therapy, and children and adolescents aged 6–17 years who receive medication, behavioral therapy, or a combination of the two26. To evaluate the success of AAP treatment guidelines and Healthy People 2020 goals and to inform effective responses, especially during the potential escalation of ADHD severity during the COVID-19 pandemic, trends related to the treatment of ADHD should be monitored10,11,12,13.

While regional patterns of ADHD diagnosis and medication prescription vary, the U.S. Southern region generally has higher rates of diagnosis and medication, while the Western region has lower rates of both9,27. National data could mask within-state variance relevant to interactions of multiple factors, including state demographics, socioeconomic status (SES) factors, and ADHD diagnosis and medication treatment rates. Targeted studies have identified SES patterns related to ADHD outcomes across counties and greater variation of ADHD outcomes have been reported among counties28,29,30. However, the U.S. lacks a comprehensive system to collect and report ADHD data at the county-level, and the existing data is subject to misrepresentation and discrepancies in sampling and data collection methodology29,30,31.

ADHD diagnosis among children has been correlated to family income, gender, race/ethnicity, primary household language, and insurance status32,33,34,35,36. Factors associated with a higher likelihood of taking ADHD medication included being non-Hispanic, living in homes with English as the primary language, living in the South, and having a currently co-occurring condition(s)25. These associations have not been reexamined collectively since trends were identified in the 2016 NSCH25. While previous studies identified variables associated with a higher likelihood of ADHD diagnosis and medication prescription, many important SES indicators from the NSCH were grouped and not individually assessed in relation to each other, potentially confounding critical trends and associations. Thus, when evaluating available data, such variation within broad SES categories should be considered and identified. This approach likely represents a more equitable survey-based epidemiological data analysis, minimizing generalization biases and potential misrepresentation of individuals’ SES designation.

In this study, we quantify more recent associations between ADHD diagnosis and a broad range of SES factors including income, insurance, area of residence, neighborhood characteristics, race/ethnicity, sex, family structure, household language, parent education, parental nativity, and school environment, utilizing data from the NSCH 2018–2019 dataset. The inclusion of available factors associated with mental health disorders, such as family structure, highlights important associations to childhood ADHD diagnosis in the U.S37,38,39. This comprehensive analysis of both national and representative state-level data differentiates national and state-specific associations.

Since the COVID-19 pandemic is projected to significantly exacerbate the effects of ADHD, specifically in children10,11, an analysis of the most recent pre-COVID-19 ADHD data and SES factors is necessary as a basis to better understand during and post COVID-19 ADHD trends10,40. Our national and state-specific hierarchical analysis, encompassing relevant SES factors, demonstrates the correlations between the various factors and ADHD outcomes. Since national data does not represent the demographic heterogeneity among states and local regions, we compared the data of two outlier states with opposing ADHD rates to national averages to better understand the effect of SES factors on ADHD diagnosis among children and adolescents in the U.S.

Results

Geographic distribution of ADHD diagnosis and treatment

Prevalence of parent-reported ADHD diagnosis and medication among noninstitutionalized children 3–17-years-old for 2018–2019 are presented by state and region in Fig. 1, with the District of Columbia treated as a state throughout this study. The nationwide prevalence of ADHD is 8.7%, and the percent of children taking medication among those diagnosed is 62.1%. States in the Southern region have the highest statewide average rates of diagnosis (10.97%) and medication (7.52%), while those in the Western region have the lowest statewide average rates of diagnosis (7.46%) and medication (4.14%) (Fig. 1a–c).

Figure 1
figure 1

Distribution of parent-reported ADHD prevalence and medication throughout the United States. United States map showing percentages by state of (a) children aged 3–17 years currently diagnosed with ADHD and (b) children with ADHD currently taking ADHD medication. Darker shades indicate states with relatively higher percentages, and lighter shades indicate states with relatively lower percentages. California (5.6%) has the lowest percentage, Nevada (5.9%) has the second-lowest percentage, and Louisiana (15.7%) has the highest percentage of children who currently have ADHD. Nevada (32.2%) has the lowest percentage, Nebraska (81.76%) has the highest percentage, and Louisiana has the third-highest percentage (76.28%) of children with ADHD currently taking ADHD medication. Percentages of diagnosis and medication were mapped onto the U.S. states using the map function of the graph builder in JMP 14.3.0. (c) Compared to the Midwest (n = 11,946) and Northeast (n = 9,102), regional state averages of children diagnosed with ADHD are highest in the South (n = 17,889) and lowest in the West (n = 12,979). Red indicates children diagnosed with ADHD who are not receiving medication. Blue indicates children diagnosed with ADHD who receive medication. Darker shades of both gradients indicate children in older age brackets. (d) Cell plot of states ordered lowest (top) to highest (bottom) by region as ratios of the percent of children aged 3–17 years diagnosed with ADHD taking medication to the percent of all children aged 3–17 years diagnosed with ADHD. Darker shades of green indicate higher ratios relative to other states. There are 12 states in the Midwest, 9 states in the Northeast, 17 states in the South, and 13 states in the West. Data were collected from the National Survey of Children’s Health (NSCH) 2018–19 combined dataset.

Within each region, the oldest age group (12–17-year-old) are the largest group of children diagnosed with ADHD and receiving medication (Midwest 2.81%, Northeast 3.04%, South 3.89%, West 2.32%) and the youngest age group (3–5-year-old) is the smallest in each region (Midwest 0.10%, Northeast 0.14%, South 0.20%, West 0.05%). This pattern is consistent among children diagnosed but not receiving medication, with 12–17-year-old being the largest group (Midwest 2.14%, Northeast 2.17%, South 1.74%, West 1.85%) and 3–5-year-old being the smallest group (Midwest 0.14%, Northeast 0.27%, South 0.30%, West 0.15%). In all regions, more children diagnosed with ADHD are receiving medication than not within the age groups of 6–11 and 12–17. The opposite trend is evident within the youngest age group of 3–5-year-old (Fig. 1c).

When statewide ratios of medication prevalence among children diagnosed with ADHD (Fig. 1b) to overall child ADHD diagnosis prevalence (Fig. 1a) is organized by region; the Midwest ranges from 5.5 in Missouri (MO) to 11.4 in Nebraska (NE), and the Northeast ranges from 5.0 in Massachusetts (MA) to 7.0 in Pennsylvania (PA). In the region with the highest average prevalence of ADHD, the South, the ratios range from 4.4 in the District of Columbia (DC) to 9.1 in Maryland (MD). In the region with the lowest average prevalence, the West, the ratios range from 5.1 in Oregon (OR) to 8.3 in Hawaii (HI; Fig. 1d).

ADHD diagnosis and treatment variation among representative states

Of all states, Louisiana (LA; 15.7%) has the highest percentage of children diagnosed with ADHD and California (CA; 5.6%) has the lowest, followed by Nevada (NV; 5.9%; Figs. 1a, 2). The percentage of children taking ADHD medication among those diagnosed (Fig. 1b) is highest in NE (80.0%), third highest in LA (76.3%), and lowest in NV (32.2%). Kansas (KS) diagnosis (9.2%) and medication (62.4%) rates are closest to the median rates of all states. NE and South Dakota (SD) have low diagnosis rates (7.0% for both), but high medication rates (80.0% and 74.3%, respectively) relative to the median rates. Massachusetts (MA) and West Virginia (WV) have high diagnosis rates (10.5% and 14.5%, respectively), but MA has a low medication rate (52.9%), while WV has a near-median medication rate (64.6%; Fig. 2).

Figure 2
figure 2

Compared to estimates of (a) the national population (n = 51,916), (bh) the represented states vary from 5.9 to 15.7% of children age 3–17 years diagnosed with ADHD and 32.2% to 76.3% of children taking ADHD medication among those diagnosed. From lowest to highest rate of diagnosis, the states represented are Nevada (n = 994), Nebraska (n = 953), South Dakota (n = 960), Kansas (n = 1,045), Massachusetts (n = 1,029), West Virginia (n = 1,082), and Louisiana (n = 1,080). Gray indicates children that are not currently diagnosed with ADHD and red indicates children that are currently diagnosed with ADHD. Blue indicates children that are currently diagnosed with ADHD and are taking ADHD medication, and orange indicates children that are currently diagnosed with ADHD but are not taking ADHD medication. (b) Nevada represents particularly low rates of diagnosis and medication, while (h) Louisiana represents particularly high rates of diagnosis and medication. Data were collected from the NSCH 2018–2019 combined dataset.

Multiple correspondence analysis of ADHD diagnosis and SES factors

Individual NSCH responses were analyzed nationally (n = 59,445), within LA (n = 1255), and within NV (n = 1154). Of all analyzed individual response data, 6.7% were missing and imputed for Multiple Correspondence Analysis (MCA; Fig. 3a–c). In MCA, inter-variable proximity indicates the strength of possible associations. The national dataset analysis resulted in 43 dimensions, with inertias ranging from 0.2% to 6.8% and a total inertia of 13.2% for the first two dimensions. LA analysis resulted in 42 dimensions, with inertias ranging from 0.2% to 7.4% and a total inertia of 13.6% for the first two dimensions. NV analysis resulted in 41 dimensions, with inertias ranging from 0.2% to 8.1%, and a total inertia of 13.9% for the first two dimensions.

Figure 3
figure 3

MCAs of ADHD diagnosis, SES factors, and demographics for samples (a) nationwide (n = 59,445), (b) from Louisiana (n = 1,255), and (c) from Nevada (n = 1,154). Data were collected from the NSCH 2018 and 2019 individual responses, and missing responses were imputed using the regularized iterative MCA algorithm of the missMDA R package. (ac) MCA plots display associations between variables, with more closely associated variables displayed closer together. For instance, in panel A, ADHD_N and CBSA_Y are closely associated, but CBSA_N and Race_P are not. Positive diagnosis (ADHD_Y) and lack of diagnosis (ADHD_N) are shown in red. The percent variance retained is indicated for each dimension. (d) Table of all variables analyzed, with their possible responses and abbreviations.

None of the national, LA, or NV MCAs indicate a relationship with ADHD diagnosis for parental nativity or Hispanic ethnicity (Fig. 4a–c). In all three analyses, unsafe school, unsafe neighborhood, and economic hardship are associated with ADHD diagnosis, while safe school and neighborhood, no economic hardship, two married parents, a household language besides English, no guardian’s death, and no parental divorce or separation are associated with a lack of ADHD diagnosis. All three MCAs display no guardian’s death and no parental divorce or separation in close proximity to a lack of ADHD diagnosis; NV data distribution does not display these adverse experiences near to the positive ADHD diagnosis, although both national and LA MCAs do.

Figure 4
figure 4figure 4

Bar plots representing MCAs of ADHD diagnosis, SES factors, and demographics for samples (a) nationwide (n = 59,445), (b) from Louisiana (n = 1,255), and (c) from Nevada (n = 1,154). Data were collected from the NSCH 2018 and 2019 individual responses, and missing responses were imputed using the regularized iterative MCA algorithm of the missMDA R package. Bars show each variable’s proximity to either diagnosis variable, calculated as (distance to ADHD_N−distance to ADHD_Y)/(distance between ADHD_N and ADHD_Y). A value closer to 1.0 indicates a variable closer to the positive diagnosis and further from the lack of diagnosis, while a value closer to − 1.0 indicates a variable closer to the lack of diagnosis and further from the positive diagnosis.

MCA reveals the lack of ADHD diagnosis associated with Asian and multi-racial children in the national dataset, Pacific Islander, White, multi-racial, and other children in LA, and White and Asian children in NV. In LA, a positive ADHD diagnosis is associated with Black children. The national MCA reveals associations between lack of ADHD diagnosis and higher urbanicity [living in a core-based statistical area (CBSA), metropolitan statistical area, or metropolitan principal city], which are not evident for LA or NV (Fig. 3a,d). Neighborhood amenities, namely parks, libraries, recreation centers, and sidewalks, are associated with a lack of ADHD diagnosis nationally, but not in LA or NV. Neighborhood maintenance, namely a lack of vandalism, rundown housing, and litter, show a clear relationship with ADHD diagnosis for LA, but a weaker relationship for NV and almost no relationship nationally. Neither the national nor NV MCAs display proximity between the household adults’ highest education and ADHD diagnosis, but the LA MCA suggests an association between a greater than high school educational attainment and a lack of ADHD diagnosis (Fig. 3a,c). National and LA MCAs display insurance coverage associated with a lack of ADHD diagnosis, while a lack of or discontinuous insurance coverage is moderately associated with ADHD diagnosis. Proximity to ADHD diagnosis is observed for grandparent households in the national dataset and for single-mother households in LA. Only the LA MCA reveals a relationship between sex and ADHD diagnosis, with a lack of diagnosis nearer to female than a positive ADHD diagnosis.

Correlation of ADHD diagnosis and SES factors

Chi-square and Cramér’s V analyses reveal levels of significance and effect sizes for each response dataset (Figs. 5, 6; Table 1). Nationally, all factors tested are significant against ADHD diagnosis after Bonferroni p-value adjustment except Women, Infants and Children (WIC) program participation, all neighborhood maintenance factors (vandalism, rundown housing, and litter), the presence of a neighborhood library, all urbanicity factors (metropolitan principal city, metropolitan statistical area, and CBSA), insurance status, and insurance gap (Fig. 5a). The largest effect sizes, namely the factors most strongly related to ADHD diagnosis for the national dataset, are sex and having experienced divorce or separation of parents. The smallest effect sizes are the presence of neighborhood sidewalks, the presence of neighborhood recreation centers, Hispanic ethnicity, and the highest educational attainment of household adults (Fig. 6a). Hierarchical clustering of variables based on national dataset p-values (Fig. 5a-II) clusters ADHD diagnosis with free school lunch, and race at the lowest level, which are all then clustered with the presence of a neighborhood library. At the highest level, sex does not cluster with ADHD and the rest of the variables.

Figure 5
figure 5figure 5

Correlation matrices and hierarchical clustering dendrograms from chi-square test adjusted p-values between ADHD diagnosis, SES factors, and demographics (a) nationwide (n = 52,065), (b) in Louisiana only (n = 1,084), and (c) in Nevada (n = 999) only. In the (I) matrices, black indicates statistical nonsignificance (p > 0.05), darker colors indicate lesser significance, and yellow indicates the greatest statistical significance. “ + ” and “−” indicate positive and negative correlations, respectively. (II) Dendrograms reflect hierarchical clustering of examined variables according to p-values from pairwise inter-variable chi-square analyses. Variables are grouped according to similarities in p-values, and greater branch length along the x-axis indicates greater cluster dissimilarity.

Figure 6
figure 6

Matrices of Cramér’s V values between ADHD diagnosis, SES factors, and demographics (a) nationwide (n = 52,065), (b) in Louisiana only (n = 1,084), and (c) in Nevada (n = 999) only. Darker colors indicate lower Cramér’s V effect sizes and yellow indicates the highest Cramér’s V effect sizes.

Table 1 Chi-square and Cramér’s V analyses values of the socioeconomic factors discussed, as compared to ADHD diagnosis.

There are no statistically significant correlations to ADHD diagnosis in the LA dataset (Fig. 5b-I). Besides sex, the largest effect sizes on ADHD diagnosis are neighborhood safety and having experienced the death of a guardian. The smallest effect sizes are WIC program participation, neighborhood vandalism, the presence of neighborhood recreation centers, the presence of neighborhood library, metropolitan statistical area, and family government cash assistance (Fig. 6b). In the LA hierarchical clustering, ADHD is grouped with guardian death at the lowest level, which both share a cluster with U.S. nativity of caregivers, household language, and metro area status.

For the NV dataset, only economic hardship is significant relative to ADHD diagnosis (Fig. 5c-I). Besides economic hardship, the largest effect sizes on ADHD diagnosis are sex and the ability to afford food (Fig. 6c). The smallest effect sizes are WIC program participation, neighborhood litter, the presence of neighborhood library, family government cash assistance, and caregiver employment status. For NV, ADHD is hierarchically clustered with safe neighborhood, food stamps or Supplemental Nutrition Assistance Program (SNAP), free school lunch, and economic hardship.

Discussion

ADHD is the most diagnosed emerging neurodevelopmental disorder in children and likely an escalating health care crisis during and post-COVID-19 U.S. pandemic. We utilized the most recent publicly available dataset (NSCH 2018–2019) to analyze a broad range of socioeconomic status (SES) factors relative to childhood ADHD and to each other. An analysis of SES factor interactions can influence future data collection strategies and help develop new hypotheses. We selected two opposing outlier states, Louisiana (LA), and Nevada (NV), based on high and low ADHD diagnosis and medication rates, respectively, to compare national to state-specific associations, since the state-level variation of ADHD diagnostic prevalence has been established27. While a more comprehensive approach is to analyze the greater variation between and within counties, the data collection systems in the U.S. that report county-specific childhood ADHD diagnosis are limited29,30. This analysis is not focused on sex-specific correlations since our chi-square results (Fig. 5a-I) are in agreement with current literature indicating a higher prevalence of ADHD among males41,42,43,44. While common trends such as higher rates of ADHD diagnosis among adolescents than children are expected, certain weak- to no-associations are surprising: non-significant correlation of recent insurance coverage status (Fig. 5) and limited direct state-specific correlations to ADHD diagnosis25.

Higher ADHD diagnosis and medication treatment rates in the Eastern U.S.

The national prevalence of ADHD in children and adolescents follows previously established regional distribution, with higher diagnosis rates in the Eastern U.S. compared to the West (Fig. 1a)9. ADHD medication treatment rates follow a similar pattern like ADHD diagnosis, but are less concentrated regionally, with evident state-specific variations (Fig. 1b): Texas (TX) and Nebraska (NE) have relatively low diagnosis rates, but higher rates of ADHD-diagnosed medicated children, while West Virginia (WV) has a relatively high diagnosis rate but low medication treatment rate. The Western U.S. has the least and Northeastern U.S. has the most variable diagnosis rates (1a); relative to ADHD medication treatment, the Midwestern U.S. is the least variable, while Northeastern U.S. is the most variable (1b). The pattern of ADHD diagnosis by age group is relatively consistent between regions (1c). Considering the ratios of the medication to diagnosis rates, which is arguably the most relevant comparison when assessing the regional distributions, the Midwestern U.S. is the most, while Northeastern U.S. is the least variable (Fig. 1d).

Regional ADHD medication treatment rates are more variable than diagnosis rates

To better understand ADHD diagnosis associations with SES factors, within-state variation is arguably more informative than either the national or regional distribution, especially since the variation among counties in ADHD outcomes is associated with the socioeconomic composition of the population27,30. LA and NV are discussed in detail because each is a national outlier on the opposite ends of the spectrum in their respective regions, marked by the highest and lowest combined average diagnosis and medications rates. Nebraska (NE), South Dakota (SD), Kansas (KS), Massachusetts (MA), and West Virginia (WV) represent in-between state-specific prevalence (Fig. 2). The order of states from NV to LA represents an increasing ADHD diagnosis rate, but the corresponding medication treatment trends do not show a clear pattern related to diagnosis rates. Moreover, while the rates of ADHD diagnosis across states in the Midwestern region (NE, SD, and KS) are consistent (90.8–93.0%), there is considerable variation in medication treatment rates (NE–80%, SD–74.3%, and KS–62.4% respectively; Fig. 2). Notably, although KS is the median state most closely reflecting the national prevalence, the state-specific ADHD data distribution varies and does not necessarily fit national patterns. Within-state variation is consistent with previously identified regional and state pattern discrepancies. Although the state variation determinants are not known, at least two factors, the characteristics of the statewide physician population and laws concerning a schools’ function in ADHD treatment recommendations, are associated with differences in ADHD diagnosis and medication treatment rates27. While local physician population and policies regarding schools’ roles in ADHD treatment have potential confounding effects on the variability of ADHD diagnosis between and within states, these data were not included in the NSCH and could not be quantified.

Older children and adolescents have higher ADHD diagnosis and medication treatment rates

As the cumulative effect of ADHD symptoms results in older children (12–17-year-old) having higher diagnosis rates, regional trends indicate that ADHD-diagnosed adolescents more often take medication than not (Fig. 1c)25. Such a tendency is anticipated because medication is more often recommended to patients in this age group as per American Academy of Pediatrics (AAP) guidelines22,24. The youngest age group (3–5-year-old) has the lowest diagnosis and medication treatment rates across all regions (Fig. 1c). This is expected, since young children are difficult to diagnose due to a lack of well-defined symptoms in social settings and diagnosis delays, particularly since there is no objective test to diagnose ADHD27,45,46,47,48. Once diagnosed, the preferred initial treatment for the younger age groups is behavioral therapy22,24. Therefore, these results are predictable and indicate a need to investigate SES factors that influence variation in ADHD diagnosis among all age groups.

Correlation of ADHD diagnosis and race

MCA analysis considers non-linear associations among variables and reveals patterning in complex data sets with categorical dependent variables, while matrices provide statistical correlations between ADHD diagnosis and 30 SES factors nationally and within two outlier states (Figs. 36)49. The relationships between ADHD and SES factors reveal four overarching categories on the national and state levels: race/ethnicity, financial status, family structure, and neighborhood characteristics. Nationally, ADHD diagnosis is significantly positively correlated with White children and negatively correlated with Asian children (Fig. 5a). Associations of race/ethnicity and SES factors to ADHD diagnosis have generated inconsistent results in the past, which may be due to confounding variables32,38. Our analysis supports the association of higher ADHD diagnosis rates among White children and highlights the direct correlation between ADHD and race, accounting for various SES factors33,34. It has been documented that race/ethnicity is associated with other SES factors, but the variation of the confounding effects of individual interactions is not as established50. The effect size of race to other SES factors indicates the strength of confounding relationships, notably stronger correlations to neighborhood sidewalks, family structure, adults’ highest education, and financial factors including food stamps or SNAP and free school lunch (Fig. 6a).

Race has a higher effect size than other factors nationally and within LA and NV for all variables, except sex and safe schools nationally (Fig. 6). The number of significant correlations and high effect size indicate any relationships between ADHD and race could be due to systemic inequities, especially considering the large racial, ethnic, and cultural diversity which interact with other SES factors51. While ADHD diagnosis and race are not significantly correlated within LA and NV (Fig. 4b,c), the implications of race and ethnicity, such as speaking a household language other than English or Spanish, are strongly associated with not being diagnosed with ADHD (Fig. 5a–c); which is expected as English-speaking households have higher rates of ADHD diagnosis9,34. Historically, children with special health care needs from non-English primary-language homes lack access to medical homes, a common source of family-oriented primary care52,53. Our data shows a significant negative correlation of having current health insurance to a household language other than English or Spanish, which is consistent with previous findings of insurance gaps in non-English primary-language homes52,54. National and state-specific overlaps of race and ADHD diagnosis associations highlight the need to better understand both direct and indirect effects of race- and ethnicity-related indicators on ADHD diagnosis and eventual medication trends among children.

Correlation of ADHD diagnosis and family financial status

We also investigated the relationship between ADHD and indicators of family financial status, including government cash assistance, free school lunches, food stamps/SNAP program, economic hardship, both caregivers being employed, and food affordability. Nationally, all indicators of financial status, except having a gap in health insurance in the last twelve months, are significantly correlated to ADHD diagnosis. This is unexpected because lack of health insurance is an identified factor related to a decreased risk of ADHD diagnosis34. Relative to other SES factors, insurance status did not have a direct association with ADHD diagnosis but is significantly correlated to other financial status factors, suggesting an indirect effect. Financial status is significantly correlated to other SES factors associated with ADHD diagnosis. Like race, financial status appears to have a systemic association with ADHD diagnosis, either directly, or indirectly via correlations to parents being divorced or separated, living in a safe neighborhood, and children attending a safe school (Fig. 5a).

Expectedly, financial status indicators are correlated to adults’ highest education and neighborhood amenities, since families with a higher income can afford to live in communities with more resources resulting in better overall health outcomes55. Nationally, adults’ highest education level is significantly correlated to every financial status indicator, following the pattern of education higher than high school resulting in better financial outcomes, which is consistent with findings of overall higher education being correlated to lower income inequality56. Although not statistically significant for every correlation, that pattern is consistent across LA and NV, as within-state associations relevant to ADHD diagnosis are more variable. Indirect association of ADHD diagnosis could be argued as four indicators of financial hardship (caregivers being unemployed, being on food stamps or SNAP, qualifying for free school lunch, and government cash assistance) are correlated to divorce or separation of parents, which is strongly associated to ADHD diagnosis39 (Figs. 4, 5). In NV, ADHD diagnosis is only significantly correlated to economic hardship, although economic hardship correlates to divorce or separation of the parents and family structure, suggesting an indirect correlation57. Family structure, specifically divorce/separation and being raised by a single parent, has been reported to increase the likelihood of childhood ADHD diagnosis, reinforcing the existence of correlations between financial status and family structure38,39. These findings support the importance of investigating the various interactions between SES factors and ADHD outcomes.

Correlation of ADHD diagnosis and family structure

Family structure factors, associated with ADHD diagnosis both nationally and within LA and NV, highlight the relationship between the household environment and ADHD diagnosis in children. Family structure, specifically living in a household without two married parents, increases the likelihood of children experiencing economic hardship and being exposed to violence and disruptive households58. Adverse childhood experiences (ACEs), including economic adversity and household violence, are reportedly significantly associated with childhood ADHD diagnosis, even when controlling for known parental and child-risk factors of ADHD59,60. Hence, family structure is important to investigate in relation to other SES factors and ADHD outcomes. Nationally, grandparent households, the death of a guardian, and divorce or separation are associated with parent-reported ADHD diagnosis (Fig. 4a). Divorce or separation, with the highest effect size (Fig. 6a), is statistically correlated to ADHD and every factor except sex and metro-principal status (Fig. 5a). This reinforces the strong association of living in a household with divorced or separated parents and a higher likelihood of ADHD diagnosis, whether through direct correlation or indirect effects on other SES factors, as well as a lower likelihood of childhood ADHD diagnosis for those living in a two-parent household (Fig. 5a). Therefore, family structure, combined with other relevant SES factors, is a strong marker of ADHD diagnosis nationally and is evident in outlier states: for example, in LA, the death of a guardian is strongly associated with parent-reported ADHD diagnosis (Fig. 4b). Living in a stable home with two married parents results in children with higher standards of living and experiencing fewer stressful events, which is evident in lower rates of ADHD diagnosis among kids58,61. The relationship of ACEs and unstable home environments are related to other negative outcomes in young adults, including poor physical and mental health, chronic disease, anxiety, and risk behaviors such as substance abuse which can create a positive feedback loop with symptoms of ADHD58,62,63,64. These associations highlight the importance to better understand the effects of ACEs and family structure on ADHD outcomes.

Correlation of ADHD diagnosis and neighborhood characteristics

Neighborhood characteristics are strongly associated with ADHD diagnosis in children65,66. Nationally, indicators of urbanicity, including living in a core-based statistical area (CBSA), metropolitan statistical area, or metropolitan principal city, are positively correlated with a lack of ADHD diagnosis, which is expected since urban areas have greater access to resources57. Neighborhood amenities, including the presence of parks, sidewalks, and recreation centers, promote physical activity which is associated with improved cognition and behavior and could benefit childhood ADHD outcomes65,67,68,69,70. Lack of neighborhood amenities is strongly associated with ADHD diagnosis (Fig. 4a) and vice versa, while both LA and NV show a distinct pattern of neighborhood maintenance, including rundown housing, litter, and vandalism, more strongly associated with ADHD diagnosis (Fig. 4b,c). These associations are evident in LA, with metropolitan principal city status significantly correlated to neighborhood library, park, recreation centers, and sidewalks access (Fig. 5b). Notably, urbanicity and neighborhood characteristic factors have low effect sizes on ADHD diagnosis highlighting a need for further investigation of their association with ADHD diagnosis (Fig. 6). However, CBSA status in LA and urbanicity factors (CBSA, metropolitan statistical area, and metropolitan principal city status) in NV were suppressed for respondent confidentiality and could not be analyzed.

National and statewide correlations should only be compared with careful consideration of differing sample sizes. Compared to the national analysis, both statewide analyses use a small sample size, resulting in lesser statistical power. Therefore, the overall greater p-values associated with the statewide analyses are expected, given that chi-square tests of smaller sample sizes can yield false negatives, and larger sample sizes can yield false positives. Small sample sizes may also result in less accuracy overall, indicating further needs for proper and comprehensive data collection.

Conclusions, limitations, and recommendations

Our study highlights important direct and indirect associations between ADHD diagnosis in children and critical SES factors, namely race, financial status, family structure, and neighborhood characteristics. While certain SES factors appear to have strong associations, this analysis is limited by the type of data sourced from the NSCH, representing a cross-section of parent-reported responses57. Consequently, we were unable to determine temporal order or causality between SES factors and ADHD diagnosis. While our analysis was limited to ADHD diagnosis and did not account for ADHD treatment differences or diagnosis severity, our focus was on evaluating associations of critical SES factors with ADHD diagnosis to establish a reference for the possible post-COVID-19 increase in ADHD severity among children and adolescents in the U.S.10,11. We also limited state-specific analyses to two opposing outlier states and did not account for county-level variation which has been shown to be more variable than among state variation in ADHD outcomes29,30. Future studies assessing relationships between SES factors and childhood ADHD diagnosis, severity, and treatment should rely on longitudinal data to establish causality, which is a critical metric for predicting trends and devising effective health care policies.

This report utilizes the most recent publicly accessible NSCH data from a nationally representative sample to highlight associations of 30 SES factors and parent-reported childhood ADHD diagnosis. Although focused on ADHD diagnosis, this report includes important interactions between SES factors, suggesting indirect correlations and confounding factors. This is timely because the COVID-19 pandemic is expected to exacerbate the severity of ADHD outcomes in the U.S.10,11. Kazda et al. reports pre-pandemic trends in ADHD overdiagnosis in children and adolescents, highlighting the potential for counterproductive, if not harmful, effects on misdiagnosed or mild-symptom-diagnosed youths16. Our results may serve as a reference for the post-COVID-19 data, help devise new hypotheses, and highlight gaps within the national data collection system. Data collection and reporting systems for county-level ADHD data are limited in the U.S.30 While there are studies addressing the variation between ADHD medication treatment among counties with identified profiles, the medication treatment prevalence data is subject to misrepresentation29,30,31. Moreover, prevalence of ADHD diagnosis is not systematically reported at the county or community level, and studies based on county variation reported discrepancies in sampling and data collection methodology28,30. To effectively model and mitigate the ADHD epidemic and similar national health crises, the U.S. should rely on comprehensive, county-specific, near real-time survey and epidemiological data reporting. Such data streams paired with longitudinal data normalized to census numbers are critical components of a primary prevention program development and implementation.

Methods

All data used represent portions of aggregate parent-reported ADHD diagnosis and treatment data for non-institutionalized children 3–17 years of age, disaggregated into three categorizations: 3–5-year-old, 6–11-year-old, and 12–17-year-old. Data are collected from the annual National Survey of Children’s Health (NSCH) 2018–2019 combined dataset, nationally and by state (https://www.nschdata.org/browse/survey/). NSCH is a United States population-based, cross-sectional survey sponsored by the Health Resources and Services Administration’s Maternal and Child Health Bureau. The NSCH is designed to collect data on health markers of noninstitutionalized children aged 0–17, which has been distributed via an online/mail-based format annually since 2016, previously conducted as a telephone survey every four years. Households that are more likely to contain children are oversampled but are otherwise randomly selected for an invitation household screen. Households with children, as indicated by the screen, receive a questionnaire regarding the health of one randomly selected child. A full explanation of the NSCH methods can be found at https://mchb.hrsa.gov/data/national-surveys/data-user.

Diagnosis and treatment rates were established using either weighted population estimates or survey sample counts, which were both collected directly from the NSCH datasets. The weights for population estimates account for nonresponse, selection probability, household size, household poverty threshold, educational attainment of the respondent, number of children with special healthcare needs in the household, and demographic characteristics of the child. Weighted population estimates and confidence intervals that were not directly provided by the NSCH were calculated from individual response datasets using the provided weights and the xlogit method from the survey package in R 3.6.371,72 following the survey design outlined by the NSCH source and accuracy statement (Table 2).

Table 2 Demographics by ADHD status, as counts of the NSCH sample and as estimated percentages of the United States 3–17-year-old population, from 2018–2019.

To determine regional ADHD diagnosis and medication rates by age group, the fifty U.S. states and the District of Columbia were grouped into four regions (Midwest, Northeast, South, and West)27. Weighted population estimates of children with ADHD by age group and medication status were calculated into a percentage of the total children aged 3–17 years per state. To examine states in the context of each region, ratios for each state were calculated as the percent of children aged 3–17 years diagnosed with ADHD that are taking medication divided by the percent of children aged 3–17 years that are diagnosed with ADHD. Weighted estimates of ADHD prevalence and medication use among children diagnosed with ADHD nationally and within individual states were used to further examine a variety of states.

Individual household survey responses from the NSCH 2018 and 2019 Topical Data and Input were combined and used to analyze national and statewide associations in ADHD diagnosis and SES factors. Each respondent was weighted to account for the survey year in addition to the same factors used to weight population estimates. The NSCH 2018 Topical SAS Codebook was used to code the questionnaire questions and parent-reported responses (https://www.childhealthdata.org/learn-about-the-nsch/nsch-codebooks). Approximately 71,000 screens were completed for the 2018 NSCH, with 30,530 of the 38,140 qualified households completing the child-specific questionnaire. Approximately 68,500 screens were completed for the 2019 NSCH, with 29,433 of the 36,196 qualified households completing the child-specific questionnaire (https://www.childhealthdata.org/learn-about-the-nsch/methods). Unrelated survey questions were removed before analysis. For multiple correspondence analysis (MCA), missing responses were imputed using the regularized iterative MCA algorithm of the missMDA R package73. Missing responses were excluded for chi-square analyses and associated Cramér’s V calculations.

Data analysis

The data analyzed in this study are from households that completed the child-specific questionnaire and recorded a valid response to the question “has a doctor or other health care provider EVER told you that this child has Attention Deficit Disorder or Attention Deficit/Hyperactivity Disorder, that is, ADD or ADHD?” (n = 9,445). If the parent or guardian responded “yes”, they were directed to answer the following questions:

  • Does this child CURRENTLY have the condition?

    • If yes, is it Mild, Moderate, or Severe?

  • Is this child CURRENTLY taking medication for ADD or ADHD?

  • At any time DURING THE PAST 12 MONTHS, did this child receive behavioral treatment for ADD or ADHD, such as training or an intervention that you or this child received to help with his or her behavior?

ADHD outcomes and SES factors were analyzed at the national, regional, and state levels. Local physician population and policies regarding schools’ roles in ADHD treatment were not included in the NSCH and could not be quantified. Percentages of children with a reported ADHD diagnosis currently taking medication in each state were calculated from weighted population estimates of parent-reported diagnosis and medication use. Percentages of diagnosis and medication were mapped onto the U.S. states using the map function of the graph builder in JMP 14.3.0.

To determine ADHD prevalence within each age group by region, weighted population estimates of each category by the state were used to calculate a percentage out of the total estimated population and then normalized to the total weighted percentage of children with ADHD in that state. Regional percentages are the averages of that category across all states in that region, normalized to the average total weighted percentage of children with ADHD in that region. States selected for further analysis reflect a combination of relatively low or high ADHD diagnosis and medication rates. The two states represent either end of both spectrums, and the states in between represent high diagnosis and low medication rates, low diagnosis and high medication rates, or median medication and diagnosis rates. These data were visualized using the bar graph (Fig. 1c) and cell plot functions (Fig. 1d) of JMP 14.3.0.

To analyze associations among ADHD diagnosis and SES factors, related items from individual household survey responses were analyzed via multiple correspondence analysis (MCA) using the FactoMineR package74, then plotted using the Factoextra package in R75. Associated bar plots were created using the ggplot2 package, with proximity values calculated as (distance to ADHD_N–distance to ADHD_Y) / (distance between ADHD_N and ADHD_Y). ADHD diagnosis was compared to qualitative survey items related to the family’s SES via Pearson’s chi-square analyses with Bonferroni correction and an adjusted p-value cut-off of p < 0.05 to quantify significance. To quantify the strength of association, Cramér’s V values were calculated for each chi-square value. For items with four possible ordinal responses (e.g. Definitely agree, Somewhat agree, Somewhat disagree, or Definitely disagree), to increase the counts in the contingency tables, the two most positive and two most negative responses were combined to yield two categories of responses (e.g. Definitely/somewhat agree or Definitely/somewhat disagree). To examine similarities in significance between variables, hierarchical clustering based on p-values using the complete linkage agglomerative method was performed in R 3.6.372. Associated matrices were created using the ggplot2 package, software version 3.3.076.

No protocol approval was needed for this study because it uses only publicly available, de-identified data. The author’s institutional review board does not consider the secondary analysis of publicly available data as research on human subjects. All methods were performed in accordance with relevant guidelines and regulations.