Skip to main content

Assessing frontline HIV service provider efficiency using data envelopment analysis: a case study of Philippine social hygiene clinics (SHCs)



Globally, local and frontline HIV service delivery units have been deployed to halt the HIV epidemic. However, with the limited resources, there is a need to understand how these units can deliver their optimum outputs/outcomes efficiently given the inputs. This study aims to determine the efficiency of the social hygiene clinics (SHC) in the Philippines as well as to determine the association of the meta-predictor to the efficiencies.


In determining efficiency, we used the variables from two data sources namely the 2012 Philippine HIV Costing study and 2011 Integrated HIV Behavioral and Serologic Surveillance, as inputs and outputs, respectively. Various data management protocols and initial assumptions in data matching, imputation and variable selection, were used to create the final dataset with 9 SHCs. We used data envelopment analysis (DEA) to analyse the efficiency, while variations in efficiencies were analysed using Tobit regression with area-specific meta-predictors.


There were potentially inefficient use of limited resources among sampled SHC in both aggregate and key populations. Tobit regression results indicated that income was positively associated with efficiency, while HIV prevalence was negatively associated with the efficiency variations among the SHCs.


We were able to determine the inefficiently performing SHCs in the Philippines. Though currently inefficient, these SHCs may adjust their inputs and outputs to become efficient in the future. While there were indications of income and HIV prevalence to be associated with the efficiency variations, the results of this case study may only be limited in generalisability, thus further studies are warranted.

Peer Review reports


HIV/AIDS is one of the gravest public health issues in the world. According to the 2014 Global Statistics [1], 1.2 million people died of AIDS-related illnesses, with about 36.9 million people infected with HIV, of which 2 million accounts for the annual newly-infected cases. The global response to date, through the prevention and treatment services for sufferers and population at risk, has succeeded in decreasing the number of People Living with HIV (PLHIV) steadily [2, 3].

In recent years, country-led frontline service delivery programs/interventions have been carried out to halt the epidemic across the globe [4,5,6]. These programs/interventions, focusing on better access to basic prevention, care and treatment (PCT) services, have catered to both the PLHIV population and key populations such as, male having sex with male (MSM), female sex workers (FSW) and people who inject drugs (PWID). Each country devised their own intervention as a response to the epidemic. In Ghana, TB integration into the HIV services through “one-stop shops” were found to have improved HIV screening [7]. While in Asia, Hong Kong and the Philippines have social hygiene clinics (SHC), which provide basic PCT services to key populations, particularly for FSW [8].

