10124165432320Clin Toxicol (Phila)Clin Toxicol (Phila)Clinical toxicology (Philadelphia, Pa.)1556-36501556-951926609895510711310.3109/15563650.2015.1113543NIHMS820053ArticleRegional Variations in Pediatric Medication Exposure: Spatial Analysis of Poison Center Utilization in Western PennsylvaniaNguyenMargaret B.MD1PizonAnthony F.MD2BranasCharles C.PhD3FabioAnthonyPhD4Division of Emergency Medicine, Department of Pediatrics, Rady Children’s Hospital, University of California San Diego Division of Toxicology, Department of Emergency Medicine, University of Pittsburgh School of Medicine Department of Biostatistics and Epidemiology, School of Medicine, University of Pennsylvania, Philadelphia, PA Department of Epidemiology, Epidemiology Data Center, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, PA Corresponding Author: Margaret Nguyen, Rady Children’s Hospital San Diego, 3020 Children’s Way San Diego, CA 92123. mbnguyen1@rchsd.org. Phone: (858) 966-8036.2310201626112015201601120175414752Context

Medication drug exposures among young children continue to rise despite current poison prevention efforts. These exposures result in increased healthcare utilization and medical costs. New tactics are needed to reduce injuries related to pediatric drug exposures.

Objective

We aimed to identify cluster patterns in 1) calls for pediatric medication drug exposures and 2) a subset of calls that resulted in medical evaluation referrals. We identified and evaluated population characteristics associated with cluster patterns.

Methods

We analyzed 26,685 pharmaceutical drug exposures involving children less than 5 years of age based on calls reported to the Pittsburgh Poison Center from January 1, 2006 to December 31, 2010. We performed spatial statistics to assess for clustering. We used logistic regression was used to estimate population characteristics associated with clustering.

Results

Spatial analysis identified 22 exposure clusters and 5 referral clusters. 65% of 89 ZIP codes in the clusters of drug exposure with HCF referral were not identified in the exposure clusters. ZIP codes in the HCF referral clusters were characterized as rural, impoverished, and with high rates of unemployment and school dropouts.

Discussion

Our principal findings demonstrate pediatric drug exposures do exist in discrete geographic clusters and with distinct socioeconomic characteristics.

Conclusion

This study offers a starting point for subsequent investigations into the geographic and social context of pediatric medication drug exposures. This is an important step in revising pediatric poison prevention strategies.

spatial analysisGISinjury preventionsocioeconomic
INTRODUCTION

Over the last century, American childhood deaths associated with accidental poisoning had been significantly reduced. Recent reports, however, indicate that not only are pediatric drug exposures on the rise, but the associated clinical outcomes are worsening in severity. These medication exposures ultimately result in increased healthcare utilization. 15 Poisoning is one of the top reasons for Emergency Department (ED) visits in the 1 to 5 year age group in the United States(US).6 Due to their exploratory nature, children in this age group are particularly vulnerable to unintentional self-exposures.7,8 New prevention strategies are required to reduce drug exposures and to decrease the associated healthcare burden.9,10 Identification of medication poisoning locations and understanding the communities at risk for injury is essential to revising prevention strategies. Geographic information systems (GIS) is a computerized platform designed to collect, manage, and analyze spatial and non-spatial data. 11.Public health campaigns have utilized GIS and spatial scan statistics to assess regional variation in substance abuse prevalence, alcohol related injury, as well as to plan and assess the impact of preventive interventions in childhood lead poisoning. 1214

A critical gap exists in the literature regarding geographic trends and population characteristics of unintentional, drug exposures in children. In this study, we use poison center data, United States (US) Census information and spatial statistics to assess regional patterns of drug exposure in children < 5 years. The purpose of this study was to assess regional pharmaceutical poisoning rates, identify locations and populations at risk for injury, and determine geographic variations in the poisoning injury resulting in health services utilization. We hypothesized that drug exposures exist in clusters and are not randomly distributed over space. We performed a spatial analysis on all calls to the poison center regarding pediatric drug exposure and on the subset of calls that resulted in referrals for medical evaluation. We then compared population characteristics of ZIP codes with and without cluster patterns.

METHODSData Collection

