Among the various healthcare facilities, community pharmacies (henceforth referred to as pharmacies) are the most easily accessible for laypeople in that individuals can enter a pharmacy without an appointment and receive professional advice immediately. Moreover, pharmacies are increasingly diversifying and already offer a number of complementary services, such as screening, disease prevention and harm-reduction interventions (Xu and Rojas-Fernandez, 2003; Benrimoj and Frommer, 2004; Stopka et al., 2013). Geographical access to healthcare services is one of the determinants with respect to their utilisation and this parameter may produce two diametrically opposite effects. Thus, people living in areas with an abundant supply of healthcare services are more prone to utilise those services (supplier-induced demand); by contrast, zones where healthcare provision is low may be subject to so-called supplier-suppressed demand (Rice and Smith, 2001). Spatial inequalities in healthcare provision occur in both developing countries (Kumar, 2004; Abolhallaje et al., 2014) and industrialised ones (Houghton, 2005; Perez and Decoupigny, 2009). Limited geographical access not only reduces the availability of pharmaceuticals but also limits recommendations regarding their use (Norris et al., 2014). Qato et al. (2014) and Amstislavski et al. (2012) have even coined the terms pharmacy deserts and medication deserts, respectively, to describe neighbourhoods where geographic access to pharmacies and the availability of commonly used drugs are poor.
Research has shown that geographical access to pharmacies depends on several factors, such as socio-economic determinants (Casey et al., 2002; Amstislavski et al., 2012; Qato et al., 2014) including levels of deprivation (Ward et al., 2014; Todd et al., 2014), community segregation (Qato et al., 2014), landscape features (Perez and Decoupigny, 2009) and degree of urbanity vis-à-vis remoteness (Lin, 2004; Todd et al., 2014; Norris et al., 2014). Among these, urban-rural inequalities are of particular relevance in Italy, since rural areas account for more than 90% of the surface area and host about 57% of the population. As in other countries, Italian rural areas are characterised by a high concentration of elderly residents and are subject to depopulation, especially in some regions (OECD, 2009). It should however be borne in mind that rural areas are not, even when located in the same region, homogeneous in terms of geomorphologic characteristics, e.g., altitude. Analogously, not all mountainous areas are rural, and many residents live in urban centres. Compared with lowland areas, however, the people living in mountainous areas are often penalised by complex physical topography, remoteness, low population density and significant population loss, all of which may negatively influence the provision of services (Schuler et al., 2004). Indeed, access to primary healthcare in the Province of Cuneo in Piedmont, Italy has been found to be lower in its mountainous areas than in the plains (Perez and Decoupigny, 2009). To guarantee more equitable access to pharmaceuticals, single nations have adopted different legislative and regulatory strategies (WHO, 2001; Vogler et al., 2014). In Italy, a new deregulatory reform became law (Italian Republic, 2012) in 2012. In terms of pharmacy licensing, the reform aimed to increase the number of pharmacies and ensure more competitive and broader pharmaceutical services by relaxing previous restrictions (Garattini et al., 2012). Indeed, it has been claimed that in Italy the regulations governing pharmacies are anti-competitive (Antitrust Authority, 1998). The new law loosened the demographic criterion and fixed the number of inhabitants per pharmacy at 3300 (instead of 4000 and 5000 per pharmacy in municipalities with more than or less than 12,500 inhabitants, respectively) also establishing that additional pharmacies might be opened if the above criterion was exceeded by 50% (1650 inhabitants). Additional pharmacies up to 5% of the existing number may be opened in crowded zones (e.g., in proximity to stations). According to this legislation, the pharmaceutical service must also be ensured in remote and under-populated areas. As a result of the new legislative degree, the number of pharmacies is destined to grow (Farma.For, 2012), and Garattini et al. (2012) have estimated that an increase of more than 20% pharmacy licenses can be expected. While the number of pharmacies per population is a well-recognised and widely used indicator (Xu and Rojas-Fernandez, 2003; Moïse and Docteur, 2007; Amstislavski et al., 2012), the number of pharmacies within a certain geographical or administrative area provides additional useful insights into the accessibility to pharmaceuticals: the fewer pharmacies there are within a zone, the farther a consumer has to travel to purchase needed medicines (Moïse and Docteur, 2007).
In both the Italian and European contexts, Liguria has clearly distinguishable features of both a demographic and geographic nature. The region’s concentration of elderly residents is among the highest in Europe (Eurostat, 2015). The fact that this hilly-mountainous region occupies a long, narrow strip of the Ligurian coast creates several logistic challenges in the face of the expected increase of the number of pharmacies, making it useful to rationalise the pharmaceutical supply in a more equal manner throughout the region. We have hypothesised some degree of geographical variation with reagrd to the density of pharmacies at the municipal level across the region; indeed it has already been shown that both health-related outcomes and healthcare supply vary spatially (Odoi and Busingye, 2014; Yao et al., 2015). The present study aimed: i) to explore the patterns of potential spatial accessibility to pharmacies in Liguria from the perspective of local characteristics, and ii) to identify disparities and their sources.
Materials and Methods
Liguria is a northwestern Italian region along the Ligurian sea with relatively high population density (approximately 290 inhabitants per km2). With the Apennines to the East and the Alps to the West, two thirds of the region’s surface is mountainous and one third hilly. Administratively, Liguria is composed of 235 municipalities (comuni) grouped into three provinces (La Spezia, Savona and Imperia) and the Metropolitan City of Genoa (previously Genoa Province). According to the geomorphologic classification (Istat, 2014), Italian municipalities are divided into totally/partially mountainous (at least 80% of the surface is above 600 m or has a difference in altitude >600 m) and non-mountainous areas; four fifths of Liguria’s municipalities belong to the former category. In accordance with the Eurostat definition (population density <100 inhabitants per km2 or proportion of employees in the agricultural sector higher than the EU average), 149 municipalities are defined as rural (Istat, 2014). Although most mountainous municipalities are rural, the two above-mentioned categories are not the same: of the 187 totally/partially mountainous municipalities, only 137 are rural. Demographically, Liguria has a higher proportion (28%) of elderly residents than any other Italian region (Istat, 2015a).
Study variables and data sources
The main geographical unit of analysis was the municipality; reference to the provincial level was occasionally made for descriptive and comparative purposes only. The outcome was pharmaceutical accessibility, as measured by two proxies, namely the number of pharmacies per 1000 inhabitants (Ph/inhab) and the number of pharmacies per km2 (Ph/km2) (Moïse and Docteur, 2007).
A set of independent variables was chosen on the basis of regional socio-economic features, previous research (Casey et al., 2002; Perez and Decoupigny, 2009; Amstislavski et al., 2012; Norris et al., 2014; Qato et al., 2014) and data availability at the municipal level. This dataset included population (1000s of inhabitants)/population density (inhabitants per km2), mean altitude (100 m units) above the mean sea level (msl), dummy variables of the status of rural according to the Eurostat classification (Istat, 2014), municipality (0=urban), coastal municipality (0=landlocked), proportion of >65-year-olds, <5-year-olds, immigrants, and a proxy of municipal income (mean declared personal income per taxpaper in 1000s of Euro).
The data on existing pharmacies were extracted (October 22, 2015) from the Italian Ministry of Health’s open dataset (Ministry of Health, 2015), which is updated on a regular basis and contains a full description of single pharmacies, including their geographic location. All types of dispensing points, including main pharmacies (farmacie ordinarie), pharmacy branches (farmacie succursali), annual dispensaries (dispensari annuali) and seasonal ones (dispensari stagionali) were eligible. Pharmacies operating in Liguria and open on January 1, 2015 were included in the study. Geographic and socio-economic data were retrieved from the Italian National Institute of Statistics (Istat, 2015a, 2015b) and Ministry of Economy and Finance (2012).
As a preliminary step, both indicators of density of pharmacies at the provincial level were calculated and described. Municipalities without their own pharmacy were identified and their frequency compared at the provincial level. This subset was additionally analysed in order to evaluate the distance from and travelling time to the nearest pharmacy – expressed as mean [standard deviation (SD)] – located in surrounding municipalities (bordering municipalities of other regions were also considered). For this purpose, Google maps (https://maps.google.com) were used, owing to their acceptable positional accuracy (Lemke et al., 2015) and a low likelihood of significant errors in travel time estimation (Delamater et al., 2012). Availability of at least one pharmacy was compared between mountainous and non-mountainous municipalities by means of Fisher’s exact test. The mean altitude of municipalities without a pharmacy and those with at least one pharmacy was compared by means of the t test and expressed as standardised mean difference (d) with 95% confidence interval (CI).
For the main study, both study outcomes were transformed by cube-root transformation in order to ensure their approximately normal distributions. Thematic choropleth mapping with single hue or bi-polar colour progressions was used to map the local data.
Global spatial autocorrelation analysis was performed using Moran’s I statistics (Moran, 1950). Negative, near zero and positive values indicate dispersed, random and cluster patterns, respectively. To compute the spatial autocorrelation analysis, k nearest-neighbour spatial-weights matrices were built. Since there is no universally accepted method to establish an optimal k (Hassanat et al., 2014), the square root of the number of observations (k=15 in this case) was used. A Monte Carlo simulation with 999 rounds was calculated in order to verify the statistical significance of Is.
To identify the presence of potential spatial clusters (with significantly positive I) or outliers (with significantly negative I), local spatial analysis was performed by means of local Moran’s I coefficients (Anselin, 1995). Choropleth maps were subsequently created to visualise and distinguish among four types of clusters/outliers, namely high-high, i.e. hotspots that signified municipalities with a high density of pharmacies surrounded by municipalities with a high density of pharmacies, while low–low (coldspots) signified municipalities with a low density of pharmacies surrounded by municipalities with a low density of pharmacies. Low-high and high-low outcomes represented municipalities with low/high density of pharmacies surrounded by municipalities with high/low density of pharmacies, respectively.
Prediction modelling of the municipal density of pharmacies involved several steps. First, univariate ordinary least-squares (OLS) regressions were computed in order to highlight global unadjusted associations between the study outcomes and independent variables. Final multiple OLS models were then constructed by minimising the corrected Akaike’s information criterion (AICc). Considering possible multicollinearity issues, if two or more nested models showed approximately equal AICc values, the most parsimonious one was preferred. Finally, given that the residual spatial autocorrelation is a well-known challenge in spatial datasets, the geographically weighted regression (GWR) approach (Fotheringham et al., 2002) was also exploited. GWR is a relatively recent modelling technique that overcomes issues related to spatial heterogeneity by allowing regression estimates to vary from location to location. In particular, the OLS model (eq. 1) can be converted to the GWR model (eq. 2) as follows:2015). An optimal bandwidth was re-calibrated each time in order to obtain a more refined model. The Monte Carlo randomisation tests with 999 permutations were used to test for significance of the spatial variability of GWR model parameter estimates. Mixed GWR modelling (Mei et al., 2004), which allows some predictors to be constant, was then applied if parameter estimates of an independent variable did not vary spatially according to the randomisation test.
In all multiple models, several interactions were assessed. In the event of a significant interaction, simple-slopes analysis was conducted and its results plotted. All models underwent the diagnostics of spatial autocorrelation of residuals by means of Moran’s I test. Since multicollinearity was likely, its presence was checked by means of global and local variance inflation factors (VIFs); VIFs <4 were considered acceptable (Miles and Shevlin, 2001). All analyses were performed in the R environment (R Core Team, 2014).
Pharmacies in Liguria
A total of 605 pharmacies were identified, corresponding to 0.38 pharmacies per 1000 inhabitants and 0.11 pharmacies per km2. At the provincial level, both indicators of density of pharmacies were unevenly distributed (Ph/Inhab: 0.36, 0.43, 0.37, 0.37 and Ph/km2: 0.07, 0.08, 0.18, 0.09 in Imperia, Savona, Genoa and La Spezia, respectively). At the municipal level, this gap was much wider (Ph/inhab: 0-11.07/1000 inhabitants; Ph/km2: 0-0.90/km2) and highly positively skewed (a skewness coefficient of 5.58 and 2.52 for Ph/Inhab and Ph/km2, respectively). Mapping of the cube-root-transformed study outcomes is shown in Figure 1.
A total of 43/235 (18.3%) municipalities did not have a pharmacy of their own; these were asymmetrically distributed among the four provinces: 7.2% (5/69), 9.4% (3/32), 14.9% (10/67) and 37.3% (25/67) in Savona, La Spezia, Genoa and Imperia, respectively. In these municipalities, the mean distance to the nearest pharmacy was 6.8 (SD 3.7) km, corresponding to 11.3 (SD 6.3) min of driving. All municipalities without a pharmacy belonged to totally/partially mountainous categories, while only 75.0% of municipalities with at least one pharmacy belonged to these categories (Fisher’s exact test: P<0.001). A more detailed analysis revealed that the mean altitude [396.2 (SD 225.6) m] of the municipalities with at least one pharmacy was significantly lower (t=6.63, P<0.001) than that [698.1 (SD 278.8) m] of municipalities without a pharmacy, and the effect size was large [d=1.12 (95% CI: 0.77-1.47)].
Spatial autocorrelation analysis
The global spatial autocorrelation analysis of the density of pharmacies revealed a modest but highly significant (P=0.001) clustered pattern: Moran’s I of Ph/Inhab and Ph/km2 were 0.082 and 0.295, respectively. The local spatial autocorrelation analysis of Ph/Inhab (Figure 2A) identified 27 municipalities with a significant clustered pattern (N=7 hot spots and N=20 cold spots), while 14 were significant outliers (N=8 hot-cold and N=6 cold-hot). Both clusters and outliers were situated in the mountainous hinterland. The hot spots were located in the Province of Savona, while the cold spots were concentrated in the Provinces of Imperia and Genoa. The local spatial autocorrelation analysis of Ph/km2 (Figure 2B) showed evident coast-hinterland and urban-rural patterns: 21 of 24 significant hot spots were coastal municipalities, including all four provincial capitals, while all significant (N=32) cold spots were municipalities located in the hinterland. Only five significant outliers (N=2 hot-cold and N=3 cold-hot) were identified.
Predicting the density of pharmacies: global models
Table 1 reports simple OLS models for establishing relationships between the density of pharmacies and the independent variables of interest. A significantly negative association emerged between Ph/inhab on the one hand and both the mean altitude and the proportion of elderly residents on the other, while a higher municipal income determined a higher Ph/inhab. However, these models explained only a small (<3%) amount of variance. All but one (% of immigrants) of the simple OLS models for predicting Ph/km2 proved highly significant (P<0.001): the association with population, coastal location, proportion of <5 year olds and income was positive, while it was negative with respect to mean altitude, rural status and percentage of elderly residents. Residuals of all univariate OLS models were significantly (P<0.01) spatially auto-correlated. Notably, the independent variable of mean altitude explained the highest proportion of variance in predicting both the indicators of density of pharmacies (Table 1).
The best-fitting multiple OLS model for predicting Ph/inhab in Liguria incorporated the independent variables of mean altitude and income level and their interaction. As shown in Table 2, higher income was associated with higher Ph/inhab, while mean altitude was no longer statistically significant (P=0.93). The altitude-by-income interaction term was also highly (P<0.001) positively significant, indicating that the impact of income level on Ph/inhab depended on municipal altitude. The model explained approximately 8% of variance (Table 2). Residuals of the model were still significantly (P=0.024) spatially auto-correlated (I value of 0.042). No signs of multicollinearity were detected (VIFs of 1.71, 1.87 and 1.53 for mean altitude, municipal income and interaction term, respectively). As shown by the simple-slopes analysis (Figure 3), the association between Ph/inhab and income was significantly (P<0.001) positive (b=0.073) among municipalities located at high (+1 SD above the mean) altitudes, while no association (P=0.37) was observed among lower-lying (-1 SD below the mean) municipalities. The selected multiple OLS model with the outcome of Ph/km2 contained more predictors and explained approximately 60% of the variance. Model coefficients were significantly positive for population, coastal location and income, while they were significantly negative for mean altitude and proportion of the elderly (Table 3). Although the altitude-by-income interaction was not significant, the analysis of simple slopes demonstrated that the relationship between Ph/km2 and income level was significantly (P=0.014) positive (b=0.019) among municipalities located at high altitudes, while it was not significant (P=0.089) among those situated at lower altitudes (Figure 4). Model diagnostics revealed no substantial multicollinearity (VIFs in a range of 1.06-2.81), although spatial autocorrelation was present (I=0.097, P=0.001).
Predicting the density of pharmacies: mixed geographically weighted regression
The best-fitting GWR model for predicting Ph/inhab contained the independent variables of mean altitude, proportion of the elderly and immigrants and income level. The model explained 24.8% of the variance, significantly outperforming the multiple OLS model (ΔadjR2: 0.08 ΔAICc: -7.22). The introduction of interaction terms worsened model fit. Residuals of the GWR model were spatially independent (I=-0.0002, P=0.41) and no concerns regarding multicollinearity were detected (maximum local VIFs of 2.69, 1.74, 1.28 and 1.79 for mean altitude, proportion of the elderly and immigrants and income level, respectively). The Monte Carlo simulation showed a significant variability of local estimates for intercept (P=0.048), mean altitude (P=0.037) and percentage of immigrants (P=0.046), while it did not for the proportion of elderly residents (P=0.16) and income level (P=0.54). The last two predictors were therefore inserted as global terms in a mixed GWR model (Table 4), which determined a slight improvement (ΔAICc: -1.1). As shown by the choropleth maps in Figure 5A, the association between Ph/inhab and mean altitude displayed a high level of spatial variability, the regression coefficients being mostly positive in the eastern part of the region (especially on the border between Genoa and Savona) and highly negative in the westernmost part, roughly corresponding to Imperia Province (especially in the inland Arroscia Valley). By contrast, the association between the relative numbers of foreign-born people and Ph/inhab was positive across almost the whole region and the estimates were relatively high in the eastern part of the Genoa (Figure 5B).
Ph/km2 was best predicted by the explanatory variables of population, mean altitude, rural status and income level, and the model showed a significant improvement over the multiple OLS model (ΔadjR2: 0.06 ΔAICc: -4.72). There were no problems with regard to model diagnostics (I=-0.023, P=0.80; maximum local VIFs in the range of 2.51-3.10).
The randomisation test showed a significantly high variation in the local coefficients for population (P<0.001), altitude (P<0.001) and rural status (P=0.006); by contrast, estimates of the model intercept and the predictor income level varied little (P=0.27 and P=0.51, respectively). The final mixed GWR model is shown in Table 5; in comparison with the GWR model containing all local variables, the mixed model showed a significant improvement (ΔAICc: -5.5). The relationship between Ph/km2 and population was highly positive (local bs of 0.10– 0.20) only in a few municipalities of the Arroscia valley, while these regression coefficients tended to zero in the overwhelming majority of the Liguria municipalities (Figure 6A). Conversely, municipal altitude was consistently negatively associated with Ph/km2 across the whole study area, with about an 11-fold variation in parameter estimates, the highest being in Imperia (Figure 6B). The status of rural municipality exerted a negative effect on Ph/km2 in some municipalities in Arroscia Valley and the eastern part of Genoa (Figure 6C).
Key findings and comparison with previous research
We show here that the distribution of pharmacies (both Ph/inhab and Ph/km2) in Liguria follows a modest but significant clustered pattern, identifying possible under-supplied areas, generally corresponding to less affluent, small mountainous municipalities. It is known that spatial disparities in accessibility to healthcare resources are almost unavoidable owing to constraints of various nature (Delamater et al., 2012), especially in areas with particular topographic features and scattered small mountain villages as typified by Liguria. However, assessing geographical accessibility undoubtedly helps to identify areas characterised by lower or higher access possibilities, which is an essential step towards mitigating such disparities (Apparicio et al., 2008). This documentation should be useful for decision-makers, healthcare planners and other stakeholders with respect to the consequences following opening, closing or relocating facilities including change of the spectrum of services offered (Delamater et al., 2012).
Our analysis shows that both socio-economic and landscape characteristics are associated with both indicators of the geographical accessibility to pharmacies. While urban-rural inequalities in pharmaceutical supply have been established in several studies carried out in different geographic settings (Lin, 2004; Norris et al., 2014; Todd et al., 2014), the present study highlights the fact that, in coastal and hilly-mountainous regions such as Liguria, location (coastal/mountainous) may be an independent predictor of healthcare supply. Indeed, differences in accessibility were observed between plains and mountainous places, and between coastal and hinterland areas. This is in line with the results of local spatial autocorrelation analysis, since most coldspots were situated in inland municipalities, mostly in the region’s highest mountain ranges: Marguareis Alps and Ligurian Apennines. Moreover, we found that the effects of economic status and altitude were not simply additive (especially in the model for predicting Ph/inhab): wealthier municipalities had a higher density of pharmacies, but this relationship depended on height, being significant only in locations at a relatively high altitude. These findings are also consistent with previous studies on healthcare services (Liu et al., 2011; Fang et al., 2014). Fang et al. (2014) have in fact established that economic determinants influence the utilisation of healthcare services to a greater extent in mountainous regions than in the plains. Analogously, Liu et al. (2011) have reported that maternal healthcare utilisation is closely related to both wealth and altitude. In our opinion, the difference in Ph/km2, observed between coastal and interior municipalities, is probably linked to migration from inland to coastal areas (Gandolfo, 2005) and to the relatively scant and constantly declining agricultural activity in Liguria’s mountainous areas (Liguria Ricerche, 2011). Another possible explanation may lie in the predominance of seaside tourism in Liguria (unlike mountain tourism in other Italian Alpine regions) and the consequent relative increase in population during the summer (Liguria Region, 2015). As a result, entrepreneurial activities may be more compromised in landlocked municipalities. Moreover, a study (Perez and Decoupigny, 2009) on access to healthcare in the southern part of the Alps (including the French Alpes-Maritime’s and the Italian provinces Cuneo and Imperia), has demonstrated coast-hinterland disparities in access to pharmacies. With regard to the tertiary sector, a gap between the hinterland and the coast has also been described by some authors (Jordan, 2000; Andriotis, 2006).
Municipalities with a higher proportion of elderly residents displayed lower values of Ph/km2, which potentially increases travelling time; it is thus not surprising that greater distance is associated with lower utilisation of healthcare facilities, including pharmacies as shown by Nemet and Bailey (2000) and Hiscock et al. (2008). Proximity to a pharmacy is of particular importance to the elderly; in one survey, Italian seniors living in both urban and rural areas ranked the availability of a pharmacy in 2nd place (after food stores) (Marcellini et al., 2007). The elderly need pharmacies more than the general population (Lin, 2004); indeed, more than 90% of elderly people in Italy take at least one medicine on a long-term basis (Pilotto et al., 2006). The elderly are also less likely to utilise health services when they feel unwell (Fang et al., 2014). It may be argued that older adults in Italy have access to cars, either as drivers or as passengers (Mollenkopf et al., 2004; Marcellini et al., 2007). However, senior citizens often give up driving (Hakamies-Blomqvist and Wahlström, 1998) or at least reduce it (Raitanen et al., 2003) for several reasons including health-related ones, and their journeys become shorter and closer to home (Whelan et al., 2006). Moreover, mountain roads may be dangerous (Fu et al., 2011), especially in the winter. In addition, more than a quarter of the elderly in Italy live alone (National Observatory on Health Status in the Italian Regions, 2009), which reduces their transport options; even those living with their families may be reluctant to ask to be driven somewhere (Adler and Rottunda, 2006). The above-mentioned factors may contribute to under-utilisation of healthcare facilities and non-adherence to medical prescriptions. In this regard, the ability to drive and the distance to pharmacies have been found to be the only determinants of non-adherence to medication among patients with epilepsy (Welty et al., 2010).
Although most mountainous inland municipalities have a relatively low population density, this was not found to be associated with Ph/inhab confirming the report by Mangano (2010) that population density is not the main predictor, at least not in Italy, of pharmacy location. Indeed, the study by Mangano (2010), which was conducted in the Province of Catania, Sicily, has established that pharmacies are often opened independently from population density. Of note, Catania Province is also characterised by a substantial between-municipality difference in altitude.
Usefulness of geographically weighted regression in modeling the density of pharmacies
The usefulness of GWR in studying issues of spatial accessibility to primary healthcare has already been demonstrated (Bagheri et al., 2009). In our hands, the GWR approach outperformed the OLS models by being capable to completely reset the residual spatial autocorrelation frequent in spatial datasets. It allowed us to highlight the fact that relationships between the density of pharmacies and some predictors varied significantly across Liguria. We found that the relationship between municipal altitude and both indicators of the density of pharmacies was highly negative in several municipalities of the westernmost province of Imperia, many of which corresponded to the coldspots revealed in the local spatial autocorrelation analysis. Similarly, the rural status of a municipality was highly and negatively associated with Ph/km2 in two distinct inland areas. Local coefficients for the percentage of immigrants showed somewhat opposite patterns, in that the eastern half of Liguria had much higher values that the western half. This could be explained by the uneven distribution of immigrants in Italy (Bratti and Conti, 2013). Indeed, foreign-born residents are less numerous in eastern provinces (Genoa and La Spezia) than in the western ones (Savona and Imperia). Interestingly, in the GWR model predicting Ph/inhab, both spatially varying explanatory variables also cut the region into halves along the longitude. It should be noted that the Ligurian coast has historically been subdivided into a eastern part (Riviera di Levante) and a western one (Riviera di Ponente) parts, which have distinct geographical features. On the other hand, while the above-described predictors varied spatially, municipal income and the proportion of elderly residents were constant across the study area, thus supporting the application of the mixed GWR approach.
Deregulation and spatial accessibility to community pharmacies
The 2012 deregulatory reform explicitly stated that, to ensure greater accessibility to pharmaceutical services and an equal distribution of pharmacies in a municipality, the Local Health Unit and the Order of Pharmacists must first be consulted in order to identify areas for new pharmacies. It was also underlined that the location of new pharmacies should guarantee accessibility to pharmaceutical services by residents of remote areas and those with low population density (Italian Republic, 2012). According to our spatial analysis, it seems that the most under-served zone in Liguria is situated in inland Imperia. This finding is in line with the analysis carried out by Farma.For (2012), which states that, as a result of deregulation, a total of 75 new pharmacies would be opened in Liguria, and that most of these (N=33 or 0.15 new pharmacies/1000 inhabitants) should be located in the Imperia Province. By contrast, these proportions are substantially lower in the other provinces (0.06, 0.02 and 0.05 new pharmacies/1,000 inhabitants in Savona, Genoa and La Spezia, respectively).
The country large study by Vogler et al. (2014) reveals that in countries having adopted deregulation policies (in primis Norway and Sweden, but also England, The Netherlands and Ireland) the number of pharmacies increased. However, the newly opened pharmacies were mostly concentrated in urban areas, which already had a high number of pharmacies in the pre-reform period, while no indications of improved accessibility in rural areas were found. In the 2.5 years following the Italian reform, only 5 new pharmacies were licensed in Liguria. Three of these were opened in small mountainous municipalities that did not previously have their own pharmacy. Paradoxically, another new pharmacy was opened in the city of Genoa, which, according to the new criterion, had previously been deemed to have an excess of 32 pharmacies (Farma.For, 2012). This pattern is consistent with the results reported by Wagner et al. (2009), who found that the deregulation reform in England had only slightly redressed the geographic imbalance in the distribution of pharmacies. Garattini et al. (2012) conclude that the Italian deregulatory reform is unlikely to benefit consumers, and that opening new small pharmacies will further aggravate financial pressure on the National Health Service in terms of reimbursements of the costs of medicines. Moreover, according to these authors, under-supplied municipalities could be treated exceptionally, for example by allowing general practitioners to dispense pharmaceuticals in more remote areas (Garattini et al., 2012). An interesting opportunity to improve access to pharmaceutical services in isolated and mountainous municipalities, especially those without their own pharmacy, could be the experimental project Farma-click launched by the Italian Ministry of Health. This project consists of installing automatic drug dispensers equipped with a camera; once the patient has described his or her symptoms directly to a pharmacist and/or scanned the doctor’s prescription, the pharmacist authorises the machine to dispense the medicines required (OECD, 2009).
It has been noted that opening a rural pharmacy in every small village that lacks its own might not be necessary, since these municipalities are usually close to one another and it would not be difficult to reach a pharmacy in the vicinity (Garattini et al., 2012). Since walking is the most common means of locomotion among the elderly (Mollenkopf et al., 2004), it is obviously advantageous if the nearest pharmacy is situated within walking distance. Our analysis reveals, however, that considering a walkable distance of either 800 m (Guerra et al., 2012) or 1600 m (Villanueva et al., 2013), only 3 centroids of 43 in Liguria municipalities had the closest pharmacy within these distances. Moreover, walkability may be further compromised by the location of all 43 municipalities in mountainous areas.
The interpretation of our results is subject to some limitations. First of all, the outcome variables used in the present study are proxies of potential accessibility and, as such, refer to hypothetical rather than documented accessibility, which measures the actual utilisation of healthcare services by residents of a given area. Potential accessibility may underestimate travelling distances in both urban and rural areas (Lin et al., 2005). Indeed, Lin et al. (2005) have estimated that the travelling distance calculated by means of revealed access is at least twice that calculated by means of potential access. This implies that many people do not access the closest pharmacy, perhaps owing to inconvenient opening hours, personal preferences or financial issues. Secondly, although it has been suggested that analysing potential accessibility to pharmacies at the municipal level and focusing on remote areas may be optimal (Moïse and Docteur, 2007), the minimum local level in Italy is, theoretically, the census tract. Indeed, dividing large municipalities with this in mind would provide useful, more detailed information. However, data at the sub-municipality level are rarely utilised in Italy and thus unrecoverable (Testi et al., 2005). Thirdly, the relatively low percentage of explained variance, especially in the models predicting Ph/inhab, suggests that other explanatory variables, e.g., density of general practitioners (Mangano, 2010) may be important in predicting the density of pharmacies. Finally, as in all observational studies, no causality can be inferred.
The results of the study presented here show that, in order to mitigate spatial access inequalities in Liguria, new drug-dispensing points should be established in mountainous inland zones, particularly in those municipalities without their own pharmacy. Considering that licensing and subsequent opening of new pharmacies in remote and sparsely populated areas may be economically disadvantageous (Garattini et al., 2012), and that such areas have relatively old populations and low rates of entrepreneurship, the implementation of alternative dispensing strategies look attractive. The density of pharmacies at the municipal level is associated with several socio-economic and landscape factors, some of which vary spatially, while others do not. In such situations, application of the mixed GWR technique is a valuable solution, and mapping its results may facilitate the work of policy-makers and healthcare planners.