With the increasing number of HIV/AIDS-related frontline health service providers, there is a growing need to assess whether these providers have been efficiently delivering the health services subject to the available resources. In literature, effectiveness/outcomes and costing studies of the frontline HIV/AIDS services have been extensively studied independently [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. However, efficiency studies, which account for both effectiveness/outcomes and costs at the same time, were relatively limited [10, 28,29,30,31,32,33,34]. Efficiency studies provide information about the efficient/inefficient use of finite resources, which could be subsequently utilised as potential suggestions in the improvement of both management and policy. As a tool to analyse efficiency, data envelopment analysis (DEA) has been used most often, since it is non-parametric, which imposes no restrictive hypothesis on the data generating process and requires minimal assumptions about the technology [35, 36]. Although DEA has been widely used in different disciplines [37,38,39,40,41,42], its utilisation in HIV/AIDS research is limited [28,29,30]. HIV/AIDS research using DEA has been instrumental to various areas across the globe as a means for resource allocation and service delivery optimisation efforts. Health centers in Rwanda were performing at an average efficiency at 78%, which suggests that there are still room for improvement in service delivery [28]. Inclusion of efficiency determination in national HIV/AIDS programmes in resources-needs estimation, using macro-level data, can narrow the spending gap in the HIV/AIDS response [29]. Lepine, et al. [35] has observed that Avahan NGOs in India could have reduced their inputs by 43% given the outputs they reached. In a more comprehensive analysis by Zanakis, et al. [43], countries with lower population density that manage to provide better health system performance and per capita support with lower GNP and better media information, tend to have lower HIV/AIDS indicators. In Kenya, Omondi Aduda, et al. [44] has observed that significant improvement in the technical efficiency among the outreach facilities were attributable to the voluntary medical male circumcision using DEA.

Most of the HIV-related DEA studies have focused on the efficiency of the delivery of health services to the general population, with little or no focus on the key populations. Key populations such as the MSM, FSW and PWID also access the same basic health services but are at higher risk of acquiring sexually transmitted diseases. Their access to these health services may vary from the general population, hence, determining the efficiency of service delivery for these populations would provide insightful suggestions to current service provider.

To the best of our knowledge, this is the first Philippine case study which aims to determine the less efficient frontline service providers, SHCs, among the aggregate and key populations, as well as to identify its possible determinants of efficiency, as a mean to improve the SHC’s operational practice compared to the efficient ones in search for better management.


Philippine HIV situation

As of July 2015, the Philippines has recorded a total of 27,138 HIV/AIDS cases, which has been increasing annually [45]. The country responded to the worsening HIV/AIDS epidemic, by enabling grassroots level primary health care service providers, such as the Rural Health Units (RHU) and City Health Office (CHO), to function as SHCs in providing frontline STI PCT services [46]. In the 1980s, there were roughly at total of 170 SHCs, however, the number went down to 70 SHCs due to the decrease in the accreditation of the RHU/CHO [47].

Data envelopment analysis

In this study, DEA was used to evaluate the efficiency of the SHCs. DEA benchmarks the performance of decision-making units (DMU) such as healthcare facilities, like the SHC, taking into consideration the inputs and outputs at the same time [43, 48]. A set of efficient DMUs were identified to form the best-practice frontier, and efficiency measures were estimated relative to the remaining DMUs. After constructing an efficient frontier line, DMUs which lie on the line were considered to be efficient, while those not on the line were inefficient [49]. By benchmarking to the best performers, this will help the inefficient DMUs to improve their functional organisation based on the best performing DMUs in the frontier line [50]. DEA imposes less assumptions on the functional shape of the relationship of the outputs and the inputs [51, 52], and can handle multiple inputs and outputs at the same time. Efficiency in DEA was modeled as the maximum ratio of weighted outputs (yrk) to weighted inputs (xik) subject to the similar ratios for every DMU be less than or equal to unity. This can be summarised in a mathematical formulation found below [36, 53,54,55] (Eq. 1) [56]:

$$ {\displaystyle \begin{array}{l}\mathit{\operatorname{Max}}\ {\theta}_k=\frac{\sum_{r=1}^S{u}_r{y}_{rk}}{\sum_{i=1}^m{v}_i{x}_{ik}}\\ {} Subject\ to\ \frac{\sum_{r=1}^S{u}_r{y}_{rk}}{\sum_{i=1}^m{v}_i{x}_{ik}}\le 1,\\ {}{u}_r\ge \varepsilon, {v}_i\ge \varepsilon \forall r,i\end{array}} $$


k ≡ SHCs;

θ ≡ efficiency coefficient;

vi, i = 1, 2, …, m, are the weights assigned to the i-th inputs;

ur, r = 1, 2, …, m, are the weights assigned to the r-th outputs

DEA was used to facilitate the input-output modeling in determining the efficiency frontier, returns to scale and inefficient-to-efficient movement in the first stage analysis. DEA explores how output variables interacted with input variables; an expansion of the basic principle of production function taking into account the multi-output, multi-input modeling.

Due to DEA’s feature of multiple input-output analysis, it has been widely used in most operational research analyses for efficiency determination [36, 57, 58]. Recent advances in DEA have further emphasised in identifying the possible determinants of efficiency by using area-specific, meta-predictors [29, 30, 59]. The efficiency coefficients in the first stage were then regressed with the area-specific parameters using Tobit regression in determining which meta-predictors can explain the variations among the efficiencies.

Inputs and outputs

The inputs and outputs used in this study were from secondary data sources, which are highlighted in the Additional file 1. Inputs can be in the form of unit costs or physical units [60], however, due to data availability, we were only able to acquire the latest input data of health service unit costs from the 2012 Philippine HIV Costing Study [61]. Output variables, i.e. number of people accessing a specific health service, were taken from the 2011 Integrated HIV Behavioral and Serologic Surveillance (IHBSS). The input and output variables were then managed using the robust protocol; as shown in Fig. 1. Both the 13 SHCs identified from the Philippine HIV Costing Study 2012 and 18 SHCs from the 2011 IHBSS were sampled from the remaining 70 functional SHCs across the country. After matching the input and output variables at the SHC level, we were left with 13 SHCs, having five inputs and five outputs each.

Fig. 1
figure 1

Schematic diagram of Data Processing and Analyses. Shows the flow of the robust protocols used from data matching until Tobit regression

Fig. 2
figure 2

Location of the nine study SHCs. Geographical location of the operational non-study SHCs (blue dots), and the study SHCs (red dots) across the Philippines. Fig. 2 was generated using R statistical programming [86] software through the "ggmaps" package

In total there were 10 variables, with five pairs of input-output correspondence. However, DEA has been observed to have varying discriminatory power with respect to the DMUs and the input-output proportion [62]. Considerable research on the field has been done to determine the ideal proportion of DMU taking into consideration its inputs and outputs. A research by Doyle and Green [63] has noted that the number of DMUs should be at least twice that of the combined total number of inputs and outputs in order for DEA to discriminate among the efficiencies. In this study, we utilised principal component analysis (PCA) in facilitating variable dimension reduction to determine the optimum combination of inputs and outputs [64]. Based on previous studies [64, 65], PCA was noted to work well with small number of samples (n < 25), which is ideal for the dataset having only 9 observation data points (9 SHCs).

We used a Matlab program, PCA-DEA, by Adler and Yazhemsky [64], to perform the variable dimension reduction. In principle, the PCA replaces the original inputs and outputs with concise groups of principal components that explain the variance structure of a matrix of data through linear combinations of variables [66]. If most of the population variance can be attributed to the first few components, then those principal components can replace the original variables with minimal loss of information [67]. From the 5-input, 5-output variable in the initial dataset, the dimension reduction favoured the 3-input, 2-output combination, which was able to explain 76–80% of the variations in the component analysis on both input and output sides. Adler and Yazhemsky [64] notes that dimensionality reduction should at least be able to explain 76% [for variable returns to scale (VRS) specification] and 80% [for constant returns to scale (CRS) specification] the variation to provide good approximation of the efficiency classification.

From the 13 input-output matched SHCs, 9 SHCs were completely matched; after excluding SHCs which have less than 80% complete data. The 9 SHCs at the aggregate level were further stratified, to determine the variations between the efficiency scores if they were to be disaggregated into key populations. Although the 9 SHCs have constraints with respect to the generalisability of the results, we believe that the dataset we built based on the robust assumptions from the established secondary data sources, can provide intuitive results for well-informed policymaking for micro-level health facility managers in a resource and data-limited setting worldwide.

Area-specific, meta-predictors

There are numerous possible candidate meta-predictors present in the secondary output data source, which may explain the variations. However, after initial diagnostics, the 55 candidate meta-predictors, which included but are not limited to education, civil status, clients, oral sex, frequency of condom use, and others, were found to have issues related to high multicollinearity, hence, we chose not to include them in the analysis. Instead, we used income and HIV prevalence, through the AMTP 5 categorisation, as meta-predictors. The AMTP 5 classification signifies the priority level of investment for a specific area, depending on how big such area contributes towards the epidemic [68]. Cities with “A” were coded with 1, while those with “B” are 0.

SHC field location

The 9 SHCs identified in this study are located across the country with four in the national capital region. Most of the SHCs are located in highly urbanised areas of the country wherein the HIV epidemic is concentrated [69].

Bootstrapping and Tobit regression

Even though DEA is computationally flexible, since there is no prior structural function imposed [70, 71], important issues, namely DEA technical efficiency scores being sensitive to sampling variation and that the efficiency estimates being serially correlated due to the small sample size [70, 72], can make DEA prone to biased results and overestimated efficiency scores; thereby erroneously classifying DMUs into either efficient or inefficient [70, 71, 73]. To overcome this issue, we used bootstrapping to determine whether the ranks, based on the efficiency scores, of the SHCs would vary or not.

Simar and Wilson [74] proposed a bootstrapping methodology for analysing the sampling variation and estimate the confidence intervals of the radial DEA measures (\( \hat{\theta} \)). Utilising homogeneity bootstrapping approach, it assumes that there is an underlying data generation process which generates the data set (x, y) from the production possibility set [64].

Bootstrapping approximates the unknown distribution using the difference between the original efficiency estimates and ‘true’ efficiency, through the distribution of the difference between the bootstrapped efficiency estimates and the original efficiency estimator, conditioned on the original data [64]. Further technical specifications of DEA-related bootstrapping are discussed in full elsewhere [64, 74, 75].

We further examined the variations with the efficiencies by using Tobit regression (in Eq. 2) to identify the possible determinants of efficiency using various area-specific meta-predictors.

$$ {\displaystyle \begin{array}{l}\mathrm{y}=\left\{\begin{array}{c}{y}^{\hbox{'}};0\le {y}^{\hbox{'}}\le 1\\ {}0;{y}^{\hbox{'}}<0\ \\ {}1;1<{y}^{\hbox{'}}\ \end{array}\right.\\ {}{\mathrm{y}}^{\hbox{'}}=\beta {x}_i+\varepsilon \end{array}} $$

Whereby: y is the efficiency score; y as the latent (unobservable) variable; β vector of parameters; xi are the meta-predictors.

In this study, both input and output variables are secondary data sources, with no individual information, hence ethical consent is not applicable. All analyses were done using R programming and STATA 13.


Table 1 shows the summary of the inputs and outputs as well as the area-specific meta-predictors. In the aggregate population, SHCs incur the greatest cost with the HIV test 380 (± 173), which constitutes the largest percentage (84%) of input costs. While for the output side, condom access and access to gram staining services were of the same percentage (at 50%). The wide standard deviation as observed in the input side with nearly 50% of the mean, is indicative of the variations among the cost of goods/services.

Table 1 Summary table of the inputs, outputs and area-specific meta-predictors

Table 2 shows both raw and bias-corrected technical efficiencies derived from DEA and bootstrapping, respectively. At the aggregate level, 3 SHCs were found to be inefficient, with Pasay as the lowest with a technical efficiency score of 0.54. When disaggregated into the key populations, MSM was observed to have the least average technical efficiency score at 0.81, which seconded FFSW (at 0.84) and has still more means to improve compared to the perfect average efficiency score of RFSW (at 1.00). Though there is an overestimation of the efficiencies, with values of the bias-corrected DEA efficiencies lower than the raw technical efficiencies, the ranking of the SHCs with respect to the efficiency scores are still the same, in both the aggregate and key populations. Consistent and evident low efficiency can be observed in Pasay SHC.

Table 2 Raw and bootstrapped DEA efficiencies in different key populations

Table 3 shows the input and output adjustments needed for the SHC to operate on the efficient frontier line. At the aggregate level, Baguio has much to improve in relation to the reduction of its input costs for gram stain (− 0.50 Php), condom (− 4.90 Php) and particularly HIV testing (− 465 Php), and likewise the increment the number of people accessing its condom services (+ 139 pax). Across the key populations, there is no clear pattern of reductions or increments among the SHCs, however, Baguio and Pasay SHCs have been observed to need more adjustments compared to the remaining SHCs. On the other hand, from among the 9 SHCs, Davao and Quezon City SHCs don’t need any change in their input or output components, which may be partially indicative, yet inconclusive signs of good management performance.

Table 3 Input reduction and output increment of the SHCs

Income and HIV prevalence, through AMTP 5 indicator, were used as the area-specific meta-predictors, instead of the output dataset’s candidate meta-predictors. HIV prevalence was observed to be negatively associated with the efficiency scores, with notable significant inverse relationship observed at the aggregate level (at − 0.60). On the other hand, income is positively associated to the efficiency scores, with a significant association observed in the FFSW (at 2.42e− 4). We have observed a positive association with income and negative association with HIV prevalence, as shown in Table 4. RFSW results were not shown due to the homogeneity of the efficiency scores, which can’t be analysed using Tobit regression.

Table 4 Tobit regression of the efficiencies with area-specific meta-predictors


The frontline service delivery is an effective response to the HIV/AIDS epidemics [4, 5]. In the Philippines, SHCs operate as one of the major frontline health service providers. In this case study, we have successfully identified the potentially inefficient use of finite resources among sampled SHCs. As exemplified by Pasay SHC, DEA was able to objectively identify its inefficiency among the SHCs, and even across the different key populations. Likewise, further analyses showed that both income and HIV prevalence were statistically related to efficiency.

Efficiency is an important index for benchmarking a DMU’s performance. Performance, according to ranking, based on either input or output variables (Table 1), would only provide limited information of the DMUs performance. Through DEA, inefficiency can be detected and be resolved by adjusting either the input or output side. At the aggregate population, in Table 2, an example of Pasay SHC, with the lowest efficiency score even after bootstrapping, can overcome this inefficiency by following the adjustments/suggestions of DEA, such as finding cheaper unit cost of inputs (e.g. 12.60 Php cheaper gram stain service and 58.30 Php cheaper HIV testing kit as shown in Table 3). Pasay SHC’s inefficiency can be attributed to the higher priced inputs in relation to the outputs (services provided) (shown in Table 3), which may be reflective of the bargaining transactions of the local government unit (LGU). The devolved functions of the government system in the Philippines is also apparent in health, particularly in the procurement of resources. LGUs have jurisdiction over the procurement of resources in their respective areas. While the procurement process is a standard protocol across LGUs, the bargaining transactions may vary, thus resulting to varied procurement prices. In achieving optimum efficiency, aside from Pasay, other inefficient SHCs in both the aggregate and key populations, such as Baguio, Cagayan de Oro, Cebu and others, can also benefit from DEA’s input reduction and output increment features (as shown in Table 3). The potential adjustments in the inputs/outputs provide an indication of future directions regarding SHC management. In Table 3, inputs can be reduced, or outputs can be increased per specific units to achieve optimum efficiency. While these DMUs are currently operating away from the efficient frontier, there are opportunities where resource use (input) and service delivery (output) can be optimised to approach the frontier, and prospectively become efficient. Specifically, the temporal changes in technology, in terms of the operational management of the respective inputs and outputs, provide an opportunity for DMUs to optimise their efficiency. In Nepal, Silwal and Ashton [76] notes the potential for hospitals to increase their outputs given the additional resources through the years. While in Greece, Polyzos [77] observed that middle-sized hospitals were able to achieve 100% technical efficiency in the span of 2009–2011 as a result of technological changes (i.e. spending cuts and constant reforms).

It should be noted that although we assume SHCs as DMUs, Philippine SHCs don’t decide independently, instead, they are under the mandate of the LGUs. In most cases, bargaining of goods/services are dependent on the LGUs, rather than the SHCs. The complicated functionality of the SHC-LGU nexus in relation to decision-making and bargaining of goods, is further affected by the uncertainty introduced by the local market. Market price dictates unit cost, which then affects the DMUs’ decision. These interrelated complexities in the functional and structural relation of SHCs, LGUs and local markets were not accounted, and brings us to one of the limitations of the study, which assumes that SHCs are DMUs, and that unit costs of goods/services from local markets have low variance in a short time horizon, thus such complexities were assumed to be null in this study.

Studies conducted in the field of HIV service facility efficiency assessment focused more on the aggregate population [28,29,30], which assumes a homogeneous population rather than the disaggregated, heterogeneous populations. Santos, et al. [48] notes that aggregation of data might undermine the performance of other DMUs providing the same services. In the same note, Gori, et al. [78] observed that disaggregated data provides better estimation of the effects than aggregated models. Disaggregation of service provision can provide better insights regarding the utilisation of resources per service, compared to aggregate ones [79]. This scenario was also observed in this case study, wherein among the key populations, efficiency estimates were varied, with some estimates consistently lower than the aggregate population; a case exemplified by the Pasay SHC. This suggests that Pasay SHC, should not only focus at the aggregate level, but most importantly at the disaggregated key populations in addressing the issue of inefficiency. We have also observed that SHCs performed better in providing services to the RFSW compared to the other two key populations. The difference in the performance (in terms of the efficiencies) may be related to the nature of the access by key populations. RFSW as the term implies are registered under the SHCs, whereby they are required to access the services in a timely manner. MSMs and FFSWs, on the other hand, are populations which access the services in a discretionary manner.

After bootstrapping the DEA efficiency scores (as shown in Table 2), we have determined that though there is consistent overestimation of results in the raw efficiencies, the ranking of the efficiencies even after being bias-corrected are relatively the same, which only shows that the results of this study are robust. We further examined the variations among the efficiencies by determining which exogenous area-specific, meta-predictors may be able to explain the difference among the efficiency scores. Both income and HIV prevalence (through the AMTP classification), were observed to be related with the efficiency scores. As shown in Table 4, income is positively associated with efficiency, an increase in the income of the local area may contribute to the increase in the efficiency of the SHCs. Similarly, if there is a larger proportion of income, this may translate to higher allocation for the specific programs, such that of HIV PCT [80, 81]. The devolved system of government of the country gives LGU exclusive allocative powers with respect to the management of its territorial resources as well as its own transactions. The crude assumption of greater income equates greater allocation is anchored on the premise of prioritisation. We can assume that everything has a proportional increase or decrease based on income and the level of prioritisation.

On the other hand, high prevalence areas have been significantly associated to less efficiency, in the aggregate population, as shown in Table 4. The same trend was observed in the statistically non-significant coefficients of MSM and FFSW. This inverse association is an eventual paradox to the usual understanding of the economies of scale. We believe that this paradox might have been due to two following limitations of the data sources: 1) the nature of the data, leaning towards a capital-intensive perspective, and 2) non-inclusion of treatment unit costs due to lack of complete data across SHCs. The capital-intensive approach did not include the number of health personnel in the SHCs. Highly prevalent areas may have greater number of service personnel, but since the available data included only the service delivery time and the salary wage per health personnel, labour, such as actual number of health service providers, were not accounted. Likewise, after thorough data management, most treatment variables 1in the input data source were not included in the study, as only few SHCs have sufficient unit cost data for treatment services. This may be attributable to the paradox, whereby highly prevalent areas may provide more treatment services, but the model was only able to capture the prevention and control intervention components.

With the various data management techniques used in this study, caution should be regarded with the limitations in the 1) specifications of the protocols, 2) generalisability of the results, and 3) data availability. The protocols in data generation, data matching, data imputation, and variable selection are based on the published literature and the observed plausibility for the analysis to be carried out. Since this study is a case study of the Philippine SHCs, the number of SHCs used in this study compared to the 70 SHCs all over the country may prove that the results acquired can’t be fully generalised with other units. We also acknowledge the possibility that the locations whereby the data are available may already be performing better than the rest of the SHCs with no available data, which, in turn may bias the results towards efficient performers. Though this may be the case with the current data source, there is also another possibility whereby the SHCs with no data may perform better than the study SHCs. While this possibility is plausible, we can’t ascertain the validity of this assumption especially with the limited data. As mentioned earlier, some results of this study are constrained by the data availability of actual number of health service personnel and the complete unit cost data of the SHCs. The inclusion of these variables into DEA warrants further investigation.


