Skip to main content

Efficiency of fine scale and spatial regression in modelling associations between healthcare service spatial accessibility and their utilization



Healthcare accessibility, a key public health issue, includes potential (spatial accessibility) and realized access (healthcare utilization) dimensions. Moreover, the assessment of healthcare service potential access and utilization should take into account the care provided by primary and secondary services. Previous studies on the relationship between healthcare spatial accessibility and utilization often used conventional statistical methods without addressing the scale effect and spatial processes. This study investigated the impact of spatial accessibility to primary and secondary healthcare services on length of hospital stay (LOS), and the efficiency of using a geospatial approach to model this relationship.


This study focused on the ≥ 75-year-old population of the Nord administrative region of France. Inpatient hospital spatial accessibility was computed with the E2SFCA method, and then the LOS was calculated from the French national hospital activity and patient discharge database. Ordinary least squares (OLS), spatial autoregressive (SAR), and geographically weighted regression (GWR) were used to analyse the relationship between LOS and spatial accessibility to inpatient hospital care and to three primary care service types (general practitioners, physiotherapists, and home-visiting nurses). Each model performance was assessed with measures of goodness of fit. Spatial statistical methods to reduce or eliminate spatial autocorrelation in the residuals were also explored.


GWR performed best (highest R2 and lowest Akaike information criterion). Depending on global model (OLS and SAR), LOS was negatively associated with spatial accessibility to general practitioners and physiotherapists. GWR highlighted local patterns of spatial variation in LOS estimates. The distribution of areas in which LOS was positively or negatively associated with spatial accessibility varied when considering accessibility to general practitioners and physiotherapists.


Our findings suggest that spatial regressions could be useful for analysing the relationship between healthcare spatial accessibility and utilization. In our case study, hospitalization of elderly people was shorter in areas with better accessibility to general practitioners and physiotherapists. This may be related to the presence of effective community healthcare services. GWR performed better than LOS and SAR. The identification by GWR of how these relationships vary spatially could bring important information for public healthcare policies, hospital decision-making, and healthcare resource allocation.


Accessibility to healthcare is widely recognized as a critical issue in public health and plays a fundamental role in improving health outcomes [1,2,3,4,5,6,7,8,9,10,11]. One of the most common conceptualizations of healthcare accessibility was developed by Andersen who characterized access along several dimensions, including potential and realized access [12, 13]. Potential access, also known as spatial access, refers to the spatial distribution of healthcare facilities, and the possibility to reach healthcare activity locations or services from a given position or by an individual [14,15,16,17,18]. Potential access represents only one dimension of healthcare accessibility because the presence of nearby healthcare facilities may not result in realized access [19,20,21], and individuals who can access such healthcare services may choose not to use them [22]. Realized access refers to the actual utilization of healthcare services [23] and the real interaction with the healthcare system [24]. In the last two decades, much geographic research on healthcare access has focused on the potential access to services by measuring the spatial accessibility of medical services, rather than on the realized access (i.e. the utilization of health services) [25, 26]. However, measuring both potential and realized access is mandatory in order to address the diversity of demands by providers, patients as well as policy makers. Moreover, healthcare service utilization could be strongly influenced by the potential access [27,28,29]. Yet, less research effort has been dedicated to understand how potential accessibility affects healthcare utilization [26].

The assessment of the spatial access to healthcare services and their utilization should also take into account the care provided by the various facility types: ambulatory/home care, which is the basis of primary care in many countries [30, 31], and inpatient hospital care, also called secondary care, including both acute and long-term care hospitals [32,33,34,35]. Spatial access to and utilization of these two healthcare types are closely linked. Previous research findings suggest that patients who face geographic barriers to primary care may use more frequently hospital services [36,37,38,39]. This is the case of patients living in medically underserved areas [38, 40, 41]. Moreover, the ability of primary healthcare services to take care of discharged patients has a significant impact on the hospitalization length [33, 42, 43]. Length of hospital stay (LOS) is a classical indicator of hospital care utilization, and is interpreted as a measure of the healthcare supply and treatment efficiency [33, 44,45,46]. Investigating the effect of spatial accessibility to primary and secondary care services on LOS is important to better understand the complex links between healthcare spatial accessibility and utilization, and between hospital and ambulatory/home healthcare [47,48,49]. This knowledge is crucial for efficient resource planning and allocation [50,51,52].

Finally, many previous studies on the relationship between healthcare spatial accessibility and utilization used conventional statistical methods without addressing the scale effect and spatial processes [37, 39, 53,54,55]. However, this relationship relies on spatial data and is influenced by the associated spatial effects [29, 56,57,58].

Healthcare utilization data have already been used at the country, regional, district or zip code levels to study the geographical and temporal dynamics of epidemics and their correlation with and environmental factors [59,60,61,62]. Moreover, it has been suggested that for public health planning, analyses at a finer spatial scale are more efficient to better identify and target critical areas that require intervention scaling up [63,64,65]. Indeed, such analyses could significantly increase the accuracy of the estimates about future healthcare utilization demands (e.g. number of hospital beds) [66, 67], and could capture the heterogeneity of spatial accessibility [68]. However, few studies have addressed the scale effect by comparing quantitatively the efficiency of fine and large spatial scale analyses, especially when analysing the relationship between healthcare spatial accessibility and utilization.

Moreover, spatial data exhibit spatial non-stationarity and spatial autocorrelation [69]. Spatial autocorrelation can be defined as the self-similarity of nearby observations [70]. Therefore, autocorrelated residuals may increase the level of uncertainty of the regression coefficients, and usually lead to larger prediction intervals [71]. Spatial autocorrelation has been described for geographical and spatial economic phenomena and may be readily dealt using spatial regression models [72, 73]. The residual spatial autocorrelation issue is commonly addressed using the spatial autoregressive (SAR) model that takes into account spatial dependency. Non-stationarity is another spatial issue. Most studies suggested generalized results across all locations [29, 37, 39, 42, 55, 74], and ignored how the relationship between healthcare spatial accessibility and utilization may take different forms in different places. Spatial non-stationarity occurs when the strength and direction of the relation in one place does not apply in another. Compared with global models in which results are assumed to be stationary and could be generalized to all locations, the outcomes of local models are location-specific. A generalized conclusion does not fully take into account specific local situations and may mislead public policy-makers [75,76,77,78]. Geographically weighted regression (GWR) allows such analysis and can help to identify relationships that remain hidden when using global models, such as ordinary least squares (OLS) and SAR. Using both SAR and GWR allows considering the global and local spatial scales, and can help to identify all factors, their impact on dependent variables, and their interdependence [79]. Moreover, to the best of our knowledge, no study has compared the efficiency of a traditional regression (OLS) with a geospatial approach (SAR and GWR) for modelling the relationship of potential and realized healthcare access.

Therefore, the main objective of this work was to assess the efficiency of fine spatial scale analyses for investigating the impact of spatial accessibility to primary and secondary healthcare services on the length of hospital stay, relative to geospatial variables. This study had two specific objectives: i) to assess the efficiency of a fine spatial scale analysis to identify hidden factors in the relationship between LOS and healthcare spatial accessibility; and ii) to explain the spatial variability of this association by comparing traditional (OLS) and geospatial modelling approaches (SAR and GWR), thus providing a case study on the application of geospatial techniques for the analysis of this relationship.

Material and methods

Study setting and population

This study was carried out in the Nord administrative region that is located in the north of France, with a surface area of 5743 km2 and a population density of 456 inhabitants per km2. Several primary care spatial accessibility indicators are available for this department [80,81,82,83]. Moreover, it has been shown that edge effects lead to minor spatial accessibility variations in this area [84]. The study focused on the ≥ 75-year-old group of the population of this region, because many different healthcare resources, such as hospital facilities and primary care professionals, are involved in their management. Moreover, as their recovery period after a hospital stay is often longer, their length of stay could be more influenced by the capacity of the primary healthcare services to manage them [47,48,49].

Statistical unit

As the first objective was to assess the efficiency of a fine spatial scale analysis, this study was carried out first at the French Geographic Code unit (FGC) level and then at the sub-municipal French census block group level (known as IRIS: “Ilots Regroupés pour l’Information Statistique”). The FGC statistical unit is defined by the French national discharge database, and is approximately equivalent to the municipality postal code. The Nord administrative region is divided into 240 FGCs with a population ranging from 1,000 to 227,000 inhabitants per FGC. IRIS is defined by the French National Institute of Statistics and Economic Studies (INSEE) [85], and is the smallest infra-urban level for which complete French census data are available. There are 1,346 IRIS units in the Nord administrative region with a population ranging from 6 to 5,414 inhabitants per IRIS. A FGC can include 1 to 110 IRIS units, but one IRIS only belongs to one FGC.

