Acessibilidade / Reportar erro

Medical centers location and specialists’ allocation: a healthcare planning case study

Abstract

Paper aims

To set the locations of new medical centers to meet the population’s secondary care needs, the additional number of specialists, equipment, and an installation sequence at municipalities.

Originality

We developed descriptive cost functions models and adopted aggregate data from official sources to set parameters of an integrated MILP model.

Research method

A case study at the Brazilian state of Minas Gerais.

Main findings

For every scenario, the recommended locations set centers dispersed over the state area, in cities with the minimum required infrastructure. We also propose a scenario of secondary care network re-design and demonstrate the reduced cost of such a strategy.

Implications for theory and practice

To automate the decision process, we developed a web-based system, providing flexibility and scientific-based results. Finally, we propose a sequence for installing 43 new medical centers and improving the capacity of 27 existing infrastructure based on equality principles.

Keywords:
Public healthcare planning; Facility location; Mixed integer linear programming

1. Introduction

In Brazil, decreasing economic progress imposes challenges to the deployment of good public healthcare planning. Setting the location of healthcare units in a hierarchical system, the size of hospitals, of medical centers or emergency units is critical since these decisions directly affect the population access to health services. Besides, cost-intensive procedures as computed tomography or nuclear magnetic resonance imaging technologies point to the promotion of inequalities in equipment access, which implies, in Brazil, a detachment from the principles of universality, integrality, and equity in health systems (Abreu, 2017Abreu, G. R. F. (2017). Acesso e utilização às tecnologias de diagnóstico por imagem pelo SUS no Estado da Bahia: o caso do tomógrafo computadorizado e equipamento de ressonância magnética (Dissertação de mestrado). Universidade Federal da Bahia, Salvador.).

The 1988’s Brazilian Federal Constitution established a national health system subsidized by social security payments and taxes (Sistema Único de Saúde - SUS) to provide free healthcare at the primary, secondary and tertiary levels to all residents (Victora et al., 2011Victora, C. G., Barreto, M. L., Leal, M., Monteiro, C. A., Schmidt, M. I., Paim, J., Bastos, F. I., Almeida, C., Bahia, L., Travassos, C., Reichenheim, M., & Barros, F. C. (2011). Health conditions and health-policy innovations in Brazil: the way forward. Lancet, 377(9782), 2042-2053. http://dx.doi.org/10.1016/S0140-6736(11)60055-X. PMid:21561659.
http://dx.doi.org/10.1016/S0140-6736(11)...
). The system follows principles of universalization, equity, and integrality. Every citizen can access primary, secondary and tertiary care in an integrated and decentralized public health system (Brasil, 2019aBrasil, Ministério da Saúde. (2019a). Sistema Único de Saúde (SUS): estrutura, princípios e como funciona. Brasília. Retrieved in 2019, July 15, from http://portalms.saude.gov.br/sistema-unico-de-saude
http://portalms.saude.gov.br/sistema-uni...
). Nonetheless, the Brazilian fast-aging community has raised services demand and increased the health system's requirements, but a recent federal law (Brasil, 2016Brasil. (2016, December 15). Altera o Ato das Disposições Constitucionais Transitórias, para instituir o Novo Regime Fiscal, e dá outras providências (Emenda constitucional nº 95, de 15 de dezembro de 2016). Diário Oficial da República Federativa do Brasil. Retrieved in 2019, October 10, from http://www.planalto.gov.br/ccivil_03/constituicao/emendas/emc/emc95.htm
http://www.planalto.gov.br/ccivil_03/con...
) has established even tighter bounds on the public health budget. Consequently, the location of health units and the allocation of resources have been a topic of both political and methodological relevance, pressuring leaders to strive on setting priorities in allocating resources.

In this study, we evaluate and analyze the secondary care of Minas Gerais (MG), a Brazilian state, where the 2000 and 2010 census demonstrate the fast aging diagnostic (see Figure 1). We acquire health needs on demographic projections (Instituto Brasileiro de Geografia e Estatística, 2019Instituto Brasileiro de Geografia e Estatística – IBGE. (2019). Projeção da população do Brasil e das Unidades da Federação. Rio de Janeiro. Retrieved in 2019, October 10, from https://www.ibge.gov.br/apps/populacao/projecao/
https://www.ibge.gov.br/apps/populacao/p...
) and the Ministerial Decree 1.631, and adopt mathematical programming models to solve an optimization problem of location and resource allocation (Brasil, 2015Brasil, Ministério da Saúde. (2015, October 2). Aprova critérios e parâmetros para o planejamento e programação de ações e serviços de saúde no âmbito do SUS (Portaria nº 1.631, de 1º de outubro de 2015). Diário Oficial da República Federativa do Brasil. Retrieved in 2017, September 28, from http://bvsms.saude.gov.br/bvs/saudelegis/gm/2015 /prt1631 01 10 2015.html
http://bvsms.saude.gov.br/bvs/saudelegis...
).

Figure 1
Demographic pyramid shows fast aging population of MG from 2000 to 2010 (Instituto Brasileiro de Geografia e Estatística, 2019)Instituto Brasileiro de Geografia e Estatística – IBGE. (2019). Projeção da população do Brasil e das Unidades da Federação. Rio de Janeiro. Retrieved in 2019, October 10, from https://www.ibge.gov.br/apps/populacao/projecao/
https://www.ibge.gov.br/apps/populacao/p...
.

Although SUS has improved access to primary and tertiary care over the past years, the provision of secondary is still problematic. The three service levels are interdependent (Macedo & Martin, 2014Macedo, L. M. D., & Martin, S. T. F. (2014). Interdependência entre os níveis de atenção do Sistema Único de Saúde (SUS): significado de integralidade apresentado por trabalhadores da Atenção Primária. Interface: Comunicação, Saúde, Educação, 18(51), 647-660. http://dx.doi.org/10.1590/1807-57622013.0597.
http://dx.doi.org/10.1590/1807-57622013....
). The primary care consists of low-cost prevention events, reaching in 2010 nearly 98 million people in 85% of Brazilian municipalities. On the extreme side, public teaching hospitals and contracted private sectors providers execute tertiary care, paid by the SUS at about market value. With intermediate technological density between primary and tertiary care, secondary care is responsible for performing medium complexity procedures (Minas Gerais, 2019Minas Gerais, Secretaria de Estado de Saúde. (2019). SUS. Belo Horizonte. Retrieved in 2019, June 10, from http://www.saude.mg.gov.br/sus
http://www.saude.mg.gov.br/sus...
) and providing specialized services to outpatient at the hospital level.

The provision of secondary care is problematic. Medium complexity procedures are often restricted to patients with private health plans. Since there is little regulation at the second level, SUS is highly reliant on the private sector. In Brazil, only 24.1% of Computed Tomography (CT) scanners and 13.4% of Magnetic Resonance Imaging (MRI) scanners are public (Paim et al., 2011Paim, J., Travassos, C., Almeida, C., Bahia, L., & Macinko, J. (2011). The Brazilian health system: history, advances, and challenges. Lancet, 377(9779), 1778-1797. http://dx.doi.org/10.1016/S0140-6736(11)60054-8. PMid:21561655.
http://dx.doi.org/10.1016/S0140-6736(11)...
). Furthermore, the medical centers, qualified for healthcare on the second level, are not geographically dispersed, which motivates this study of medical centers' location and equipment allocation for health care.

Healthcare location problems have been an active research area (Rais & Viana, 2011Rais, A., & Viana, A. (2011). Operations research in healthcare: a survey. International Transactions in Operational Research, 18(1), 1-31. http://dx.doi.org/10.1111/j.1475-3995.2010.00767.x.
http://dx.doi.org/10.1111/j.1475-3995.20...
). The models are an extension of the classic p-median, p-centers (Hakimi, 1964Hakimi, S. L. (1964). Optimum locations of switching centers and the absolute centers and medians of a graph. Operations Research, 12(3), 450-459. http://dx.doi.org/10.1287/opre.12.3.450.
http://dx.doi.org/10.1287/opre.12.3.450...
), set-covering (Toregas et al., 1971Toregas, C., Swain, R., ReVelle, C., & Bergman, L. (1971). The location of emergency service facilities. Operations Research, 19(6), 1363-1373. http://dx.doi.org/10.1287/opre.19.6.1363.
http://dx.doi.org/10.1287/opre.19.6.1363...
), and the max-covering model (Church & ReVelle, 1974Church, R., & ReVelle, C. (1974). The maximal covering location problem. Papers in Regional Science, 32(1), 101-118. http://dx.doi.org/10.1007/BF01942293.
http://dx.doi.org/10.1007/BF01942293...
), and the location-allocation model (Schilling et al., 1979Schilling, D., Elzinga, D. J., Cohon, J., Church, R., & ReVelle, C. (1979). The TEAM/FLEET models for simultaneous facility and equipment siting. Transportation Science, 13(2), 163-175. http://dx.doi.org/10.1287/trsc.13.2.163.
http://dx.doi.org/10.1287/trsc.13.2.163...
). The p-centers and the location-allocation models are of particular importance to this work because, collectively, they address issues of equity, minimizing the distance between remote patients, and the equipment allocation to facilities.

From a historical perspective, we refer the interested readers to a comprehensive review (Hale & Moberg, 2003Hale, T. S., & Moberg, C. R. (2003). Location science research: a review. Annals of Operations Research, 123(1), 21-35. http://dx.doi.org/10.1023/A:1026110926707.
http://dx.doi.org/10.1023/A:102611092670...
), successful case studies (Brandeau, 2016Brandeau, M. L. (2016). Creating impact with operations research in health: making room for practice in academia. Health Care Management Science, 19(4), 305-312. http://dx.doi.org/10.1007/s10729-015-9328-0. PMid:26003321.
http://dx.doi.org/10.1007/s10729-015-932...
) and a recent survey (Ahmadi-Javid et al., 2017Ahmadi-Javid, A., Seyedi, P., & Syam, S. S. (2017). A survey of healthcare facility location. Computers & Operations Research, 79, 223-263. http://dx.doi.org/10.1016/j.cor.2016.05.018.
http://dx.doi.org/10.1016/j.cor.2016.05....
). Applications include modeling a 3-level location system of perinatal facilities (Galvão et al., 2002Galvão, R. D., Acosta Espejo, L. G., & Boey, B. (2002). A hierarchical model for the location of perinatal facilities in the municipality of Rio de Janeiro. European Journal of Operational Research, 138(3), 495-517. http://dx.doi.org/10.1016/S0377-2217(01)00172-2.
http://dx.doi.org/10.1016/S0377-2217(01)...
); the use of a genetic algorithm for a maximal covering location for health care (Shariff et al., 2012Shariff, S. S. R., Moin, N. H., & Omar, M. (2012). Location allocation modeling for healthcare facility planning in Malaysia. Computers & Industrial Engineering, 62(4), 1000-1010. http://dx.doi.org/10.1016/j.cie.2011.12.026.
http://dx.doi.org/10.1016/j.cie.2011.12....
); and a facility location model for primary care system approaching physicians’ preferences (Güneş et al., 2014Güneş, E. D., Yaman, H., Çekyay, B., & Verter, V. (2014). Matching patient and physician preferences in designing a primary care facility network. The Journal of the Operational Research Society, 65(4), 483-496. http://dx.doi.org/10.1057/jors.2012.71.
http://dx.doi.org/10.1057/jors.2012.71...
). Following, Khodaparasti et al. (2016)Khodaparasti, S., Maleki, H., Jahedi, S., Bruni, M., & Beraldi, P. (2016). Enhancing community-based health programs in Iran: a multi-objective location-allocation model. Health Care Management Science, 20(4), 485-499. PMid:27115695. address a hierarchical location model for community-based organizations with health care providers by applying multi-objective programming. Additional applications consists on the development of a system to optimize the equipment allocation (Treurnicht & Van Dyk, 2014Treurnicht, M. J., & Van Dyk, L. (2014). A decision support system for equipment allocation in a telemedicine referral network. South African Journal of Industrial Engineering, 25(1), 29-38. http://dx.doi.org/10.7166/25-1-641.
http://dx.doi.org/10.7166/25-1-641...
); the tactical workforce planning and capacity allocation minimizing salary costs considering vacation and subcontracting opportunities (Van der Veen et al., 2015Van der Veen, E., Hans, E. W., Veltman, B., Berrevoets, L. M., & Berden, H. J. (2015). A case study of cost-efficient staffing under annualized hours. Health Care Management Science, 18(3), 279-288. http://dx.doi.org/10.1007/s10729-014-9292-0. PMid:24997580.
http://dx.doi.org/10.1007/s10729-014-929...
); and a performance evaluation of 87 private and public primary care health units (Novignon & Nonvignon, 2017Novignon, J., & Nonvignon, J. (2017). Improving primary health care facility performance in Ghana: efficiency analysis and scal space implications. BMC Health Services Research, 17(1), 399. http://dx.doi.org/10.1186/s12913-017-2347-4. PMid:28606131.
http://dx.doi.org/10.1186/s12913-017-234...
). On health care location-allocation problem, examples include the location of preventive health care units (Gu et al., 2010Gu, W., Wang, X., & McGregor, S. E. (2010). Optimization of preventive health care facility locations. International Journal of Health Geographics, 9(1), 17. http://dx.doi.org/10.1186/1476-072X-9-17. PMid:20298608.
http://dx.doi.org/10.1186/1476-072X-9-17...
), two-phase procedures to set the location and size of medical departments (Stummer et al., 2004Stummer, C., Doerner, K., Focke, A., & Heidenberger, K. (2004). Determining location and size of medical departments in a hospital network: a multiobjective decision support approach. Health Care Management Science, 7(1), 63-71. http://dx.doi.org/10.1023/B:HCMS.0000005399.23600.69. PMid:14977095.
http://dx.doi.org/10.1023/B:HCMS.0000005...
), and for improving spatial accessibility (Luo et al., 2017Luo, J., Tian, L., Luo, L., Yi, H., & Wang, F. (2017). Two-step optimization for spatial accessibility improvement: a case study of health care planning in rural China. BioMed Research International, 2017, 2094654. http://dx.doi.org/10.1155/2017/2094654. PMid:28484707.
http://dx.doi.org/10.1155/2017/2094654...
), and an integrated approach aiming inequality reduction (Sathler et al., 2019Sathler, T. M., Almeida, J. F., Conceição, S. V., Pinto, L. R., & Campos, F. C. (2019). Integration of facility location and equipment allocation in health care management. Brazilian Journal of Operations & Production Management, 16(3), 513-527. http://dx.doi.org/10.14488/BJOPM.2019.v16.n3.a13.
http://dx.doi.org/10.14488/BJOPM.2019.v1...
).

Despite the rich literature, most studies address primary care on the municipality level. Little attention, however, has been paid to location-allocation problems for health care on the secondary level, so there are potential gaps and open issues yet to investigate, besides, the literature presents few studies that address realistic problems (Güneş et al., 2014Güneş, E. D., Yaman, H., Çekyay, B., & Verter, V. (2014). Matching patient and physician preferences in designing a primary care facility network. The Journal of the Operational Research Society, 65(4), 483-496. http://dx.doi.org/10.1057/jors.2012.71.
http://dx.doi.org/10.1057/jors.2012.71...
; Ahmadi-Javid et al., 2017Ahmadi-Javid, A., Seyedi, P., & Syam, S. S. (2017). A survey of healthcare facility location. Computers & Operations Research, 79, 223-263. http://dx.doi.org/10.1016/j.cor.2016.05.018.
http://dx.doi.org/10.1016/j.cor.2016.05....
). Our model: (i) adopts descriptive models of cost functions and dada aggregation to elaborate parameters; (ii) evaluates the existing equipment idle capacity on a hierarchical health care system; (iii) considers both the population requirements for local multiple medical services and physicians’ preferences on living on municipalities with minimum infrastructure; and (iv) the government financial condition, as stated in objective function.

The authors developed the method following flexibility principles (De Neufville & Scholtes, 2011De Neufville, R., & Scholtes, S. (2011). Flexibility in engineering design. Cambridge: MIT Press. http://dx.doi.org/10.7551/mitpress/8292.001.0001.
http://dx.doi.org/10.7551/mitpress/8292....
) for dealing with the problem faced in a project. Since we recognized that the context was of a generic style, representing an unstudied component of real health care location-allocation problems (Almeida et al., 2018Almeida, J. F. F., Pinto, L. R., Conceição, S. V., Campos, F. C. C., & Miranda Júnior, G. (2018, December). A two-step optimization process for medical centers location and capacity allocation. In A. Leiras, C. A. González-Calderón, I. Brito Junior, S. Villa & H. T. Y. Yoshizaki (Eds.), Operations Management for Social Good: 2018 POMS International Conference in Rio (Springer Proceedings in Business and Economics). USA: Springer International Publishing.), this work can directly contribute to general health care planning by supporting policymaking and providing scenario analyses.

In the following sections, we present a Mixed Integer Linear Programming (MILP) location-allocation model, a descriptive model based on cost functions for setting health care costs, and a strategy for installing health care units based on equality principles. Next, we provide a case study applying the proposed method to solve the problem in MG. We evaluate scenarios by distance, by budget availability are discuss a sequence for health units’ installation. We finish the paper by providing conclusions and suggestions for future research.

2. Method

This paper presents a MILP model for healthcare planning at the state level, particularly, secondary care. Cost functions descriptive models and aggregate data are the base for the MILP model’s parameters. The mathematical model aids the planner to select secondary care facilities locations and designate medical specialists, and equipment to them. The planner, here represented by the government health care management team, must combine patients' and physicians’ inclinations, adequately size capacity, and guarantee service quality to justify the public infrastructure financing. As medical procedures often depend on equipment, we include these recourses and its’ demands into the model considering official standards’ amount of equipment per inhabitants and the actual available capacity.

The aim is minimizing the annual overall cost of a strategy that includes (i) the maintenance costs of medical centers; (ii) the cost of procurement of new equipment for existing infrastructure and new medical centers, and balance with (iii) a social opportunity cost of society, estimated by cost functions. We also evaluate available specialists and equipment in a hierarchical system. Secondary care specialists working in primary or tertiary care includes the available capacity. Their production is decreased from the estimated demand, set on Ministerial Decree 1.631. Datasus Outpatient and Hospital Information System and Applications (Departamento de Informática do SUS, 2018aDepartamento de Informática do SUS – DATASUS. (2018a). Cadastro Nacional de Estabelecimento de Saúde. Brasília. Retrieved in 2018, September 10, from http://cnes.datasus.gov.br/
http://cnes.datasus.gov.br/...
, bDepartamento de Informática do SUS – DATASUS. (2018b). Sistemas e aplicativos. Brasília. Retrieved in 2018, September 10, from http://datasus.saude.gov.br/sistemas-e-aplicativos
http://datasus.saude.gov.br/sistemas-e-a...
) provides a history of outpatient and hospital production.

We denote a state with I municipalities, offering E medical specialties services and Q equipment. Patients move within their municipality or via available paths i,jKI×I between pairs of municipalities to meet E medical specialties (see Table 1).

Table 1
Input parameters of the facility location and resource allocation model.

The number of inhabitants Bi per municipality i is sourced from Instituto Brasileiro de Geografia e Estatística (2019)Instituto Brasileiro de Geografia e Estatística – IBGE. (2019). Projeção da população do Brasil e das Unidades da Federação. Rio de Janeiro. Retrieved in 2019, October 10, from https://www.ibge.gov.br/apps/populacao/projecao/
https://www.ibge.gov.br/apps/populacao/p...
projection to 2019. The observed physicians’ preference for municipalities with a population higher than a minimum amount Bmin and better infrastructure is considered. A basic measure for selecting acceptable distances is defining a lower bound parameter of the maximal allowed patients’ displacement as maximinijK|ijDij for critical medical specialties services, like cardiology. Open Route online Services provides the distance between municipalities. The data helps setting candidates based on the maximum displacement DMeij of patients, which results from the combination of medical specialty and distances (or time) to be established according to government parameters for service levels and budget availability. Therefore, the candidate municipalities J=iNi are a set of locations that meet the following requirements Ni=j|i,jK,DijmaxeDMeijBjBmin, and |Ni| must not be empty for any municipality. Services often don’t meet nominal capacities based on projections, therefore, the global efficiency of services Ej provided on municipality j can be adjusted based on observations combined with data envelopment analysis (DEA), moreover, the decision-maker may have the option to consider Pj municipalities that are provider of specialties services with available capacity on healthcare services.

The demand of patients for medical specialties services HEie and exams on equipment HQiq is available on official government publications; however, the decision-maker can estimate demand based on benchmarking indicators of OECD countries (Organisation for Economic Co-operation and Development, 2017Organisation for Economic Co-operation and Development – OECD. (2017), Health at a Glance 2017: OECD indicators. Paris: OECD Publishing. https://doi.org/10.1787/health_glance-2017-en.
https://doi.org/10.1787/health_glance-20...
), for instance. The demand does not remain constant through every geographical region but decays proportionally to the distance. We propose a simplified decay function αeij=maxDMeijDij,0/DMeij based on distance; however, this function can be improved by a floating catchment area (FCA) method based on real data (Bauer & Groneberg, 2016Bauer, J., & Groneberg, D. A. (2016). Measuring spatial accessibility of health care providers–introduction of a variable distance decay function within the Floating Catchment Area (FCA) method. PLoS One, 11(7), e0159148. http://dx.doi.org/10.1371/journal.pone.0159148. PMid:27391649.
http://dx.doi.org/10.1371/journal.pone.0...
).

The (i) nominal capacity of the equipment Kq is different from the (ii) available capacity of the equipment KQiq. The first consists in the Original Equipment Manufacturer (OEM) information of annual production capacity of each equipment, individually, while the second is related to the municipality capacity on providing the service of exams on equipment, therefore, the information is obtained by the number of available equipment on healthcare units. The availability of physicians KEie of each specialty is based on historical data of similar periods in the past and consists of registries of the historical production of physicians on each municipality. A matrix Meq assigns each type of medical specialty to the equipment used.

The government sets parameters of a minimum amount of secondary health care facilities F and equipment Pq according to its financial conditions and to the desired service level, it aims at providing to its population. The cost parameters of the average annual cost of Medical Centers operations, acquisition of equipment, and patients’ displacement between municipalities guide the objective function that aims at minimizing the sum of government and population costs.

We propose the use of cost functions (Horngren et al., 2002Horngren, C. T., Bhimani, A., Datar, S. M., Foster, G., & Horngren, C. T. (2002). Management and cost accounting. Harlow: Financial Times/Prentice Hall.) to estimate fixed and variable costs of the annual cost of medical center operations. Let S be a set of support team that enables the operations of health centers, Ns, the number of professionals of each team S, and Ss, the salary of each professional of support (Brasil, 2019bBrasil, Ministério do Trabalho. (2019b). PDET/CAGED: Programa de Dissseminação das Estatísticas do Trabalho. Brasília. Retrieved in 2019, June 10, from pdet.mte.gov.br/caged
pdet.mte.gov.br/caged...
). The fixed costs FC include a periodic time t payment of salaries of multifunctional teams, as nurses, pharmacists, psychologists, social workers, and, administrative support, their corresponding social taxesT, and a General Costs (GC), therefore,FC=stTNsSs+GC. The variable costs VC consist of the average number of hours provide by each medical specialty e and their category’s average salary Se (Brasil, 2019bBrasil, Ministério do Trabalho. (2019b). PDET/CAGED: Programa de Dissseminação das Estatísticas do Trabalho. Brasília. Retrieved in 2019, June 10, from pdet.mte.gov.br/caged
pdet.mte.gov.br/caged...
), and divided by 40 to meet a full-time equivalent professional (FTE), therefore, VCe=tweeSe, and the average annual cost of a medical center facility operationsCI=FC+eVCe.

The additional annual cost of new equipment is equivalent to the cost of the equipment divided by the depreciation period. General costs may include average operational equipment costs. The values of each type of equipment follow average market prices. To estimate the annual patients’ displacement costs, let r be the average rate per kilometer of vehicle, s, the average speed vehicle, t, the annual number of trips, h, the time on queue and meeting physician, and m the average income of patient per hour. A medical care meet requires two trips (go, and return). The annual estimated number of meetings is base for setting the number of trips. For example, a physician that meets one patient per hour may meet 1×8×5×4.25×12=2,040 meets per year, and t=2,040×2=4,080 trips per year. The opportunity cost spent by patients considers the time on round trips2Dij/s, plus the average time h. This time is multiplied by the average hourly income of a patient (m), therefore, the displacement costs are represented byCDij=2rDij+2Dij/s+hmt.

The data for parameters above are base for answering some questions, translated into model variables: Which municipality should receive a secondary level healthcare unit? How many physicians of specialty e should be hired to such municipalities and how many equipment of type q should be acquired for each of them? The demand of patients from municipality i for medical specialty e and their demand for exams on equipment q should be satisfied in which municipality? Table 2 presents the variables from the current questions.

Table 2
Input variables of the facility location and resource allocation model.
m i n i m i z e j C I z j + S e w e e j + q j C Q q w q q j + e i j α e i j C D i j H E i e x e i j (1.0)

The objective function 1.0 minimizes the total annual cost composed of three classes of costs shared by the government and the society, that is (i) the maintenance of medical centers, (ii) the depreciation costs of newly acquired equipment of different types for each unit, and a social cost of patients’ displacement for having medical care. This function is subject to constraints 1.1-1.17.

j x e i j = 1 e E , i I (1.1)
j y q i j = 1 q Q , i I (1.2)

In Equation 1.1 and Equation 1.2, the demand for medical services and the demand for exams based on equipment, respectively, of each district must be satisfied by only one medical center. Thus, the patient municipality does not have to assign the demand of every medical specialty nor every equipment to a single destination municipality, but it can be assigned to different municipalities, however, the split of demand of a patient municipality is not allowed for each medical specialty or equipment.

x e i j z j e E , i I , j J (1.3)
y q i j z j q Q , i I , j J (1.4)

Following, the demand for medical specialty, on the Constraint 1.3, or for medical equipment, on the Constraint 1.4, can only be assigned to a selected municipality, that is, the demand assignment is conditioned to the installation of a Medical Center on municipality j.

x e j j = z j e E , j J (1.5)

The Equation 1.5 establishes that the demand of a candidate municipality that was selected to receive a Medical Center must be locally satisfied since patients would not displace to different municipalities if its own municipality has healthcare units. Even if its health unit is overcrowded, the patient may not know nor may not accept being moved to different municipalities for healthcare services.

x e i j = y q i j e E , q Q , i I , j J | M e q = 1 (1.6)
w q q j z j F P q q Q , j J (1.7)

Constraint 1.6 link the medical specialty demand to the type of equipment each specialty depends on, such that the patient does not have to move to different municipalities to take exams. On Constraint 1.7, the number of equipment of type q is limited by an upper bound on selected municipalities or equal zero, if the municipality is not selected to receive a healthcare unit.

j z j F j J (1.8)
j w q q j P q q Q (1.9)

The healthcare service level desired by the government is presented on Constraints 1.8 and 1.9 which defines the minimum number of Medical Centers and equipment, respectively. The more healthcare units, the higher is the healthcare service level, but the higher is the operating cost with administration, physicians, and equipment. Besides, the number of medical centers and equipment is dependent on the state’s long-term financial conditions.

z j P j j J (1.10)

Sometimes, an optimal solution may yield to municipalities with no infrastructure, which may be very close to other municipalities with available infrastructure and idle capacity. For such cases, Constraint 1.10 is interesting for decision-makers to force the selection of some candidates with better infrastructure, which may reduce the healthcare unit installation cost.

K q E j w q q j i max H Q i q K Q i q , 0 y q i j q Q , j J (1.11)
w e e j i max H E i e K E i e , 0 x e i j e E , j J (1.12)

Constraints 1.11 and 1.12 project the minimum capacity to satisfy the demand of patients for exams and medical specialties, respectively. The above constraints consider the existing capacity of the available infrastructure of healthcare on the state to acquire only the extra recourses required to satisfy the demand. Finally, Constraints 1.13 to 1.17 set the domain of the model’s variables.

x e i j 0,1 e E , i I , j J (1.13)
y q i j 0,1 q Q , i I , j J (1.14)
z j 0,1 j J (1.15)
w q q j + q Q , j J (1.16)
w e e j + e E , j J (1.17)

On solving the optimization problem, we obtain the number of municipalities, medical specialists, and equipment that satisfy all described constraints at a minimum overall cost. This solution, however, is static; therefore, we suggest a sequence for installing medical centers guided by inequality indexes. The inequality indexes are the Social Vulnerability Index (SVI) (Costa & Marguti, 2015Costa, M. A., & Marguti, B. O. E. (2015). Atlas da vulnerabilidade social nos municípios brasileiros. Brasília: IPEA.), the Inequality-adjusted Human Development Index (IHDI) from the United Nations Development Programme (2018)United Nations Development Programme – UNDP. (2018). Human development indices and indicators: 2018 statistical update. New York., or the Slope Index of Inequality (SII), and the Concentration Index (CIX), based on demographic health survey (Silva et al., 2018Silva, I. C. M. D., Restrepo-Mendez, M. C., Costa, J. C., Ewerling, F., Hellwig, F., Ferreira, L. Z., Ruas, L. P. V., Joseph, G., & Barros, A. J. D. (2018). Measurement of social inequalities in health: concepts and methodological approaches in the Brazilian context. Epidemiologia e Serviços de Saúde, 27(1), e000100017. http://dx.doi.org/10.5123/S1679-49742018000100017. PMid:29513856.
http://dx.doi.org/10.5123/S1679-49742018...
).

In this study, the authors propose the above methodology to aid Minas Gerais (MG) health care policy. MG is the Brazilian state with the largest number of municipalities, 853, the second most populous, with more than 20 million inhabitants, and fourth biggest, with 586,528 km2. Solving this health care problem to MG indicates that the method can solve the health care problem for the remaining states of the country.

3. Results and discussion

Minas Gerais proposed in 2015 advances in health care including a decentralization process in the second level to increase the population access to specialized centers. The government inaugurated the first medical center, in this new arrangement, in 2016 (Jornal Estado de Minas, 2016Jornal Estado de Minas. (2016, 12 abril). Primeiro Centro de Especialidades Médicas de Minas é inaugurado em Pirapora. Belo Horizonte. Retrieved in 2017, September 28, from https://www.em.com.br/app/noticia/gerais/2016/04/12/interna_gerais,752198/primeiro-centro-de-especialidades-medicas-de-minas-e-inaugurado-em-pir.shtml
https://www.em.com.br/app/noticia/gerais...
). Some constraints were established. Candidate municipalities should satisfy a combination of requirements of least 30 thousand inhabitants to receive nine medical specialties services and equipment, which reveals physicians’ preference on minimal infrastructure and on working with similar professionals. The patients’ preference for minimal displacement establishes a maximal allowed distance from patients’ municipalities to medical centers’ municipalities.

The authors adopted open route services (Neis & Zipf, 2007Neis, P., & Zipf, A. (2007). A web accessibility analysis service based on the OpenLS route service. In International Conference on Geographic Information Science of the Association of Geograpic Information Laboratories for Europe (AGILE 2007). Aalborg, Denmark.) to obtain time and distance between the 853 municipalities (i,j) yielding 727,609 registries for each parameter. Following, we defined the lower bound parameter maximinijK|ijDij=77km of the maximal allowed distance; however, 4.2% of municipalities do not meet these both requirements (more than 30 thousand inhabitants, and at most 77 km from candidate municipalities). Those municipalities that do not meet these requirements are an exception to a general rule and require a practical solution with special conditions. In general, we established a displacement ranging from one to three hours according to the accessibility need to health care services, setting 80 km for cardiology, and 120 km to 180 km to the remaining specialties resulting in the selection of 121 candidates.

We set the demand of patients from 853 municipalities for physicians of specialties HEie and exams HQiq on equipment as stochastic parameters Ua,b. Demands are uniformly distributed ranging from a to b. The parameters a,b range from a central value (central value – range, central value + range) which is selected on Ministerial Degree 1,631 (Brasil, 2015Brasil, Ministério da Saúde. (2015, October 2). Aprova critérios e parâmetros para o planejamento e programação de ações e serviços de saúde no âmbito do SUS (Portaria nº 1.631, de 1º de outubro de 2015). Diário Oficial da República Federativa do Brasil. Retrieved in 2017, September 28, from http://bvsms.saude.gov.br/bvs/saudelegis/gm/2015 /prt1631 01 10 2015.html
http://bvsms.saude.gov.br/bvs/saudelegis...
). We adopt a range variability of 0.26 based on indicators of Mexico physicians’ density, which is a country included in Organisation for Economic Co-operation and Development (2017)Organisation for Economic Co-operation and Development – OECD. (2017), Health at a Glance 2017: OECD indicators. Paris: OECD Publishing. https://doi.org/10.1787/health_glance-2017-en.
https://doi.org/10.1787/health_glance-20...
and comparable to Brazil in indicators, like GDP.

The OEM information of standard equipment set base parameters of the nominal capacity of MRI scanners, Mammograph, CT scanners, Ultrasound, and Doppler. DATASUS database provides the municipality capacity on offering the service of exams for each type of equipment and the number of machines available on 853 municipalities (Departamento de Informática do SUS 2018aDepartamento de Informática do SUS – DATASUS. (2018a). Cadastro Nacional de Estabelecimento de Saúde. Brasília. Retrieved in 2018, September 10, from http://cnes.datasus.gov.br/
http://cnes.datasus.gov.br/...
, bDepartamento de Informática do SUS – DATASUS. (2018b). Sistemas e aplicativos. Brasília. Retrieved in 2018, September 10, from http://datasus.saude.gov.br/sistemas-e-aplicativos
http://datasus.saude.gov.br/sistemas-e-a...
). We obtained 10 years of the historical production of physicians. We reached the information of availability of FTE physicians’ specialists in angiology, cardiology, endocrinology, gynecology, mastology, nephrology, ophthalmology, pediatrics, and urology by dividing the average values of physicians’ historical production from 2008 to 2017 on each municipality by a standard period of work of 40 hours of FTE physician.

The government sets the minimum number of secondary health care facilities and equipment according to its financial conditions and the desired service level it aims at providing to its population. This study evaluates the optimal strategy with minimum annual overall costs of heal care on the secondary level. The method considers simultaneously the installation of medical centers in different municipalities of MG, the provision of nine medical specialties services, and exams on five types of equipment.

The annual cost function of medical center operations considers the salary of multifunctional teams of nurses, pharmacists, psychologists, social workers, and, administrative support include fixed costs. The variable costs consist of the average number of hours provided by each medical specialty and their category’s average salary, available at PDET (Brasil, 2019bBrasil, Ministério do Trabalho. (2019b). PDET/CAGED: Programa de Dissseminação das Estatísticas do Trabalho. Brasília. Retrieved in 2019, June 10, from pdet.mte.gov.br/caged
pdet.mte.gov.br/caged...
). Following, we include the annual cost of new equipment, traditionally depreciated in 10 years. The values are based on average prices of catalogs of the standard equipment of each type, therefore, we set depreciation costs of R$300,000 for MRI scanner, R$41,000 for CT scanners, R$5,600 for mammograph, R$4,500 for ultrasound, and R$850 for doppler. Finally, the annual patients’ displacement costs between municipalities of MG considered: (i) the average rate of Reals per kilometer (R$0.50/km) of bus companies for two (round trip) inter-municipal trips; (ii) the opportunity cost spent by patients; and (iii) the average number of trips, for health care with physicians, in a year. The number of trips is a stochastic parameter Nμ,σ based on the estimated number of meetings. Therefore, medical specialties meeting from one to two patients per hour yields about N2000,300 trips per year. The opportunity cost considers the time on round trips at 60 km/h, plus an average time of 5 hours to arrive, wait on queues, receive the medical care, and return to the bus. This time is multiplied by the average hourly income of an MG’s citizen (R$9.45), available at Instituto Brasileiro de Geografia e Estatística (2010aInstituto Brasileiro de Geografia e Estatística – IBGE. (2010a). Censo demográfico 2010: Rendimento mensal domiciliary per capita segundo as grandes regiões e unidades da Federação. Rio de Janeiro. Retrieved in 2019, October 10, from https://censo2010.ibge.gov.br/sinopse/index.php?dados=P21&uf=00
https://censo2010.ibge.gov.br/sinopse/in...
, bInstituto Brasileiro de Geografia e Estatística – IBGE. (2010b). Censo demográfico 2010. Rio de Janeiro. Retrieved in 2019, September 10, from https://www.ibge.gov.br/estatisticas-novoportal/sociais/populacao/9662-censodemogra_co-2010.html
https://www.ibge.gov.br/estatisticas-nov...
), therefore, displacement costs CDij depend on the distance and are represented by R$20.50Dij+2Dij/60+59.45N2000,300.

We evaluated eight optimized scenarios. The scenarios [1-3] consider the maximal patients’ displacement, highlighting municipalities that need special treatment. The scenarios [4-6] takes the government budget availability into account. Following, the scenario [7] proposes a re-design of MG’s network of medical centers, and finally, on scenario [8] we suggest a sequence for setting medical centers, based on equity principles. We implemented in MathProg (GLPK) and run the MILP models in a Linux Mint 17.3 64-bit, RAM of 8 GB, Intel Core I5 2.50 GHz x 2 processor. The mathematical programming problem comprises 1,445,574 constraints and 709,193 variables, all of which are integers, being 683,586 of them, binary. Instances were optimally solved within two hours.

The recent disclosure of MG State Secretary of Health shows that the health care policy of installing new medical specialties centers units has not been implemented yet (Ricardo, 2019Ricardo, P. (2019). Centro de Especialidades Médicas amplia acesso aos serviços de saúde no Norte de Minas. Belo Horizonte: Secretaria de Estado de Saúde de Minas Gerais. Retrieved in 2019, October 10, from http://www.saude.mg.gov.br/ngc/story/8735-centro-de-especialidades-medicas-amplia-acesso-aos-servicos-de-saude-no-norte-de-minas
http://www.saude.mg.gov.br/ngc/story/873...
). The state government has recently faced a serious financial crisis (Agência Minas, 2019Agência Minas. (2019, February 22). Nota de esclarecimento. Belo Horizonte: Secretaria de Estado de Saúde de Minas Gerais. Retrieved in 2019, September 12, from http://www.agenciaminas.mg.gov.br/noticia/nota-de-esclarecimento-2019-02-22-14-50-00-0300
http://www.agenciaminas.mg.gov.br/notici...
), and much of the health resource has been used for working capital expenditures rather than investments (Assembleia Legislativa do Estado de Minas Gerais, 2018Assembleia Legislativa do Estado de Minas Gerais – ALMG. (2018). Relatório comprova baixo investimento na saúde em Minas. Belo Horizonte. Retrieved in 2019, September 29, from https://www.almg.gov.br/acompanhe/noticias/arquivos/2018/07/11_audiencia_saude_relatorio_SUS.html
https://www.almg.gov.br/acompanhe/notici...
). The State government has transferred financial resources (R$19,350,970.03) to the existing 26 Specialized Care Centers (with Brazilian Portuguese initials CEAE) and 1 Medical Center (with Brazilian Portuguese initials CEM) (Minas Gerais, 2018Minas Gerais. (2018, December 29). Altera o Art. 2° e o Anexo Único da Resolução SES/MG nº 6.495, de 13 de novembro de 2018, que autoriza o repasse de incentivo financeiro, em caráter excepcional, para apoio e fortalecimento dos Centros Estaduais de Atenção Especializada (CEAE) e do Centro de Especialidades Médicas (CEM) (Resolução SES/MG nº 6.563, de 21 de dezembro de 2018). Diário Oficial do Estado de Minas Gerais. Retrieved in 2019, August 15, from http://saude.mg.gov.br/images/documentos/RESOLUÇÃO_6563.pdf
http://saude.mg.gov.br/images/documentos...
). Since the 27 health units are dispersed over the state (see Figure 2), we adopt the plausible assumption that the model must capacitate the existing health care centers before creating new units.

Figure 2
The 27 health units (Specialized Care Centers and a Medical Center) dispersed over the state area. Source: Software Tabwin 2019.

3.1. Scenarios [1-3] based on patients’ displacement

Scenarios based on patients’ displacement provides the minimum number of municipalities, equipment, and specialists that satisfy the maximum patients’ displacement constraints. Therefore, we deactivate constraints 1.8 and 1.9 since parameter DMeij imposes the limitation. We present three scenarios (see Figure 3) based on patients that move distances ranging from 180 km to 80 km (or take from 3 hours to approximate 1 hour). We describe the third scenario in more detail; evaluating the number of selected municipalities, the additional amount of each type of equipment, and the suggestion of hiring extra FTE specialists. Besides, we evaluate optimal overall costs and compare the state’s costs with social costs.

Figure 3
Municipalities selected (43, 45, and 64) for patients’ displacement of at most 180, 120, and 80 km, respectively, to meet health care on medical centers. Municipalities in purple (scenarios 2 a 3) are exception requiring special treatment.

In the scenario [3], patients should move at most 80 km, or approximate 1 hour and 20 minutes to reach a health unit. In this scenario, 64 municipalities (in black and green, see Figure 4) offer secondary care at medical centers. Since 27 municipalities (in black, see Figure 4) already have a basic infrastructure (health units), the government should install 37 additional medical centers (in green, see Figure 4).

Figure 4
Municipalities selected in black and green enable the patients’ displacement of at most 80 km to meet health care on medical centers. Municipalities in purple are exception requiring special treatment.

In this strategy, the government should hire 73 angiologists, 137 cardiologists, 69 endocrinologists, 183 gynecologists, 72 mastologists, 102 nephrologists, 160 ophthalmologists, 211 pediatrics, and 91 urologists (see Figure 5). Besides, it requires the procurement of 133 MRIs, 75 CTs, 142 mammograms, 957 ultrasounds, and 63 ultrasound Doppler (see Figure 6). The strategy satisfies the demand of 816 municipalities, 95.8% of the state demand. The 37 municipalities in purple (see Figure 4) are Águas Formosas, Águas Vermelhas, Aimorés, Cachoeira de Pajeú, Arantina, Arinos, Bertópolis, Bom Jardim de Minas, Bonfinópolis de Minas, Brasilândia de Minas, Buritis, Chapada Gaúcha, Crisólita, Cuparaque, Divisa Alegre, Formoso, Fronteira dos Vales, Itabirinha, Itueta, Juvenília, Machacalis, Mantena, Medina, Miravânia, Montalvânia, Monte Formoso, Nova Belém, Passa-Vinte, Resplendor, Riachinho, Salto da Divisa, Santa Fé de Minas, Santa Helena de Minas, Santa Rita de Jacutinga, Santa Rita do Itueto, São João do Manteninha, and Urucuia. They represent 4.2% of municipalities that deviates the general rule, where patients move 110±20 km to reach destination. To follow equality principles, these municipalities require special treatment, as dedicated public transport, for example.

Figure 5
Additional hiring of medical specialties for patients' maximum displacement of 80 km.
Figure 6
Additional equipment for patients' maximum displacement of 80 km.

The analysis revealed each medical center serving an average of 14 districts patients traveling on average 57 km and a standard deviation of 9 km from their municipality to the nearest medical center. The maximal distance is 80 km, limited by the scenario, while the shortest distance, excepting patients with medical center on their municipality, is 33 km. The Appendix A Appendix A Equipment assignment and number of specialists for each municipality of MG. This appendix presents the scenario [3] selection of medical centers and equipment considering the 853 municipalities of MG. We compute the gap of each medical specialty for each municipality, considering both, its demand and the assignment of specialists and equipment to the 63 medical centers. Figure A1 to Figure A5 propose assignment of equipment, while Figure A6 to Figure A14 suggest the number of specialists for each municipality. Figure A1 MRIs procurement for scenario [3]. Figure A2 CTs procurement for scenario [3]. Figure A3 Mammograms procurement for scenario [3]. Figure A4 Ultrasound procurement for scenario [3]. Figure A5 Ultrasound Doppler procurement for scenario [3]. Specialists hiring suggestion for scenario [3]. Figure A6 Andiology hiring suggestion for scenario [3]. Figure A7 Cardiology hiring suggestion for scenario [3]. Figure A8 Endocrinology hiring suggestion for scenario [3]. Figure A9 Gynecology hiring suggestion for scenario [3]. Figure A10 Mastology hiring suggestion for scenario [3]. Figure A11 Nephology hiring suggestion for scenario [3]. Figure A12 Ophthalmology hiring suggestion for scenario [3]. Figure A13 Pediatrics hiring suggestion for scenario [3]. Figure A14 Urology hiring suggestion for scenario [3]. provides the number of equipment and specialist per municipality for this scenario in detail.

Figure 7 presents the financial results of the overall annual costs of 64 medical centers, equipment maintenance, and a comparison with social cost. The medical centers’ costs include a maintenance cost for existing 27 health care units (SES Resolution 6563, 2018), the costs for expanding its capacity to meet patients’ requirements, and the costs of additional 37 medical centers. Following, the results present the annual depreciation costs of equipment. The overall state annual costs are approximate R$250 M, which is comparable to the overall annual social opportunity cost, of R$259 M.

Figure 7
Financial results of overall annual costs of 64 medical centers and equipment maintenance, and a comparison with social cost.

3.2. Scenarios [4-6] based on budget availability

Scenarios based on budget availability focus on service offers. The government sets the number of municipalities and the minimum number of equipment it desires to provide according to its financial conditions. The model satisfies the patients’ demand and government offer constraints. Therefore, we activate constraints 1.8 and 1.9. We present three scenarios (see Figure 8) based on requirements of 70, 60, and 50 medical centers and at least the number of equipment selected for the first scenario (maximum 180 km), that is 125 MRIs, 62 CTs, 134 mammograms, 952 ultrasounds, and 43 ultrasound Doppler. We describe the financial implications for scenario 4 [70 CEM]; evaluating the average patients’ displacement, and the optimal overall costs and compare the state’s costs with social costs.

Figure 8
Municipalities selected to operate 70, 60, and 50 medical centers. Since 27 health units (in black) are fixed and available, the state would have to install 43, 33, and 23 new CEM.

With 70 medical centers of scenario 4, MG patients’ displacement follows a normal distribution with N156,50 km. As presented in Figure 9 and Figure 10, the number and cost of additional equipment are relatively similar to scenario 1, with 131 MRIs, 78 CTs, 147 mammograms, 961 ultrasounds, and 70 ultrasound Doppler. Figure 10 also demonstrates that when the MG state offers 70 medical centers, its costs are higher than the social opportunity cost.

Figure 9
Municipalities (70) selected in black (27 CEAE) and blue (43 CEM). Patients’ average displacement of 156 km. All municipalities demand are covered.
Figure 10
Financial results of overall annual costs of 70 medical centers and equipment maintenance, and a comparison with social cost.

3.3. Scenario [7] based on patients’ displacement with no fixed infrastructure

So far, we have evaluated scenarios considering the available infrastructure. In the following scenario, we propose a maximum patient’s displacement of 120 km (or two hours) and we relax the constraint 1.10 that fixes medical centers to existing infrastructure. The goal in this scenario is to evaluate a re-design option of the MG secondary care network. Under this strategy, the MG state requires 41 municipalities with medical centers (see Figure 11). The solution requires the procurement of 119 MRIs, 61 CTs, 136 mammograms, 946 ultrasounds, and 41 ultrasound Doppler. Besides, the government should hire 60 angiologists, 124 cardiologists, 59 endocrinologists, 176 gynecologists, 61 mastologists, 90 nephrologists, 152 ophthalmologists, 199 pediatrics, and 80 urologists. MG patients’ average displacement takes 1.5 hours, following a normal distribution with N95,31 km.

Figure 11
Financial results of overall annual costs of 41 medical centers, equipment maintenance, and social cost.

From the 41 municipalities selected to receive medical centers, 12 municipalities already have the available infrastructure, that is Capelinha, Frutal, Governador Valadares, Janaúba, Januária, Leopoldina, Patos de Minas, Patrocínio, Pirapora, São João del Rei, São Lourenço, and Teófilo Otoni. Municipalities without services are assigned to one of the 41 municipalities. The model generates this strategy for all scenarios. The strategy is aligned with the Intermunicipal Health Consortiums and provided for in the organic law of the SUS (Law 8,080/1990) to assure integral care to the population of the connected municipalities (Brasil, 1990Brasil. (1990, September 20). Dispõe sobre as condições para a promoção, proteção e recuperação da saúde, a organização e o funcionamento dos serviços correspondentes e dá outras providências (Lei nº 8.080, de 19 de setembro de 1990). Diário Oficial da República Federativa do Brasil. Retrieved in 2019, October 10, from http://conselho.saude.gov.br/legislacao/lei8080_ 190990.htm
http://conselho.saude.gov.br/legislacao/...
).

3.4. Priority-base schedule of medical centers installation

A time-based schedule for creating new medical specialties centers depends on the government budget and its priorities, therefore, for the future health units, we suggest the priority-based schedule adopting the SVI index (Instituto de Pesquisa Econômica Aplicada, 2019Instituto de Pesquisa Econômica Aplicada – IPEA. (2019). IVS: Índice de Vulnerabilidade Social. Brasília. Retrieved in 2019, June 10, from http://ivs.ipea.gov.br/index.php/pt/planilha
http://ivs.ipea.gov.br/index.php/pt/plan...
) for the selected municipalities. We designated 70 municipalities of the previous scenarios from which 43 are new medical centers, and 27 are available health units. Table 3 presents the proposed sequence based on SVI for installing new medical centers and improving the existing health units’ capacity.

Table 3
Proposed sequence for installing medical centers in MG, and improve CEAE’s capacity.

The authors and the IT lab team developed a system to provide flexibility and automate the analysis of the secondary care planning process. The system provides the results of optimizations run, with the location of medical centers, the number of medical specialties, the assignment of patients' demand to medical centers on municipalities. Simultaneously, the analyst can evaluate the assignment of patients' demand for exams to medical centers and the number of the equipment for procurement within a maximum distance. The system, presented in Figure 12, enables scenario investigations by changing the quantity of new equipment and a maximum distance of coverage.

Figure 12
General view of the web-based system interface (text in Portuguese). Location and demand coverage of new 20 MRI scanners considering available MRI scanners on MG for a maximum distance of 150 km. (the connections show municipalities and services of the northern region of the state). Source: Nescon (2018)Nescon. (2018). Retrieved in 2018, June 18, from https://www.nescon.medicina.ufmg.br/
https://www.nescon.medicina.ufmg.br/...
.

4. Conclusion

This study proposed an optimization method for simultaneously setting out the location of medical centers and the equipment allocation for secondary care taking into account the tradeoff between patients’ desire for minimum displacement, and physicians’ inclination for working on settled MG municipalities. The proposed method considered existing infrastructure on 27 municipalities of MG and developed analysis based on patients’ displacement, government budget availability.

For patients’ displacement scenarios, we observed that it is possible to implement a solution where a patient may meet secondary care at most one hour and a half from its residence. For this scenario, the social opportunity cost is equivalent to government costs. For budget availability scenarios, we observed that the service level can increase in quality (with numerous medical centers and equipment), however, with 70 medical centers, for example, patients would take longer to reach such services, besides, the government costs would be higher than the citizens’ opportunity costs. Following, we presented a scenario that re-designs the secondary care network of MG and demonstrated that 41 municipalities would be enough to satisfy the state’s demand. Finally, the authors presented a table with 70 municipalities and a sequence for installing medical centers and increasing the capacity of the available infrastructure based on the SVI index. The table took into account 43 similar municipalities of optimization runs and the 27 municipalities of MG that already offer specialized care.

In general, the outcomes suggested that this integrated optimization modeling approach is a prospective method to endorse the integration of professionals of public health, medicine, and engineering, to support the decision-making process. We consulted specialists offer for decades and provided a deep financial evaluation of a wide system projected to ensure universality, equity, and long-term sustainability.

The developed system is a potential tool to provide knowledge-based policymaking. It is also fundamental to highlight that both descriptive models, as cost functions, as the mathematical model do not replace managers. Its’ contribution resides on aiding decision-makers interacting to find the best solution by consensus by evaluating scenarios.

Some limitations of the study are worth mentioning. First, physicians move among municipalities, therefore, the offer of medical specialties is not accurate when we evaluate the values by municipalities. Micro or macro-regions describes better medical specialties’ offer. Second, we did not consider the possibility of physicians of a municipality meet the requirements of other municipalities. This strategy would increase the dynamics of patients’ and physicians’ movement increasing considerably the model complexity. Such an analysis would have provided answers to questions of political aspects. Finally, a common adoption is to consider population coverage in static computations. Although it is desirable, we did not consider such a parameter since it trades-off maximum distance, but both parameters can be evaluated simultaneously in future works.

For future works, we recommend (i) a study about municipalities’ efficiency on providing secondary care using Data Envelopment Analysis (DEA) scores (Mitropoulos et al., 2013Mitropoulos, P., Mitropoulos, I., & Giannikos, I. (2013). Combining DEA with location analysis for the effective consolidation of services in the health sector. Computers & Operations Research, 40(9), 2241-2250. http://dx.doi.org/10.1016/j.cor.2012.01.008.
http://dx.doi.org/10.1016/j.cor.2012.01....
); (ii) evaluating further decades’ projection of both, population (Instituto Brasileiro de Geografia e Estatística, 2019Instituto Brasileiro de Geografia e Estatística – IBGE. (2019). Projeção da população do Brasil e das Unidades da Federação. Rio de Janeiro. Retrieved in 2019, October 10, from https://www.ibge.gov.br/apps/populacao/projecao/
https://www.ibge.gov.br/apps/populacao/p...
), and specialists’ offer. For this last, we recommend systems dynamics, since it is dependent on several causal variables with data not readily available. Finally, (iii) the decision of installing medical centers can be modelled by two-stage stochastic programming (Errarhout et al., 2016Errarhout, A., Kharraja, S., & Corbier, C. (2016). Two-stage stochastic assignment problem in the home health care. IFAC-PapersOnLine, 49(12), 1152-1157. http://dx.doi.org/10.1016/j.ifacol.2016.07.659.
http://dx.doi.org/10.1016/j.ifacol.2016....
; Zarrinpoor et al., 2017Zarrinpoor, N., Fallahnezhad, M. S., & Pishvaee, M. S. (2017). Design of a reliable hierarchical location-allocation model under disruptions for health service networks: a two-stage robust approach. Computers & Industrial Engineering, 109, 130-150. http://dx.doi.org/10.1016/j.cie.2017.04.036.
http://dx.doi.org/10.1016/j.cie.2017.04....
), where the first stage decision is the installation of medical centers, and on the second stage, the uncertain demand and offer of medical specialties are revealed.

Appendix A Equipment assignment and number of specialists for each municipality of MG.

This appendix presents the scenario [3] selection of medical centers and equipment considering the 853 municipalities of MG. We compute the gap of each medical specialty for each municipality, considering both, its demand and the assignment of specialists and equipment to the 63 medical centers. Figure A1 to Figure A5 propose assignment of equipment, while Figure A6 to Figure A14 suggest the number of specialists for each municipality.

Figure A1
MRIs procurement for scenario [3].
Figure A2
CTs procurement for scenario [3].
Figure A3
Mammograms procurement for scenario [3].
Figure A4
Ultrasound procurement for scenario [3].
Figure A5
Ultrasound Doppler procurement for scenario [3]. Specialists hiring suggestion for scenario [3].
Figure A6
Andiology hiring suggestion for scenario [3].
Figure A7
Cardiology hiring suggestion for scenario [3].
Figure A8
Endocrinology hiring suggestion for scenario [3].
Figure A9
Gynecology hiring suggestion for scenario [3].
Figure A10
Mastology hiring suggestion for scenario [3].
Figure A11
Nephology hiring suggestion for scenario [3].
Figure A12
Ophthalmology hiring suggestion for scenario [3].
Figure A13
Pediatrics hiring suggestion for scenario [3].
Figure A14
Urology hiring suggestion for scenario [3].
  • How to cite this article: Almeida, J. F. F., Pinto, L. R., Conceição, S. V., & Campos, F. C. C. (2019). Medical centers location and specialists’ allocation: a healthcare planning case study. Production, 29, e20190006. https://doi.org/10.1590/0103-6513.20190006.

References

  • Abreu, G. R. F. (2017). Acesso e utilização às tecnologias de diagnóstico por imagem pelo SUS no Estado da Bahia: o caso do tomógrafo computadorizado e equipamento de ressonância magnética (Dissertação de mestrado). Universidade Federal da Bahia, Salvador.
  • Agência Minas. (2019, February 22). Nota de esclarecimento Belo Horizonte: Secretaria de Estado de Saúde de Minas Gerais. Retrieved in 2019, September 12, from http://www.agenciaminas.mg.gov.br/noticia/nota-de-esclarecimento-2019-02-22-14-50-00-0300
    » http://www.agenciaminas.mg.gov.br/noticia/nota-de-esclarecimento-2019-02-22-14-50-00-0300
  • Ahmadi-Javid, A., Seyedi, P., & Syam, S. S. (2017). A survey of healthcare facility location. Computers & Operations Research, 79, 223-263. http://dx.doi.org/10.1016/j.cor.2016.05.018
    » http://dx.doi.org/10.1016/j.cor.2016.05.018
  • Almeida, J. F. F., Pinto, L. R., Conceição, S. V., Campos, F. C. C., & Miranda Júnior, G. (2018, December). A two-step optimization process for medical centers location and capacity allocation. In A. Leiras, C. A. González-Calderón, I. Brito Junior, S. Villa & H. T. Y. Yoshizaki (Eds.), Operations Management for Social Good: 2018 POMS International Conference in Rio (Springer Proceedings in Business and Economics). USA: Springer International Publishing.
  • Assembleia Legislativa do Estado de Minas Gerais – ALMG. (2018). Relatório comprova baixo investimento na saúde em Minas Belo Horizonte. Retrieved in 2019, September 29, from https://www.almg.gov.br/acompanhe/noticias/arquivos/2018/07/11_audiencia_saude_relatorio_SUS.html
    » https://www.almg.gov.br/acompanhe/noticias/arquivos/2018/07/11_audiencia_saude_relatorio_SUS.html
  • Bauer, J., & Groneberg, D. A. (2016). Measuring spatial accessibility of health care providers–introduction of a variable distance decay function within the Floating Catchment Area (FCA) method. PLoS One, 11(7), e0159148. http://dx.doi.org/10.1371/journal.pone.0159148 PMid:27391649.
    » http://dx.doi.org/10.1371/journal.pone.0159148
  • Brandeau, M. L. (2016). Creating impact with operations research in health: making room for practice in academia. Health Care Management Science, 19(4), 305-312. http://dx.doi.org/10.1007/s10729-015-9328-0 PMid:26003321.
    » http://dx.doi.org/10.1007/s10729-015-9328-0
  • Brasil. (1990, September 20). Dispõe sobre as condições para a promoção, proteção e recuperação da saúde, a organização e o funcionamento dos serviços correspondentes e dá outras providências (Lei nº 8.080, de 19 de setembro de 1990) Diário Oficial da República Federativa do Brasil. Retrieved in 2019, October 10, from http://conselho.saude.gov.br/legislacao/lei8080_ 190990.htm
    » http://conselho.saude.gov.br/legislacao/lei8080_
  • Brasil, Ministério da Saúde. (2015, October 2). Aprova critérios e parâmetros para o planejamento e programação de ações e serviços de saúde no âmbito do SUS (Portaria nº 1.631, de 1º de outubro de 2015) Diário Oficial da República Federativa do Brasil. Retrieved in 2017, September 28, from http://bvsms.saude.gov.br/bvs/saudelegis/gm/2015 /prt1631 01 10 2015.html
    » http://bvsms.saude.gov.br/bvs/saudelegis/gm/2015
  • Brasil. (2016, December 15). Altera o Ato das Disposições Constitucionais Transitórias, para instituir o Novo Regime Fiscal, e dá outras providências (Emenda constitucional nº 95, de 15 de dezembro de 2016) Diário Oficial da República Federativa do Brasil. Retrieved in 2019, October 10, from http://www.planalto.gov.br/ccivil_03/constituicao/emendas/emc/emc95.htm
    » http://www.planalto.gov.br/ccivil_03/constituicao/emendas/emc/emc95.htm
  • Brasil, Ministério da Saúde. (2019a). Sistema Único de Saúde (SUS): estrutura, princípios e como funciona Brasília. Retrieved in 2019, July 15, from http://portalms.saude.gov.br/sistema-unico-de-saude
    » http://portalms.saude.gov.br/sistema-unico-de-saude
  • Brasil, Ministério do Trabalho. (2019b). PDET/CAGED: Programa de Dissseminação das Estatísticas do Trabalho Brasília. Retrieved in 2019, June 10, from pdet.mte.gov.br/caged
    » pdet.mte.gov.br/caged
  • Church, R., & ReVelle, C. (1974). The maximal covering location problem. Papers in Regional Science, 32(1), 101-118. http://dx.doi.org/10.1007/BF01942293
    » http://dx.doi.org/10.1007/BF01942293
  • Costa, M. A., & Marguti, B. O. E. (2015). Atlas da vulnerabilidade social nos municípios brasileiros Brasília: IPEA.
  • De Neufville, R., & Scholtes, S. (2011). Flexibility in engineering design Cambridge: MIT Press. http://dx.doi.org/10.7551/mitpress/8292.001.0001
    » http://dx.doi.org/10.7551/mitpress/8292.001.0001
  • Departamento de Informática do SUS – DATASUS. (2018a). Cadastro Nacional de Estabelecimento de Saúde Brasília. Retrieved in 2018, September 10, from http://cnes.datasus.gov.br/
    » http://cnes.datasus.gov.br/
  • Departamento de Informática do SUS – DATASUS. (2018b). Sistemas e aplicativos Brasília. Retrieved in 2018, September 10, from http://datasus.saude.gov.br/sistemas-e-aplicativos
    » http://datasus.saude.gov.br/sistemas-e-aplicativos
  • Errarhout, A., Kharraja, S., & Corbier, C. (2016). Two-stage stochastic assignment problem in the home health care. IFAC-PapersOnLine, 49(12), 1152-1157. http://dx.doi.org/10.1016/j.ifacol.2016.07.659
    » http://dx.doi.org/10.1016/j.ifacol.2016.07.659
  • Galvão, R. D., Acosta Espejo, L. G., & Boey, B. (2002). A hierarchical model for the location of perinatal facilities in the municipality of Rio de Janeiro. European Journal of Operational Research, 138(3), 495-517. http://dx.doi.org/10.1016/S0377-2217(01)00172-2
    » http://dx.doi.org/10.1016/S0377-2217(01)00172-2
  • Gu, W., Wang, X., & McGregor, S. E. (2010). Optimization of preventive health care facility locations. International Journal of Health Geographics, 9(1), 17. http://dx.doi.org/10.1186/1476-072X-9-17 PMid:20298608.
    » http://dx.doi.org/10.1186/1476-072X-9-17
  • Güneş, E. D., Yaman, H., Çekyay, B., & Verter, V. (2014). Matching patient and physician preferences in designing a primary care facility network. The Journal of the Operational Research Society, 65(4), 483-496. http://dx.doi.org/10.1057/jors.2012.71
    » http://dx.doi.org/10.1057/jors.2012.71
  • Hakimi, S. L. (1964). Optimum locations of switching centers and the absolute centers and medians of a graph. Operations Research, 12(3), 450-459. http://dx.doi.org/10.1287/opre.12.3.450
    » http://dx.doi.org/10.1287/opre.12.3.450
  • Hale, T. S., & Moberg, C. R. (2003). Location science research: a review. Annals of Operations Research, 123(1), 21-35. http://dx.doi.org/10.1023/A:1026110926707
    » http://dx.doi.org/10.1023/A:1026110926707
  • Horngren, C. T., Bhimani, A., Datar, S. M., Foster, G., & Horngren, C. T. (2002). Management and cost accounting Harlow: Financial Times/Prentice Hall.
  • Instituto Brasileiro de Geografia e Estatística – IBGE. (2010a). Censo demográfico 2010: Rendimento mensal domiciliary per capita segundo as grandes regiões e unidades da Federação Rio de Janeiro. Retrieved in 2019, October 10, from https://censo2010.ibge.gov.br/sinopse/index.php?dados=P21&uf=00
    » https://censo2010.ibge.gov.br/sinopse/index.php?dados=P21&uf=00
  • Instituto Brasileiro de Geografia e Estatística – IBGE. (2010b). Censo demográfico 2010 Rio de Janeiro. Retrieved in 2019, September 10, from https://www.ibge.gov.br/estatisticas-novoportal/sociais/populacao/9662-censodemogra_co-2010.html
    » https://www.ibge.gov.br/estatisticas-novoportal/sociais/populacao/9662-censodemogra_co-2010.html
  • Instituto Brasileiro de Geografia e Estatística – IBGE. (2019). Projeção da população do Brasil e das Unidades da Federação Rio de Janeiro. Retrieved in 2019, October 10, from https://www.ibge.gov.br/apps/populacao/projecao/
    » https://www.ibge.gov.br/apps/populacao/projecao/
  • Instituto de Pesquisa Econômica Aplicada – IPEA. (2019). IVS: Índice de Vulnerabilidade Social Brasília. Retrieved in 2019, June 10, from http://ivs.ipea.gov.br/index.php/pt/planilha
    » http://ivs.ipea.gov.br/index.php/pt/planilha
  • Jornal Estado de Minas. (2016, 12 abril). Primeiro Centro de Especialidades Médicas de Minas é inaugurado em Pirapora Belo Horizonte. Retrieved in 2017, September 28, from https://www.em.com.br/app/noticia/gerais/2016/04/12/interna_gerais,752198/primeiro-centro-de-especialidades-medicas-de-minas-e-inaugurado-em-pir.shtml
    » https://www.em.com.br/app/noticia/gerais/2016/04/12/interna_gerais,752198/primeiro-centro-de-especialidades-medicas-de-minas-e-inaugurado-em-pir.shtml
  • Khodaparasti, S., Maleki, H., Jahedi, S., Bruni, M., & Beraldi, P. (2016). Enhancing community-based health programs in Iran: a multi-objective location-allocation model. Health Care Management Science, 20(4), 485-499. PMid:27115695.
  • Luo, J., Tian, L., Luo, L., Yi, H., & Wang, F. (2017). Two-step optimization for spatial accessibility improvement: a case study of health care planning in rural China. BioMed Research International, 2017, 2094654. http://dx.doi.org/10.1155/2017/2094654 PMid:28484707.
    » http://dx.doi.org/10.1155/2017/2094654
  • Macedo, L. M. D., & Martin, S. T. F. (2014). Interdependência entre os níveis de atenção do Sistema Único de Saúde (SUS): significado de integralidade apresentado por trabalhadores da Atenção Primária. Interface: Comunicação, Saúde, Educação, 18(51), 647-660. http://dx.doi.org/10.1590/1807-57622013.0597
    » http://dx.doi.org/10.1590/1807-57622013.0597
  • Minas Gerais. (2018, December 29). Altera o Art. 2° e o Anexo Único da Resolução SES/MG nº 6.495, de 13 de novembro de 2018, que autoriza o repasse de incentivo financeiro, em caráter excepcional, para apoio e fortalecimento dos Centros Estaduais de Atenção Especializada (CEAE) e do Centro de Especialidades Médicas (CEM) (Resolução SES/MG nº 6.563, de 21 de dezembro de 2018) Diário Oficial do Estado de Minas Gerais. Retrieved in 2019, August 15, from http://saude.mg.gov.br/images/documentos/RESOLUÇÃO_6563.pdf
    » http://saude.mg.gov.br/images/documentos/RESOLUÇÃO_6563.pdf
  • Minas Gerais, Secretaria de Estado de Saúde. (2019). SUS Belo Horizonte. Retrieved in 2019, June 10, from http://www.saude.mg.gov.br/sus
    » http://www.saude.mg.gov.br/sus
  • Mitropoulos, P., Mitropoulos, I., & Giannikos, I. (2013). Combining DEA with location analysis for the effective consolidation of services in the health sector. Computers & Operations Research, 40(9), 2241-2250. http://dx.doi.org/10.1016/j.cor.2012.01.008
    » http://dx.doi.org/10.1016/j.cor.2012.01.008
  • Neis, P., & Zipf, A. (2007). A web accessibility analysis service based on the OpenLS route service. In International Conference on Geographic Information Science of the Association of Geograpic Information Laboratories for Europe (AGILE 2007) Aalborg, Denmark.
  • Nescon. (2018). Retrieved in 2018, June 18, from https://www.nescon.medicina.ufmg.br/
    » https://www.nescon.medicina.ufmg.br/
  • Novignon, J., & Nonvignon, J. (2017). Improving primary health care facility performance in Ghana: efficiency analysis and scal space implications. BMC Health Services Research, 17(1), 399. http://dx.doi.org/10.1186/s12913-017-2347-4 PMid:28606131.
    » http://dx.doi.org/10.1186/s12913-017-2347-4
  • Organisation for Economic Co-operation and Development – OECD. (2017), Health at a Glance 2017: OECD indicators Paris: OECD Publishing. https://doi.org/10.1787/health_glance-2017-en
    » https://doi.org/10.1787/health_glance-2017-en
  • Paim, J., Travassos, C., Almeida, C., Bahia, L., & Macinko, J. (2011). The Brazilian health system: history, advances, and challenges. Lancet, 377(9779), 1778-1797. http://dx.doi.org/10.1016/S0140-6736(11)60054-8 PMid:21561655.
    » http://dx.doi.org/10.1016/S0140-6736(11)60054-8
  • Rais, A., & Viana, A. (2011). Operations research in healthcare: a survey. International Transactions in Operational Research, 18(1), 1-31. http://dx.doi.org/10.1111/j.1475-3995.2010.00767.x
    » http://dx.doi.org/10.1111/j.1475-3995.2010.00767.x
  • Ricardo, P. (2019). Centro de Especialidades Médicas amplia acesso aos serviços de saúde no Norte de Minas Belo Horizonte: Secretaria de Estado de Saúde de Minas Gerais. Retrieved in 2019, October 10, from http://www.saude.mg.gov.br/ngc/story/8735-centro-de-especialidades-medicas-amplia-acesso-aos-servicos-de-saude-no-norte-de-minas
    » http://www.saude.mg.gov.br/ngc/story/8735-centro-de-especialidades-medicas-amplia-acesso-aos-servicos-de-saude-no-norte-de-minas
  • Sathler, T. M., Almeida, J. F., Conceição, S. V., Pinto, L. R., & Campos, F. C. (2019). Integration of facility location and equipment allocation in health care management. Brazilian Journal of Operations & Production Management, 16(3), 513-527. http://dx.doi.org/10.14488/BJOPM.2019.v16.n3.a13
    » http://dx.doi.org/10.14488/BJOPM.2019.v16.n3.a13
  • Schilling, D., Elzinga, D. J., Cohon, J., Church, R., & ReVelle, C. (1979). The TEAM/FLEET models for simultaneous facility and equipment siting. Transportation Science, 13(2), 163-175. http://dx.doi.org/10.1287/trsc.13.2.163
    » http://dx.doi.org/10.1287/trsc.13.2.163
  • Shariff, S. S. R., Moin, N. H., & Omar, M. (2012). Location allocation modeling for healthcare facility planning in Malaysia. Computers & Industrial Engineering, 62(4), 1000-1010. http://dx.doi.org/10.1016/j.cie.2011.12.026
    » http://dx.doi.org/10.1016/j.cie.2011.12.026
  • Silva, I. C. M. D., Restrepo-Mendez, M. C., Costa, J. C., Ewerling, F., Hellwig, F., Ferreira, L. Z., Ruas, L. P. V., Joseph, G., & Barros, A. J. D. (2018). Measurement of social inequalities in health: concepts and methodological approaches in the Brazilian context. Epidemiologia e Serviços de Saúde, 27(1), e000100017. http://dx.doi.org/10.5123/S1679-49742018000100017 PMid:29513856.
    » http://dx.doi.org/10.5123/S1679-49742018000100017
  • Stummer, C., Doerner, K., Focke, A., & Heidenberger, K. (2004). Determining location and size of medical departments in a hospital network: a multiobjective decision support approach. Health Care Management Science, 7(1), 63-71. http://dx.doi.org/10.1023/B:HCMS.0000005399.23600.69 PMid:14977095.
    » http://dx.doi.org/10.1023/B:HCMS.0000005399.23600.69
  • Toregas, C., Swain, R., ReVelle, C., & Bergman, L. (1971). The location of emergency service facilities. Operations Research, 19(6), 1363-1373. http://dx.doi.org/10.1287/opre.19.6.1363
    » http://dx.doi.org/10.1287/opre.19.6.1363
  • Treurnicht, M. J., & Van Dyk, L. (2014). A decision support system for equipment allocation in a telemedicine referral network. South African Journal of Industrial Engineering, 25(1), 29-38. http://dx.doi.org/10.7166/25-1-641
    » http://dx.doi.org/10.7166/25-1-641
  • United Nations Development Programme – UNDP. (2018). Human development indices and indicators: 2018 statistical update New York.
  • Van der Veen, E., Hans, E. W., Veltman, B., Berrevoets, L. M., & Berden, H. J. (2015). A case study of cost-efficient staffing under annualized hours. Health Care Management Science, 18(3), 279-288. http://dx.doi.org/10.1007/s10729-014-9292-0 PMid:24997580.
    » http://dx.doi.org/10.1007/s10729-014-9292-0
  • Victora, C. G., Barreto, M. L., Leal, M., Monteiro, C. A., Schmidt, M. I., Paim, J., Bastos, F. I., Almeida, C., Bahia, L., Travassos, C., Reichenheim, M., & Barros, F. C. (2011). Health conditions and health-policy innovations in Brazil: the way forward. Lancet, 377(9782), 2042-2053. http://dx.doi.org/10.1016/S0140-6736(11)60055-X PMid:21561659.
    » http://dx.doi.org/10.1016/S0140-6736(11)60055-X
  • Zarrinpoor, N., Fallahnezhad, M. S., & Pishvaee, M. S. (2017). Design of a reliable hierarchical location-allocation model under disruptions for health service networks: a two-stage robust approach. Computers & Industrial Engineering, 109, 130-150. http://dx.doi.org/10.1016/j.cie.2017.04.036
    » http://dx.doi.org/10.1016/j.cie.2017.04.036

Publication Dates

  • Publication in this collection
    2 Dec 2019
  • Date of issue
    2019

History

  • Received
    02 Feb 2019
  • Accepted
    23 Oct 2019
Associação Brasileira de Engenharia de Produção Av. Prof. Almeida Prado, Travessa 2, 128 - 2º andar - Room 231, 05508-900 São Paulo - SP - São Paulo - SP - Brazil
E-mail: production@editoracubo.com.br