We were able to determine the inefficiently performing SHCs in the Philippines. Though currently inefficient, these SHCs may adjust their inputs and outputs to become efficient in the future. While there were indications of income and HIV prevalence to be associated with the efficiency variations, the results of this case study may only be limited in generalisability, thus further studies are warranted.

Operational definitions

  • Meta-predictor – variables which include income and HIV prevalence, which serve as area-specific characteristics.

  • Registered Female Sex Workers (RFSW) – Born female, 15 years or older, who accepted payment (cash or kind) in exchange for sex in the past 1 month and is based in an entertainment establishment registered at the local social hygiene clinic [82].

  • Freelance Female Sex Workers (FFSW) – Born female, 15 years or older, who accepted payment (cash or kind) in exchange for sex in the past 1 month and is street-based or based in an entertainment establishment not registered at the local social hygiene clinic. This includes those found in all kinds of cruising sites, and those who practice indirect sex work [82].

  • Males who have Sex with Males (MSM) – Born male, 15 years or older, who reported oral or anal sex with another male in the past 12 months [82].

  • (Best practice) frontier – represents the maximum output that a DMU can produce from all inputs combined. DMUs on the best practice frontier are considered to be technically efficient, while those distant ones are otherwise considered to be inefficiently performing [83].

  • Local Government Unit (LGU) – “are institutional units whose fiscal, legislative and executive authority extends over the smallest geographical areas distinguished for administrative and political purposes” [84].

  • Constant return to scale (CRS) – assumes that the average productivity, in terms of the ratio of output (y)/input (x) is not dependent on the scale of production [85].

  • Variable return to scale (VRS) – assumes a constant but also increasing or decreasing returns to scale at different scale sizes [85].



