 Research
 Open Access
 Published:
Horizontal equity in access to public GP services by socioeconomic group: potential bias due to a compartmentalised approach
International Journal for Equity in Health volume 18, Article number: 187 (2019)
Abstract
Background
Horizontal equity in access to public general practitioner (GP) services by socioeconomic group has been addressed econometrically by testing the statement “equal probability of using public GP services for equal health care needs, regardless of socioeconomic status”. Based on survey data, the conventional approach has been to estimate binomial econometric models in which when the respondent reports having visited a public GP, it counts as 1, otherwise it counts as 0. This is what we call a compartmentalised approach. Those respondents who did not visit a public GP but visited instead another doctor (specialist or private GP) would count as 0 (despite having used instead other modes of health care), thus conclusions of the compartmentalised approach might be biased. In such cases, a multinomial econometric model that we called comprehensive approach would be more appropriate to analyse horizontal equity in access to public GP services. The objective of this paper is to test for this potential bias by comparing a compartmentalised and a comprehensive approach, when analysing horizontal equity in access to public GP.
Methods
Using data from the 2016/17 Spanish National Health Survey, we estimate the probability of visiting a public GP as determined by socioeconomic status, health care need and demographic characteristics. We use binomial and multinomial logit and probit models in order to highlight the potential differences in the conclusions regarding socioeconomic inequities in access to public GP services. Socioeconomic status is proxied by education level, social class and employment situation.
Results
Our results show that conclusions are sensitive to the approach selected. Particularly, the horizontal inequity favouring individuals with lower education that resulted from the compartmentalised approach disappears under a comprehensive approach and only a social class effect remains.
Conclusion
An analysis of horizontal equity in access to a particular health care service (like public GP services) undertaken following a compartmentalised approach should be compared with a comprehensive approach in order to test that there is no bias as a consequence of considering as zeros the utilisation of other types of health care.
Introduction
There is ample literature in health economics that analyses horizontal inequities in access to public sector general practitioner (GP) services by socioeconomic groups following an econometric approach. In this paper we refer to those studies that using survey data test the principle of “equal probability of visiting a GP for equal health care needs” by means of econometric models. Generally, these studies estimate the probability of visiting a GP in a particular period, as a function of health care need and other socioeconomic and demographic characteristics of the individual. Socioeconomic status is generally approached by education, social class or/and income.
Most of these studies show evidence of a gradient in the probability of using GP services favouring those in lower socioeconomic status. For example, Abásolo et al. [1] analyse horizontal equity in access and utilisation of public GP services in Spain applying probit models to data from the Spanish national health survey of 1993. They find evidence of horizontal inequity that favours men and women in lower socioeconomic groups, women with fewer educational qualifications and men who are not working. Urbanos [2] uses the Spanish national health surveys of 1987, 1993 and 1995 and probit models to explain the inequality in the use of public health care services among individuals with similar medical needs. For the case of public GP services, she concludes that in general, more highly educated individuals are more likely to underconsume compared to the population whose education is minimal. In addition, her results indicate that the propensity towards excessive consumption increases as social class decreases, except for the year 1995. Morris et al. [3] investigate inequity in the use of general practitioner consultations and other health care services of the National Health Service in England. Their analysis is based on pooled data from three rounds (1998, 1999, 2000) of the English health survey. They use linear probability models to estimate the probability of visiting a GP finding evidence that those with lower education attainment are more likely to visit their GP. They also find that income (and being unemployed) is negatively correlated with use. Regidor et al. [4] study horizontal equity in the utilisation of public GP services in Spain using data from the 2003 Spanish national health survey. They estimate the probability of visiting a GP through a binomial regression considering education level and social class variables among covariates. They show the existence of inequity in public GP visits favouring individuals with lower education levels and also those belonging to lower social classes. Abásolo et al. [5] using the national health survey of 2006 and through multilevel binary logistic functions, analyse horizontal equity in the probability of utilisation of public health care services in Spain by socioeconomic groups. They find that education level is negatively related to the probability of visiting a public GP. A negative correlation is also found for income and for being employed. González and Clavero [6] analyse, through dynamic models, the explanatory factors of income related inequality in health care utilisation in Spain (they consider only those patients covered by the public health care system, excluding those individuals with private insurance). They use the European Community Household Panel. They find that individuals with higher income have a lower probability to visit a GP; in addition, education is the main determinant of the inequality of utilisation as those with primary studies or less have a higher probability to visit a GP. Similar conclusions are reached by Gonzalez and Clavero [7] in another study with the same data base in which they study the inequalities in health care utilisation caused by the double health insurance coverage in Spain.
Other studies show no evidence of horizontal inequity by socioeconomic status or show horizontal inequity favouring the most advantaged. For example, Van Der Heyden et al. [8], with data from the 1997 Belgian National Health Interview Survey, used logistic regressions to adjust by demographic and health status variables (in addition to socioeconomic variables), finding no evidence of inequity in the probability to use GP services by education or income variables. San Sebastian et al. [9] use the 2006, 2010 and 2014 cross sectional household surveys (“health on equal terms surveys”) for the four northernmost counties in Sweden. Using probit regressions and analysing and computing the horizontal inequity indexes find a prorich inequality in access to GPs, particularly for 2010 and 2014, being the most important contributor income and education, attributable they argue to underutilisation among those with lowest education. Wagenius et al. [10] use the same data base than [9] but just for 2014, selecting young respondents between 16 and 24 years. They undertake binomial regressions to explain the probability of visiting the GP, youth clinics and nurses. They found propoor inequity in access to GP services both for men and women. However, they also found a higher probability of using GP services for those with higher education levels for men, but not for women for whom there is no evidence of horizontal inequity.
A common characteristic in all of the above studies is that they use binomial models (eg. logit, probit, linear probability models). That is, the dependent variable is a binary variable, so that if a respondent reports that she visited a public GP, then it counts as 1, otherwise it counts as zero. In addition, when these studies analyse also the equity in the utilisation of other health care services (specialist, GP, emergency services, etc.), these are considered separately, in different binomial models, independent of each other. That is, all the above studies take a compartmentalized approach to the problem of horizontal equity in access to health care services. However, in the analysis of equity in access to a particular service (eg. public GP service) it may be the case that the respondent reports to have used the health care system, but instead of visiting a public GP, she was in another point of the public health care system (i.e. visiting a specialist doctor) or was voluntarily visiting a private doctor. Binomial econometric models for equity in access to public GP services would also consider such visits as zeros (i.e. nonutilisation), thus potentially obtaining biased conclusions regarding horizontal equity in access. With the aim to test for this potential bias, the objective of this research is to study whether the compartmentalized approach to the problem of horizontal equity in access to public GP services reaches the same conclusions as a comprehensive approach in which, in addition to adjusting for variables of health care need and other demographic and socioeconomic characteristics, we also consider the use of other public and private health care services, thus avoiding an inappropriate (we believe) treatment of zeros. We undertake this analysis using the same data base (the Spanish National Health Survey 2016–17) in order to maximise comparability of the two approaches. To the best of our knowledge, this exercise has not been carried out in the related literature so far. In what follows, Section 2 presents the data, variables and methods used. The results are presented in Section 3 and discussed in Section 4. Conclusions are presented in Section 5.
Data, variables and methods
Data and variables
The data used in this research come from the 2016–17 Spanish National Health Survey. This is a facetoface, crosssectional populationbased survey that employs a threestage, stratifiedrandom design to identify samples of adults (aged 15 or over) within the target number of survey interviews (23,089). The firststage units are the census sections, which are stratified according to the size of the municipality to which the section belongs. The secondstage units are the main households. The third stage units are chosen from a list of persons within the household who can be interviewed and asked to fill in the questionnaire at the time the survey is carried out (for more details, see [11]). Data on health care utilisation, selfreported morbidity and other demographic and socioeconomic characteristics are also collected.
On the one hand, for the compartmentalised approach, the dependent variable of interest is a dummy variable representing whether or not the individual has visited a public sector GP in the past four weeks. On the other hand, for those who report not having visited a public GP, we have also considered whether they have visited instead a public specialist doctor or a private doctor (specialist or general practitioner); thus for the comprehensive approach the dependent variable refers to whether or not the individual has visited one of the four different health services. The type of visit public or private sector visit, and whether it was a GP visit or a specialist visit is only known for the last visit of the respondent and we focus our study on this last visit. Therefore, the different outcomes of the dependent variables have been built upon the combination of four survey questions. Firstly, respondents are asked whether they have visited a GP in the past four weeks. Secondly, respondents are asked whether they have visited a specialist doctor in the past four weeks. Thirdly those who have visited both sorts of doctors in the same period are asked the nature (GP or specialist) of the very last visit. Finally, individuals are asked whether the doctor they had consulted the last time was in the public health system, was from a private insurance company or was in a private consultation. We considered the former case as a public visit, whilst the second and third cases were considered as private visits. Unfortunately, information about other health care services like emergency services is referred to a different period (12 months) so it has not been possible to take them into account in the analysis.
Regarding the explanatory variables, socioeconomic status (main focus of this paper) is proxied through three variables: education, social class and employment situation. Education is measured by a categorical variable indicating the highest level of education achieved by the respondent: no studies, primary studies, secondary studies and university studies. Social class of the reference person in the household is also measured by a categorical variable with four categories based on the National Classification of Occupations: high social class (Directors and managers with university degrees), mediumhigh social class (intermediate professions and selfemployed), mediumlow social class (skilled and partlyskilled occupations), and low social class (unskilled workers). Employment situation of the respondent is measured by a categorical variable with three possible activity statuses: employed (the individual is currently employed), unemployed (the individual is currently unemployed) or inactive (the individual is retired, disabled or housemaker). We have not included an income variable due to the number of missing cases (24% in the Spanish National Health Survey of 2017), so for the purpose of this paper, we considered sufficient the three proxies for socioeconomic status included in the analysis.
In order to adjust by health care need, we have considered selfreported measures of individuals’ health state. These include a categorical indicator that records whether individuals considered their general health during the twelve months prior to the survey to be ‘very good’, ‘good’, ‘fair’, ‘poor’ or ‘very poor’. We have also considered a dummy variable indicating whether the respondent reports the presence of any chronic conditions. Regarding demographic characteristics, we have included age and sex of the respondents.
Methods