Data sources

Multiple data sources were combined for the present study.

Index of spatial accessibility (ISA)

ISA is a previously developed indicator to measure spatial accessibility to hospital care (inpatients). Based on the enhanced two-step floating catchment area (E2SFCA) method, ISA takes into account the number of beds in Medical, Surgical and Obstetrics facilities (Médecine, Chirurgie, Obstétrique: MCO) and Postoperative and Rehabilitation Care facilities (Soins de Suite et de Réadaptation: SSR) facilities, the car travel time, and the population distribution. ISA construction steps have been previously described [83, 86]. It is a special form of physician-to-population ratio, expressed as the number (N) of beds per 10,000 inhabitants. Higher scores indicate higher accessibility. This indicator was initially developed at the census block (IRIS) scale, and was then summarized at the French Geographic Code (FGC) scale for MCO and SSR facilities. On average, there are 22.69 beds in MCO and 5.49 beds in SSR for 10,000 inhabitants in the Nord administrative region. ISA spatial distribution revealed important variations within IRIS units (Fig. 1). Specifically, the highest ISA values for MCO were observed in urban areas located in the northern part of the studied territory and also in the centre. Conversely, the lowest values were observed mostly in the southern part. The highest ISA scores for SSR were concentrated in the middle part of the region, whereas access was lower in the North and South.

Fig. 1

Index of Spatial Accessibility distribution at the IRIS level. a Index of Spatial Accessibility (ISA) for Medical, Surgical and Obstetrics (MCO) and b Postoperative and Rehabilitation Care (SSR) centres. For each map, the Nord administrative region is represented using a graduated colour scale to highlight the ISA score variability among IRIS units

Length of stay

The Length Of Stay (LOS) was defined as the mean hospital stay length of elderly people (≥ 75 years of age) relative to the total ≥ 75-year-old population in a given FCG (Eq. 1).

$${LOS}_{i}=\sum_{g>75}\frac{{Average\, length\, of\, stay}_{gi}}{{P}_{gi}}$$

where gi represents the three groups of ≥ 75-year-old people (75–84, 85–94 and > 95 years) for a given spatial unit i, and Pgi the corresponding total population for that age group. The numerator represents the average length of stay of each age group, standardized to the total population for that age group (inpatients or not) in the denominator. The mean LOS was chosen instead of the median LOS in order not to underweight extreme LOS values. Indeed, in this context of healthcare resource allocation, extremely long LOS do exist and need to be taken into account.

The LOS was calculated from the French national discharge database that collects hospital activity and patient discharge data [87,88,89,90] using MCO and SSR inpatient care data for the year 2014. The number of ≥ 75-year-old people in the Nord administrative region was obtained from the 2016 French national Census [85].

Localized potential accessibility (APL)

The spatial accessibility to primary (ambulatory/home) care services was described using the Localized Potential Accessibility (Accessibilité potentielle localisée: APL) database [80]. Based on the E2SFCA method, APL indices are available at the national level but only for eight types of healthcare professionals: general practitioners (GPs), physiotherapists, home-visiting nurses, paediatricians, dental surgeons, midwives, gynaecologists, and ophthalmologists. For this study, the APL indices for GPs, physiotherapists, and home-visiting nurses were used because these three types of health professionals could contribute substantially to the primary care services to older adults. Furthermore, the healthcare provided by them might interact with inpatient hospital care. APL indices are available at the FGC level and are expressed for 10,000 inhabitants.

Socioeconomic variables

Although this was not the principal aim of the study, several neighbouring socioeconomic variables were introduced because they also could affect healthcare utilization behaviours, as previously reported [36, 91,92,93,94]. From the initial database of economic variables [85] that might influence care utilization, data pre-processing was performed using the variance inflation factor (VIF) to identify the presence of multicollinearity among covariates. VIF values lower than 10 are considered acceptable [95, 96]. The final variables included in the model for further analysis and the scale at which they are available are summarized in Table 1.

Table 1 Variable description


Ordinary least squares model

Multiple linear regression models are frequently used for predictive and explanatory analyses and the Ordinary Least Squares (OLS) method allows estimating the best fit. An OLS model treats data as independent, assuming that any variation due to spatial variability within the area units is not captured by the method.

In our study, the OLS model was expressed by the following function:

$$ln\left({LOS}^{*}\right)={\beta }_{0}+ {\beta }_{1}.ln\left({ISA}^{*}\right)+{\beta }_{2}.precarious +{\beta }_{3}.nonOwner+\sum {\beta }_{n}.\mathrm{ln}\left({APL}_{n}\right)+ \varepsilon ;with \varepsilon \approx iid\left(0;{\sigma }^{2}\right)$$


  • signifies that the LOS for MCO (with the corresponding ISA as independent variable) and for SSR (with the corresponding ISA as independent variable) were modelled separately.

  • n defines the number of different types of primary healthcare professionals considered in the analysis (n = 3: GPs, physiotherapists, and home-visiting nurses).

  • The other covariates, namely Precarious, Non-Owner and APL, are defined in Table 1.

Spatial non-stationarity and spatial autocorrelation: towards spatial statistical methods

Traditional regression models, such as OLS, consider only the stationary relationship between dependent and independent variables, ignoring the spatial autocorrelation and non-stationarity [97]. However, spatial effects, such as spatial heterogeneity and spatial autocorrelation, can affect the accuracy of OLS estimates by increasing the regression errors and uncertainty [98]. Therefore, two alternative spatial modelling methods were assessed to overcome these limitatons.

Spatial AutoRegressive model
$$ln\left({{LOS}^{*}}_{i}\right)=\rho \sum {w}_{ij}ln\left({{LOS}^{*}}_{j}\right)+{\beta }_{0}+ {\beta }_{1}.ln\left({{ISA}^{*}}_{i}\right)+{\beta }_{2}.{precarious }_{i}+{\beta }_{3}.{nonOwner}_{i}+\sum {\beta }_{n}.\mathrm{ln}{\left({APL}_{n}\right)}_{i}+ {\varepsilon }_{i};with {\varepsilon }_{i}\approx iid\left(0;{\sigma }^{2}\right).$$

The Spatial AutoRegressive (SAR) model takes into account spatial dependency and addresses the spatial autocorrelation problem by adding a spatial lag term, \(\rho WY\), in Eq. 2 that defined the OLS model. The SAR model is based on the assumption that a dependent variable at a location is affected by the dependent variable of neighbouring locations in addition to the effects of independent variables [99]. The values of the dependent variable of neighbouring geographic units are averaged and become a term on the independent side of the equation. The model equation is as follows: where:

  • \(\rho \) is the spatial autoregressive coefficient. It quantifies the effect of neighbour observations of \(ln\left({{LOS}^{*}}_{j}\right)\) and the direction of that effect [100].

  • \(ln\left({{LOS}^{*}}_{i}\right)\) is partly explained by the values taken by \(ln\left({LOS}^{*}\right)\) in the neighbouring geographic units.

  • \(\sum {w}_{ij}\) is the spatial lag for the n x n weight matrix \(W\), and represents the proximity between each pair of geographic units \(\left(i,j\right)\).

  • The spatialreg R package was used for this model [101]. To construct the spatial weight matrix, both rook and queen contiguity weights were analysed and the rook contiguity, which defines adjacent polygons as those sharing edges, was selected due to its slightly higher R2.

Geographically weighted regression model

Global models, such as OLS and SAR, may mask potential spatial non-stationarity. To explore LOS spatial non-stationarity, the GWR model was implemented. This model considers that coefficients vary across space, by estimating different relationships between the dependent and independent variables for each geographic location. The purpose of GWR is to embed the geographical location into parameters based on a traditional regression, in order to establish a spatial-weighting matrix and run a local weighted regression for each area unit to allow the analysis of the spatial variation and related driving factors of the research object at a specific scale. Regression coefficients are then defined as local coefficients. Therefore, the GWR model is more suitable for studying the spatial heterogeneity and local effects, and truly depicts the local influence of independent variables on dependent variables. The model is described by the following equation:

$$ln({{LOS}^{*}}_{i})={\beta }_{0}\left({u}_{i},{v}_{i}\right)+ {\beta }_{1}\left({u}_{i},{v}_{i}\right).ln\left({ISA}^{*}\right)+{\beta }_{2}\left({u}_{i},{v}_{i}\right).precarious +{\beta }_{3}\left({u}_{i},{v}_{i}\right).nonOwner+\sum {\beta }_{n}\left({u}_{i},{v}_{i}\right).\mathrm{ln}\left({APL}_{n}\right)+ \varepsilon ;with \varepsilon \approx iid\left(0;{\sigma }^{2}\right)$$


  • \({(u}_{i},{v}_{i})\) denotes the coordinates of each data point where local regressions are calculated.

  • At location \(i,\) \(ln({{LOS}^{*}}_{i})\) is the local dependent variable; \({\beta }_{0}\) and \({\beta }_{k}\) represent the local estimate intercept and coefficient of factor \(k,\) respectively [96, 102, 103].

At each regression area unit, local coefficients are estimated according to a spatial weighting scheme, characterized by two elements: (1) the neighbouring window of the area unit, and (2) the distance decay function used to calculate the spatial weights for each area unit in the window.

The spatial weighting function employed in a GWR model assumes that neighbouring area units in a window have more similar characteristics compared with those distant from each other. Thus, the area unit parameters \({\beta }_{k}\left({u}_{i},{v}_{i}\right)\) are more strongly influenced by closer observations than by units further away. The neighbouring window of each area unit can be determined using fixed or adaptive kernel types. The fixed type selects an optimal global bandwidth for the whole area; all area units that fall within the bandwidth are included in the regression. The adaptive type adjusts the size of the spatial window by choosing a specified number of nearest neighbours. The kernel function modifies the weights given to each neighbouring unit according to its distance from the regression unit. In this study, the adaptive kernel approach was used because the geographical unit size in the study area varied greatly, and therefore the fixed bandwidth approach appeared irrelevant. The Gaussian function was selected to minimize the Akaike information criterion (AIC). The spatial weight matrix was constructed with the rook contiguity method and the GWR models were implemented using the gwr R package [104].

Model assessment and comparison

Global Moran’s I

The Moran test was first employed to assess the global autocorrelation of the dependent variable LOS. The Global Moran’s I index was estimated to quantify the level of spatial autocorrelation (or clustering) of LOS. Moran’s I statistic values range from − 1 to + 1, depending on the degree and direction of the association. Values significantly close to 1 indicate positive spatial autocorrelation, and values significantly close to − 1 indicate negative spatial autocorrelation. A number approaching 0 indicates the absence or a small spatial autocorrelation level, which means that data have a random spatial relationship.

Monte-Carlo test for the GWR model

It is important to note that the GWR model does not assume that relationships vary across space, but can identify whether they do or not. If the relationships do not vary across space, the global model could be appropriate for the data. Therefore, the Monte-Carlo test is used to evaluate the spatial variability of individual parameters or coefficients in a GWR model [102, 105,106,107]. The main idea of the Monte-Carlo test is that in the absence of relevant spatial phenomena, the geographical coordinates of the observations could be permuted randomly for a certain number of times. Consequently, the variance remains unchanged and the p-values can be estimated.

Measures of goodness of fit

Two goodness of fit measures were implemented to compare the OLS, SAR and GWR approaches: the adjusted coefficient of determination (adjusted R2) and the AIC [108, 109]. R2, one of the most popular measures, indicates the variance proportion of a dependent variable that is explained by an independent variable in a regression model. The model performance is proportional to the R2 value. AIC estimates the robustness of each model, relative to each of the other models. The best model is the one with the lowest AIC index value.

Analysis of residuals: Moran’s I and Local Indicators of Spatial Association

A posteriori, the global Moran’s I and local indicators of spatial association (LISA) were used to determine whether the three models (OLS, SAR and GWR) could efficiently eliminate spatial autocorrelation in the residuals of their estimates.

The LISA indicator, also known as Local Moran’s I, was developed by Anselin in 2005 [110]. While the global Moran’s I statistic shows whether overall, spatial autocorrelation is present or not, LISA identifies the presence of spatial autocorrelation clusters at specific locations. Four types of local spatial associations between an observation point and its neighbours can be detected: High-High (HH), Low-Low (LL), High-Low (HL) and Low–High (LH). A first-order rook contiguity relationship was chosen as the spatial weight for the LISA clustering analysis. The Geoda 0.9.5 software was used for the LISA (Local Moran’s I) analyses [111].

Statistical strategy

The main objective was to assess the efficiency of fine spatial scale analyses.

First, the OLS regression model was used to examine the global linear relationship between LOS and spatial accessibility to inpatient hospital care and to the three types of primary care services (GPs, physiotherapists, and home-visiting nurses). This part of the study was carried out first at the FGC level and then at the IRIS level. The OLS model was considered as the reference model in this study.

Then, the efficiency of a geospatial modelling approach to enhance our understanding of the relationship between LOS and healthcare spatial accessibility was tested. After assessment of the spatial autocorrelation with the Global Moran’s I test, the SAR and GWR models were applied and compared with the results of the OLS model previously carried out. The performance of each model (OLS, SAR, and GWR) was assessed using measures of goodness of fit.

Finally, the spatial autocorrelation and spatial non-stationarity of the residuals of these models were assessed, and the efficacy of the spatial statistical methods to reduce or eliminate spatial autocorrelation in the residuals was also explored [112].

In summary, by using the OLS, SAR and GWR models, the LOS model can be expressed as follows:

OLS: Length of Hospital Stays = F (inpatient hospital spatial accessibility, primary care service spatial accessibility, Precarious, Non-Owner)

SAR: Length of Hospital Stays = F (Spatially Adjacent Neighbour’s length of hospital stays, inpatient hospital spatial accessibility, primary care service spatial accessibility, Precarious, Non-Owner)

GWR: Length of Hospital Stays = f (inpatient hospital spatial accessibility, primary care service spatial accessibility, Precarious, Non-Owner),

where f is the regression equation for each observation in a dataset, influenced or weighted to a greater degree by the variables of other observations nearer to it.


LOS basic statistics

Analysis of the LOS for the ≥ 75-year-old population in the Nord administrative region (Table 2) showed that in 2014, the mean LOS values were 0.26 and 0.85 for MCO and SSR, respectively. LOS was initially constructed at the FGC level, and then disaggregated at the IRIS level. The spatial variation was important, with standard deviations of 0.20 and 0.92 for MCO and SSR, respectively. The LOS_MCO and LOS_SSR score distribution highlighted a non-homogeneous repartition with higher values close to the border with other French regions, especially in the southern part (Fig. 2). Lower scores were observed from Dunkerque to Bailleul, and also around Lille, Roubaix, Orchies and Tourcoing.

Table 2 LOS of ≥ 75-year-old people in MCO and SSR facilities at the IRIS scale – Nord administrative region
Fig. 2

Length Of Stay spatial distribution for ≥ 75-year-old people at the IRIS level. a Length Of Stay in Medical, Surgical and Obstetrics and b Post-operative and Rehabilitation Care centres

Efficiency of the fine scale analysis

OLS regression analysis was used to compare the efficiency of the spatial analysis at the FGC and IRIS scales. For analyses at the FGC level, independent variables available at the IRIS level (ISA, Non-Owner and Precarious) were summarized by taking into account the population density. Conversely, for analyses at the IRIS level, variables available only at the FGC level (APL and LOS) were disaggregated by giving the same value to each IRIS included in the same FGC.

The OLS regression analysis (regression coefficient estimates and the corresponding p-values are in Table 3) showed that for LOS at MCO facilities, three variables were significant at the FGC level (p-value < 0.01): spatial accessibility to GPs (APL_GPs), percentage of inhabitants who do not own their main property (Non-Owner), and percentage of inhabitants with precarious employment (Precarious). APL_GPs and Non-Owner values were negatively correlated with LOS, suggesting that a better spatial accessibility to GPs might decrease LOS, and that hospital stays at MCO facilities tend to be shorter for elderly people living in disadvantaged area units where more people do not own their main property.

Table 3 OLS regression analysis of length of stay (LOS) at MCO and SSR facilities

When the OLS model was run at the IRIS scale, all independent variables showed the same direction of association (positive or negative) with LOS as observed at the FCG level, but for the Precarious variable that was no longer significant at an α-risk level of 0.05. The number of variables showing significant relationships with LOS increased from three to five. Specifically, spatial accessibility to home-visiting nurses and physiotherapists, and to MCO facilities became significant. LOS was negatively correlated with APL_Physiotherapists and APL_GPs. However, the regression results at the IRIS level suggested that better spatial accessibility to home-visiting nurses corresponded to longer LOS. Furthermore, LOS was shorter for patients with easier access to a MCO facility (ISA_ MCO).

Similar results were obtained for the analysis of LOS at SSR facilities, with only two and five significant variables at the FGC and IRIS level, respectively. Among the significant variables at the IRIS level, only ISA_SSR displayed an opposite sign compared with the analysis for MCO facilities, showing longer LOS for patients with easier access to SSR facilities.

The adjusted R2 for the OLS regression analyses increased from 0.389 and 0.352 at the FGC level to 0.508 and 0.485 at the IRIS level for MCO and SSR facilities, respectively. The standard errors decreased accordingly (from 44 to 77% for all coefficients).

In the next steps of this study, the OLS model at the IRIS scale was considered as the reference model and was compared with the other two spatial models.

Efficiency of the geospatial modelling approach

SAR and GWR analyses were carried out to address the issues of spatial autocorrelation and non-stationarity, respectively.

Coefficient estimates

Table 4 summarizes the results obtained with the three different models. In the OLS regression model, LOS in MCO was negatively correlated with spatial accessibility to MCO facilities (ISA_MCO), spatial accessibility to GPs and physiotherapists (APL_GPs and APL_Physiotherapists), and percentage of inhabitants who do not own their main property (Non-Owner). It was positively correlated with spatial accessibility to home-visiting nurses (APL_Nurses).

Table 4 Comparison of the OLS, SAR and GWR models for LOS in MCO and SSR facilities

Compared with the results obtained with OLS, in the SAR model that controls the global spatial effect, the number of variables significantly associated with LOS decreased from five to three (APL_GPs, APL_Physiotherapists and Non-Owner). These three independent variables remained negatively associated with LOS. The absolute value of the coefficients and standard errors of each variable decreased. Furthermore, the magnitude of the regression coefficients suggested that spatial accessibility to physiotherapists had a greater global effect on LOS than spatial accessibility to GPs. Other factors, such as spatial accessibility to MCO facilities and to home-visiting nurses and percentage of inhabitants with a precarious employment (Precarious), did not significantly affect LOS in the SAR model.

By incorporating the mean LOS in neighbouring IRIS units into the equation, SAR took into account an estimate of the potential effect of neighbouring IRIS on LOS. For example, the spatial lag coefficient \(\uprho \) for MCO was 0.767 (p < 0.001) (Table 4). This means that for each 1 unit increase in LOS at a given IRIS, their neighbours’ LOS increased by 0.767 units. However, this spatial effect was assumed to be the same across all IRIS units in the study area.

GWR was developed to highlight local effects. It adds additional information by estimating the regression coefficients for each area unit. Unlike OLS and SAR that give a global estimate, GWR produces local regression coefficients to highlight non-stationarity. The regression coefficients of GWR (Table 4) varied in the study area. The sign of the coefficients for different IRIS changed for all independent variables. This indicated that the effects of independent variables on LOS were not constant throughout the studied region and differed among IRIS units. For example, the local regression coefficients for APL_GPs varied from − 37.882 to 37.084, showing that the influence of spatial accessibility to GPs on LOS had various magnitudes and even different directions.

Spatially varying parameters given by GWR

From the GWR model, local \(\beta s\) were estimated for each independent variable. Figures 3 and 4 shows the spatial variations of the regression coefficients \(\beta \).

Fig. 3

Significant GWR coefficient estimates \((\alpha \)-risk level of 0.05), local R2 and bandwidth (MCO facilities) (part 1). Coefficients for a ISA to MCO, b APL to GPs, c APL to physiotherapists, d APL to home-visiting nurses

Fig. 4

Significant GWR coefficient estimates \((\alpha \)-risk level of 0.05), local R2 and bandwidth (MCO facilities) (part 2). Coefficients for a Non-Owners, and b Precarious, c Local R2. d Bandwidth

The results of the GWR model showed the spatial heterogeneity of the associations with LOS. For instance, the spatial variation in the association between LOS in MCO and spatial accessibility to GPs was apparent (Fig. 3b), with 38.11% and 23.11% of IRIS units having a significant negative (mainly in the south or around Dunkerque and Bailleul) and positive association (mainly in urban areas in the centre, around Lille, Roubaix and Valenciennes), respectively. The associations between LOS in MCO and spatial accessibility to physiotherapists (Fig. 3c) displayed a completely different spatial distribution. For 55.35% of IRIS units, this relationship was negative, as expected (i.e. when spatial accessibility to physiotherapists increased, LOS number decreased); however, for 13.97% of IRIS units, particularly in the urban areas of the centre and in the southern part, the relationship was positive. These findings indicate that these two indexes of spatial accessibility to primary care services had completely different local influence on LOS. The local adjusted R2 values varied across the studied region, and were higher in the urban areas in the centre (Fig. 4c). The adaptive bandwidth used in the GWR model to take into account the hugely variable size of the geographical units was smaller for small IRIS (Fig. 4d).

Coloured areas show areas where the indicated parameter is significantly associated with length of hospital stay (LOS). Red, positive correlation; blue, negative correlation.

Only \(\beta \) coefficient values that were significant at an α-risk level of 0.05 are presented. To be continued in Fig. 4.

Coloured areas show areas where the indicated parameter is significantly associated with length of hospital stay (LOS). Red, positive correlation; blue, negative correlation.

Only \(\beta \) coefficient values that were significant at an α-risk level of 0.05 are presented.

Spatial model assessment & comparison

AIC and adjusted R2

The performance of each model (OLS, SAR, and GWR) used in this study was assessed by measuring the coefficient of determination (adjusted R2) and AIC. The very high R2 for the GWR model (0.955 for MCO and 0.931 for SSR) indicated that in this model, the independent variables explained the largest part of the variance of the dependent variable compared with the other two. R2 was 0.775 and 0.508 for MCO, and 0.769 and 0.485 for SSR, with the SAR and OLS models, respectively. The AIC value for the GWR model was 539 for MCO, lower than the AIC values for the SAR (2238) and OLS (3285) models.

Monte-Carlo test for GWR

The Monte-Carlo test was used to verify the significant spatial variability of the GWR model individual coefficients. The p-values < 0.05 of the independent variables, except for Non-Owner and Precarious, confirmed the significant spatial variation in the local coefficient estimates for all four indexes of spatial accessibility to primary and secondary care services. This result confirmed that most of the independent variables varied significantly across space (for both MCO and SSR), and stressed the importance of using models that take into account the spatial non-stationarity in our study.

Residual analysis using Moran’s I and local indicators of spatial association

The LISA indicator maps showed that in the OLS, SAR, and GWR models (for MCO), spatial clustering of LOS residuals occurred in statistically significant patterns (Fig. 5). In the OLS model, residuals exhibited spatial autocorrelation in many areas across the studied region (Fig. 5a). Autocorrelation clusters were smaller with the SAR model (Fig. 5b). The spatial autocorrelation of residuals was reduced and residuals were spatially clustered only in 204 IRIS units (68HH, 74LL, 34LH, and 28HL). The number of spatially associated clusters was lowest in the GWR model (Fig. 5c): only 126 IRIS units (47HH, 48LL 32LH, and 28HL).

Fig. 5

LISA map of OLS, SAR and GWR Model Residuals (MCO analysis). The red areas indicate areas with high LOS residuals surrounded by other areas with high LOS residuals (HH). The blue areas indicate areas of low LOS residuals surrounded by other areas with low LOS residuals (LL). The light blue areas indicate areas with low LOS residuals surrounded by areas with high LOS residuals (LH). The pink areas indicate areas with high LOS residuals surrounded by areas of low LOS residuals (HL). White areas, no significant spatial clustering. p < 0.01

Moreover, the Moran’s I values for the residuals varied similarly in the OLS, SAR and GWR models (Table 4). Moran’s I value significantly decreased between the OLS and GWR models (0.609 to 0.086). The lowest Moran’s I (0.031) was obtained with the SAR model. This value was very close to 0, suggesting an absence of spatial autocorrelation in LOS residuals.


To our knowledge, this is the first study that investigated the efficiency of fine scale and spatial regression in modelling associations between healthcare service spatial accessibility and length of stay. Different spatial scales and regression models were examined.

Comparison with the international literature

Previous studies explored the association between healthcare spatial accessibility and utilization by different perspectives. In most cases, researchers focused on the spatial accessibility and utilization of only one type of healthcare service (primary or secondary). For instance, Jones et al. found that due to longer travel time, patients living further away from the provider were less likely to make doctor’s appointments compared with those living closer [113]. Similarly, Arcury et al. [107] demonstrated that greater distance resulted in fewer regular check-up visits. Kim et al. found that a reduction of the travel time to hospital increases their utilization by patients with chronic obstructive pulmonary disease [74]. Another study did not observe any significant association between spatial accessibility by car or public transport and utilization of gynaecologists and GPs [29]. In these studies, the interaction between ambulatory/home care and inpatient hospital care was not considered. Other authors studied the relationship between primary care spatial accessibility and inpatient hospital care utilization. Bindman et al. [37] and Daly et al. [41] demonstrated significant associations between preventable hospitalizations and access to primary care. Several groups examined the role of primary care spatial accessibility in inpatient care/emergency department utilization. Specifically, Bindman et al. showed that access to primary care was inversely associated with the hospitalization rates for five chronic medical conditions. Daly et al., Fishman et al., and Huang et al. found that the odds of preventable emergency department use are higher in patients living in medically underserved areas. Kjekshus et al. analysed the interaction between primary and inpatient hospital healthcare services in Norway and the effect on LOS. They highlighted that LOS was influenced by the primary healthcare provider capacity and also by the percentage of elderly in the hospital catchment area, type of patients, coordination procedure, and intrinsic features of the hospital [47]. In many of these studies, hospital care utilization was measured using the admission rate or LOS. The evaluation of spatial accessibility was based on the self-rated access level [37], network distance/travel time [29, 74, 113, 114], or the E2SFCA method that takes into account both healthcare density and proximity [39, 41, 55]. The statistical methods used to investigate these associations were based on classical regression models [29, 37, 39, 47, 55, 74], or global spatial lag models [41].

Interpreting the study results at the global scale

The main strengths of our study are the concomitant assessment of 1) primary care spatial accessibility, 2) inpatient hospital spatial accessibility, and 3) hospital care utilization. Moreover, the OLS, SAR and GWR models allowed obtaining results both at the global and local scale. At the global level, the analysis revealed a significant and negative association between LOS of elderly people (≥ 75 years) and spatial accessibility to GPs and physiotherapists. In other words, better spatial accessibility to these two primary care services corresponded to shorter hospital stays. This finding could be explained by the hypothesis that inpatient hospital care facilities and primary care services may interact in a complementary way. Specifically, in areas with better spatial accessibility to primary care, hospital stays could be shorter thank to the presence of effective primary care services (e.g. outpatient care and neighbourhood healthcare services). This result appears coherent with the findings of previous studies [115,116,117]. As mentioned by Kjekshus et al., in most countries, after hospital treatment, a patient should be followed by primary healthcare services, and the ability of the primary healthcare services to take care of such patients is believed to have a significant impact on LOS, particularly for elderly people in whom the recovery period is often longer [47].

Local scale results: Efficiency of the geospatial modelling approach

As no previous study compared the efficiency of traditional regression analyses with a geospatial approach for modelling the relationship of potential and realized healthcare access, this study compared the results obtained with the OLS, SAR, and GWR models. The main weakness of the OLS approach is that it ignores spatial non-stationarity and spatial autocorrelation. SAR, which addressed the spatial autocorrelation issue, gave more robust results than OLS. However, one limitation of the SAR model is that all observations share a common spatial effect ρ. Yet, spatial variation in geographic data is seldom constant across a study area. This issue was addressed by the GWR model that can highlight regional patterns of spatially varying parameter estimates, thus revealing additional insights at local scales. This allowed determining whether the relationship between potential and realized access was intrinsically different across space. Whatever the model, our study found a negative association between LOS and access to GPs and to physiotherapists at the global scale. However, the GWR model revealed that for some areas (i.e. local analysis), this relationship was positive, and that the distribution of areas positively or negatively correlated with LOS varied when considering the access to GPs and to physiotherapists. Only the GWR model allowed highlighting this level of information.

Statistically speaking, SAR and GWR improved the model performance. However, GWR ranked first (highest adjusted R2 and lowest AIC). After applying spatial models, the spatial autocorrelation of residuals was significantly reduced in the SAR model. This reduction was lower in the GWR model, as indicated by its Moran’s I values.

In terms of public health, the ability to identify how these relationships vary in space could bring important information to engage discussions about healthcare public policy, hospital decision-making and healthcare resource allocation: 1) for areas where primary care spatial accessibility is negatively correlated with LOS, it would be important to make sure that there is already a coordinated approach between primary and secondary care services, and to allocate more resources to community care (especially specific healthcare professional types in function of the population needs); 2) for areas with high spatial accessibility to primary healthcare services and long LOS (i.e. positive correlation), it should be determined whether the primary services can properly follow patients after hospital discharge. Moreover, a consolidated approach should be developed to facilitate the care pathway coordination. The objective is to contribute to the sustainability of inpatient and outpatient care services, to complement inpatient hospital care with primary care, and to increase healthcare efficiency.

Finally, although spatial analyses provide a tool for exploring the impact of potential accessibility to primary and secondary healthcare services on LOS, our preliminary quantitative results should be supplemented by a qualitative approach to better understand them. Future investigative fieldwork and analyses of the spatial accessibility effect on LOS should include different dimensions related to care pathway coordination, such as patient and health professional types and their behaviours, hospital characteristics, procedure performed, financial issues and service quality. These dimensions should provide a more comprehensive analysis of healthcare service access and utilization.

Contributions and limitations

Our study investigated the impact of spatial accessibility to hospital and primary healthcare services and socio-economic factors on LOS. The use of multiple spatial scales and regression models allowed assessing the spatial heterogeneity and spatial dependence of different factors. Moreover, it assessed the efficiency of fine scale and spatial regression model analysis. However, some limitations need to be acknowledged. First, we chose the LOS to describe healthcare utilization. Besides spatial accessibility, this indicator could be influenced by other factors, such as disease severity, disease type, treatments reimbursed and chosen by the patients. These factors were not included in the present case study on the relationship between LOS and spatial accessibility. Then, in our potential and realized access framework, we included hospital and primary healthcare service spatial accessibility, and a hospital service utilization indicator. We did not have information about primary care service utilization. The socio-economic variables involved in this study were limited in number, and therefore their analysis was preliminary. As more data become available, more influential factors could be added to the spatial models to obtain a better explanation of the relationships between potential and realized healthcare access. In addition, as our models are based on E2SFCA metrics and data aggregated at the FGC or IRIS level, there may be an ecological fallacy and data may not reflect the associations at the individual level within the territory. Furthermore, some data were only available at the FGC level and were disaggregated to the IRIS level in a homogeneous way. For future studies, we want to construct LOS and APL indicators at a finer scale using more sophisticated disaggregation techniques.


Examining the association of LOS with primary and secondary healthcare service spatial accessibility provides information that may be exploited for public health planning. The comparison of the OLS, SAR and GWR models showed that spatial regressions are useful for these analysis. GWR performed better and could provide additional insights by revealing the hidden spatial distribution patterns of coefficient estimates and their statistical significance.

The finding of “spatially varying relationships” between healthcare spatial accessibility and LOS is important because it can help to better understand the complicated links between healthcare spatial accessibility and utilization, and between primary and secondary healthcare. It should be noted that the GWR method only revealed the spatially varying patterns, and additional research is needed to obtain more insights into the causal effects of primary and secondary healthcare service spatial accessibility and length of stay.

Availability of data and materials

All data generated or analyzed during this study are included in this published article. If readers need supplementary information, they can contact me (



Localized Potential Accessibility


Enhanced two-Step Floating Catchment Area


National File on Health and Social Institutions


French Geographic Code unit


Spatial Accessibility Index


Medical, Surgical and Obstetrics


Annual Statistical Survey of Healthcare Facilities


Postoperative and Rehabilitation Care


Length Of Stay


  1. 1.

    Commission on the Social Determinants of Health (CSDH). Closing the gap in a generation: Health equity through action on the social determinants of health. World Health Organization. CSDH 2008

  2. 2.

    Paterson I, Judge K. Equality of Access to Healthcare. London, UK: Routledge; 2002.

    Google Scholar 

  3. 3.

    Korda RJ, Butler JR, Clements MS, Kunitz SJ. Differential impacts of health care in Australia: trend analysis of socioeconomic inequalities in avoidable mortality. Int J Epidemiol. 2007;36(1):157–65.

    Article  Google Scholar 

  4. 4.

    Boerma T, Eozenou P, Evans D, Evans T, Kieny MP, Wagstaff A. Monitoring progress towards universal health coverage at country and global levels. PLoS Med. 2014; 11(9): e1001731.

  5. 5.

    Horrill T, McMillan DE, Schultz ASH, Thompson G. Understanding access to healthcare among Indigenous peoples: a comparative analysis of biomedical and postcolonial perspectives. Nurs Inq. 2018;25(3):e12237.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Shengelia B, Murray CJL, Adams OB. In: Health Systems Performance Assessment. Debates, methods and empiricism. Murray CJL, Evans DB, editor. Geneva: World Health Organization; 2003. Beyond access and utilization: defining and measuring health system coverage; pp. 221–234.

  7. 7.

    Levesque JF, Harris MF, Russell G. Patient-centred access to health care: conceptualising access at the interface of health systems and populations. Int J Equity Health. 2013;12:18.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Condelius A, Andersson M. Exploring access to care among older people in the last phase of life using the behavioural model of health services use: a qualitative study from the perspective of the next of kin of older persons who had died in a nursing home. BMC Geriatr. 2015;15:138.

    Article  Google Scholar 

  9. 9.

    Cookson R, Asaria M, Ali S, Ferguson B, Fleetcroft R, Goddard M, et al. Health Equity Indicators for the English NHS: a longitudinal whole-population study at the small-area level. Southampton: NIHR Journals Library: University of Southampton Science Park, Southampton SO16 7NS, UK; 2016. Health Services and Delivery Research.

  10. 10.

    Mbanya VN, Terragni L, Gele AA, Diaz E, Kumar BN. Access to Norwegian healthcare system - challenges for sub-Saharan African immigrants. Int J Equity Health. 2019;18(1):125.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Zhang T, Liu J, Liu C. Changes in perceived accessibility to healthcare from the elderly between 2005 and 2014 in China: an oaxaca-blinder decomposition analysis. Int J Environ Res Public Health. 2019;16(20):3824.

    Article  PubMed Central  Google Scholar 

  12. 12.

    Andersen RM. Revisiting the behavioral model andaccess to medical care: does it matter? J Health SocBehav. 1995;36(1):1–10.

    CAS  Google Scholar 

  13. 13.

    Andersen RM, Davidson P. Improving access to care in America: Individual and contextual indicators. Changing the US Health Care System: Key Issues in Health Services Policy and Management. 2014.

  14. 14.

    Hansen WG. How Accessibility Shapes Land Use. J Am Inst Plan. 1959;25:73–6.

    Article  Google Scholar 

  15. 15.

    Kwan MP. Space-time and integral measures of individual accessibility: a comparative analysis using a point-based framework. Geogr Anal. 1998;30:191–216.

    Article  Google Scholar 

  16. 16.

    Luo W, Wang F. Measures of spatial accessibility to health care in a GIS environment: synthesis and a case study in the Chicago region. Environ Plan B Plan Des. 2003;30:865–84.

    Article  Google Scholar 

  17. 17.

    Kwan MP, Weber J. Scale and accessibility: Implications for the analysis of land use-travel interaction. Appl Geogr. 2008;28:110–23.

    Article  Google Scholar 

  18. 18.

    Pan X, Kwan MP, Yang L, Zhou S, Zuo Z, Wan B. Evaluating the Accessibility of Healthcare Facilities Using an Integrated Catchment Area Approach. Int J Environ Res Public Health. 2018;15(9):2051.

    Article  PubMed Central  Google Scholar 

  19. 19.

    Donabedian A. Aspects of Medical Care Administration: Specifying Requirements for Health Care. Harvard University Press; 1973.

  20. 20.

    Aday LA, Andersen RM. Equity of access to medical care: a conceptual and empirical overview. Med Care. 1981;19(12 Suppl):4–27.

    CAS  Article  Google Scholar 

  21. 21.

    Strand de Oliveira, J. Access to care. Journal of the American Academy of Physician Assistants, 26(11), 40–45. 2013.

  22. 22.

    Andersen R, Aday LA. Access to medical care in the U.S.: realized and potential. Med Care. 1978;16(7):533–46.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Borders TF, Booth BM. Research on rural residence and access to drug abuse services: where are we and where do we go? J Rural Health. 2007;23(s1):79–83.

    Article  PubMed  Google Scholar 

  24. 24.

    Mahmood A, Mosalpuria K, Wyant DK, Bhuyan SS. Association between Having a Regular Health Provider and Access to Services Linked to Electronic Health Records. Hospital Topics. 2018.

    Article  PubMed  Google Scholar 

  25. 25.

    Higgs G. The role of GIS for health utilization studies: literature review. Health Serv Outcomes Res Methodol. 2009;9:84–99.

    Article  Google Scholar 

  26. 26.

    Delamater PL, Shortridge AM, Kilcoyne RC. Using floating catchment area (FCA) metrics to predict health care utilization patterns. BMC Health Serv Res. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Guagliardo MF. Spatial accessibility of primary care: concepts, methods and challenges. Int J Health Geogr. 2004;3(1):3.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Focke A. Regionale Leistungs und Krankenhausplanung. EinSimulationsmodell auf Basis eines Ameisenalgorithmus. Wiesbaden: Deutscher Universitäts-Verlag/GWV Fachverlage; 2006.

    Book  Google Scholar 

  29. 29.

    Stentzel U, Bahr J, Fredrich D, Piegsa J, Hoffmann W, van Berg N. Is there an association between spatial accessibility of outpatient care and utilization? Analysis of gynecological and general care. BMC Health Serv Res. 2018;18:1.

    Article  Google Scholar 

  30. 30.

    Karpiel MS. Using patient classification systems to identify ambulatory care costs. CBS Business Network. Retrieved 25 July 2011.

  31. 31.

    Berman P. Organization of ambulatory care provision: a critical determinant of health system performance in developing countries. Bull World Health Organ. 2000;78(6):791–802.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32. Accessed 22 May 2020.

  33. 33.

    Gandré C, Gervaix J, Thillard J, Macé JM, Roelandt JL, Chevreul K. The development of psychiatric services providing an alternative to full-time hospitalization is associated with shorter length of stay in french public psychiatry. Int J Environ Res Public Health. 2017;14(3):325.

    Article  PubMed Central  Google Scholar 

  34. 34.

    Bourgueil Y, Marek A. Three Models of Primary Care Organisation in Europe, Canada, Australia and New-Zealand. 2009. Questions d’économie de la santé. N°141.

  35. 35.

    Health Care System. The Free Medical Dictionary. Retrieved December 21, 2020.

  36. 36.

    Billings J, Zeitel L, Lukomnik J, Carey TS, Blank AE, Newman L. Impact of socioeconomic status on Hospital Use in New York City. Health Aff. 1993;12(1):162–73.

    CAS  Article  Google Scholar 

  37. 37.

    Bindman AB, Grumbach K, Osmond D, et al. Preventable hospitalizations and access to health care. JAMA. 1995;274(4):305–11.

    CAS  Article  Google Scholar 

  38. 38.

    Mathison DJ, Chamberlain JM, Cowan NM, Engstrom RN, Fu LY, Shoo A, Teach SJ. Primary care spatial density and nonurgent emergencydepartment utilization: a new methodology for evaluating access to care. Acad Pediatr. 2013;13(3):278–85.

    Article  PubMed  Google Scholar 

  39. 39.

    Fishman J, McLafferty S, Galanter W. Does spatial access to primary care affect emergency department utilization for nonemergent conditions? Health Serv Res. 2018;53(1):489–508.

    Article  PubMed  Google Scholar 

  40. 40.

    Bauder AR, Gross CP, Killelea BK, Butler PD, Kovach SJ, Fox JP. The relationship between geographic access to plastic surgeons and breast reconstruction rates among women undergoing mastectomy for cancer. Ann Plast Surg. 2017;78(3):324–9.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Daly MR, Mellor JM, Millones M. Do avoidable hospitalization rates among older adults differ by geographic access to primary care physicians? Health Serv Res. 2017;53:3245–64.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Gautun H, Syse A. Earlier hospital discharge: a challenge for Norwegian municipalities. Nordic J Soc Res. 2017.

    Article  Google Scholar 

  43. 43.

    Barra M, Dahl FA, Faiz KW, Lurås H. Associations between length of stay and municipal characteristics for a large Norwegian stroke unit. Res Square. 2020.

  44. 44.

    Horwood CR, Moffatt-Bruce SD, Rayo MF. Continuous cardiac monitoring policy implementation: three-year sustained decrease of hospital resource utilization. Adv Health Care Manag. 2019;24:18.

    Article  Google Scholar 

  45. 45.

    Sly KA, Lewin TJ, Frost BG, Tirupati S, Turrell M, Conrad AM. Care pathways, engagement and outcomes associated with a recovery-oriented intermediate stay mental health program. Psychiatry Res. 2020;20(286):112889.

    Article  Google Scholar 

  46. 46.

    Shiyovich A, Gilutz H, Arbelle JE, Greenberg D, Plakht Y. Temporal trends in healthcare resource utilization and costs following acute myocardial infarction. Isr J Health Policy Res. 2020;9(1):6.

    Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Kjekshus LE. Primary health care and hospital interactions: Effects for hospital length of stay. Scand J Public Health. 2005;33(2):114–22.

    Article  PubMed  Google Scholar 

  48. 48.

    Deraas TS, Berntsen GR, Hasvold T, et al. Does long-term care use within primary health care reduce hospital use among older people in Norway? A national five-year population-based observational study. BMC Health Serv Res. 2011;11:287.

    Article  Google Scholar 

  49. 49.

    Levin KA, Crighton EM. The Reshaping Care for Older People programme and changes in unscheduled hospital care: Analysis of routinely collected hospital data. Maturitas. 2017.

    Article  PubMed  Google Scholar 

  50. 50.

    Bullen N, Moon G, Jones K. Defining localities for health planning: a GIS approach. Soc Sci Med. 1996;42(6):801–16. (PMID: 8778994).

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Shortt NK, Moore A, Coombes M, Wymer C. Defining regions for locality health care planning: a multidimensional approach. Soc Sci Med. 2005;60(12):2715–27. (Epub 2005 Jan 8 PMID: 15820582).

    Article  PubMed  Google Scholar 

  52. 52.

    Alegana VA, Wright JA, Pentrina U, Noor AM, Snow RW, Atkinson PM. Spatial modelling of healthcare utilisation for treatment of fever in Namibia. Int J Health Geogr. 2012;11:6.

    Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Fortney JC, Steffick DE, Burgess JF Jr, Maciejewski ML, Petersen LA. Are primary care services a substitute or complement for specialty and inpatient services? Health Serv Res. 2005;40(5 Pt 1):1422–42.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Zhou D, Lin Z, Lim SH. Spatial characteristics and risk factor identification for land use spatial conflicts in a rapid urbanization region in China. Environ Monit Assess. 2019;191(11):677. (PMID: 31654141).

    Article  PubMed  Google Scholar 

  55. 55.

    Huang Y, P Meyer, L Jin. Neighborhood socioeconomic characteristics, healthcare spatial access, and emergency department visits for ambulatory care sensitive conditions for elderly. Preventive Medicine Reports, Vol. 12, 2018, pp. 101–105.

  56. 56.

    Konerding U, Bowen T, Elkhuizen SG, et al. The impact of accessibility and service quality on the frequency of patient visits to the primary diabetes care provider: results from a cross-sectional survey performed in six European countries. BMC Health Serv Res. 2020;20:800.

    Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Cui B, Boisjoly G, Wasfi R, et al. Spatial access by public transport and likelihood of healthcare consultations at hospitals. Transp Res Rec. 2020;2674(12):188–98.

    Article  Google Scholar 

  58. 58.

    Mohamad Fadli K, Muhamad Hanafiah J, Rosliza AM. Spatial analysis in primary health care utilization: a systematic review. Malays J Med Health Sci. 2019;15(203):89–96.

    Google Scholar 

  59. 59.

    Philippon S, Broutin H, Constantin de Magny G, Toure K, Diakite CH, et al. Meningococcal meningitis in Mali: a long-term study of persistence andspread. Int J Infect Dis. 2009;13:103–9.

    Article  Google Scholar 

  60. 60.

    Khairat S, Haithcoat T, Liu S, Zaman T, Edson B, Gianforcaro R, Shyu CR. Advancing health equity and access using telemedicine: a geospatial assessment. J Am Med Inform Assoc. 2019;26(8–9):796–805.;PMCID:PMC6696489.

    Article  PubMed  PubMed Central  Google Scholar 

  61. 61.

    Werner RM, Coe NB, Qi M, Konetzka RT. Patient outcomes after hospital discharge to home with home health care vs to a skilled nursing facility. JAMA Intern Med. 2019;179(5):617–23.;PMCID:PMC6503560.

    Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Lieberman-Cribbin W, Tuminello S, Flores RM, Taioli E. Disparities in COVID-19 Testing and Positivity in New York City. Am J Prev Med. 2020;59(3):326–32.

    Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Curtis AJ, Lee WA. Spatial patterns of diabetes related health problems for vulnerable populations in Los Angeles. Int J Health Geogr. 2010;9:43.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Paireau J, Girond F, Collard JM, Maïnassara HB, Jusot JF. Analysing spatio-temporal clustering of meningococcal meningitis outbreaks in Niger reveals opportunities for improved disease control. PLoS Negl Trop Dis. 2012;6(3):e1577.

    Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Kabaghe AN, Chipeta MG, McCann RS, et al. Access and adequate utilization of malaria control interventions in rural Malawi: a descriptive quantitative study. Malar J. 2018;17(1):104.

    Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Miller P. Medical center uses desktop mapping to cut costs and improve efficiency. Geogr Inform Syst. 1994;4:40–1.

    Google Scholar 

  67. 67.

    Jia P, Shi X, Xierali IM. Teaming up census and patient data to delineate fine-scale hospital service areas and identify geographic disparities in hospital accessibility. Environ Monit Assess. 2019;191(Suppl 2):303.

    Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Joyce KM, Burke RC, Veldman TJ, Beeson MM, Simon EL. Use of fine-scale geospatial units and population data to evaluate access to emergency care. West J Emerg Med. 2018;19(6):1043–8.

    Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Anselin L. Spatial econometrics: Methods and models. Dordrecht. The Netherlands: Kluwer; 1988.

    Book  Google Scholar 

  70. 70.

    Getis A. Reflections on spatial autocorrelation. Reg Sci Urban Econ. 2007;37:491–6.

    Article  Google Scholar 

  71. 71.

    Anselin L. Under the hood: Issues in the specification and interpretation of spatial regression models. Agric Econ. 2002;27(3):247–67.

    Article  Google Scholar 

  72. 72.

    Ge Y, Zhao X. Regional economy and development: a viewpoint and application of spatial statistics. Spatial Statistics. 2017;21:337–8.

    Article  Google Scholar 

  73. 73.

    Havard S, Deguen S, Zmirou-Navier D, Schillinger C, Bard D. Traffic-related air pollution and socioeconomic status: a spatial autocorrelation study to assess environmental equity on a small-area scale. Epidemiology. 2009;20(2):223–30.

    Article  PubMed  Google Scholar 

  74. 74.

    Kim M, Ren J, Tillis W, Asche CV, Kim IK, Kirkness CS. Explaining the link between access-to-care factors and health care resource utilization among individuals with COPD. Int J Chron Obstruct Pulmon Dis. 2016;11:357–67.

    CAS  Article  Google Scholar 

  75. 75.

    Jargowsky PA. Poverty and Place. New York: Russell Sage Foundation; 1997.

    Google Scholar 

  76. 76.

    Cooke TJ. Geographic context and concentrated urban poverty within the United States. Urban Geogr. 1999;20(6):552–66.

    Article  Google Scholar 

  77. 77.

    Benson T, Chamberlin J, Rhinehart I. An investigation of the spatial determinants of the local prevalence of poverty in rural Malawi. Food Policy. 2005;30(5–6):532–50.

    Article  Google Scholar 

  78. 78.

    Curtis KJ, Voss PR, Long DD. Spatial variation in poverty-generating processes: child poverty in the United States. Soc Sci Res. 2012;41(1):146–59.

    Article  PubMed  Google Scholar 

  79. 79.

    Guo Y, Chang SS, Sha F, Yip PSF. Poverty concentration in an affluent city: Geographic variation and correlates of neighborhood poverty rates in Hong Kong. PLoS One. 2018;13(2):e0190566.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Localized Potential Accessibility (APL) database:,490,530. Accessed 22 Mar 2019.

  81. 81.

    Lucas-Gabrielli V, Mangeney C. Comment améliorer les mesures d’accessibilité aux soins pour mieux cibler les zones sous-denses ? [How can accessibility measures be improved to better target underserved areas?]. Rev Epidemiol Sante Publique. 2019;67 Suppl 1:S25-S32. French.

  82. 82.

    Lucas-Gabrielli V, Chevillard G. Déserts médicaux » et accessibilité aux soins: de quoi parle-t-on ? ["Medical deserts" and accessibility to care: what are we talking about? Med Sci (Paris). 2018;34(6–7):599–603.

    Article  Google Scholar 

  83. 83.

    Gao F, Kihal W, Le Meur N, Souris M, Deguen S. Assessment of the spatial accessibility to health professionals at French census block level. Int J Equity Health. 2016;15(1):125.

    Article  PubMed  PubMed Central  Google Scholar 

  84. 84.

    Gao F, Kihal W, Le Meur N, Souris M, Deguen S. Does the edge effect impact on the measure of spatial accessibility to healthcare providers? Int J Health Geogr. 2017;16:46.

    Article  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Institut national de la statistique et des études économiques: Accessed 2 may 2018.

  86. 86.

    Gao F, Jaffrelot M, Deguen S. Potential inpatient hospital care accessibility scores: association with healthcare utilization and non-hospital care service accessibility. BMC Health Serv Res. Revisions being processed.

  87. 87.

    Fourquet F, Demont F, Lecuyer AI, Rogers MA, Bloc DH. French medical hospital information system and cross infection surveillance: theory and feasibility. Médecine Mal Infect. 2003;33(2):110–3.

    Article  Google Scholar 

  88. 88.

    Lamarsalle L, Hunt B, Schauf M, Szwarcensztein K, Valentine WJ. Evaluating the clinical and economic burden of healthcar e-associate d infections during hospitaliza tion for surgery in France. Epidemiol Infect. 2013;141(12):2473–82.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  89. 89.

    Exposto F, Petrică N, Davies E, Beaudet A. Identification of a pulmonary arterial hypertension (PAH) patient cohort and study of its burden of illness in Programme de Médicalisation des Systèmes d’information (PMSI). Int J Cardiol. 2020.

    Article  PubMed  Google Scholar 

  90. 90.

    Charles-Nelson A, Lazzati A, Katsahian S. Obes Surg. Analysis of Trajectories of Care After Bariatric Surgery Using Data Mining Method and Health Administrative Information Systems. 2020.

    Article  Google Scholar 

  91. 91.

    Lambrew JM, DeFriese GH, Carey TS, Ricketts TC, Biddle AK. The effects of having a regular doctor on access to primary care. Med Care. 1996;34(2):138–51. (PMID: 8632688).

    CAS  Article  PubMed  Google Scholar 

  92. 92.

    Salit SA, Kuhn EM, Hartz AJ, Vu JM, Mosso AL. Hospitalization costs associated with homelessness in New York City. N Engl J Med. 1998;338(24):1734–40.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Vahedi S, Rezapour A, Mohammadbeigi A, Khosravi A. Economic inequality in outpatient healthcare utilization: the case of Iran. J Res Health Sci. 2018;18(3):e00424.

    PubMed  PubMed Central  Google Scholar 

  94. 94.

    Gupta N, Crouse DL. Social disparities in the risk of potentially avoidable hospitalization for diabetes mellitus: an analysis with linked census and hospital data. Can Stud Popul. 2019;46:145–59.

    Article  Google Scholar 

  95. 95.

    Belsley DA, Kuh E, Welsch RE. Regression Diagnostics: Identifying Influential Data and Sources of Collinearity. New York: Wiley; 2005.

    Google Scholar 

  96. 96.

    Rybarczyk G, Taylor D, Brines S, Wetzel R. A geospatial analysis of access to ethnic food retailers in two michigan cities: investigating the importance of outlet type within active travel neighborhoods. Int J Environ Res Public Health. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  97. 97.

    Mohammadinia A, Alimohammadi A, Saeidian B. Efficiency of geographically weighted regression in modeling human leptospirosis based on environmental factors in Gilan Province Iran. Geosciences (Switzerland). 2017;7(4):136.

    Article  Google Scholar 

  98. 98.

    Bertazzon S, Johnson M, Eccles K, Kaplan GG. Accounting for spatial effects in land use regression for urban air pollution modeling. Spat Spatiotemporal Epidemiol. 2015;14–15:9–21.

    Article  PubMed  Google Scholar 

  99. 99.

    Lichstein JW, Simons TR, Shriner SA, Franzreb KE. Spatial autocorrelation and autoregressive models in ecology. Ecol Monogr. 2002;72:445–63.

    Article  Google Scholar 

  100. 100.

    Duzgun HS, Kemec S. Spatial and geographically weighted regression. Encyclopedia of GIS. 2008;1:1073–7.

    Article  Google Scholar 

  101. 101.

    Bivand R, Piras G. Comparing implementations of estimation methods for spatial econometrics. J Stat Softw. 2015;63(18):1–36.

    Article  Google Scholar 

  102. 102.

    Fotheringham S, Brundson C, Charlton M. Geographically Weighted Regression: The Analysis of Spatially Varying Relationships; Wiley: New York, NY, USA, 2002.

  103. 103.

    Matthews SA, Yang TC. Mapping the results of local statistics: using geographically weighted regression. Demogr Res. 2012;26:151–66.

    Article  Google Scholar 

  104. 104.

    Bivand R, Yu D, Nakaya T, Garcia-Lopez, M-A. Package ‘spgwr’. 2015. Accessed 2 Apr 2020.

  105. 105.

    Hope ACA. A simplified Monte Carlo significance test procedure. J Royal Stat Soc. 1968;30(3):582–98.

    Google Scholar 

  106. 106.

    Brunsdon C, Fotheringham AS, Charlton M. Geographically weighted regression modeling spatial nonstationarity. Statistician. 1998;47:431–43.

    Google Scholar 

  107. 107.

    Brunsdon C, Fotheringham AS, Charlton M. Geographically weighted summary statistics - a framework for localized exploratory data analysis. Comput Environ Urban Syst. 2002;26:501–24.

    Article  Google Scholar 

  108. 108.

    Lin CH, Wen TH. Using geographically weighted regression (GWR) to explore spatial varying relationships of immature mosquitoes and human densities with the incidence of dengue. Int J Environ Res Public Health. 2011.

    Article  PubMed  PubMed Central  Google Scholar 

  109. 109.

    de Assis I, Arcoverde MAM, Ramos ACV, Alves LS, Berra TZ, et al. Social determinants, their relationship with leprosy risk and temporal trends in a tri-border region in Latin America Social determinants, their relationship with leprosy risk and temporal trends in a tri-border region in Latin America. PLOS Negl Trop Dis. 2018;12(4):e0006407.

    Article  Google Scholar 

  110. 110.

    Anselin L. Exploring Spatial Data With GeoDa. Spatial Analysis Laboratory: University of Illinois and the Center for Spatially Integrated Social Science; 2005.

    Google Scholar 

  111. 111.

    Anselin L, Ibnu S, Youngihn K. GeoDa: an introduction to spatial data analysis. Geogr Anal. 2006;38(1):5–22.

    Article  Google Scholar 

  112. 112.

    Slagle M. A comparison of spatial statistical methods in a school finance policy context. J Educ Finance. 2010;35:199–216.

    Article  Google Scholar 

  113. 113.

    Jones AP, Haynes R, Sauerzapf V, Crawford SM, Zhao H, Forman D. Travel times to health care and survival from cancers in northern England. Eur J Cancer. 2008;44(2):269–74.

    CAS  Article  Google Scholar 

  114. 114.

    Arcury TA, Gesler WM, Preisser JS, Sherman J, Spencer J, Perin J. The effects of geography and spatial behavior on health care utilization among the residents of a rural region. Health Serv Res. 2005;40(1):135–55.

    Article  Google Scholar 

  115. 115.

    Barra M, Dahl FA, Faiz KW, Lurås H. Associations between length of stay and municipal characteristics for a large Norwegian stroke unit. Health Econ Outcomes Res. 2020.

    Article  Google Scholar 

  116. 116.

    Gautun H, Bratt C, Billings J. Nurses’ experiences of transitions of older patients from hospitals to community care A nation-wide survey in Norway’. Health Sci Rep. 2020;3(3):e174.

    Article  PubMed  PubMed Central  Google Scholar 

  117. 117.

    Mathias B, Fredrik AD, Kashif WF, Hilde L. Associations between length of stay and municipal characteristics for a large Norwegian stroke unit. Health Econ Outcomes Res. 2020.

    Article  Google Scholar 

Download references


This research is supported by EHESP Rennes, REPERES—EA 7449, IPLESP and ENSAI. Points of view or opinions in this article are those of the authors and do not necessarily represent the official position or policies of the EHESP Rennes, REPERES—EA 7449, IPLESP and ENSAI.


Not applicable.

Author information




FG designed the project, performed the statistical analysis and drafted the manuscript. SD supervised the overall project, oversaw the statistical analysis, and helped to draft and revised the manuscript. CL, KK and MG performed the statistical analysis with FG and SD. BB gave important suggestions to this study. All authors interpreted the data and reviewed the manuscript for important intellectual content. All authors have read and approved the final version of the manuscript.

Corresponding author

Correspondence to Fei Gao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gao, F., languille, C., karzazi, K. et al. Efficiency of fine scale and spatial regression in modelling associations between healthcare service spatial accessibility and their utilization. Int J Health Geogr 20, 22 (2021).

Download citation


  • Spatial accessibility
  • Utilization of healthcare
  • Hospital care
  • Primary care
  • Length of stay
  • Administrative data utilization
  • OLS
  • GWR
  • SAR