Email updates

Keep up to date with the latest news and content from IJHG and BioMed Central.

Open Access Highly Accessed Methodology

Comparing alternative approaches to measuring the geographical accessibility of urban health services: Distance types and aggregation-error issues

Philippe Apparicio1*, Mohamed Abdelmajid2, Mylène Riva13 and Richard Shearmur1

Author Affiliations

1 Spatial Analysis and Regional Economics Laboratory, Université du Québec, Institut national de la recherche scientifique, Urbanisation, Culture et Société, 385 rue Sherbrooke est, Montréal (Québec), H2X 1E3, Canada

2 Department of Geography, Université du Québec à Montréal, Pavillon Hubert-Aquin, 1255 rue Saint-Denis, Montréal (Québec), H2X 3R9, Canada

3 Department of Social and Preventive Medicine, Faculty of Medicine, University of Montréal, P.O. Box 6128, Downtown Station, Montréal (Québec), H3C 3J7, Canada

For all author emails, please log on.

International Journal of Health Geographics 2008, 7:7  doi:10.1186/1476-072X-7-7


The electronic version of this article is the complete one and can be found online at: http://www.ij-healthgeographics.com/content/7/1/7


Received:31 August 2007
Accepted:18 February 2008
Published:18 February 2008

© 2008 Apparicio et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Over the past two decades, geographical accessibility of urban resources for population living in residential areas has received an increased focus in urban health studies. Operationalising and computing geographical accessibility measures depend on a set of four parameters, namely definition of residential areas, a method of aggregation, a measure of accessibility, and a type of distance. Yet, the choice of these parameters may potentially generate different results leading to significant measurement errors.

The aim of this paper is to compare discrepancies in results for geographical accessibility of selected health care services for residential areas (i.e. census tracts) computed using different distance types and aggregation methods.

Results

First, the comparison of distance types demonstrates that Cartesian distances (Euclidean and Manhattan distances) are strongly correlated with more accurate network distances (shortest network and shortest network time distances) across the metropolitan area (Pearson correlation greater than 0.95). However, important local variations in correlation between Cartesian and network distances were observed notably in suburban areas where Cartesian distances were less precise.

Second, the choice of the aggregation method is also important: in comparison to the most accurate aggregation method (population-weighted mean of the accessibility measure for census blocks within census tracts), accessibility measures computed from census tract centroids, though not inaccurate, yield important measurement errors for 5% to 10% of census tracts.

Conclusion

Although errors associated to the choice of distance types and aggregation method are only important for about 10% of census tracts located mainly in suburban areas, we should not avoid using the best estimation method possible for evaluating geographical accessibility. This is especially so if these measures are to be included as a dimension of the built environment in studies investigating residential area effects on health. If these measures are not sufficiently precise, this could lead to errors or lack of precision in the estimation of residential area effects on health.

Background

In recent urban health studies, an increased focus has been directed to evaluating the accessibility of urban resources, either for individuals or for populations living in residential areas. Although the concept of "accessibility" is multidimensional (accessibility may be defined in terms of affordability, acceptability, availability and spatial accessibility [1]) evaluating geographical accessibility in residential areas offers critical information for public policy in planning and service provision as it allows for the identification of areas with lower (or higher) access to urban resources and the assessment of spatial and social inequalities in access [2,3].

Geographical accessibility refers to the ease with which residents of a given area can reach services and facilities [2]. Most common approaches for defining geographical accessibility are based on distance or travel time to a resource (for a review, please refer to [4]). These measures assume that every member of the population is a potential user of the service; the pattern of spatial accessibility will depend on the relative location of the population and services [5,6]. Table 1 synthesises approaches for conceptualizing and measuring different dimensions of geographical accessibility.

Table 1. Approaches for conceptualizing and measuring the geographical accessibility of services and facilities for residential areas

Several studies have measured the geographical accessibility in residential areas of services and facilities that have the potential to contribute to the population's well-being and health such as health care services [2,7-18], recreational facilities [2,16,18,19], and food supermarkets [16,18,20-23]. Accessibility to these types of resources is especially important for populations with limited mobility and revenue since more direct and easier access confers opportunities by reducing the time and financial costs of access, and by potentially influencing life choices [24]. Other studies have measured geographical accessibility of resources potentially associated with more negative health outcomes, such as waste facilities, fast food restaurants, and pollution from large motorways [25-27].

Over the past two decades, the operationalization of geographical accessibility measures in urban and health studies has become easier, largely due to developments in GIS transportation softwares and modules. These measures require the specification of a set of four parameters, namely 1) a spatial unit of reference for the population, i.e. a definition of residential areas; 2) an aggregation method, i.e. to account for the distribution of population in the residential area; 3) a measure of accessibility; and 4) a type of distance for computing the accessibility measures selected. The choice of these parameters is likely to generate different results, potentially leading to significant measurement errors [2,28,29].