Binary models for a compartmentalized approach.
In a compartmentalized approach, the dependent variable (Y) is a binary variable with only two possible values: 1 if the individual has used the public GP service, and 0 otherwise. The binary response models most commonly used are the logit and probit models [12].
Both models assume that the conditional probability of Y = 1 given a vector of covariates x is:
where β is the vector of parameters of the model. A binary response model is referred to as a logit model if F is the cumulative logistic distribution function, and it is called a probit model if F is the cumulative normal distribution function. Both logit and probit models can also be described as latent variable models, where there is an auxiliary random variable Y^{∗} = x ′ β, where the dependent variable Y can be viewed as an indicator for whether this latent variable is positive (Y = 1 if Y^{∗} > 0 and 0 otherwise):
\( \mathrm{logit}\ \mathrm{model}:\frac{P\left(Y=1x\right)}{1P\left(Y=1x\right)}={e}^{x\prime \beta } \)
probit model : P(Y = 1 x) = P(Y^{∗} > 0) = Φ(x ′ β)
The logistic and normal distributions are both symmetrical around zero, but the logistic distribution has fatter tails. Both methods yield similar inferences, although logit models have become more popular in social sciences because coefficients can be interpreted in terms of odds ratios.

Multinomial models for a comprehensive approach.
In a comprehensive approach, the dependent variable is a nominal variable with five categories, indicating the health service used on the last visit (public or private GP or specialist), in addition to the possibility of not having used any service in the last four weeks. In this case, it is necessary to consider multinomial (logistic or probit) regressions [13, 14].
A multinomial logistic model for k possible outcomes can be described as a combination of (k1) independent binary logit models, in which one outcome is chosen as a reference.
\( \frac{P\left(Y=1x\right)}{P\left(Y=kx\right)}={e}^{x\prime {\beta}_1},\frac{P\left(Y=2x\right)}{P\left(Y=kx\right)}={e}^{x\prime {\beta}_2},\dots, \frac{P\left(Y=k1x\right)}{P\left(Y=kx\right)}={e}^{x\prime {\beta}_{k1}} \)
Since the k probabilities must sum to one, it is straightforward to obtain their expressions:
\( P\left(Y=1x\right)=\frac{e^{x\prime {\beta}_1}}{1+\sum \limits_{j=1}^{k1}{e}^{x\prime {\beta}_j}},P\left(Y=2x\right)=\frac{e^{x\prime {\beta}_2}}{1+\sum \limits_{j=1}^{k1}{e}^{x\prime {\beta}_j}},\dots, \)
\( P\left(Y=k1x\right)=\frac{e^{x\prime {\beta}_{k1}}}{1+\sum \limits_{j=1}^{k1}{e}^{x\prime {\beta}_j}},P\left(Y=1x\right)=\frac{1}{1+\sum \limits_{j=1}^{k1}{e}^{x\prime {\beta}_j}} \)
The multinomial logit relies on the assumption of independence of irrelevant alternatives (IIA), which states that the odds of preferring one alternative over another do not depend on the presence of other alternatives. Multinomial probit can be used as an alternative when IIA assumption is violated.
The multinomial probit model can be described in terms of a latent variable model. The latent variables for each alternative are \( {Y}_1^{\ast }=x^{\prime }{\beta}_1+{\varepsilon}_1 \), \( {Y}_2^{\ast }=x^{\prime }{\beta}_2+{\varepsilon}_2 \), …, \( {Y}_k^{\ast }=x^{\prime }{\beta}_k+{\varepsilon}_k \), where ε = (ε_{1}, ε_{2}, …, ε_{k}) follows a multivariate normal distribution with mean vector zero and a specified covariance matrix. Then
This model does not respect IIA, since error terms can be arbitrary correlated.