The Pittsburgh Poison Center is the referral center for western Pennsylvania and covers 51 counties and a population of 5.9 million people. Data from phone calls are recorded in real time with a standardized electronic format and can serve as a surveillance system for pediatric drug exposure. We queried original data from the Pittsburgh Poison Center and searched for records meeting the following criteria: exposures in children < 5 years, call-types originating from non-healthcare facilities, and events occurring between January 1, 2006 and December 31st, 2010. Analysis of these records was limited to children with an unintentional exposure to a pharmaceutical agent. Unintentional exposures that resulted from a therapeutic error (i.e. caregivers giving a medication too frequently or in excess of recommended dosing) and adverse effects were excluded. We chose the 2006-2010 period as it corresponded to US Census collection of demographic data, in particular, the number of children by age group. We excluded records if the reported ZIP codes were located outside of Pennsylvania or if the US Census reported no children < 5 years. The primary outcome was the classification of a ZIP code in a cluster of drug exposure, (exposure cluster).

When a call regarding an exposure was placed to the poison center, a poison center specialist assessed the clinical signs, symptoms, and substance with a pre-determined practice guideline. If the exposure was deemed concerning, the specialist referred the child to a healthcare facility (HCF) for a medical evaluation. Occasionally, the poison center received a call when the child’s caregivers self-referred for a medical evaluation and were already en-route to a HCF. We grouped both of these referral types in the category of drug exposure with HCF referral (referral cluster).

Spatial Data Preparation

We used ArcGIS 10.2 Software (Esri Inc., Redlands, CA) to geocode (map) all calls according to the five digit ZIP codes reported to the Pittsburgh Poison Center by the caller at the time of the exposure. We created an incidence rate map for 1) all drug exposures and 2) the subset of drug exposures resulting in HCF referrals. The denominator for these incidence rates was the total population of children < 5 years in the respective ZIP code as reported by the US Census. We downloaded ZIP code shapefiles from ArcGIS Online and obtained demographic variables from 2006-2010 US Census summary files. We apportioned demographic variables to the ZIP code level. To characterize a ZIP code as urban or rural, we superimposed US Census urban shapefiles were superimposed on ZIP code centroids (centers). ZIP codes with centroids that fell within the urban area were designated as urban, while those that did not were characterized as a rural.

Spatial Analysis

We applied the spatial scan statistic in SatScan v9.3 (Kulldorff M. and Information Management Services, Inc, www.satscan.org) to identify nonrandom clusters. For this study, the null hypothesis assumed drug exposures were randomly distributed in geographic space and expected exposure counts were proportional to the population at risk. Using a Poisson-based model and adjusting for the underlying spatial inhomogeneity of the background population, the statistic placed a circular scanning window at each location and calculated the number of observed and expected observations both in and outside the window. The area with the maximum likelihood was designated the primary cluster; in other words, the most likely cluster and the least likely to have occurred by chance. A cluster with p-value <0.05 was considered to be statistically significant. 15,16 A circular scanning window was placed at different coordinates with maximum size set to contain 50%, 40%, 30%, 20%, 10% and 5% of the population at risk. The statistic performed 999 simulations and did not allow overlapping clusters.

Population Characteristics Analysis

We performed descriptive statistics to characterize the demographics of the sample region. We applied the chi-squared test to compare categorical variables and the Mann-Whitney U test to compare non-parametric continuous variables. We performed a backwards stepwise logistic regression in SPSS (Armonk, NY) to calculate odds ratios (ORs) for the association of population characteristics with clusters of drug exposure and clusters of drug exposure associated with HCF. The variables used in these models included education, poverty, unemployment, race, household size, income, women head of household, population density and the classification of a ZIP code as urban. The education variable is reported as the percentage of the population without a high school diploma and is a reflection of the high school dropout rate. Selection of these variables was based on plausible associations with drug exposure. Variables associated with the outcome in unadjusted bivariate analyses were selected by using a backward elimination logistic regression model. An iterative process was used to assess for multicollinearity. Education, poverty, unemployment, race, household size, income, women head of household, and population density were ranked in tertiles. Insignificant variables (p>0.2) were eliminated from the model. The Hosmer-Lemeshow test was used to assess model fit. University of Pittsburgh’s Institutional Review Board granted approval for this study.

RESULTS

From January 1st, 2006 to December 31st, 2010, there were 94,081 unique records involving children < 5 years reported to the Pittsburgh Poison Center from a non-healthcare facility (Figure 1). There were 390 calls from outside the service area; 99 calls were located in 20 ZIP codes where no children < 5 years were reported by 2010 US Census. The final sample included 26,685 exposures with 877 drug exposures resulting in HCF referral (Table 1). The study area included 280,855 children < 5 years across 873 ZIP Codes. A scanning window set at 50% identified a large primary cluster occupying a large portion of southwest Pennsylvania. By restricting the window size to 5%, smaller clusters were identified with minimum risk of missing significant regional-level clusters. We identified 22 exposure clusters and 5 referral cluster (Figure 2). The primary cluster from each of these samples was identified in distinctly different geographic regions; the primary exposure cluster was located in Warren, PA (RR=1.7), while the primary HCF referral cluster was identified in Erie, PA (RR=37.5). While 9.5% of the 324 ZIP codes in the exposure cluster were also identified in a referral cluster (n = 31), 65% of 89 ZIP codes in referral clusters were not initially identified in the exposure clusters, demonstrating minimally overlapping geographic pattern (n =58).

ZIP codes within exposure clusters were urban and more densely populated compared to ZIP codes not in an exposure cluster (Table 2). These ZIP codes were also statistically different in terms of median household size, educational attainment, and women head of household; the results from the regression analysis are depicted in Table 3.

DISCUSSION

We hypothesized pediatric drug exposures exist in clusters and are not randomly distributed across western Pennsylvania. Our principal findings demonstrate pediatric drug exposures do, in fact, exist in distinct geographic patterns. Moreover, ZIP codes in referral clusters were geographically distinct from the exposure clusters and were characterized as largely rural, impoverished, and with high unemployment and school dropout rates.

These results have significant implications as they provide a starting point for further investigation into the relationship between location and the social environment of children with unintentional drug poisoning. While previous studies have applied spatial analysis to understand patterns in adult prescription drug use, we found a dearth of studies using spatial techniques to investigate geographic patterns of medication poisoning in children.1720 Moreover, there is a critical gap in our understanding of the socioeconomic context of American childhood poisoning.2124 This investigation’s demographic characterization offers some first insights to the socio-economic context of American childhood poisoning. The regression model indicates population density had a moderate effect on a ZIP code falling within an exposure cluster while education status appears to have a small effect. These findings are consistent with previous studies identifying population density as a positive predictor of poison center utilization, while less educated populations are associated with a decrease in poison center utilization. Poison center utilization is influenced by poison center awarenesss.30 In more densely populated and more educated communities, knowledge of poison center service may proliferate more easily, resulting in the increase of a population’s awareness of poison center services. This may translate into more calls regarding drug exposures.

International literature has investigated the relationship of socioeconomic disparity and childhood poisoning.25,26 In the United Kingdom, children from lower socioeconomic status were at greater risk of death associated with poisoning when compared to children from more affluent status. Lower socioeconomic status may be a surrogate marker of suboptimal child supervision which might result in neglect and increased child injury.2729 Our findings suggest a similar relationship may exist in impoverished, rural American communities with high unemployment and school dropout rates.

This study offers a novel strategy in closing this gap in the poison prevention literature. Understanding where clusters exist at a local level offers target areas for revised prevention strategies. For example, within these high risk areas, poison injury resources may be dispensed at the point-of-sale of prescription medications. Professional associations might create promotional campaigns to educate pharmacists in directing customers to keep medications in original child-resistant containers; pediatricians can be prompted to discuss poison prevention tactics with families of young children living in identified cluster areas. Neighborhood and community social networks in high risk areas can be encouraged to organize “Take Back Medications” campaigns to safely and properly dispose of prescription medications. Drug exposure and drug poisoning cluster identification can also be useful to health departments by providing a means to measure the impact and performance of any poison intervention programs.

An important concept to consider is that the data sampled are from phone calls placed to the poison center; these calls may not necessarily reflect true exposure incidence, but rather a population’s utilization of poison center services. Decreased poison center utilization has been associated with increased healthcare utilization. Moreover, poison center underutilization in rural areas is believed to be associated with increased hospitalization.31,32 Our findings support a similar pattern; if exposure clusters serve as a proxy for poison center utilization, then there exists a distinct pattern of referral clusters in rural communities, outside of these areas of high poison center use. Poison Center specialists do not consider proximity to a health care facility in the decision to refer a child for a medical evaluation. Specialists gather information based on pre-determined criteria and refer when clinical concern is high. In order to gain insight on this geographic disparity, future research will investigate potential differences in the type of medications that urban and rural children are exposed to.