In this paper, we investigate differences in results when geographical accessibility of residential areas (census tracts) to health care services is computed using different distances types and different aggregation methods. In the next section, we describe methods for defining the four parameters for computing accessibility measures. With these methods in mind, we then provide an overview of methodological issues in measuring accessibility.

Evaluating geographical accessibility of services and facilities in residential areas: specifying a set of parameters

Spatial unit of reference and aggregation methods

Selecting the appropriate spatial unit of analysis, i.e. the operational definition for residential areas, is critical for minimizing aggregation errors [2,30]. Aggregation error arises from the distribution of individuals around the centroid of spatial units [2]. As spatial units vary in size from smaller areas, e.g. census blocks, to larger ones, e.g. census tracts, accessibility measured for smaller units is less subject to aggregation error than that measured for larger spatial units [2].

To evaluate the geographical accessibility of a service for a population living in a residential area, e.g. a census tract, three methods can be used [2]; they are illustrated in Figure 1. The first method consists in computing the distance between the centroid of the census tract and the service (Figure 1.a). This method shows the inappropriateness of ignoring the spatial distribution of the population inside the census tract [2].

thumbnailFigure 1. Choosing the spatial unit of reference for calculating distances and error aggregation.

The second method consists of calculating the population-weighted mean centre of the census tracts (Equation 1) and then evaluating the distance between this new location and the service. Toward this end, smaller spatial units entirely contained by the census tracts can be used, such as dissemination areas, census blocks, or postal codes. This method accounts for the spatial distribution of the population inside the census tract in order to minimize aggregation error.

( x i ¯ , y i ¯ ) = ( b i w b x b b i w b , k i w b y b b i w b ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaWaaeWaaeaadaqdaaqaaiabdIha4naaBaaaleaacqWGPbqAaeqaaaaakiabcYcaSmaanaaabaGaemyEaK3aaSbaaSqaaiabdMgaPbqabaaaaaGccaGLOaGaayzkaaGaeyypa0ZaaeWaaeaajuaGdaWcaaqaamaaqafabaGaem4DaC3aaSbaaeaacqWGIbGyaeqaaiabdIha4naaBaaabaGaemOyaigabeaaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqacqGHris5aaqaamaaqafabaGaem4DaC3aaSbaaeaacqWGIbGyaeqaaaqaaiabdkgaIjabgIGiolabdMgaPbqabiabggHiLdaaaOGaeiilaWscfa4aaSaaaeaadaaeqbqaaiabdEha3naaBaaabaGaemOyaigabeaacqWG5bqEdaWgaaqaaiabdkgaIbqabaaabaGaem4AaSMaeyicI4SaemyAaKgabeGaeyyeIuoaaeaadaaeqbqaaiabdEha3naaBaaabaGaemOyaigabeaaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqacqGHris5aaaaaOGaayjkaiaawMcaaaaa@63D3@ (1)

Where:

wb = total population of spatial unit b completely within census tract i (i.e. dissemination area or census block or postal code).

xb and yb = X and Y coordinates of spatial unit b.

Finally, the third method consists of computing the distance between the services and each centroid of spatial units completely within census tracts, and then calculating the average of these distances weighted by the total population of each unit (Figure 1.b and 1c). In comparison with the previous methods, this one is more accurate because it more exactly accounts for the distribution of the population inside the census tract.

Accessibility measures

The five most commonly used measures of accessibility are: 1) the distance to the closest service, 2) the number of services within n metres or minutes, 3) the mean distance to all services, 4) the mean distance to n closest services, and 5) the gravity model. If the more accurate aggregation method detailed previously is selected, these accessibility measures can be written as:

Z i a = b i w b ( min | d b s | ) b i w b , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdggaHbaakiabg2da9KqbaoaalaaabaWaaabuaeaacqWG3bWDdaWgaaqaaiabdkgaIbqabaGaeiikaGIagiyBa0MaeiyAaKMaeiOBa42aaqWaaeaacqWGKbazdaWgaaqaaiabdkgaIjabdohaZbqabaaacaGLhWUaayjcSdGaeiykaKcabaGaemOyaiMaeyicI4SaemyAaKgabeGaeyyeIuoaaeaadaaeqbqaaiabdEha3naaBaaabaGaemOyaigabeaaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqacqGHris5aaaakiabcYcaSaaa@5225@ (2)

Where:

Z i a MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdggaHbaaaaa@2FE3@ = mean distance between census tract i and closest service.

wb = total population of spatial unit b completely within census tract i.

dbs = distance between spatial unit b and service s.

Z i b = b i W b j S S j b i W b , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaacbmGae8NwaO1aa0baaSqaaiabdMgaPbqaaiabdkgaIbaakiabg2da9maalaaabaWaaabuaeaacqWFxbWvdaWgaaWcbaGaemOyaigabeaaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqaniabggHiLdGcdaaeqbqaaiab=nfatnaaBaaaleaacqWGQbGAaeqaaaqaaiabdQgaQjabgIGiolabdofatbqab0GaeyyeIuoaaOqaamaaqafabaGae83vaC1aaSbaaSqaaiabdkgaIbqabaaabaGaemOyaiMaeyicI4SaemyAaKgabeqdcqGHris5aaaakiabcYcaSaaa@4D15@ (3)

Where:

Z i b MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdkgaIbaaaaa@2FE5@ = mean number of services within n metres or minutes of census tract i.

wb = total population of spatial unit b completely within census tract i.

S = all services.

Sj = number of services within n metres or minutes of spatial unit centroid b with Sj = 1 where dbs n and Sj = 0 where dbs > n.

Z i c = b i w b d b s b i w b , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdogaJbaakiabg2da9KqbaoaalaaabaWaaabuaeaacqWG3bWDdaWgaaqaaiabdkgaIbqabaGaemizaq2aaSbaaeaacqWGIbGycqWGZbWCaeqaaaqaaiabdkgaIjabgIGiolabdMgaPbqabiabggHiLdaabaWaaabuaeaacqWG3bWDdaWgaaqaaiabdkgaIbqabaaabaGaemOyaiMaeyicI4SaemyAaKgabeGaeyyeIuoaaaGccqGGSaalaaa@4933@ (4)

Z i c MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdogaJbaaaaa@2FE7@ = mean distance between census tract i population and all services.

wb = total population of spatial unit b completely within census tract i.

dbs = distance between spatial unit centroid b and service s.

Z i d = b i w b s d b s n b i w b , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdsgaKbaakiabg2da9KqbaoaalaaabaWaaabuaeaacqWG3bWDdaWgaaqaaiabdkgaIbqabaWaaabuaeaadaWcaaqaaiabdsgaKnaaBaaabaGaemOyaiMaem4CamhabeaaaeaacqWGUbGBaaaabaGaem4CamhabeGaeyyeIuoaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqacqGHris5aaqaamaaqafabaGaem4DaC3aaSbaaeaacqWGIbGyaeqaaaqaaiabdkgaIjabgIGiolabdMgaPbqabiabggHiLdaaaOGaeiilaWcaaa@4E24@ (5)

Z i d MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdsgaKbaaaaa@2FE9@ = mean distance between census tract i and n closest services.

wb = total population of spatial unit b completely within census tract i.

dbs = distance between spatial unit centroid b and service s, dbs is sorted in ascending order.

n = number of closest services to be included in measure.

Z i e = b i w b S S w s d b s α b i w b , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdwgaLbaakiabg2da9KqbaoaalaaabaWaaabuaeaacqWG3bWDdaWgaaqaaiabdkgaIbqabaWaaabuaeaacqWGtbWudaWgaaqaaiabdEha3jabdohaZbqabaGaemizaq2aa0baaeaacqWGIbGycqWGZbWCaeaacqGHsisliiGacqWFXoqyaaaabaGaem4uamfabeGaeyyeIuoaaeaacqWGIbGycqGHiiIZcqWGPbqAaeqacqGHris5aaqaamaaqafabaGaem4DaC3aaSbaaeaacqWGIbGyaeqaaaqaaiabdkgaIjabgIGiolabdMgaPbqabiabggHiLdaaaOGaeiilaWcaaa@533B@ (6)

Z i e MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemOwaO1aa0baaSqaaiabdMgaPbqaaiabdwgaLbaaaaa@2FEB@ = mean value of potential gravity.

wb = total population of spatial unit b completely within census tract i.

S = number of services in study area.

dbs = distance between spatial unit centroid b and service s.

α = friction parameter (usually 1, 1.5 or 2).

Sws = weight given to the service s such as its size (for example, number of beds for a hospital).

Types of distance

Four types of distance are typically used for calculating accessibility measures: Euclidean distance (straight-line), Manhattan distance (distance along two sides of a right-angled triangle opposed to the hypotenuse), shortest network distance and shortest network time (Figure 2) [28,31]. Euclidean and Manhattan distances can easily be computed using geographic coordinates:

thumbnailFigure 2. Several types of distance.

d i j = ( x i x j ) 2 + ( y i y j ) 2 , MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaemizaq2aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpdaGcaaqaaiabcIcaOiabdIha4naaBaaaleaacqWGPbqAaeqaaOGaeyOeI0IaemiEaG3aaSbaaSqaaiabdQgaQbqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaakiabgUcaRiabcIcaOiabdMha5naaBaaaleaacqWGPbqAaeqaaOGaeyOeI0IaemyEaK3aaSbaaSqaaiabdQgaQbqabaGccqGGPaqkdaahaaWcbeqaaiabikdaYaaaaeqaaOGaeiilaWcaaa@46F8@ (7)

dij = |xi - xj| + |yi - yj|,(8)

Where:

Xi and Yi = X and Y coordinates of point i with a plane projection.

On the other hand, calculation of network distance and network time distance – which represent respectively the shortest and fastest paths between two points using a street network – is more complex. Indeed, the computation of these two distances necessitates geometric network files – with directions, speed limits, turning restrictions, and delays available for each street segment – integrated into GIS, and a GIS module or GIS software specialized in transportation analysis (ESRI Network Analyst Extension or TransCad software, for example).

Shortest network and network time distances represent two different objectives. Shortest network distance is useful for evaluating the path between two points as if taken on foot; consequently, it is frequently used in studies on the accessibility of "proximal" services and facilities [4,18,20,32-34]. Shortest time distance is more accurate for evaluating distances by car or public transportation.

Methodological issues and accuracy in measuring geographical accessibility

When evaluating geographical accessibility, the choice of these parameters is likely to generate different results, potentially leading to significant measurement errors [2,28,29].

Most studies examining the geographical accessibility of health care and health-related services have been concerned with measuring the accessibility of the closest facility using Euclidean distance [2,13], shortest network distance [8,18,32], shortest network time distance [9,11,16,17,35], or a combination of distance types [10,12,15]. Others have also examined proximity to diversity by measuring the average shortest network distance to selected services [32], and the offer provided by the immediate surroundings, i.e. of the number of services within a given distance [7,14,18,32]. Few studies have conceptualized different dimensions of geographical accessibility within one investigation (for exceptions, see [14,15,32]), although this would be useful in order to describe the complexity of geographical accessibility of a given service [32]. Furthermore, within a given set of data, the choice of the accessibility measure is fundamental since accessibility varies depending on the indicator used [3,4,36].

Some studies have compared discrepancies in results when geographical accessibility was measured using different types of distances. In a study exploring trade-offs between various types of distance, Apparicio and colleagues [28] compared distance matrices based on simple (Euclidean, Manhattan) and network (road, time) measures of distance between all census tracts in Canada's eight largest metropolitan areas. They examined whether the Euclidean and Manhattan approximations are correlated with a more accurate measure of distance, i.e. travel time along the road network, at the metropolitan and census tract levels. The authors concluded that, at the metropolitan level, the use of Euclidean or Manhattan distances to estimate shortest network time does not introduce major errors. However, in sub-metropolitan areas, or areas located away from the central business district, the use of Euclidean or Manhattan approximations of shortest network time may lead to substantial errors. In measuring accessibility in these areas, network-based distance/time matrices may provide more appropriate results. Similar results were also observed [10,12,15].

With respect to operational definition of residential areas, a wide variety of area units has been used, ranging from smaller units such as census meshblocks [9,10,16,18], enumeration districts [11,13] and census output areas [12] to larger units such as census tracts [7,17,32,35], communities and city-defined neighbourhoods [2,8], and wards [13-15]. Some studies have controlled for the location of the population within the spatial unit by calculating the population-weighted mean centre of the spatial unit [12,14,16,17,35] or by calculating the population-weighted mean accessibility of smaller spatial units located within the spatial unit of interest [2,13,32]. Nonetheless, a considerable number of studies do not employ methods for minimizing aggregation errors, i.e. they compute accessibility for the geometric centroid of the spatial unit.

For public policy and planning, measuring geographical accessibility of urban resources and facilities is of interest as it conveys information on the presence of enabling resources [37] or opportunity structures [38] in the residential environment. Geographical accessibility measures are however prone to a variety of methodological problems, among which is error induced by using different distance types [28] and aggregation methods [2,29].

Study objectives

In this paper, we investigate differences in results when geographical accessibility of residential areas (census tracts) to selected health care services is computed using different distances types and different aggregation methods. Specific objectives are to: 1) Compare measures of distances; and 2) Estimate aggregation errors for several accessibility measures.

Data and methods

Study area and health services