An alternative Binomial model for a comprehensive approach.
A different way to address the existence of multiple alternatives of health service utilisation but without including them in the regression model would be to drop from the sample those individuals who in the last visit reported to use other health services different from public GP services (private GP services or private or public specialist services). Then a binomial probit or logit model would be estimated for the remaining sample of individuals. However, if individuals finally selected are not at random, this would lead to biased estimates. The probit model with sample selection [15] allows to test for selection bias. It considers two probit models: one for the probability of visiting a public GP (probit equation) and another one for the probability of being selected for the probit equation which, in our case is 1 minus the probability of visiting a doctor other than the public GP (selection equation). It is assumed that the error terms for both equations follow a bivariate normal distribution with a rho correlation. When rho is nonzero, the standard probit model applied to the probit equation yields biased results.
For the model to be well identified, the selection equation should have at least one variable that is not in the probit equation. For this reason, the selection equation includes a dummy variable for having double coverage, that is, individuals who have private health insurance in addition to public health care insurance, on the basis that having private insurance may explain participation (i.e. through the demand of private GP or specialist visits) but probably does not affect the probability to visit a public GP.
Results
Table 1 shows the descriptive statistics for the dependent and independent variables. Regarding the dependent variable, 24.8% of individuals reported that the last visit they had in the past four weeks was made to a public GP. According to the compartmentalised approach, the remaining 75.2% would be considered as zeros in the binomial models. There were 12.1% of individuals who reported that the last visit they had in the past four weeks was to another professional, different from a public GP. Specifically, 7.7% visited a public specialist, 2.9% visited a private specialist and 1.5% visited a private GP (all of them referred to the last visit in the same 4week period). According to the comprehensive approach, these 12.1% of individuals would not count as zeros (zeros would be the case in the remaining 63.1%, who did not visit a doctor at all).
To begin with the compartmentalised approach, estimates for binomial logit and probit models are presented in Table 2. Sign and statistical significance of coefficients for both models are very similar. With respect to health care need variables, as expected, individuals who have a worse health state have a higher propensity of reporting visits to a public GP within the past four weeks, with a clear gradient as the state of health worsens. The probability of visiting a public GP is also greater for those individuals who suffer a chronic disease.
As for socioeconomic factors, regarding education level, the results show that in comparison with the reference category (primary studies or less), individuals with secondary studies have a significantly lower probability of visiting a public GP only for the probit model (p < 0.10). However, individuals with occupational training and individuals with university studies have a significantly lower probability of visiting a public GP for both models (p < 0.01). A comparison of the predicted probabilities fixing values for the relevant categories of the socioeconomic variables is presented in Table 3. Predicted probability for individuals with primary studies or less is 25.9% for binomial logit and 26.0% for binomial probit. Predicted probability for occupational training is 23.5% (for both binomial logit and probit) and predicted probabilities for those with university studies are 22.5% for binomial logit and 22.7% for binomial probit. Coefficients and predicted margins show, therefore, a clear gradient as education level decreases.
Regarding social class, as compared with high social class, the lower the social class, the higher the probability of visiting a public GP (p < 0.01 for all coefficients in both binomial models), showing a clear gradient. While individuals in the highest social class have a predicted probability to visit a public GP of 20.3% (for binomial logit) and 20.4% (for binomial probit), individuals belonging to the low social class have 27.5 and 27.6%, respectively (Table 3). Finally, being employed (as compared with being inactive) is also negatively related to the probability of visiting a public GP both for the binomial logit and probit models (p < 0.01). Regarding demographic variables, older adults (aged more than 64) and females are more likely to visit a public GP.
With respect to the comprehensive approach, estimates for the multinomial logit and probit models of the public GP equation are presented in Table 4 (results for the equations of the rest of health care services can be seen in Table 6 of the Appendix). As with the binomial case, sign and statistical significance of coefficients for both models are similar. Regarding health care need, individuals with worse selfreported health states are more likely to visit a public GP with a gradient as health state worsens. Also, those with a chronic disease have a higher probability of visiting a public GP.
Regarding socioeconomic factors, social class is also negatively related to the probability of visiting a public GP. Results for the multinomial models (Table 3) indicate that individuals belonging to the high social class have a 20.5% predicted probability to visit a public GP for the multinomial logit (20.7% for the multinomial probit) whilst individuals belonging to the low social class have on average a 27.6% probability for the multinomial logit (27.7% for the multinomial probit). This difference of seven percentage points that results to be statistically significant for both models is quite similar to that obtained through the binomial models. However, as can be seen in Table 4, education level no longer influences the probability of visiting a public GP for both multinomial logit and probit models. In addition, the negative effect of being employed on the probability of visiting a GP also remains in this comprehensive approach.
To test the null hypothesis of independence of irrelevant alternatives (IIA), we conducted the Small–Hsiao tests [16]. None of the five tests undertaken by omitting each alternative was significant (P > 0.1), indicating no evidence of violation of the independence of irrelevant alternatives assumption. So, both the multinomial logit (that assumes IIA) and the multinomial probit model (that does not require the fulfilment of the IIA assumption) are valid to undertake this comprehensive approach.
One alternative way to proceed would have been to drop the 12.1% of individuals who reported having visited other doctors different from a public GP in their last visit and follow a compartmentalised approach with the remaining individuals of the sample. Estimates for the probit model with sample selection are shown in Table 5. The correlation coefficient (rho) is statistically different from zero (p < 0.01), indicating the presence of selection bias. In addition, sign, magnitude and tratios of coefficients of the probit model with selection are quite different to those of the simple probit estimation. Thus, we can reject the null hypothesis that there is no selection bias if we drop those observations that instead of using public GP services, have used specialist services or private GP services.
Discussion
Horizontal equity in access to public general practitioner (GP) services by socioeconomic group has been addressed econometrically in the related literature by testing the statement “equal probability of using public GP services for equal health care needs, regardless of socioeconomic status” using healthcare surveys that have information on utilisation, health status and other socioeconomic and demographic characteristics of individuals. In this research, we have claimed that there are two different approaches that could lead to different conclusions.
One is the compartmentalised approach that has been the common way to address this issue in the literature. Under this alternative, if the respondent reports that she visited a public GP, then it counts as 1, otherwise it counts as zero. Our results to this approach applying binomial logit and probit models to data for the Spanish National Health Survey 2016–17, indicate that the lower the education level and the lower the social class, the higher the probability of visiting a public GP, thus showing a prolower socioeconomic group horizontal inequity in access, in line with a great part of the related literature [1,2,3,4,5,6,7]. Other related studies show no evidence of horizontal inequity by socioeconomic status or sometimes inequity favouring the most advantaged [8,9,10]. The latter studies consider overall (public and private) GP services. The extent to which this difference can be partly explained by the difference in the definition of the dependent variable (public versus overall GP services), by the consideration of different set of socioeconomic covariates, by other methodological issues, by the differences in the organization and provision of health care in the countries analysed (Spain, Belgium, Norway, Sweden or the UK) or simply by a different empirical evidence, is not known and would involve further research. With this caveat in mind, what is common to all these studies is that they share the compartmentalised approach.
As it has been highlighted, it may be the case that the respondent’s last visit was not a public GP but was say a specialist doctor or a private GP. Binomial econometric models would wrongly consider such visits as zeros (i.e. nonutilisation), thus potentially obtaining biased conclusions regarding horizontal equity in access. In order to test for this potential bias, we undertook a comprehensive approach through multinomial logit and probit models to account for utilisation different from the public GP in the last visit. On the one hand, conclusions regarding social class remain very similar to those of the binomial logit and probit models: the lower the social class, the higher the probability of visiting a public GP, also in line with some of the revised previous studies [1, 2, 4]. In addition, the lower propensity to visit a public GP of those employed (with respect to those inactive individuals) found in the binomial approach, remains under the multinomial approach (a relatively higher opportunity cost of time may help to explain this negative effect). However, interestingly, we find that the education effect found in the compartmentalised approach, by which individuals with lower education level have a higher propensity to visit a public GP (a finding which would be consistent with one of the predictions of the Grossman model, in which individuals with higher education level are more efficient at producing health, and therefore they need to use health care services less [17]), disappears under the comprehensive approach. The lack of evidence to reject the independence of irrelevant alternative hypothesis indicates that both multinomial logit and probit models are appropriate. The fact that both models give similar results reinforces the strength of our conclusions. This evidence contradicts previous findings regarding the effect of education on access to public GP services [1,2,3,4,5,6,7] but is in line with results obtained in other studies [8,9,10]. One possible explanation is that individuals with higher education levels (as compared with those with primary studies or less) instead of using public GP services, are using relatively more public specialist services (and/or other private specialist care or GP services) either as substitute or complementary, thus, not losing (if not increasing) access to the public health care system. For example, Negrín et al. show evidence that suggests that the relatively higher propensity to combine public and private specialist visits of those with higher education levels, is related to a complementary use of private health care services in order to reduce waiting times in the public sector specialist care [18]. And our results regarding the multinomial logit and probit models seem to reinforce such hypothesis, according to the positive and significant coefficients of the education categories for the corresponding three equations (see Table 6 in the Appendix).
Another possible approach would have been to drop those individuals who reported to be in another part of the health system and then undertake a binomial logit/probit model to estimate the probability of visiting a public GP. However, if as a consequence of dropping that part of the sample, a problem of selection bias arises, then this approach is not reliable. This is what actually happened with our database.
To conclude, although we have undertaken our analysis with the same survey under the same health care system (so any bias encountered can be attributed to the differences in approaches) we must be cautious when generalising our results to other contexts. However, we strongly suggest that it is appropriate to undertake a comparison between both approaches. If the binomial/compartmentalised approach which, as compared with the multinomial approach, is more straightforward to interpret and allows for other outputs like the construction of inequity indices gives results that are not different from those of the multinomial/comprehensive approach, then this comparison exercise adds robustness to their results. On the other hand, when a significantly different result is obtained, a comprehensive approach seems to be more appropriate. If that is the case, the conclusion reached by the binomial approach about horizontal inequity would be incorrect due to the fact that individuals using other health care services are wrongly considered as nonusers and this situation should not be regarded as inequitable.
Regarding the paper’s limitations, four points must be made. First, it should be noted that our main database (Spanish National Health Survey) does not have enough information to analyse horizontal equity taking into account the intensity of use (i.e. number of GP visits), as we only have the relevant information for the very last visit. A comprehensive approach that considers the number of visits/contacts would be much more demanding informationwise and probably very difficult to undertake. Second, conclusions on equity in access to public GP visits must be taken with caution as the health survey used only has information on doctors’ visits in the past four weeks. A longer reference period might change the conclusions, although it must also be said that we would then face a greater risk of recall bias. Third, health care need is not the same for different health services (general practitioner vs. specialist). However, the available data do not allow for a distinction in the health needs of different services. In any case, we consider that selfassessed health and selfreported chronic diseases are good predictors of demand for public GP visits (whilst for specialist visit, particularly those of the public health care system, a good indicator of health care need would be the specialist appointment itself, insofar as it must have been indicated previously by a doctor (GP or specialist). Finally, the comprehensive approach assumes that the five options are available for all the sample. Civil servants in Spain have the right to choose annually between the public provider and a private insurance provider. Thus, those civil servants who have chosen a private health care provider should not have access to public services. However, they actually do have access to public services, possibly because of the little control over this fraudulent double coverage (SanchezBayle and Beiras [19] and Rodríguez and Stoyanova [20]). In addition, it must be noted that this group represents the 3.05% of the whole sample.
Conclusion
The conclusion supported by a considerable part of the compartmentalised approach literature that access to public GP services favours the lower socioeconomic groups is debatable, at least regarding the effect of education level. The multinomial analysis suggests that the social class effect remains, but the education level effect that resulted from the binomial analysis disappears when the utilisation of other health care services (apart from public GP health care services) is considered within the analysis. An analysis of horizontal equity in access to a particular health care service undertaken following a compartmentalised approach should be compared with a comprehensive approach in order to test that there is no bias as a consequence of considering as zeros the utilisation of other types of health care.
Availability of data and materials
The data base is available in open in the following website:
http://www.mscbs.gob.es/estadEstudios/estadisticas/encuestaNacional/encuesta2017.htm
References
 1.
Abásolo I, Manning R, Jones A. Equity in utilization of and access to publicsector GPs in Spain. Appl Econ. 2001;33:349–64 https://doi.org/10.1080/00036840122511.
 2.
Urbanos R. Explaining inequality in the use of public health care services: evidence from Spain. Health Care Management Science. 2001;4:143–57 https://doi.org/10.1023/A:1011461728666.
 3.
Morris S, Sutton M, Gravelle H. Inequity and inequality in the use of health care in England: an empirical investigation. Soc Sci Med. 2005;60(6):1251–66 https://doi.org/10.1016/j.socscimed.2004.07.016.
 4.
Regidor E, Martínez D, Calle ME, Astasio P, Ortega P, Domínguez V. Socioeconomic patterns in the use of public and private health services and equity in health care. BMC Health Serv Res. 2008;8(183):1–9 https://doi.org/10.1186/147269638183.
 5.
Abásolo I, Pinilla J, Negrín M. Equidad en la utilización de servicios sanitarios públicos por Comunidades Autónomas en España: un análisis multinivel. Hacienda Pública Española. 2008;187(4):87–106 https://dialnet.unirioja.es/servlet/articulo?codigo=2841797.
 6.
González ML, Clavero A. Análisis de las desigualdades socioeconómicas en la utilización de asistencia sanitaria mediante modelos dinámicos. 2008. Hacienda Pública Española. 2008;186:9–42 https://dialnet.unirioja.es/servlet/articulo?codigo=2762719.
 7.
González ML, Clavero A. Inequalities in health care utilization in Spain due to double insurance coverage: an Oaxacaransom decomposition. Soc Sci Med. 2009;69(5):793–801 https://doi.org/10.1016/j.socscimed.2009.06.037.
 8.
Van Der Heyden JHA, Demarest S, Tafforeau J, Van Oyen H. Socioeconomic differences in the utilisation of health services in Belgium. Health Policy. 2003;65(2):153–65 https://doi.org/10.1016/S01688510(02)002130.
 9.
San Sebastián M, Mosquera PA, Ng N, Gustafsson PE. Health care on equal terms? Assessing horizontal equity in health care use in northern Sweden. Eur J Pub Health. 2017;4:637–43 https://doi.org/10.1093/eurpub/ckx031.
 10.
Wagenius CM, San Sebastián M, Gustafsson PE, Goicolea I. Access for all? Assessing vertical and horizontal inequities in health care utilisation among young people in northern Sweden. Scandinavian Journal of Public Health. 2019;47:1–8 https://doi.org/10.1177/1403494818774965.
 11.
Ministry of Health, Social Services and Equity. National Health Survey in Spain, 2006–07, 2011–12 and 2016–17. Householder’s questionnaire and adult’s questionnaire [in Spanish] [Available at: http://www.mscbs.gob.es/estadEstudios/estadisticas/encuestaNacional/encuesta2017.htm
 12.
Greene W. Econometric Analysis. Fifth ed: Prentice Hall; 2003.
 13.
Engel J. Polytomous logistic regression. Statistica Neerlandica. 1988;42(4):233–52 https://doi.org/10.1111/j.14679574.1988.tb01238.x.
 14.
Elrod T, Keane MP. A factoranalytic probit model for representing the market structure in panel data. J Mark Res. 1995;32(1):1–16 https://doi.org/10.1177/002224379503200103.
 15.
Van de Ven WPMM, Van Pragg BMS. The demand for deductibles in private health insurance: a probit model with sample selection. J Econ. 1981;17:229–52 https://doi.org/10.1016/03044076(81)900282.
 16.
BenAkiva M, Lerman SR. Discrete Choice Analysis: Theory and Application to Travel Demand. Cambridge: MIT Press; 1985.
 17.
Wagstaff A. The demand for health: some new empirical evidence. Health Econ. 1986;5(3):195–233 https://doi.org/10.1016/01676296(86)900159.
 18.
Negrín MA, Pinilla J, Abásolo I. Simultaneous use of public and private specialist services as a potential strategic behaviour in access to the National Health Service: is there a pattern by socioeconomic level? Soc Justice Res. 2019;32(3):255–76 https://doi.org/10.1007/s1121101900330y.
 19.
SanchezBayle M, Beiras H. The trouble with conservative healthcare counterreforms in Spain. In: Sen K, editor. Restructuring health services: changing contexts and comparative perspectives. New York: ZedBooks; 2003.
 20.
Rodríguez M, Stoyanova A. The effect of private insurance access on the choice of GP/specialist and public/private provider in Spain. Health Econ. 2004;13:689–703 https://doi.org/10.1002/hec.832.
Acknowledgements
Not applicable.
Funding
This study was funded by the Spanish Ministry of Economy and Competitiveness (project ECO2016–79588R). We are grateful for such financial support.
Author information
Affiliations
Contributions
MAN, JP and IA had the idea for the paper, requested the data, performed the literature search, wrote the introduction, methods, results, analysis, discussion and conclusions, and approved the final manuscript.
Corresponding author
Correspondence to I. Abásolo.
Ethics declarations
Ethics approval and consent to participate
Not applicable
Consent for publication
Not applicable
Competing interests
There are no conflicts of interest for any of the authors.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Negrín, M.A., Pinilla, J. & Abásolo, I. Horizontal equity in access to public GP services by socioeconomic group: potential bias due to a compartmentalised approach. Int J Equity Health 18, 187 (2019). https://doi.org/10.1186/s1293901910912
Received:
Accepted:
Published:
Keywords
 Public GP visits
 Horizontal equity in access
 Bias
 Health care services
 Specialist visits
 National health surveys