Indicators
In terms of input-output indicator selection, most studies select the corresponding indicators from three types of health resources: human, financial, and material [18–19]. In contrast, the service efficiency of PHC institutions is a complex system that requires consideration of many factors, which includes multiple input and multiple output indicators. The results of different index combinations may be significantly different [22–23]. The inputs of PHC institutions are resources in several areas such as human, material, and financial resources, mainly including the number of employees, the number of health technicians, equipment and fixed assets, management costs and so on. The output indicators Usually include the number of emergency visits, hospitalizations, surgeries, financial income, and overhead costs [29–34].
Based on the existing literature and the availability of the corresponding indicators, this study selected the number of PHC institutions (X1), the actual number of beds (X2) and the number of health technicians (X3) as input indicators, the number of consultations (Y1) and the number of discharges (Y2) as output indicators from two dimensions—service volume and economic efficiency (Table 1).
The input indicator data of this study from 2011–2020 were obtained from the statistical indicators of 31 provinces (municipalities directly under the central government and autonomous regions) in the China Health and Health Statistical Yearbook, because of inconsistent statistical standards, Hong Kong, Macao and Taiwan provinces were not included. The map data were based on public geographic data downloaded from the Institute of Geographical Sciences and Natural Resources, Chinese Academy of Sciences, and the shp files containing only 31 provinces (municipalities directly under the central government and autonomous regions) were self-produced after processing with ArcGIS software.
Table 1
PHC institution service efficiency indicators system
Item
|
Specific Indicators
|
Explanation of Indicators
|
Significance of Indicators
|
Input
|
X1: Number of Health Institutions (unit)
|
PHC institutions in China mainly composed of community health service centers (stations), street health centers, township health centers, village health offices, outpatient clinics, and clinics (infirmaries)
|
Reflecting the input of medical and health material resources in each year in the study area.
|
X2: Number of Beds (unit)
|
The number of beds is the sum of the number of beds in various medical and health institutions each year.
|
X3: Number of health technicians (person)
|
Health technicians mainly include registered nurses, licensed (assistant) physicians, pharmacists, health supervisors, technicians, and other health technicians.
|
Reflecting the investment scale of medical and health human resources in each year in the study area.
|
Output
|
Y1: Number of Outpatients and Emergency Visits (10,000 person)
|
Outpatient and emergency visits are the sum of outpatient visits and emergency visits.
|
Reflecting the supply of medical and health services of the research subjects in that year.
|
Y2: Number of hospital discharges(10,000 person)
|
Refers to all discharges after hospitalization during the reporting period
|
Research Methods
DEA
1.DEA- SE - SBM model
The DEA method is one of the more mature methods of nonparametric efficiency analysis methods [35–36]. Its traditional models include CCR, BCC and other models, which have the following two drawbacks: firstly, they are all measured based on radial angles and do not consider the slack variables of input and output components, the number of discharges, bed utilization and other outputs do not vary proportionally in the PHC institution service efficiency which made the measurement results not accurately; secondly, the measured efficiency values range from 0 to 1; at the same time, the units that achieve efficiency (efficiency value of 1) cannot be compared with.
Tone K proposed a non-radial, non-angle slack-based measure (SBM) model based on the traditional DEA model, which can effectively compensate for the shortcomings of some traditional models that do not incorporate slack variables when measuring low efficiency and improve the accuracy [37]. Based on this, Tone K combined the SBM model with the super-efficiency model [38] and proposed the super-efficiency slack measurement model (SE-SBM), which can further compare decision units (DMU) that simultaneously reach an efficiency value of 1 to compensate for the second deficiency mentioned above.
Therefore, in this study, the DEA-SE-SBM model based on input volume measurement was used in order to accurately reflect the efficiency of primary health care institution service in 31 provinces in China. The formula is shown in (1).
$$s.t.\left\{\begin{array}{c}{x}_{ik}=\sum _{j=1,j\ne k}^{n}{x}_{ij}{\lambda }_{j}+{s}_{i}^{-}\\ {y}_{rk}=\sum _{j=1,j\ne k}^{n}{y}_{rj}{\lambda }_{j}-{s}_{r}^{+}\\ {\lambda }_{j},{s}_{i}^{-},{s}_{r}^{+}>0\\ i=\text{1,2},\dots .,m\\ r=\text{1,2},\dots .,z\\ j=\text{1,2},\dots ,n;j\ne k\end{array}\right.$$
1
n is the number of measured regions (n = 31 in this study), j is the j-th province (j = 1, 2, ..., n), and each region has m inputs and z outputs. \({x}_{ij}\)and \({y}_{rj}\) denote the i-th input and r-th output of the j-th region. \({x}_{ik}\) and \({y}_{rk}\)denote the i-th input and r-th output of the k-th region. ρ* denotes the efficiency value of each region. λ denotes the weight vector; \({s}_{i}^{-}\) and \({s}_{i}^{+}\)are the slack variables of the input vector and output vector, which denote the savings of health care resource inputs and the increase of health care service outputs in the pilot regions, respectively. In addition, according to model (1), the service efficiency of PHC institution in 31 provinces in China from 2011–2020 can be calculated. As shown in Eq. (2), \({Score}_{kt}\)denotes the efficiency of health care services in city k in year t.
$${Score}_{kt}={\rho }_{kt}^{*}\left(k=\text{1,2},3\dots ..31;t=\text{2011,2012},\dots .,2020\right)$$
2
2.DEA-Malmquist model
In terms of efficiency evaluation, the DEA-SE-SBM model can only measure the efficiency within a certain time. It cannot continuously compare different time periods and cannot reflect the dynamic evolutionary characteristics of efficiency. Therefore, this study uses the DEA-Malmquist model to further calculate the total factor productivity change (TFPCH) of health services [39–41] to reveal its dynamic evolution.
Let (\({x}^{t}\),\({y}^{t}\),) and (\({x}^{t+1}\), \({y}^{t+1}\)) denote the input-output vectors of the health care system in periods t and t+1, respectively, the TFPCH of the two adjacent periods are:
$$TFPCH=M\left({x}^{t},{y}^{t},{x}^{t+1},{y}^{t+1}\right)={\left[\frac{{D}^{t}\left({x}^{t+1},{y}^{t+1}\right)}{{D}^{t}\left({x}^{t},{y}^{t}\right)}\times \frac{{D}^{t+1}\left({x}^{t+1},{y}^{t+1}\right)}{{D}^{t+1}\left({x}^{t},{y}^{t}\right)}\right]}^{\frac{1}{2}}$$
3
\({ D}^{t}\) (\({x}^{t}\),\({y}^{t}\)) and \({D}^{t}\)(\({x}^{t+1}\), \({y}^{t+1}\)) denote the distance function of the desired region in period t and period t+1, respectively, with period t serving as the technical reference. \({D}^{t+1}\)(\({x}^{t}\), \({y}^{t}\)) and \({D}^{t+1}\)(\({x}^{t+1}\),\({y}^{t+1}\)) have similar meanings. According to the decomposition of the Malmquist index by Fare [42], it can be further decomposed into two components: the index of technical efficiency change (EFFCH) and the index of technical change (TECHCH). In addition, with variable returns to scale (VRS) as an assumption, EFFCH can be decomposed into pure efficiency change index (PECH) and scale efficiency change index (SECH). If TFPCH>1, it means that the total factor productivity of the test subject has improved during t ~ t + 1. If TFPCH < 1, it indicates a decrease. If TFPCH = 1, it indicates that it remains unchanged. TECHCH > 1 represents technological progress and vice versa. PECH represents change in efficiency due to change in management level, and PECH > 1 represents improvement in management level and vice versa. SECH > 1 represents optimization of scale efficiency and vice versa.
Exploratory Spatial Data Analysis
ESDA is a research method in spatial econometrics that is mainly used to study the correlation between a phenomenon and the attribute values of its neighboring units in geographic space [43], to measure the aggregation or dispersion of the spatial element attributes of the research object [44]. It mainly includes global spatial autocorrelation analysis and local spatial autocorrelation analysis [45]. In this study, global Moran’s I and local Moran's I were used to reveal the overall spatial correlation and internal correlation characteristics of the service efficiency of PHC in China.
1. Global spatial autocorrelation
Global spatial autocorrelation analysis measures the degree of spatial correlation and variation among the study units from the overall level. GeoDa 1.18 software was used to calculate the global Moran's I to quantify the overall spatial correlation of the service efficiency of PHC. The calculation formula is as follows:
$$I=\frac{n\sum _{j=1}^{n}\sum _{k=1}^{n}{w}_{jk}\left({\rho }_{j}^{*}-\stackrel{-}{{\rho }^{*}}\right)\left({\rho }_{k}^{*}-\stackrel{-}{{\rho }^{*}}\right)}{\sum _{j=1}^{n}\sum _{k=1}^{n}{W}_{jk}\sum _{j=1}^{n}{\left({\rho }_{j}^{*}-\stackrel{-}{{\rho }^{*}}\right)}^{2}}$$
4
I is the global Moran's I; \({\rho }_{j}^{*}\) and \({\rho }_{k}^{*}\) are the PHC institution service efficiency values of each region j and k. The \({W}_{jk}\) spatial weight matrix, which measures the spatial relationship between regions j and k; the Rook spatial adjacency method was used in this study. \(\stackrel{-}{{\rho }^{*}}\) is the average value of health service efficiency. The global Moran's I take values in the range [-1,1][46]. When I > 0, it indicates that there is a spatial positive correlation in efficiency, the greater the spatial correlation; when I < 0, it indicates a spatial negative correlation, the more obvious the spatial difference. When I = 0, it indicates that there is no spatial correlation and exhibits a random distribution [47]. The z-value is also usually needed to assist in determining the significance of the Moran index.
2.Local spatial autocorrelation
Local spatial autocorrelation analysis is used to analyze the degree of spatial correlation between each spatial object and its neighboring units in a certain region to reflect the local characteristic differences in the distribution of spatial objects. It can be used as a supplement to the global spatial correlation analysis to compensate for the possible local spatial potential instability [48]. This study will reveal the degree of correlation between the PHC institution service efficiency in a province and its neighboring regions by calculating the local Moran's I. For province j, the local Moran's I index can be expressed as:
$${I}_{j}=\frac{n({\rho }_{j}^{*}-\stackrel{-}{{\rho }^{*})}}{\sum _{j=1}^{n}{({\rho }_{j}^{*}-\stackrel{-}{{\rho }^{*}})}^{2}}\sum _{k=1}^{n}{W}_{jk}\left({\rho }_{k}^{*}-\stackrel{-}{{\rho }^{*}}\right)$$
5
The variables in the above equation have the same meaning as those in Eq. (4), when \({I}_{j}\) is positive, it means that the province is similar to its neighboring provinces in terms of attributes, and when \({I}_{j}\) is negative, it means that it is not similar. Meanwhile, four different spatial aggregation patterns were shown by plotting the Moran scatter plot and the Local Spatially Associated Aggregation Indicator (LISA) [49–50].
The first and third quadrants in the Moran scatter plot represent spatial agglomeration effects, and the second and fourth quadrants represent spatial divergence effects. The first quadrant is the high agglomeration region (H-H agglomeration), which is used to represent the relationship of high efficiency provinces surrounded by high efficiency provinces; the third quadrant is the low agglomeration region (L-L agglomeration), which indicates that low efficiency provinces are surrounded by low efficiency provinces; the second quadrant (L-H agglomeration) is used to represent the relationship of low efficiency provinces surrounded by high efficiency provinces; and the fourth quadrant (H-L agglomeration) is used to represent the relationship of high efficiency provinces surrounded by low efficiency provinces The fourth quadrant (H-L aggregation) is used to represent the relationship of high-efficiency provinces surrounded by low-efficiency provinces.
Medical And Health Input-output In China
From 2011 to 2020, China's overall investment in PHC increased significantly (Table 2). In terms of the number of health technicians, as China continues to strengthen the construction of health care personnel, the number of health technicians has increased significantly, with Guangdong and Jiangsu having the largest increases, from 210,938 and 192,143 in 2011 to 307,240 and 285,750 in 2020, respectively, an increase of 45.65% and 48.72%. From the point of view of the number of PHC institution, there are nine regions for negative growth, other regions for growth, including Shandong, the largest increase in the number of 15,172. In terms of the number of beds, Hainan had the largest increase (64.12%), followed by Guangxi and Jiangxi, up 63.87% and 63.30% respectively. In terms of the output of PHC institution service efficiency, Heilongjiang and Gansu saw the most significant decline in the number of visits, with growth rates of -42.17% and − 23.88%, respectively. In terms of the number of hospital discharges, most regions showed a decline, with Tianjin showing a greater decline (-81.82%). Overall, most Chinese provinces have continued to increase their investment in PHC in the last decade, which is conducive to improving the accessibility of healthcare services. However, the output of PHC services in some provinces did not improve accordingly, especially in the number of hospital discharges, which showed a decreasing trend in most regions, and the efficiency of PHC resources utilization needs to be further studied.
Table 2
Comparison of input and output indicators of primary health care in 31 provinces (municipalities directly under the Central Government and autonomous regions) in China
Region
|
2011
|
2020
|
X1
|
X2
|
X3
|
Y1
|
Y2
|
X1
|
X2
|
X3
|
Y1
|
Y2
|
Beijing
|
8 718
|
4 423
|
54 468
|
5 259
|
4
|
9 675
|
5 145
|
86 400
|
6 919
|
1
|
Tianjin
|
3 981
|
6 883
|
21 690
|
3 253
|
11
|
5 258
|
6 007
|
36 689
|
4 043
|
2
|
Hebei
|
78 246
|
67 942
|
190 213
|
24 485
|
166
|
83 972
|
78 970
|
225 307
|
21 778
|
112
|
Shanxi
|
38 587
|
40 348
|
108 266
|
6 975
|
61
|
39 224
|
37 414
|
110 820
|
5 634
|
33
|
Neimenggu
|
21 905
|
23 992
|
66 747
|
5 364
|
45
|
23 277
|
26 769
|
78 094
|
4 311
|
24
|
Liaoning
|
33 712
|
34 236
|
97 272
|
8 110
|
66
|
32 174
|
38 346
|
103 145
|
6 891
|
34
|
Jilin
|
18 882
|
21 040
|
65 190
|
4 974
|
31
|
24 424
|
19 627
|
86 059
|
4 242
|
11
|
Heilongjiang
|
20 142
|
27 081
|
83 111
|
5 495
|
53
|
18 653
|
31 168
|
81 040
|
3 178
|
27
|
Shanghai
|
4 289
|
17 955
|
46 122
|
9 266
|
11
|
5 291
|
15 612
|
73 182
|
8 952
|
4
|
Jiangsu
|
29 659
|
68 097
|
192 143
|
23 055
|
169
|
32 703
|
101 630
|
285 750
|
28 069
|
222
|
Zhejiang
|
29 207
|
23 967
|
121 834
|
20 990
|
28
|
32 376
|
30 718
|
196 621
|
32 374
|
33
|
Anhui
|
21 434
|
56 291
|
133 737
|
13 525
|
143
|
27 400
|
79 745
|
169 973
|
21 689
|
117
|
Fujian
|
26 287
|
27 626
|
87 868
|
9 437
|
112
|
26 902
|
37 278
|
123 695
|
13 552
|
62
|
Jiangxi
|
38 063
|
38 211
|
111 183
|
12 539
|
218
|
35 214
|
62 398
|
122 870
|
13 322
|
169
|
Shandong
|
65 954
|
112 951
|
313 328
|
36 813
|
303
|
81 126
|
118 285
|
356 927
|
36 748
|
250
|
Henan
|
74 208
|
94 056
|
282 443
|
33 981
|
310
|
71 339
|
137 822
|
311 634
|
35 090
|
316
|
Hubei
|
34 509
|
62 330
|
154 855
|
17 444
|
192
|
33 852
|
98 093
|
178 432
|
16 292
|
265
|
Hunan
|
58 214
|
77 251
|
166 550
|
14 537
|
316
|
53 793
|
124 217
|
204 248
|
14 987
|
405
|
Guangdong
|
44 034
|
59 128
|
210 938
|
32 462
|
196
|
53 069
|
73 595
|
307 240
|
34 827
|
184
|
Guangxi
|
33 132
|
46 534
|
121 940
|
12 793
|
216
|
32 149
|
76 256
|
171 149
|
11 283
|
286
|
Hainan
|
4 513
|
5 831
|
18 626
|
2 241
|
11
|
5 733
|
9 570
|
31 009
|
3 079
|
6
|
Chongqing
|
17 037
|
37 196
|
77 760
|
8 264
|
145
|
19 838
|
55 901
|
102 150
|
8 995
|
199
|
Sichuan
|
73 646
|
113 909
|
228 917
|
26 802
|
436
|
79 491
|
149 620
|
285 911
|
29 333
|
458
|
Guizhou
|
24 957
|
34 256
|
73 043
|
7 663
|
185
|
27 138
|
54 113
|
124 823
|
8 459
|
130
|
Yunnan
|
21 800
|
40 488
|
84 633
|
11 178
|
121
|
24 592
|
61 932
|
154 144
|
14 868
|
162
|
Xizang
|
6 356
|
2 861
|
13 051
|
628
|
3
|
6 632
|
3 867
|
20 485
|
930
|
1
|
Shannxi
|
35 033
|
31 260
|
102 082
|
9 021
|
62
|
33 208
|
39 906
|
130 527
|
8 287
|
54
|
Gansu
|
25 884
|
25 154
|
67 831
|
7 680
|
59
|
24 597
|
32 896
|
74 988
|
5 846
|
72
|
Qinghai
|
5 608
|
4 325
|
15 610
|
1 172
|
16
|
6 018
|
5 398
|
20 138
|
1 072
|
9
|
Ningxia
|
3 886
|
2 781
|
11 660
|
1 458
|
7
|
4 247
|
4 094
|
19 915
|
1 793
|
4
|
Xinjiang
|
16 120
|
25 315
|
51 882
|
3 696
|
80
|
16 671
|
32 992
|
66 380
|
4 771
|
56
|
Data analysis
The study used Excel 2019 software to summarize the data; DEA- Malmquist index method was applied to the initial data by using DEAP 2.1 software to obtain the results of service efficiency for the years 2011–2020; MATLAB 2021B was used to conduct DEA-SE-SBM programming modeling and process the original data[51], to obtain the PHC service efficiency values of each region in the past ten years; The GeoDa 1.18.0 was used to establish the adjacency Rook matrix and bring in the service efficiency values of independent year regions to calculate the global Moran index and local Moran index to determine whether they have spatial autocorrelation. Finally, ArcGIS 10.8 software was used to visualize the data.