Eberhart, Michael G. MPH*; Yehia, Baligh R. MD, MPP, MSHP†,‡; Hillier, Amy PhD§; Voytek, Chelsea D. MPH‖; Blank, Michael B. PhD‖; Frank, Ian MD¶; Metzger, David S. PhD‖,#; Brady, Kathleen A. MD*,†
Early initiation of antiretroviral therapy (ART) for people living with HIV (PLWH) is beneficial to both the individual's health and for reducing the likelihood of HIV transmission to others.1,2 Successful HIV treatment as prevention requires individuals to be tested, receive their test results, be linked and retained in care, initiate HIV therapy, and achieve viral load suppression.3–5 Understanding the dynamics of this treatment cascade is essential to control HIV transmission on a community level.3–5
Local and federal agencies have used surveillance data to evaluate the treatment cascade by age, sex, race/ethnicity, and HIV transmission risk.6–9 However, limited data exist on how geographic factors impact access to and retention in ARV treatment and suppression of viral load. Geographic information system (GIS) technology allows for mapping and analyses able to identify geographic foci or hot spots of disease and have been effectively used to map the burden of tuberculosis, syphilis, and HIV infection in communities.10–12 In addition, GIS technologies have been used to analyze measures of proximity and access to resources.13–15 This methodology may also be useful for monitoring the HIV treatment cascade by identifying geographic areas with poor engagement in care and inadequate viral suppression.
As part of a Centers for AIDS Research supplement to support the Enhanced Comprehensive HIV Prevention and Planning initiative, investigators from the Philadelphia Department of Public Health (PDPH) and the University of Pennsylvania used GIS technologies to examine the current HIV treatment cascade (linkage to care, retention in care, and viral suppression) status for people originally diagnosed with HIV infection between 2008 and 2009 in Philadelphia, PA. To assist in public health planning of HIV care and case management services, we sought to identify geographic areas associated with not linking to care, not linking to care in 90 days, not retaining in care, and not achieving viral suppression after HIV diagnosis by comparing spatial patterns of a cohort of recently diagnosed cases.
Data Source and Study Population
Data for analyses reported here were extracted from the City of Philadelphia's Enhanced HIV/AIDS Reporting System (eHARS), a database containing information on all HIV cases reported to the PDPH AIDS Activities Coordinating Office Surveillance Unit. Philadelphia has mandatory name-based case reporting of all new HIV infections in the City. Additionally, local mandates require reporting of all CD4 cell counts <350/mL (or CD4 percent <25%) and all HIV-1 RNA levels to the PDPH.16 Thus, eHARS contains records and laboratory results of all PLWH who were diagnosed with HIV in Philadelphia, were a resident of Philadelphia at any time after their HIV diagnosis, and all PLWH who received care in Philadelphia after their HIV diagnosis.
The eHARS database contains information collected through medical record abstraction including identifiers, such as name, address, date of birth, and address at diagnosis, as well as laboratory, pharmacy, and health service utilization information. Death data from the Pennsylvania Bureau of Vital Statistics, Social Security Death Master Index, and the National Death Index are routinely matched to eHARS data to identify deceased persons and document cause of death when available. The eHARS data are routinely monitored to identify duplicate cases and undergo quality control and verification to ensure that abstracted data are correctly assigned to unique case records.
As this study required assigning a spatial location to each person, individuals were included if they had a: (1) new HIV diagnosis date in 2008 or 2009 and (2) Philadelphia address at the time of diagnosis. Those with invalid or insufficient address data and with prison addresses at the time of diagnosis were excluded from analyses (N = 157).
Predictor and Outcome Variables
For each person, we defined age, sex at birth, race/ethnicity, HIV transmission risk, and insurance status at the time of diagnosis. Age was divided into 3 groups: <25, 25–44, and ≥45 years. Race/ethnicity was categorized as white, non-Hispanic black, Hispanic, or other. Transmission risk was grouped into heterosexual, men who had sex with men (MSM), injection drug use (IDU), and other/unknown. Patients who had IDU in combination with another risk factor (eg, MSM, heterosexual transmission) were classified as IDU.
In addition, we collected information on imprisonment (dichotomously defined as 1 or more prison stays during the observation period), visits to multiple care sites (dichotomously defined as accessing 2 or more clinics during the observation period), and proximity to medical care (dichotomously defined as distance to the nearest care site less than the population average distance to the nearest care site). In Philadelphia, laboratory results are assigned a unique identifier indicating the facility associated with the requesting medical provider. This unique identifier was used to classify imprisoned patients and those accessing 2 or more clinics during the observation period.
Outcomes of interest were based on 4 steps of the HIV treatment cascade and included: (1) not linked to care, (2) not linked to care within 90 days, (3) not retained in care, and (4) not achieving viral suppression. Linkage to care was defined as documentation of 1 or more CD4 or viral load test results after the date of diagnosis. Linkage to care within 90 days described patients with 1 or more CD4 or viral load test results in the 90-day period after diagnosis (calculated as the difference between date of diagnosis and date of first laboratory test). Retention in care was defined using the National Quality Forum Medical Visit Frequency Measure.17 This measure defines retention in care as completing at least 1 medical visit with a provider with prescribing privileges in each 6-month interval of the 24-month measurement period, with a minimum of 60 days between medical visits. Date of first linkage to care defined the start of the 24-month measurement period. We used laboratory reports of CD4 and/or viral load testing as a proxy for HIV medical care visits. Previous studies have shown high correlation between laboratory test and medical visit data and retention in care.18 Viral suppression was classified as evidence of HIV-1 RNA <200 copies per milliliter closest to the end of the 24-month measurement period ±120 days.
Mapping and Spatial Analyses
Residential address data for persons meeting inclusion and exclusion criteria were imported to ArcGIS 9.3 for geocoding using ArcMap and address locator data provided by the PDPH Department of Technology. Persons were assigned spatial coordinates for subsequent analyses based on geographic location of residence at the time of diagnosis.
In its most straightforward application, the K-function provides an estimate of spatial dependence and is based on all the distances between events in a given study area.19 This function compares the actual spatial location of points with a simulated random distribution of points using 1 or more distance bands. The value of the K-function determines whether a given point pattern is more clustered or dispersed than a point pattern of complete spatial randomness. Unlike other tests that examine variations in the intensity of phenomena over a given area, either by dividing an area into segments or calculating densities, the K-function is able to use multiple distance bands to detect clustering or dispersion at different scales. The “local” K-function, which calculates a value at each point, rather than for the entire region, also identifies whether a point pattern is more clustered or dispersed but can additionally identify specific locations in the region (“hot spots” and/or “cool spots”) where the distribution of points differs significantly from the random distribution.20–22 The cross K-function uses a marked point process and can be used to compare subgroups (eg, persons linked vs. not linked to care) of a single population to the distribution of all points in the region. This test can be used to compare points marked as “not linked to care” with all points in the region to determine if the distributions are significantly different or spatially indistinguishable.20–22 The local version of the cross K-function can identify specific locations in the region where points marked as not linked to care are more clustered than all points in the region.
We used the local cross K-function to determine if: (1) the spatial patterns exhibited by persons that did not meet specific steps along the HIV treatment cascade differed significantly from the overall pattern of people diagnosed with HIV; and if so, (2) where in the City those differences were observed.
The spatial location of each case was analyzed using 3 radial distance bands (1000, 2500, and 5000 feet) to determine if points in marked pattern 1 (eg, persons not linked to care) were significantly more clustered than points in marked pattern 2 (eg, persons linked to care) for each of the 4 outcomes. Feet were used as the distance unit because the coordinates were calculated in feet. The primary distance radius of 5000 feet was selected based on one-half of the maximum distance between each point and its 5 nearest neighbors. The distance of 5000 feet also encompasses the nearest neighbor distance for 99.9% of the cohort and represents a rough approximation of 1 mile.
Cross K-function P-values were calculated using MATLAB and publicly available programs (designed by Tony E. Smith: http://www.seas.upenn.edu/∼ese502/NOTEBOOK/Part_I/5_Comparative_Analyses.pdf) to detect local variations of a marked point process. The P-values at each spatial location were then imported to ArcMap and interpolated to a continuous raster surface file using the spline method in ArcToolbox. The interpolated surface raster was then converted to contour lines to highlight specific geographic areas where persons classified with negative outcomes along the HIV treatment cascade exhibited significant clustering compared to the full cohort of patients.
Our a priori hypothesis was that each step in the cascade is a distinct outcome with independent predictors that include both individual and community factors. As a result, the chosen denominator for each outcome was the numerator of the previous step in the cascade; all eligible diagnoses for linked to care, all persons linked to care for linked within 90 days and retention in care, and retained in care for viral suppression.
Univariate statistics were used to describe the dataset. Multivariate logistic regression models were used to assess relationships between predictors and the outcomes. Models were adjusted for age, sex at birth, race/ethnicity, HIV transmission risk, and insurance status at the time of diagnosis, as well as imprisonment status, visits to multiple care sites, and proximity to nearest HIV medical provider. Additionally, we included geographic hot spot areas in the final model. To do this, contour hot spots were spatially joined to each person to assign a value indicating the distance between each person and the nearest contour. Persons were considered to be within the contour area if the calculated distance value was less than or equal to 5000 feet. Persons were assigned as either residing inside or outside of a particular area for each of the 4 outcomes. Adjusted odds ratios (AORs) with 95% confidence intervals (CIs) are presented. Relationships were considered statistically significant at P < 0.05.
Between 2008 and 2009, 1861 people newly diagnosed with HIV were identified (Table 1). Of these, 157 (8%) had invalid address data or were imprisoned at the time of diagnosis and were excluded from subsequent analyses. Excluded persons were less likely to be black or Hispanic; and more likely to be older adults, people with IDU as an HIV risk factor, and people with private insurance compared to included individuals. In addition, excluded patients were more likely to have had a prison stay during the subsequent observation period. Among the 1704 persons included in the analyses, 70% were male, 63% were black, and 30% were 45 years or older at the time of diagnosis. The 2 most common HIV transmission risk factors were heterosexual (40%) and MSM (36%).
Overall, 1404 persons (82% of 1704) linked to care. Among those linked, 1059 (75%) linked within 90 days and 526 (37%) were retained in care. Of the 526 individuals retained in care, 379 (72%) achieved viral suppression. Correspondingly, 18%, 25%, 63%, and 28% of eligible individuals were not linked to care, not linked to care within 90 days, not retained in care, and not virally suppressed, respectively (Table 2).
As shown in Figure 1, Philadelphia is a city of neighborhoods bordered on the east by the Delaware River, which separates Pennsylvania and New Jersey, and to the north and west by 3 Pennsylvania counties: Montgomery, Delaware, and Bucks. The maps presented in Figure 2 show the geographic distribution of new HIV diagnoses (A1), HIV care sites (A2), and the areas identified as hot spots for each of the 4 outcomes (A3–A6). Overall, 12–19 census tracts were included in the geographic clusters for each outcome, with a total of 59 unique tracts (15.4% of all census tracts in Philadelphia).
New HIV cases tended to cluster in areas of high population including sections of North Philadelphia, Center City, South Philadelphia, and West Philadelphia. Geographic areas associated with not linking to care included North Philadelphia, as well as smaller segments of West Philadelphia and the Lower Northeast, and a long strip along the City's western border with Montgomery County in the Roxborough section. Areas associated with not linking to care within 90 days include a small area of North Philadelphia, 2 small areas in West Philadelphia, and a long strip along the river in the Kensington section. Areas associated with not retaining in care are both more numerous and more geographically varied, including a large area in northern part of West Philadelphia near Fairmount Park, a smaller section in the neighborhood of West Philadelphia, an area at the southern end of South Philadelphia, a strip along the river in Center City, 2 sections of Oak Lane along the border with Montgomery County, and an area in the central section of the Lower Northeast. Three areas of the City were associated with not achieving viral suppression, including the Kensington section of North Philadelphia, a small area of Southwest Philadelphia, and the tip of Northwest Philadelphia in the Olney/Oak Lane neighborhood.
Logistic regression models were fit to evaluate the relationship between geographic hot spots (individually identified for each of the 4 steps evaluated in the HIV treatment cascade) and the outcomes of interest (Table 3).
Model 1—Not Linked to Care
Persons residing inside of the geographic areas identified by the local K-function analysis were more likely to not link to care compared to those residing outside of these areas [AOR 1.76 (95% CI: 1.30 to 2.40)]. Additionally, blacks (vs. whites), persons with IDU transmission risk (vs. heterosexual), and those with Medicare or no insurance (vs. private) had higher odds of nonlinkage. Males were more likely to be not linked to care than females. No significant differences in linkage to care were detected between whites and Hispanics.
Model 2—Not Linked to Care Within 90 Days
Among persons linked to care, individuals residing inside of the geographic areas identified by the local K-function analysis were more likely to not link to care within 90 days compared to those residing outside of these areas [AOR 1.49 (95% CI: 1.12 to 1.99)]. Compared to persons with private insurance, individuals with no insurance had higher odds of nonlinkage within 90 days (ie, less likely to link within 90 days). No differences were observed by sex, age, race/ethnicity, transmission risk, prison stay during the observation period, or proximity to medical care sites.
Model 3—Not Retained in Care
Among persons linked to care, individuals residing inside of the geographic areas identified by the local K-function analysis were more likely to not be retained in care compared to those residing outside of these areas [AOR 1.84 (95% CI: 1.39 to 2.43)]. Blacks and Hispanics were more likely to be nonretained (ie, less likely to be retained) than whites, while younger individuals were more likely to be non retained compared to those 45 years and older.
Model 4—Not Virally Suppressed
Among all persons retained in care, those residing inside of the geographic areas identified by the local K-function analysis were more likely to not be virally suppressed compared to those residing inside of these areas [AOR 3.23 (95% CI: 1.87 to 5.59)]. Persons whose nearest care site was more than the average distance to the nearest care site were more likely to not achieve viral suppression (ie, less likely to be virally suppressed) compared to their counterparts who resided closer to their care sites. No differences were observed by sex, age, race/ethnicity, transmission risk, insurance, or prison stay during the observation period.
In this cohort, we identified 1704 newly diagnosed persons with HIV who met inclusion criteria. Among these, 82% were linked to care; of whom 75% linked to care within 3 months of diagnosis. Only 37% of patients linked to care were retained in care over 24 months, and 72% of those who were retained in care achieved viral suppression. Overall, 22% of persons recently diagnosed were linked to care, retained in care, and virally suppressed at the end of the follow-up period. This estimate is consistent with other published estimates of the treatment cascade and viral suppression in the United States.4 Although predictors of failing to be linked to care, retained in care, and virally suppressed have been previously evaluated,3,9,18,23–26 this study is unique in that it used community factors including geography and proximity to care as predictors.
Geographic clustering was independently associated with poor outcomes at each step along the HIV treatment cascade. We identified between 12 and 19 census tracts included in the geographic clusters for each outcome, with a total of 59 unique tracts (15.4% of all census tracts in Philadelphia) across all outcomes. Interestingly, the geographic clusters identified for each step of the cascade were unique, with no geographic overlap between cascade steps. This suggests that the community-level factors responsible for worse outcomes may differ depending on the step of the HIV treatment cascade being evaluated.
Community factors related to poverty, such as crime, housing stability, and poor access to transportation and social services, may impact HIV treatment outcomes for individuals residing in geographic hot spot areas.27–29 We identified several characteristics differentiating these hot spots from the rest of the City of Philadelphia. First, hot spots for 3 of the 4 outcomes evaluated had narcotic arrest rates at least 2 times greater than the City average.27 High narcotic arrest rates may indicate high drug use and crime in these communities. In addition, hot spots associated with not remaining in care had higher rates of aggravated assault compared to other parts of the City. Surprisingly, other markers of low socioeconomic status (poverty and education level) did not differ between hot spots and the remainder of the City. This finding may be limited by the high rate of poverty in Philadelphia overall, with nearly 45% of the general population living below 200% of the federal poverty level.28 Of the 10 largest cities in the United States, Philadelphia has the highest rate of deep poverty, defined as people with incomes below half of the federal poverty level.29 Evaluation of whether there are differences in the prevalence of deep poverty in these hot spots warrants further investigation.
Previous studies have shown that low neighborhood/community socioeconomic status is associated with a lower likelihood of having a usual source of health care, obtaining preventive services, and an increased likelihood of having unmet medical need, even after controlling for patient characteristics and supply of health care providers.30 The impact of neighborhood environment and residents' reaction to and perception of that environment has been linked to both risk of chronic disease and rates of self-management behaviors necessary for managing chronic disease.31–34 Our study did not examine access to public transportation in these areas, access to social services or pharmacies, measures of social disorder such as crime and community response to crime, and community networks and norms (especially related to health and healthcare). We hypothesize that community norms and social disorder may have a greater effect on linkage to care; access to public transportation and social services may have a greater effect on retention in care; and access to pharmacies may have a greater effect on viral suppression. Differences in community factors that influence each step of the cascade may explain the lack of overlap in hot spots.
Although, the 59 unique census tracts associated with the geographic hotspot areas represent 14.5% of the population of Philadelphia, they represent 22.2% of all PLWH in the City. A higher burden of HIV disease in these communities with poor outcomes could indicate a need for additional services in these areas. However, proximity to HIV medical care site was not associated with linkage to care, linkage to care within 90 days, and retention in care, suggesting that inaccessibility to care is not the primary driver for the geographic clustering observed. Furthermore, our previous work on geography and access to HIV care in Philadelphia has shown that individuals do not tend to use HIV care at the facility nearest to their home residence.35
Similar to previous studies, linkage to care, retention in care, and viral suppression rates differed for different sociodemographic groups.18,25,26,31,32 After adjusting for geographic hot spot areas, persons with IDU, heterosexuals, males, those with Medicare, and those without insurance were less likely to be linked to care. However, only those with no insurance coverage were less likely to link to care within 90 days. Our analysis of linkage to care differentiates between individuals who were ever linked to care from those who never linked to care (the denominator includes all people diagnosed with HIV). Whereas our analysis of linkage to care within 90 days is in essence an evaluation of early versus late linkage to care, all persons in the denominator were those linked to care. Several previous studies examining predictors of delayed linkage to care did not identify any associations with demographic factors and mode of transmission, after persons who never linked to care were excluded.36–38 One study, however, did find that persons of non-white race and IDU as mode of transmission were independently associated with delayed linkage, but the authors excluded persons with AIDS.39 The exclusion of persons with AIDS may explain the associations with non-white race and IDU since these groups are more likely to present with AIDS than other demographic groups. Several studies have shown that the strongest predictor of delayed linkage to care is the type of setting in which an individual is tested and found to be HIV positive with nonmedical sites, inpatient medical sites, or while incarcerated being associated with delayed entry.38–40 Our study did not look at the type of facility where the HIV diagnosis was made but this will certainly be a focus of future work.
Blacks and Hispanics were less likely to be retained in care, as were persons younger than 25 years at the time of their HIV diagnosis. Persons who attended multiple clinics were more likely to meet the criteria for retention. Previous studies have demonstrated that patients using multiple sites of care attend more outpatient visits compared to those receiving care at only 1 care site.41 Proximity to care and geographic area were the only variables associated with viral suppression. Unlike previous studies that have identified disparities in rates of viral suppression for racial/ethnic minorities, we did not identify any demographic differences.26,42,43 Although further research is warranted, the lack of observed racial/ethnic disparities in viral suppression in our study may in part be explained by greater demographic similarity among those retained in care combined with a smaller sample size. Among those retained in care, environmental factors may be more salient predictors of treatment success.
Our study has several limitations that should be acknowledged. First, there were statistical differences in the demographic characteristics of age, mode of transmission, and insurance status between those in the cohort and those who were excluded. We believe that the impact of these differences on our results is minimal because a high proportion of the potential cases (92%) were included in the analysis. Second, we excluded patients diagnosed while in the jail system. The incarcerated population represents a vulnerable population that should be studied independently. Their exclusion here was because of their diagnosis during incarceration, which did not reflect the geographic predictors we intended to explore. Third, our findings come from routine HIV surveillance data; as such we could not assess ART coverage, may have incompletely accounted for migration out of Philadelphia among the patients in the cohort, and may have had underreporting of CD4 counts. Although, only CD4 results >350 cells per microliter or >25% of the total T-lymphocyte count were not reportable to the PDPH during the cohort follow-up period. All HIV viral loads were reportable as mandated by law, including those that were undetectable. Therefore, linkage to and retention in care would have been accurately ascertained by receipt of a viral load. Fourth, not all persons diagnosed with HIV during 2008 and 2009 would have been eligible for ART based on DHHS treatment guidelines during the follow-up period, which may in part explain the lower-than-expected viral suppression rate. Fifth, we did not assess the impact of the density of general medical facilities, hospitals, and pharmacies, access to public transportation or other social services, and housing stability as a mechanism for the geographic clustering that we observed. Future studies should examine the impact of these community services on the HIV treatment cascade.
Finally, our study cohort included individuals in a single large urban area who tested positive for HIV in one continuous 24-month period. Although the results may not be generalizable to other areas, the methods used for the spatial analyses could easily be replicated in other jurisdictions.
Our findings have very clear prevention implications. Previous studies have identified neighborhoods with high community viral loads as areas where uninfected individuals are at greater risk of acquiring HIV, compared to neighborhoods with lower community viral load.44,45 Measuring and using community viral load is a goal of the National HIV/AIDS Strategy (NHAS) for the United States.46 NHAS also calls for innovative solutions for reducing community viral load that may help reduce the number of new HIV infections in specific communities. Our study is innovative in this regard, providing practical information on both when and where the breakdown in the care cascade occurs (which likely results in higher community viral load). Furthermore, the Centers for Disease Control and Prevention's High Impact Prevention program calls for using combinations of scientifically proven, cost-effective, and scalable interventions targeted to the right populations in the right geographic areas to increase the impact of HIV prevention efforts and meet the NHAS goals.47 The methods used in our study can help identify the right populations (those not linked to care, not retained in care, and not virally suppressed) in the right geographic areas to scale-up scientifically proven, cost-effective interventions. We believe that a jurisdiction could specifically target separate linkage, retention, and adherence interventions in the areas identified with the greatest need.
This study is the first to identify spatial patterns as a strong independent predictor of linkage to care, retention in care, and viral suppression. Moving forward, it will be important to assess other, ecological, community-level, and neighborhood infrastructure factors that may influence access to HIV medical care and treatment outcomes within the geographic clusters identified. Thus, spatial analyses are a valuable tool for characterizing the HIV epidemic and treatment cascade. Geography-based tailoring of interventions to improve each aspect of the cascade will be critical to controlling the HIV epidemic in the United States.
1. Cohen MS, Chen YQ, McCauley M, et al.. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365:493–505.
2. Attia S, Egger M, Muller M, et al.. Sexual transmission of HIV according to viral load and antiretroviral therapy: systematic review and meta-analysis. AIDS. 2009;23:1397–1404.
3. Gardner EM, McLees MP, Steiner JF, et al.. The spectrum of engagement in HIV care and its relevance to test-and-treat strategies for prevention of HIV infection. Clin Infect Dis. 2011;52:793–800.
4. Centers for Disease Control and Prevention. Vital signs: HIV prevention through care and treatment—United States. MMWR Morb Mortal Wkly Rep. 2011;60:1618–1623.
5. Kilmarx PH, Mutasa-Apollo T. Patching a leaky pipe: the cascade of HIV care. Curr Opin HIV AIDS. 2013;8:59–64.
6. Hall ELF HI, Rhodes P, Holtgrave DR, et al.. Continuum of HIV care: differences in care and treatment by sex and race/ethnicity in the United States. Paper presented at: AIDS 2012; 2012; Washington, DC.
8. Julio Montaner BN, Colley G, Samji H, et al.. The evolution of the cascade of HIV care: British Columbia, Canada, 1996-2010. Paper presented at: CROI2013; Atlanta, GA.
9. Michael Horberg LH, Towner W, Gambatese R, et al.. HIV spectrum of engagement cascade in a large integrated care system by gender, age, and methodologies. Paper presented at: CROI2013; Atlanta, GA.
10. Shepard CW, Gortakowski HW, Nasrallah H, et al.. Using GIS-based density maps of HIV surveillance data to identify previously unrecognized geographic foci of HIV burden in an urban epidemic. Public Health Rep. 2011;126:741–749.
11. Goswami ND, Hecker EJ, Vickery C, et al.. Geographic information system-based screening for TB, HIV, and syphilis (GIS-THIS): a cross-sectional study. PLoS One. 2012;7:e46029.
12. Oppong JR, Tiwari C, Ruckthongsook W, et al.. Mapping late testers for HIV in Texas. Health Place. 2012;18:568–575.
13. Leibowitz AA, Taylor SL. Distance to public test sites and HIV testing. Med Care Res Rev. 2007;64:568–584.
14. Pierce SJ, Miller RL, Morales MM, et al.. Identifying HIV prevention service needs of African American men who have sex with men: an application of spatial analysis techniques to service planning. J Public Health Manag Pract. 2007(suppl):S72–S79.
15. Bliss RL, Katz JN, Wright EA, et al.. Estimating proximity to care: are straight line and zipcode centroid distances acceptable proxy measures? Med Care. 2012;50:99–106.
18. Rebeiro P, Althoff KN, Buchacz K, et al.. Retention among North American HIV-infected persons in clinical care, 2000-2008. J Acquir Immune Defic Syndr. 2013;62(3):356–362.
19. Ripley BD. The second-order analysis of stationary point processes. J Appl Probab. 1976;13:255–266.
20. Cressie NAC. Statistics for Spatial Data. New York, NY: Wiley; 1991.
21. Diggle PJ, Chetwynd AG. Second-order analysis of spatial clustering for inhomogeneous populations. Biometrics. 1991;47:1155–1163.
22. Dixon PM. Ripley's K function. In: Encyclopedia of Environmetrics. Eds. El-Shaarawi AH, Piegorsch WW. New York, NY: Wiley; 2002:1796–1803.
23. Keller SC, Yehia BR, Eberhart MG, et al.. Accuracy of definitions for linkage to care in persons living with HIV. J Acquir Immune Defic Syndr. 2013;63(5):622–630.
24. Yehia BR, Fleishman JA, Metlay JP, et al.. Comparing different measures of retention in outpatient HIV care. AIDS. 2012;26:1131–1139.
25. Fleishman JA, Yehia BR, Moore RD, et al.. Establishment, retention, and loss to follow-up in outpatient HIV care. J Acquir Immune Defic Syndr. 2012;60:249–259.
26. Yehia BR, Fleishman JA, Metlay JP, et al.. Sustained viral suppression in HIV-infected patients receiving antiretroviral therapy. JAMA. 2012;308:339–342.
27. Cartographic Modeling Lab. University of Pennsylvania. Philadelphia Neighborhood Information System. CrimeBase. Available at http://www.cml.upenn.edu
. Accessed May 6, 2013.
28. United States Census Bureau. 2006-2010 American Community Survey. US Census Bureau's American Community Survey Office, 2013. Available at http://factfinder2.census.gov
. Accessed July 1, 2013.
29. Lubrano A. Of big cities, Phila. Worst for people in deep poverty. Philadelphia Inquirer. 2013;2013:A1.
30. Kirby JB, Kaneda T. Neighborhood socioeconomic disadvantage and access to health care. J Health Soc Behav. 2005;46:15–31.
31. Brown AF, Ang A, Pebley AR. The relationship between neighborhood characteristics and self-rated health for adults with chronic conditions. Am J Public Health. 2007;97:926–932.
32. Gary TL, Safford MM, Gerzoff RB, et al.. Perception of neighborhood problems, health behaviors and diabetes outcomes among adults with diabetes in managed care. Diabetes Care. 2008;31:273–278.
33. Wilson DK, Kliewer W, Sica DA. The relationship between exposure to violence and blood pressure mechanisms. Curr Hypertens Rep. 2004;6:321–326.
34. Engelman KK, Hawley DB, Gazaway R, et al.. Impact of geographic barriers on the utilization of mammograms by older rural women. J Amer Geriatr Soc. 2002;50:62–68.
35. Eberhart MG. Travel Distance to Medical Care: A Geographic Analysis of Weighted Survey Data from the Medical Monitoring Project in Philadelphia, PA. First National CFAR/APC ECHPP Conference, November 19–20, 2012, 2012; Washington, DC.
36. Ickovics JR, Forsyth B, Ethier KA, et al.. Delayed entry into health care for women with HIV disease. AIDS Patient Care STDS. 1996;10:21–24.
37. Muthulingam D, Chin J, Hsu L, et al.. Dispartities in engagement in care and viral suppression among persons with HIV. J Acquir Immune Defic Syndr. 2013;63:112–119.
38. Jennes SM, Myers JE, Neaigus A, et al.. Delayed entry into HIV medical care after HIV diagnosis: risk factor and research methods. AIDS Care. 2012;24:1240–1248.
39. Torian LV, Wiewel EW, Liu K, et al.. Risk factors for delayed initiation of medical care after diagnosis of human immunodeficiency virus. Arch Intern Med. 2008;168:1181–1187.
40. Bamford LP, Ehrenkranz PD, Eberhart MG, et al.. Factors associated with delayed entry into primary HIV medical care after HIV diagnosis. AIDS. 2010;24:928–930.
41. Yehia B, Schranz A, Momplaisir F, et al.. Outcomes for HIV+ patients receiving care at multiple HIV clinics. Conference on Retroviruses and Opportunistic Infections; March 3–6, 2013, 2013; Atlanta, GA.
42. Benbow N, Scheer S, Wohl A, et al.. Linkage, access, ART use and viral suppression in four large cities in the United States, 2009. AIDS. 2012;2012; Washington, DC.
43. Hall HI, Frazier EL, Rhodes P, et al.. Continuum of care: differences in care and treatment by sex and race/ethnicity in the United States. AIDS. 2012;2012; Washington, DC.
44. Das M, Chu PL, Santos G-M, et al.. Decreases in community viral load are accompanied by reductions in new HIV infections in San Francisco. PLoS One. 2010;5:e11068.
45. Castel AD, Befus M, Willis S, et al.. Use of the community viral load as a population-based biomarker of HIV burden, AIDS. 2012;26:345–353.
46. White House Office of National AIDS Policy. National HIV/AIDS Strategy for the United States. Washington, DC: White House Office of National AIDS Policy; 2010.