While identifying locations of high rates for exposure allows important characterization of population attributes associated with increased poison center awareness, the locations with low rates are also of interest. Studies have indicated increased poison center utilization is associated with reductions in ED visits and hospital admissions as poison center specialists may facilitate safe medical management of low-risk poison exposures at home.32,33 Therefore, raising poison center awareness in the ZIP codes in low rate cluster locations may lead to reductions in healthcare costs associated with unnecessary medical evaluations. Past efforts to raise poison center awareness such as mass media mailings, however, have been challenging and expensive.3436 A cost-effective strategy may be to use spatial analysis to perform a community-based assessment to identify communities within low utilization clusters.

Limitations

These results offer evidence of geographic patterns in calls placed to a single regional poison center. While it is unknown whether these findings are generalizable across populations, this technique may be easily applied to other regional centers. It is also important to note the geographic scale in this analysis. Due to privacy concerns, poison centers do not record street addresses. Inferences should be interpreted with caution as ZIP codes contain inherent socioeconomic heterogeneity, suggesting an ecologic bias. Nonetheless, this information remains relevant as it provides insight to the location and environmental risk factors surrounding pediatric drug exposures. Moreover, injury prevention strategies are often aimed at aggregate areas of at-risk individuals, such as ZIP codes. The impact of these interventions is often measured in geopolitical units (e.g. rate reductions per ZIP code). Therefore, in this case, the ecological bias becomes less concerning.

Another caveat to consider is the poison center data does not capture all drug exposures. This reporting bias may be the product of limitations in access (e.g., no telephone), lack of awareness of the poison center, poisoning severity (emergency calls are placed to 911 rather than the poison center) and language barriers. In this study, we aimed to characterize the population at the origin of exposure; thus, reports from healthcare facilities were excluded from this analysis since the ZIP code associated with these calls do not necessarily reflect the location of the exposure. Therefore, our data may underestimate the true occurrence of drug exposures and the associated healthcare utilization.

For time period of this study, approximately 2.6% calls PPC accepted were from outside the PPC service area. These calls were largely cellphones from all over the country. These calls were random and not necessarily areas adjacent to the service area. 1.5% of calls were misrouted to another poison center. While a limitation, these numbers are small and are unlikely to significantly impact our results.

In terms of cluster identification, there exist several techniques and no single technique is identified as the gold standard. The SatScan technique is popular as it is freely available and is relatively easy to use. It is limited in that it may not detect “holes”, or areas that may not be clusters, within the clusters. It is also dependent on the user’s designation of window size. By designating a smaller window size, users may be too restrictive and miss potential clusters. While our data demonstrate certain variables increase the odds for potential clustering, the analysis does not imply causality.

CONCLUSIONS

We conclude that drug exposure and drug exposure with HCF referral exist in distinct geographical clusters. ZIP codes identified within these clusters have specific population characteristics. This study offers a starting point for subsequent investigations into the geographical and social patterns of pediatric drug exposures and associated healthcare utilization.

Funding Source: The project was supported by the University of Pittsburgh Clinical and Translational Science Institute (CTSI) through the National Institutes of Health through Grant Number UL1TR000005.

Financial Disclosures: The authors have no financial relationships to disclose relevant to this article.

Conflicts of Interests: The authors have no conflicts of interests to disclose relevant to this article.

Contributor’s Statement Page:

Margaret Nguyen, MD: conceptualized and designed the study, performed data collection, data processing, geospatial analysis and interpretation, drafted the initial manuscript, made revisions to the manuscript and approved the final manuscript as submitted.

Anthony Pizon, MD: Dr. Pizon contributed to the data collection, participated in data analysis and interpretation and made critical revisions to the manuscript.

Charles Branas, PhD: Dr. Branas contributed to the design of the spatial analysis, participated in data analysis and interpretation, made critical revisions to the manuscript, and approved the final manuscript as submitted.

Anthony Fabio, PhD: Dr. Fabio contributed to the design of the study, participated in data analysis and interpretation, made critical revisions to the manuscript, and approved the final manuscript as submitted.

DECLARATION OF INTEREST

The authors report no declarations of interest.

ReferencesBondGRWoodwardRWHoMThe Growing Impact of Pediatric Pharmaceutical PoisoningJ Pediatr22012160226570e121920539SetlikJHoMBondGREmergency department use after pediatric pharmaceutical ingestion: comparison of two national databasesClin Toxicol Phila Pa12010481647SchillieSFShehabNThomasKEBudnitzDSMedication Overdoses Leading to Emergency Department Visits Among ChildrenAm J Prev Med92009373181719666156LovegroveMCMathewJHamppCGovernaleLWysowskiDKBudnitzDSEmergency Hospitalizations for Unsupervised Prescription Medication Ingestions by Young ChildrenPediatrics11020141344e10091625225137SpillerHABeuhlerMCRyanMLBorysDJAleguasABosseGMEvaluation of changes in poisoning in young children: 2000 to 2010Pediatr Emerg Care520132956354023603655National Center for Injury Prevention and ControlCenters for Disease Control and Prevention (producer)National Estimates of the 10 Leading Causes of Nonfatal Injuries Treated in Hospital Emergency Departments, United States – 2011Centers for Disease Control and Prevention. Web-based Injury Statistics Query and Reporting System (WISQARS)[Online]. [Internet]. 2011 [cited 2013 Jan 15]. Available from: www.cdc.gov/ncipc/wisqarsFranklinRLRodgersGBUnintentional Child Poisonings Treated in United States Hospital Emergency Departments: National Estimates of Incident Cases, Population-Based Poisoning Rates, and Product InvolvementPediatrics1122008122612445119047241SchoenewaldSRossSBloomLLinLPatelMShahMLynchJKuffnerENew Insights into Root Causes of Accidental Unsupervised Ingestions (AUIs) of Over-the-Counter (OTC) MedicationsClin Toxicol9200947702652009 47 702–765BudnitzDSLovegroveMCThe Last Mile: Taking the Final Steps in Preventing Pediatric Pharmaceutical PoisoningsJ Pediatr220121602190222056349BudnitzDSSalisSPreventing Medication Overdoses in Young Children: An Opportunity for Harm EliminationPediatrics1620111276e1597921555494BellNSchuurmanNGIS and Injury Prevention and Control: History, Challenges, and OpportunitiesInt J Environ Res Public Health320107310021720617015SudakinDPowerLERegional and temporal variation in methamphetamine-related incidents: applications of spatial and temporal scan statisticsClin Toxicol320094732437GruenewaldPJFreisthlerBRemerLLaScalaEATrenoAEcological models of alcohol outlets and violent assaults: crime potentials and geospatial analysisAddiction200610156667716669900The Childhood Lead Poisoning Prevention Program Geographic Information System WorkgroupUsing GIS to Assess and Direct Childhood Lead Poisoning PreventionGuidance for State and Local Childhood Lead Poisoning Prevention Programs [Internet]2004[cited 2013 Feb 3]. Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1497528/pdf/12766217.pdfKulldorffMA spatial scan statisticCommunications in Statistics: Theory and MethodsKulldorffMHuangLPickleLDuczmalLAn elliptic spatial scan statisticStat Med30112006252239294316435334DasguptaNDavisJJonsson FunkMDartRUsing Poison Center Exposure Calls to Predict Methadone Poisoning DeathsPLoS ONE [Internet]197201277cited 2014 Mar 5Available from: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3400615/BrownsteinJSGreenTCCassidyTAButlerSFGeographic Informations Systems and Pharmacoepidemiology: Using spatial cluster detection to monitor local patterns of prescription opioid abusePharmacoepidemiol Drug Saf620101966273720535759GreenTCGrauLECarverHWKinzlyMHeimerREpidemiologic trends and geographic patterns of fatal opioid intoxications in Connecticut, USA: 1997 - 2007Drug Alcohol Depend1620111153221821131140NkhomaETEd HsuCHuntVIHarrisAMDetecting spatiotemporal clusters of accidental poisoning mortality among Texas counties, U.S., 1980 - 2001Int J Health Geogr2710200432515509301UzielYAdlerAAharonowitzGFrancoSFainmesserPWolachBUnintentional childhood poisoning in the Sharon area in Israel: a prospective 5-year studyPediatr Emerg Care420052142485115824684FlanaganRJRooneyCGriffithsCFatal poisoning in childhood, England & Wales 1968–2000Forensic Sci Int320051482-3121915639606RobertsICause specific social class mortality differentials for child injury and poisoning in England and WalesJ Epidemiol Community Health16199751333459229067BaileyJECampagnaEDartRCRADARS System Poison Center InvestigatorsThe underrecognized toll of prescription opioid abuse on young childrenAnn Emerg Med420095344192418774623AlwashRMcCarthyMAccidents in the home among children under 5: ethnic differences or social disadvantage?BMJ21519882966634145033132289SibertRStress in families of children who have ingested poisonsBr Med J1271975359758791139245ZielinskiDSChild maltreatment and adult socioeconomic well-beingChild Abuse Negl10200933106667819811826GillhamBTannerGCheyneBFreemanIRooneyMLambieAUnemployment rates, single parent density, and indices of child poverty: Their relationship to different categories of child abuse and neglectChild Abuse Negl2199822279909504211HapgoodRKendrickDMarshPHow well do socio-demographic characteristics explain variation in childhood safety practices?J Public Health19200022330711LitovitzTBensonBEYounissJMetzEDeterminants of U.S. poison center utilizationClin Toxicol Phila Pa6201048544957LoVecchioFCurrySCWaszolekKKlemensJHovsethKGloganDPoison control centers decrease emergency healthcare utilization costsJ Med Toxicol12200844221419031372ZaloshnjaEMillerTJonesPLitovitzTCobenJSteinerCThe Potential Impact of Poison Control Centers on Rural Hospitalization Rates for PoisoningPediatrics11120061185209410017079583ZaloshnjaEMillerTJonesPLitovitzTCobenJSteinerCThe impact of poison control centers on poisoning-related visits to EDs—United States, 2003Am J Emerg Med32008263310518358942KrenzelokEPMrvosRIs mass-mailing an effective form of passive poison center awareness enhancement?Vet Hum Toxicol62004463155615171496EversonGRondeauESKendrickMGarzaIIneffectiveness of a mass mailing campaign to improve poison center awareness in a rural populationVet Hum Toxicol4199335216578470364YudizkyMGrisemerPShepherdGRayMGarrisonJCan textbook covers be used to increase poison center utilization?Vet Hum Toxicol102004465285615487657

Flowchart of record selection process.

Cluster Patterns of Pediatric Drug Exposures, Pittsburgh Poison Center, 2006-2010

Spatial analysis demonstrates distinct exposure clusters and referral clusters. p<0.05 for all clusters.

Unintentional Drug Exposures for Children < 5 Years, Pittsburgh Poison Center, 2006-2010
Total Exposures26,685
Ages n, (%)
 0-6 months218 (0.8)
 7-12 months2253 (8.4)
 13-24 months17,536 (65.7)
 25 to <60 months6,678 (25.0)
Sex
 Female12,850 (48.2)
 Male13,833 (51.8)
 Unknown2 (<0.1)
Substance Number
 Single drug exposure25,757 (96.5)
 Multiple drug exposure928 (3.5)
Case Management
 Managed on site (no medical evaluation necessary)25,788 (96.6)
 Poison center specialist referred to health care facility772 (2.9)
 Self-referred to health care facility105 (0.4)
 Other/Unknown20 (0.1)
Population Characteristics of ZIP Codes With and Without Geographic Clustering
Calls for Drug ExposureCalls for Drug ExposureWith Healthcare FacilityReferral

ClusterNoClusterpClusterNoClusterp
ZIP Codes, n32454989784
Household size2.362.45<.0012.422.410.25
Income$33,846$34,6530.28$29,545$35,002<.001
No high school diploma, %1213<.0011512<.001
Non White, %2.22.00.421.62.1<.001
Population per square mile200.679.8<.00181.2110.10.12
Poverty, %10100.191310<.001
Unemployed, %770.7197<.001
Urban*, %66.752.7<.00143.859.40.005
Women head of household, %1413<.00114130.76
Association Between Population Characteristics and Clusters
Population CharacteristicExposure ClusterReferral Cluster

OR (95%CI) P value OR (95%CI) P value
Population density1.67 (1.34-2.01)<.001 ...
No high school diploma0.78 (0.64-0.94)0.01 ...
Household size0.56 (0.46-0.68)<.0011.4 (1.01-1.96)0.04
Non-White0.64 (0.52-0.79)<.0010.72 (0.53-0.99)0.04
Unemployment...3.2 (2.2-4.7)<.001
Income...0.34 (0.23-0.51)<.001

Notes. CI, confidence interval; OR, odds ratio. Both models included poverty, women head of household and urban classification as covariates. Ellipses indicate covariates not included in final model. Hosmer-Lemeshow goodness-of-fit p values were 0.58 and 0.03 respectively.