This study focuses on the Montréal census metropolitan area (CMA) which has a population of about 3.4 million inhabitants. The territory of the Montréal CMA is divided into 852 census tracts, 5,829 dissemination areas and 25,767 blocks with respective average population sizes of 4,022, 588 and 133 inhabitants, as defined by Statistics Canada. A total of 642 health services grouped into eight categories were integrated into geographic information systems (ArcGis) by geocoding addresses (Figure 3). These health services were inventoried from the website of the Ministère de la santé et des services sociaux du Québec (Quebec Ministry of Health and Social Services). Of the 642 services, 65 are located in a 10 km buffer zone around the boundaries of the Montréal CMA. These 65 services were included in order to prevent underestimation of accessibility in the bordering zones of the Montréal CMA.

thumbnailFigure 3. Categories of health services for the Montréal CMA, 2006.

Comparing distance types

To explore variations in results according to distance type, we calculate the four distance types – Euclidean, Manhattan, shortest network path and shortest network time distances – between the 642 health services and the centroids of census tracts, dissemination areas and blocks. In total, more than 83 million distances are computed (Table 2), with SAS software for Euclidean and Manhattan distances, and with the Network Analyst Extension of ArcView 3.3 [39] by using CanMap Streetfiles from DMTI [40] for shortest network and shortest network time distances.

Table 2. Distances calculated between health services and spatial units

Once these four distance types are computed, correlation analyses are performed globally and locally across entire census tract, dissemination area and block matrices. First, the global analysis, which yields one value for the CMA as a whole, allows us to assess the degree of correlation between the four distance types. Then, we examine correlations between the four distances for each spatial unit centroids and the 642 health services. This local analysis stage yields one mappable value for each census tract, dissemination area and block and allows to identify spatial variation in the degree of correlation between the four distance types.

Evaluating aggregation errors when measuring geographical accessibility

The same approach, i.e. global and local analyses, was used to evaluate aggregation errors for several accessibility measures at the census tract level. The global analysis involves calculating correlations between three aggregation methods: 1) the census tract centroid; 2) the population-weighted mean of the accessibility measure for dissemination areas within census tracts; and 3) the population-weighted mean of the accessibility measure for blocks within census tracts, the most accurate method. Although accessibility was computed for each of the eight categories of health services, for purposes of conciseness, results are reported only for accessibility of general and specialized care i.e. hospitals (n = 56) for census tracts. It is worth noting, that similar patterns of correlation were observed for other health services.

Results

Correlations between the four types of distances

Global correlations

Table 3 presents results for global correlation coefficients between the four types of distances computed for the entire sample of health services (n = 642). Three observations can be made. First, at the metropolitan scale, independently of the type of distance used, results are globally similar as indicated by high correlation coefficient values (greater than 0.95). Second, in comparison with Manhattan distance, Euclidean distance is most strongly correlated with the more accurate network path and time distances. Thus, if it is impossible to compute network distance in a study focussing on geographical accessibility in the Montréal CMA, Euclidean distance seems preferable to Manhattan distance. These results are in line with those of Apparicio et al. [28] for eight Canadian metropolitan areas (Toronto, Montréal, Vancouver, Ottawa-Hull, Calgary, Edmonton, Québec and Winnipeg), and with those of Fone et al. [12] for Caerphilly county borough in Wales. Finally, as expected, correlations between both network distances are almost perfect (0.992). This means that if directions and speed limits are unknown for computing the shortest network time, the shortest network distance is a very reliable alternative.

Table 3. Global Pearson correlations between alternative types of distance

Local correlations

Although global correlations are high, they are not perfect (values differ from one). For this reason, local variations at the intra-metropolitan scale must exist and should be examined in detail. Figure 4 presents local Pearson coefficients between Euclidean distance and shortest network time, and between Euclidean and Manhattan distances for the geographical accessibility of the 642 health services computed from the centroids of census tracts, dissemination areas, and blocks.

thumbnailFigure 4. Comparing alternative types of distance between spatial units and health services using local Pearson correlations.

Results show similar spatial patterns for the three spatial scales (census tract, dissemination area and block): with increasing distance from the central business district, correlations are reduced between Euclidean distance and shortest network time, and between Euclidean and Manhattan distances. For all spatial units in the centre of the Island of Montréal and those located on the south shore, correlations are higher. For those located on the periphery of the CMA, notably on the north shore, characterized by suburban areas, correlations are weaker, often lower than 0.9.

These results illustrate that for the Island of Montréal, integrating Euclidean distances at the census tract, dissemination area and block levels into statistical analysis, e.g. in regression or multilevel analysis, would yield similar results to those obtained if network distances were computed. However, if the focus is on the CMA as a whole, or on specific parts of the CMA, namely, those located in the northern suburbs, then results are likely to vary as a function of the distance type used to compute geographical accessibility.

Aggregation errors

Global errors

The global analysis of aggregation errors is performed by means of Spearman correlations between the three methods of aggregation used to calculated 20 accessibility measures at the census tract level using the more accurate distances (network distances). Results are shown in Table 4 for hospitals only, although similar patterns of correlation were observed for other health services.

Table 4. Spearman rank correlations between measures of the accessibility of hospitals by aggregation method

Correlations between the three aggregation methods are high (>0.9) for all measures of accessibility except for the number of services within 500, 1000 and 2000 metres. For example, correlation between the least exact aggregation method (census tract centroid) and the most exact based on blocks within census tracts is 0.588 for the number of hospitals within 500 metres, 0.776 for those within 1000 metres, and 0.898 for those within 2000 metres. This means that if we want to assess service provision in a close-proximity area around a census tract, it is preferable to use an aggregation method that precisely accounts for the distribution of population within it; if not, the risk of error might be considerable.

Local errors

A second stage of comparison of aggregation methods consists in assessing the absolute difference between the geographical accessibility results obtained from the methods based on census tract and dissemination areas centroids in relation to the most accurate method based on blocks within census tracts. The descriptive statistics for local errors are reported in Table 5 for hospitals.

Table 5. Aggregation errors in measures of the accessibility of hospitals at the census tract level

Not surprisingly, the local errors are on the whole quite small, though not insignificant: for example, compared with the most accurate method, the census tract centroid method misestimates the distance to the closest hospital by an average of 365 m, and the dissemination area method by an average of 134 m. Up to the third quartile (75%), the local errors are still quite small: for 75% of census tracts, the error associated with the census tract centroid approach is less than 365 m. However, in 10% of cases, the error is greater than 948 m, and in 5% of census tracts the error is greater than 1.5 km (Table 5). Despite the high correlations, significant errors in the measurement of geographical accessibility can occur in a small number of cases.

Absolute differences between aggregation methods for the closest hospital computed using shortest network distance and shortest distance time are further mapped in Figure 5. Again, stronger absolute aggregation errors are observed in suburban areas on the south and north shores of the CMA; errors remain smaller in central areas of the Island of Montréal.

thumbnailFigure 5. Evaluating local aggregation errors for hospitals.

For the purposes of statistical studies at a general level, the least precise aggregation method – based on census tract centroids – is adequate: it enables the broad identification of areas in Montréal that have the least access to health services. However, if one wishes to reach more precise conclusions for specific neighbourhoods, major errors appear for 5% to 10% of census tracts.

Conclusion

Over the past two decades, an increasing number of health studies have integrated the geographical accessibility of services and facilities as an important dimension of the built urban environment. The development of GIS with transportation module (ESRI Network Analyst Extension, for example) has largely fostered this increase. However, the results reported in this paper show that measures of geographical accessibility of urban health services may vary according to the distance type and aggregation method selected.

Although Euclidean and Manhattan distances are strongly correlated with network distances, local variations are nonetheless observed, notably in suburban areas. The choice of the aggregation method is also important: accessibility measures computed from census tract centroids, though not inaccurate, yield important measurement errors for 5% to 10% of census tracts. This is especially so in census tracts with lower population density and in those where the land use is largely non-residential. Because the accessibility of health services may be more problematic in suburban areas than in more central urban areas, geographical accessibility studies should be based on the most accurate measures. Thus, using the smallest area unit possible included in the spatial unit of interest appears to be a relevant alternative for minimizing aggregation errors.

Results obtained for Montreal – comparison of the four types of distances and evaluation of aggregation errors – may be generalisable to other North American cities where urban forms are similar. For example, in a study aiming at comparing types of distances, Apparicio and colleagues observed similar results for the eight largest Canadian metropolitan areas [28]. Moreover, results can also be extended to other services and amenities (not only for health services) although the magnitude of aggregation errors may likely vary [2].

Computing accessibility measures in GIS using network distances and more precise aggregation methods is no longer a daunting task. Nowadays, software and modules dedicated to network analysis are effective and user-friendly, notably the Network Analyst Extension of ArcGIS or the TransCad software. Moreover, street network data are easily accessible (for example, Statistics Canada or DMTI data). Because the calculation speed of computers is exponential, the time required for the computation of numerous network distances is no longer a limitation. Consequently, although errors associated with the choice of distance types are important for about 10% of census tracts, we should not avoid using the best estimation method possible for evaluating geographical accessibility. This is especially so if accessibility measures to health services or health-related resources are to be included as a dimension of the built environment in studies investigating residential area effects on health outcomes. Imprecision of accessibility measures could lead to errors or lack of precision in the estimation of area effects on health.

Future studies should investigate the extent of aggregation errors occurring when measuring accessibility to other types of services and amenities, but also in other cities where urban form may differ from that of North American cities, and finally in rural context where geographical accessibility is an important issue.

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

PA and MA are the principal investigators of the study. They carried out the GIS, statistical and mapping analyses. MR reviewed the literature and wrote parts of the paper. All authors jointly drafted and critically revised the paper, and read and approved the final manuscript.

Acknowledgements

The authors would like to thank the anonymous reviewers for their helpful comments on the original version of this paper. The study has been funded by the Social Sciences and Humanities Research Council of Canada (SSHRC).

References

  1. Penchansky R, Thomas JW: The concept of access. Definition and relationship to consumer satisfaction.

    Medical Care 1981, 19(2):127-140. PubMed Abstract | Publisher Full Text OpenURL

  2. Hewko J, Smoyer-Tomic KE, Hodgson MJ: Measuring neighbourhood spatial accessibility to urban amenities: Does aggregation error matter?

    Environment and Planning A 2002, 34(7):1185-1206. Publisher Full Text OpenURL

  3. Talen E: Visualizing fairness: Equity maps for planners.

    Journal of the American Planning Association 1998, 64(1):22-38. Publisher Full Text OpenURL

  4. Talen E, Anselin L: Assessing spatial equity: an evaluation of measures of accessibility to public playgrounds.

    Environment and Planning A 1998, 30(4):595-613. Publisher Full Text OpenURL

  5. Joseph A, Phillips D: Accessibility and Utilization-geographical perspectives on health care delivery. New York: Harper and Row; 1984.

  6. Joseph AE, Bantock PR: Measuring potential physical accessibility to general practitioners in rural areas: a method and case study.

    Social Science and Medicine 1982, 16(1):85-90. PubMed Abstract | Publisher Full Text OpenURL

  7. Allard SW, Tolman RM, Rosen D: Proximity to service providers and service utilization among Welfare recipients: The interaction of place and race.

    Journal of Policy Analysis and Management 2003, 22(4):599-613. PubMed Abstract | Publisher Full Text OpenURL

  8. Bamford EJ, Dunne L, Taylor DS, Symon BG, Hugo GJ, Wilkinson D: Accessibility to general practitioners in rural South Australia.

    Medical Journal of Australia 1999, 171(11–12):614-616. PubMed Abstract | Publisher Full Text OpenURL

  9. Brabyn L, Beere P: Population access to hospital emergency departments and the impacts of health reform in New Zealand.

    Health Informatics Journal 2006, 12(3):227-237. PubMed Abstract | Publisher Full Text OpenURL

  10. Brabyn L, Skelly C: Modeling population access to New Zealand public hospitals.

    International Journal of Health Geographics 2002., 1 PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Christie S, Fone D: Equity of access to tertiary hospitals in Wales: A travel time analysis.

    Journal of Public Health Medicine 2003, 25(4):344-350. PubMed Abstract | Publisher Full Text OpenURL

  12. Fone D, Christie S, Lester N: Comparison of perceived and modelled geographical access to accident and emergency departments: a cross-sectional analysis from the Caerphilly Health and Social Needs Study.

    International Journal of Health Geographics 2006, 5(1):16. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Haynes R, Bentham G, Lovett A, Gale S: Effects of distances to hospital and GP surgery on hospital inpatient episodes, controlling for needs and provision.

    Social Science and Medicine 1999, 49(3):425-433. PubMed Abstract | Publisher Full Text OpenURL

  14. Haynes R, Lovett A, Sünnenberg G: Potential accessibility, travel time, and consumer choice: Geographical variations in general medical practice registrations in Eastern England.

    Environment and Planning A 2003, 35(10):1733-1750. Publisher Full Text OpenURL

  15. Martin D, Roderick P, Diamond I, Clements S, Stone N: Geographical aspects of the uptake of renal replacement therapy in England.

    International Journal of Population Geography 1998, 4(3):227-242. Publisher Full Text OpenURL

  16. Pearce J, Witten K, Bartie P: Neighbourhoods and health: A GIS approach to measuring community resource accessibility.

    Journal of Epidemiology and Community Health 2006, 60(5):389-395. PubMed Abstract | Publisher Full Text OpenURL

  17. Wang F, Luo W: Assessing spatial and nonspatial factors for healthcare access: Towards an integrated approach to defining health professional shortage areas.

    Health and Place 2005, 11(2):131-146. PubMed Abstract | Publisher Full Text OpenURL

  18. Witten K, Exeter D, Field A: The quality of urban environments: Mapping variation in access to community resources.

    Urban Studies 2003, 40(1):161-177. Publisher Full Text OpenURL

  19. Smoyer-Tomic KE, Hewko JN, Hodgson MJ: Spatial accessibility and equity of playgrounds in Edmonton, Canada.

    The Canadian Geographer 2004, 48(3):287-302. Publisher Full Text OpenURL

  20. Apparicio P, Séguin A-M: L'accessibilité aux services et aux équipements: Un enjeu d'équité pour les personnes âgées résidant en HLM à Montréal.

    Cahiers de géographie du Québec 2006, 50(139):23-44. OpenURL

  21. Cummins S, Macintyre S: The location of food stores in urban areas: A case study in Glasgow.

    British Food Journal 1999, 101(7):545-553. Publisher Full Text OpenURL

  22. Donkin AJ, Dowler EA, Stevenson SJ, Turner SA: Mapping access to food at a local level.

    British Food Journal 1999, 101(7):554-564. Publisher Full Text OpenURL

  23. Smoyer-Tomic KE, Spence JC, Amrhein C: Food deserts in the prairies? Supermarket accessibility and neighborhood need in Edmonton, Canada.

    The Professional Geographer 2006, 58(3):307-326. Publisher Full Text OpenURL

  24. Pearce J, Witten K, Hiscock R, Blakely T: Are socially disadvantaged neighbourhoods deprived of health-related community resources?

    International Journal of Epidemiology 2007, 36(2):348-355. PubMed Abstract | Publisher Full Text OpenURL

  25. Pearce J, Blakely T, Witten K, Bartie P: Neighborhood Deprivation and Access to Fast-Food Retailing. A National Study.

    American Journal of Preventive Medicine 2007, 32(5):375-382. PubMed Abstract | Publisher Full Text OpenURL

  26. Block JP, Scribner RA, Desalvo KB: Fast food, race/ethnicity, and income: A geographic analysis.

    American Journal of Preventive Medicine 2004, 27(3):211-217. PubMed Abstract | Publisher Full Text OpenURL

  27. Reidpath DD, Burns C, Garrard J, Mahoney M, Townsend M: An ecological study of the relationship between social and environmental determinants of obesity.

    Health and Place 2002, 8(2):141-145. PubMed Abstract | Publisher Full Text OpenURL

  28. Apparicio P, Shearmur R, Brochu M, Dussault G: The measure of distance in a social science policy context: Advantages and costs of using network distances in eight Canadian metropolitan areas.

    Journal of Geographic Information and Decision Analysis 2003, 7(2):105-131. OpenURL

  29. Handy SL, Niemeier DA: Measuring accessibility: An exploration of issues and alternatives.

    Environment and Planning A 1997, 29(7):1175-1194. Publisher Full Text OpenURL

  30. Openshaw S: The modifiable areal unit problem.

    CATMOG (Concepts & Techniques in Modern Geography) 1983., 38 OpenURL

  31. Fotheringham AS, Brunsdon C, M C: Quantitative geography: Perspectives on spatial data analysis. London: Sage Publications; 2000.

  32. Apparicio P, Cloutier M-S, Shearmur R: The case of Montréal's missing food deserts: Evaluation of accessibility to food supermarkets.

    International Journal of Health Geographics 2007, 6(1):4. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  33. Apparicio P, Séguin AM: Measuring the accessibility of services and facilities for residents of public housing in Montréal.

    Urban Studies 2006, 43(1):187-211. Publisher Full Text OpenURL

  34. Ottensmann JR: Evaluating equity in service delivery in library branches.

    Journal of Urban Affairs 1994, 16(2):109-123. Publisher Full Text OpenURL

  35. Luo W, Wang F: Measures of spatial accessibility to health care in a GIS environment: Synthesis and a case study in the Chicago region.

    Environment and Planning B: Planning and Design 2003, 30(6):865-884. Publisher Full Text OpenURL

  36. Talen E: School, community, and spatial equity: An empirical investigation of access to elementary schools in West Virginia.

    Annals of the Association of American Geographers 2001, 91(3):465-486. Publisher Full Text OpenURL

  37. Andersen RM: Revisiting the behavioral model and access to medical care: does it matter?

    Journal of Health and Social Behavior 1995, 36(1):1-10. Publisher Full Text OpenURL

  38. Macintyre S, Ellaway A, Cummins S: Place effects on health: How can we conceptualise, operationalise and measure them?

    Social Science and Medicine 2002, 55(1):125-139. Publisher Full Text OpenURL

  39. ESRI: ArcView network analyst: Optimum routing, closest facility and service area analysis.

    Redlands: Environmental Systems Research Institute inc 1996. OpenURL

  40. DMTI Spatial: CanMap® Streetfiles.

    Markham: DMTI Spatial inc 2005. OpenURL

  41. Talen E: Neighborhoods as service providers: A methodology for evaluating pedestrian access.

    Environment and Planning B: Planning and Design 2003, 30(2):181-200. Publisher Full Text OpenURL