Data envelopment analysis


Decision making unit


Freelance Female Sex Workers


Female sex workers


Local government unit


Male having sex with males


Prevention, care and treatment


Registered Female Sex Workers


Social hygiene clinic


  1. UNAIDS. How AIDS changed everything. Geneva: Joint United Nations Programme on HIV/AIDS; 2014.

  2. Fettig J, Swaminathan M, Murrill CS, Kaplan JE. Global epidemiology of HIV. Infect Dis Clin North Am. 2014;28(3):323–37.

    Article  PubMed  PubMed Central  Google Scholar 

  3. UNAIDS. Global report: UNAIDS report on the global AIDS epidemic 2013. Geneva: Joint United Nations Programme on HIV/AIDS; 2013.

  4. Gupta GR, Parkhurst JO, Ogden JA, Aggleton P, Mahal A. Structural approaches to HIV prevention. Lancet. 2008;372(9640):764–75.

    Article  PubMed  Google Scholar 

  5. Bertozzi SM, Laga M, Bautista-Arredondo S, Coutinho A. Making HIV prevention programmes work. Lancet. 2008;372(9641):831–44.

    Article  PubMed  Google Scholar 

  6. De Cock KM, Fowler MG, Mercier E, de Vincenzi I, Saba J, Hoff E, Alnwick DJ, Rogers M, Shaffer N. Prevention of mother-to-child HIV transmission in resource-poor countries: translating research into policy and practice. JAMA. 2000;283(9):1175–82.

    Article  PubMed  Google Scholar 

  7. Ansa GA, Walley JD, Siddiqi K, Wei XL. Delivering TB/HIV services in Ghana: a comparative study of service delivery models. Trans R Soc Trop Med Hyg. 2014;108(9):560–7.

    Article  CAS  PubMed  Google Scholar 

  8. Morisky DE, Ang A, Coly A, Tiglao TV. A model HIV/AIDS risk reduction programme in the Philippines: a comprehensive community-based approach through participatory action research. Health Promot Int. 2004;19(1):69–76.

    Article  PubMed  Google Scholar 

  9. Zulliger R, Black S, Holtgrave DR, Ciaranello AL, Bekker LG, Myer L. Cost-effectiveness of a package of interventions for expedited antiretroviral therapy initiation during pregnancy in Cape Town, South Africa. AIDS Behav. 2014;18(4):697–705.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Holtgrave DR, Qualls NL, Curran JW, Valdiserri RO, Guinan ME, Parra WC. An overview of the effectiveness and efficiency of HIV prevention programs. Public Health Rep. 1995;110(2):134–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Shade SB, Kevany S, Onono M, Ochieng G, Steinfeld RL, Grossman D, Newmann SJ, Blat C, Bukusi EA, Cohen CR. Cost, cost-efficiency and cost-effectiveness of integrated family planning and HIV services. AIDS. 2013;27(Suppl 1):S87–92.

    Article  PubMed  Google Scholar 

  12. Larsson EC, Ekstrom AM, Pariyo G, Tomson G, Sarowar M, Baluka R, Galiwango E, Thorson AE. Prevention of mother-to-child transmission of HIV in rural Uganda: modelling effectiveness and impact of scaling-up PMTCT services. Glob Health Action. 2015;8:26308.

    Article  PubMed  Google Scholar 

  13. Zang X, Tang H, Min JE, Gu D, Montaner JS, Wu Z, Nosyk B. Cost-effectiveness of the ‘One4All’ HIV linkage intervention in Guangxi Zhuang autonomous region, China. PLoS One. 2016;11(11):e0167308.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Sharma M, Farquhar C, Ying R, Krakowiak D, Kinuthia J, Osoti A, Asila V, Gone M, Mark J, Barnabas RV. Modeling the cost-effectiveness of home-based HIV testing and education (HOPE) for pregnant women and their male Partners in Nyanza Province, Kenya. J Acquir Immune Defic Syndr. 2016;72(Suppl 2):S174–80.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Sanchez-Gonzalez G. The cost-effectiveness of treating triple coinfection with HIV, tuberculosis and hepatitis C virus. HIV Med. 2016;17(9):674–82.

    Article  CAS  PubMed  Google Scholar 

  16. Price JT, Wheeler SB, Stranix-Chibanda L, Hosek SG, Watts DH, Siberry GK, Spiegel HM, Stringer JS, Chi BH. Cost-effectiveness of pre-exposure HIV prophylaxis during pregnancy and breastfeeding in sub-Saharan Africa. J Acquir Immune Defic Syndr. 2016;72(Suppl 2):S145–53.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Pham QD, Wilson DP, Nguyen TV, Do NT, Truong LX, Nguyen LT, Zhang L. Projecting the epidemiological effect, cost-effectiveness and transmission of HIV drug resistance in Vietnam associated with viral load monitoring strategies. J Antimicrob Chemother. 2016;71(5):1367–79.

    Article  PubMed  Google Scholar 

  18. Oddershede L, Walker S, Stohr W, Dunn DT, Arenas-Pinto A, Paton NI, Sculpher M. Protease inhibitor monotherapy versus ongoing triple therapy trial T. cost effectiveness of protease inhibitor monotherapy versus standard triple therapy in the long-term management of HIV patients: analysis using evidence from the PIVOT trial. Pharmacoeconomics. 2016;34(8):795–804.

    Article  PubMed  Google Scholar 

  19. Moodley N, Gray G, Bertram M. The case for adolescent HIV vaccination in South Africa: a cost-effectiveness analysis. Medicine (Baltimore). 2016;95(4):e2528.

    Article  CAS  Google Scholar 

  20. Lin F, Farnham PG, Shrestha RK, Mermin J, Sansom SL. Cost effectiveness of HIV prevention interventions in the U.S. Am J Prev Med. 2016;50(6):699–708.

    Article  PubMed  Google Scholar 

  21. Kapoor S, Gupta A, Shah M. Cost-effectiveness of isoniazid preventive therapy for HIV-infected pregnant women in India. Int J Tuberc Lung Dis. 2016;20(1):85–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Jacobsen MM, Walensky RP. Modeling and cost-effectiveness in HIV prevention. Curr HIV/AIDS Rep. 2016;13(1):64–75.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Hutchinson AB, Farnham PG, Sansom SL, Yaylali E, Mermin JH. Cost-effectiveness of frequent HIV testing of high-risk populations in the United States. J Acquir Immune Defic Syndr. 2016;71(3):323–30.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Hontelez JA, Chang AY, Ogbuoji O, de Vlas SJ, Barnighausen T, Atun R. Changing HIV treatment eligibility under health system constraints in sub-Saharan Africa: investment needs, population health gains, and cost-effectiveness. AIDS. 2016;30(15):2341–50.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Hoenigl M, Chaillon A, Mehta SR, Smith DM, Graff-Zivin J, Little SJ. Screening for acute HIV infection in community-based settings: cost-effectiveness and impact on transmissions. J Inf Secur. 2016;73(5):476–84.

    Google Scholar 

  26. Heffernan A, Barber E, Thomas R, Fraser C, Pickles M, Cori A. Impact and cost-effectiveness of point-of-care CD4 testing on the HIV epidemic in South Africa. PLoS One. 2016;11(7):e0158303.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Bristow CC, Larson E, Anderson LJ, Klausner JD. Cost-effectiveness of HIV and syphilis antenatal screening: a modelling study. Sex Transm Infect. 2016;92(5):340–6.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zeng W, Rwiyereka AK, Amico PR, Avila-Figueroa C, Shepard DS. Efficiency of HIV/AIDS health centers and effect of community-based health insurance and performance-based financing on HIV/AIDS service delivery in Rwanda. Am J Trop Med Hyg. 2014;90(4):740–6.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zeng W, Shepard DS, Avila-Figueroa C, Ahn H. Resource needs and gap analysis in achieving universal access to HIV/AIDS services: a data envelopment analysis of 45 countries. health policy plan; 2015.

    Google Scholar 

  30. Zeng W, Shepard DS, Chilingerian J, Avila-Figueroa C. How much can we gain from improved efficiency? An examination of performance of national HIV/AIDS programs and its determinants in low- and middle-income countries. BMC Health Serv Res. 2012;12:74.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Wanyenze RK, Wagner G, Alamo S, Amanyire G, Ouma J, Kwarisima D, Sunday P, Wabwire-Mangen F, Kamya M. Evaluation of the efficiency of patient flow at three HIV clinics in Uganda. AIDS Patient Care STDs. 2010;24(7):441–6.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Marseille E, Dandona L, Saba J, McConnel C, Rollins B, Gaist P, Lundberg M, Over M, Bertozzi S, Kahn JG. Assessing the efficiency of HIV prevention around the world: methods of the PANCEA project. Health Serv Res. 2004;39(6 Pt 2):1993–2012.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Holmes CB, Blandford JM, Sangrujee N, Stewart SR, DuBois A, Smith TR, Martin JC, Gavaghan A, Ryan CA, Goosby EP. PEPFAR’s past and future efforts to cut costs, improve efficiency, and increase the impact of global HIV programs. Health Aff (Millwood). 2012;31(7):1553–60.

    Article  Google Scholar 

  34. de Mendoza C, Rodriguez C, Corral A, del Romero J, Gallego O, Soriano V. Evidence for differences in the sexual transmission efficiency of HIV strains with distinct drug resistance genotypes. Clin Infect Dis. 2004;39(8):1231–8.

    Article  PubMed  Google Scholar 

  35. Lepine A, Vassall A, Chandrashekar S. The determinants of technical efficiency of a large scale hiv prevention project: application of the dea double bootstrap using panel data from the indian avahan. Cost Eff Resour Alloc. 2015;13:5.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Charnes A, Cooper WW, Rhodes E. Measuring the efficiency of decision making units; 1978.

    Book  Google Scholar 

  37. Han Y, Long C, Geng Z, Zhang K. Carbon emission analysis and evaluation of industrial departments in China: an improved environmental DEA cross model based on information entropy. J Environ Manage. 2018;205:298–307.

    Article  PubMed  Google Scholar 

  38. Halkos GE, Polemis ML. The impact of economic growth on environmental efficiency of the electricity sector: a hybrid window DEA methodology for the USA. J Environ Manage. 2018;211:334–46.

    Article  PubMed  Google Scholar 

  39. Flokou A, Aletras V, Niakas D. A window-DEA based efficiency evaluation of the public hospital sector in Greece during the 5-year economic crisis. PLoS One. 2017;12(5):e0177946.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Cavalieri M, Guccio C, Rizzo I. On the role of environmental corruption in healthcare infrastructures: an empirical assessment for Italy using DEA with truncated regression approach. Health Policy. 2017;121(5):515–24.

    Article  PubMed  Google Scholar 

  41. Egilmez G, McAvoy D. Benchmarking road safety of U.S. states: a DEA-based Malmquist productivity index approach. Accid Anal Prev. 2013;53:55–64.

    Article  PubMed  Google Scholar 

  42. Biorn E, Hagen TP, Iversen T, Magnussen J. The effect of activity-based financing on hospital efficiency: a panel data analysis of DEA efficiency scores 1992-2000. Health Care Manag Sci. 2003;6(4):271–83.

    Article  PubMed  Google Scholar 

  43. Zanakis SH, Alvarez C, Li V. Socio-economic determinants of HIV/AIDS pandemic and nations efficiencies. Eur J Oper Res. 2007;176(3):1811–38.

    Article  Google Scholar 

  44. Omondi Aduda DS, Ouma C, Onyango R, Onyango M, Bertrand J. Voluntary medical male circumcision scale-up in Nyanza, Kenya: evaluating technical efficiency and productivity of service delivery. PLoS One. 2015;10(2):e0118152.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  45. NEC. HIV/AIDS and ART registry of the Philippines (HARP): National Epidemiology Center; Manila, Philippines. 2015.

  46. Tenorio A. Supply side analysis of social hygiene clinics and treatment hubs as providers of human immunodeficiency virus (hiv) infection-related care and services in the philippine setting: revised final report; 2012.

    Google Scholar 

  47. Belimac G. Personal communication; 2015.

    Google Scholar 

  48. Santos SP, Amado CAE, Santos MF. Assessing the efficiency of mother-to-child HIV prevention in low- and middle-income countries using data envelopment analysis. Health Care Manag Sci. 2012;15(3):206–22.

    Article  PubMed  Google Scholar 

  49. Nakanishi YJ, Falcocchio JC. Performance assessment of intelligent transportation systems using data envelopment analysis. Res Transp Econ. 2004;8:181–97.

    Article  Google Scholar 

  50. Zere E, Mbeeli T, Shangula K, Mandlhate C, Mutirua K, Tjivambi B, Kapenambili W. Technical efficiency of district hospitals: evidence from Namibia using data envelopment analysis. Cost Eff Resour Alloc. 2006;4:5.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Jacobs R. Alternative methods to examine hospital efficiency: data envelopment analysis and stochastic frontier analysis. Health Care Manag Sci. 2001;4(2):103–15.

    Article  CAS  PubMed  Google Scholar 

  52. Lin LC, Tseng LA. Application of DEA and SFA on the measurement of operating efficiencies for 27 international container ports. Proc E Asia Soc Tran. 2005;5:592–607.

    Google Scholar 

  53. Po R-W, Guh Y-Y, Yang M-S. A new clustering approach using data envelopment analysis. Eur J Oper Res. 2009;199(1):276–84.

    Article  Google Scholar 

  54. Cotte Poveda A. Estimating effectiveness of the control of violence and socioeconomic development in Colombia: an application of dynamic data envelopment analysis and data panel approach. Soc Indic Res. 2012;105(3):343–66.

    Article  PubMed  Google Scholar 

  55. Akazili J, Adjuik M, Jehu-Appiah C, Zere E. Using data envelopment analysis to measure the extent of technical efficiency of public health centres in Ghana. BMC Int Health Hum Rights. 2008;8:11.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Flokou A, Kontodimopoulos N, Niakas D. Employing post-DEA Cross-evaluation and Cluster Analysis in a Sample of Greek NHS Hospitals. J Med Syst. 2011;35(5):1001–14.

    Article  PubMed  Google Scholar 

  57. Baker KR. Linear programming: data envelopment analysis. Optimization modeling with spreadsheets: Wiley; New Jersey, US. 2011. p. 175–210.

  58. Gale D. Linear programming and the simplex method. Not AMS. 2007;54(3):364–69.

  59. Huang Z, Lin H, Liu Y, Zhou M, Liu T, Xiao J, Zeng W, Li X, Zhang Y, Ebi KL, Tong S, Ma W, Wang L. Individual-level and community-level effect modifiers of the temperature-mortality relationship in 66 Chinese communities. BMJ Open. 2015;5(9):e009172.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Coelli T. A guide to frontier version 4.1: a computer program for stochastic frontier production and cost function estimation; 1996.

    Google Scholar 

  61. UNAIDS-PHL. Philippine HIV costing study: selected HIV prevention and treatment services. Manila: Joint United Nations Programme on HIV/AIDS; 2013.

    Google Scholar 

  62. Cook WD, Tone K, Zhu J. Data envelopment analysis: prior to choosing a model. Omega. 2014;44:1–4.

    Article  Google Scholar 

  63. Doyle J, Green R. Data envelopment analysis and multiple criteria decision-making. Omega-Int J Manag Sci. 1993;21(6):713–5.

    Article  Google Scholar 

  64. Adler N, Yazhemsky E. Improving discrimination in data envelopment analysis: PCA-DEA or variable reduction. Eur J Oper Res. 2010;202(1):273–84.

    Article  Google Scholar 

  65. Nataraja NR, Johnson AL. Guidelines for using variable selection techniques in data envelopment analysis. Eur J Oper Res. 2011;215(3):662–9.

    Article  Google Scholar 

  66. Põldaru R, Roots J. A PCA–DEA approach to measure the quality of life in Estonian counties. Socio Econ Plan Sci. 2014;48(1):65–73.

    Article  Google Scholar 

  67. Adler N, Golany B. PCA-DEA. In: Zhu J, Cook WD, editors. Modeling data irregularities and structural complexities in data envelopment analysis. Boston: Springer; 2007.

    Google Scholar 

  68. PNAC. Fifth AIDS Medium Term Plan - PHAREA. Manila: Philippine National AIDS Council; 2012.

    Google Scholar 

  69. UNAIDS-PHL. 2014 Global AIDS response progress reporting: UNAIDS-PHL; Manila, Philippines. 2014.

  70. Song ML, Zhang LL, Liu W, Fisher R. Bootstrap-DEA analysis of BRICS' energy efficiency based on small sample data. Appl Energ. 2013;112:1049–55.

    Article  Google Scholar 

  71. Hawdon D. Efficiency, performance and regulation of the international gas industry - a bootstrap DEA approach. Energ Policy. 2003;31(11):1167–78.

    Article  Google Scholar 

  72. Obure CD, Jacobs R, Guinness L, Mayhew S, Vassall A, Initiative I. Does integration of HIV and sexual and reproductive health services improve technical efficiency in Kenya and Swaziland? An application of a two-stage semi parametric approach incorporating quality measures. Soc Sci Med. 2016;151:147–56.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Halkos GE, Tzeremes NG. Industry performance evaluation with the use of financial ratios: an application of bootstrapped DEA. Expert Syst Appl. 2012;39(5):5872–80.

    Article  Google Scholar 

  74. Simar L, Wilson PW. Statistical inference in nonparametric frontier models: the state of the art. J Prod Anal. 2000;13(1):49–78.

    Article  Google Scholar 

  75. Assaf A, Matawie KM. Improving the accuracy of DEA efficiency analysis: a bootstrap application to the health care food service industry. Appl Econ. 2010;42(27):3547–58.

    Article  Google Scholar 

  76. Silwal PR, Ashton T. Productivity of public hospitals in Nepal: a data envelopment analysis. BMJ Open. 2017;7(7):e015327.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Polyzos N. A three-year performance evaluation of the NHS hospitals in Greece. Hippokratia. 2012;16(4):350–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  78. Gori E, Grassetti L, Rossi C. Linear mixed models in efficiency analysis: evidence from validation procedures. Stat Appl. 2002;14:247–75.

    Google Scholar 

  79. Harling G, Wood R, Beck EJ. Efficiency of interventions in hiv infection, 1994-2004. Dis Manag Health Out. 2005;13(6):371–94.

    Article  Google Scholar 

  80. Vassall A, Remme M, Watts C, Hallett T, Siapka M, Vickerman P, Terris-Prestholt F, Haacker M, Heise L, Haines A, Atun R, Piot P. Financing essential HIV services: a new economic agenda. PLoS Med. 2013;10(12):e1001567.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Piot P, Quinn TC. Response to the AIDS pandemic--a global health model. N Engl J Med. 2013;368(23):2210–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. DOH. Philippine MARP and PLHIV Estimates. 2011.

    Google Scholar 

  83. Hollingsworth B. The measurement of efficiency and productivity of health care delivery. Health Econ. 2008;17(10):1107–28.

    Article  PubMed  Google Scholar 

  84. OECD. National accounts 2011. Available from:

    Google Scholar 

  85. Thanassoulis E. Introduction to the theory and application of data envelopment analysis. Boston: Springer; 2001.

    Book  Google Scholar 

  86. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2010.

    Google Scholar 

Download references


We would like to acknowledge Mr. Zimmbodilion Y. Mosende, UNAIDS-Philippines, and Dr. Gerard Belimac, NASPCP – NCDPC (DOH-Philippines), for their invaluable comments and suggestions in the inception of this manuscript.


There is no funding to be declared for this study.

Availability of data and materials

Data used in the analyses can be accessed directly from the corresponding author via appropriate requests.

Author information

Authors and Affiliations



XTS and MK conceptualised the research topic. XTS carried out the data analysis and the writing of the manuscript. MK and IO provided extensive comments in the data analysis, manuscript writing and proofreading. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xerxes T. Seposo.

Ethics declarations

Ethics approval and consent to participate

No necessary ethics approval was necessary for utilization of the secondary data sources.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Detailed data management and analyses protocols. (DOCX 49 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Seposo, X.T., Okubo, I. & Kondo, M. Assessing frontline HIV service provider efficiency using data envelopment analysis: a case study of Philippine social hygiene clinics (SHCs). BMC Health Serv Res 19, 415 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: