1 Introduction

Changes are part of daily life in business organizations. Competitiveness, pressure by stakeholders, sustainability issues, and supply chain variables, all force businesses to continually search for improvements (Arena et al., 2020; Jiang et al., 2018; Li & Wang, 2020). Several strategies have been adopted to keep or to improve business. Industrial relocation (IR) is a business strategy consisting of moving operations locations. The choice for an IR strategy must consider benefits and opportunities, and also risks and threats of such a decision (Chen et al., 2018). Reduction in costs is usually an outstanding opportunity with IR. Moreover, IR may lead to better infrastructure, lower land costs, and other possible benefits granted by some cities, regions, or countries (Li et al., 2019; Tang et al., 2020). IR reinforces complementary advantages among different regions and promotes coordinated regional economic development (Zhang et al., 2020). However, many other aspects should be analyzed, regarding the complexity and the socio-economic impacts of moving a factory location (Pan et al., 2019; Yang & Gallagher, 2017). IR plays a vital role in the promotion of economic growth, then IR has become a notable feature for the world’s economic development (Wu et al., 2014).

The purpose of this paper is to guide decision-makers how to assess alternatives for IR by using multi-attribute decision-making (MADM). With MADM, IR strategies may be assessed not only based on a single attribute, as costs, or profits, but based on several attributes (Ortiz-Barrios & Alfaro-Saiz, 2020). This paper presents the application of MADM methods in a real case of IR. The main motivation for this research was to verify the adherence of MADM to support decision-making on IR. The case of a metalworking plant located in Brazil, is a representative example, as presented in Sect. 4. Nevertheless, IR is not limited to a specific industry.

Four leading methods of MADM (Zavadskas et al., 2014) were applied: analytic hierarchy process (AHP), multi-attribute utility theory (MAUT), multi-attribute value theory (MAVT), and technique of order preference by similarity to the ideal solution (TOPSIS). AHP, MAUT, MAVT, and TOPSIS were applied with the same data, collected from a multinational company in the metal-mechanic industry. There is a trend for methods combination, in MADM literature (Ortiz-Barrios et al., 2020a). In this paper, MADM methods were not combined but fully applied and compared. Results from AHP, MAUT, MAVT, and TOPSIS applications were compared with compatibility indices (Garuti & Salomon, 2012) and correlation coefficients (Salomon & Rangel, 2015). Therefore, this paper presents more than a case study, combining mathematical modeling within a mixed research approach. With compatibility indices and correlation coefficients, it was expected to verify similarity among MADM methods. Compatibility is a measure of closeness (Lipovetsky, 2020), or similarity (Yoon, 2020), of decision vectors. Besides an emerging subject in AHP and MADM theories, literature on compatibility is deeply sparse. The joint analysis of compatibility and correlation of results from a real case is the major novelty of this paper.

The remainder of this paper is as follows. Section 2 presents a brief literature review on IR, mainly as a decision problem. Section 3 presents the research methodology, introducing MADM methods, compatibility indices, and correlation coefficients. Section 4 presents the case of MADM application in IR. Section 5 presents results from AHP, MAUT, MAVT, and TOPSIS, including their comparisons. Section 6 is dedicated to conclusions and an overview of the main findings.

2 Literature review

Industrial diversification (ID) and industrial relocation (IR) are two dominant types of strategic responses to industrial decline (Grillitsch & Asheim, 2018; Schamp, 2005). In ID, firms stay in the location but leave the industry. In IR, firms stay in the industry but leave the location. Therefore, IR is a strategy to increase business competitiveness by moving the location of its operations.

There are different tactics for IR, such as pseudo relocation, stratified relocation, and total relocation (Atakhan-Kenneweg et al., 2021; Brouwer et al., 2004). Pseudo relocation occurs by acquisitions or outsourcing of operations. Stratified relocation has been applied by manufacturers involving some segments of production (Johansson & Olhager, 2017; Liao & Chan, 2011). Total relocation is the focus of this paper. This IR tactic involves main driven forces including the need for more suitable premises, cost-saving, and policies incentives.

The main reasons for IR include environmental impacts (Aus Dem Moore et al., 2019; Koch & Basse Mama, 2019; Li & Wang, 2020; Li et al., 2019; Wang et al., 2019; Zhang et al., 2017), increasing land prices (Chen et al., 2018; Tang et al., 2020), manufacturing improvements (Jiang et al., 2018; Liu et al., 2021b), industrial restructuring (Kurata et al., 2020; Lundberg et al., 2016), shared value creation (Arena et al., 2020), and social-economic impacts (Indraprahasta et al., 2019; Pan et al., 2019; Pappas et al., 2018; Yang & Gallagher, 2017).

Several methods have been applied to IR, like linear regression (Pappas et al., 2018; Tang et al., 2020), Monte Carlo simulation (Li & Wang, 2020), non-linear regression (Wang et al., 2019), and percolation theory (Jiang et al., 2018). Those methods are pure mathematical modeling. MADM, instead requires human input, for instance to weight attributes. Human interaction is a major strength of MADM over regression and simulation models since this allows subjectiveness to be incorporated in problem-solving. MADM is gaining attention thanks to the possibility of tackling conflicting interests in very complex problems (Mardani et al., 2015b; Ortiz-Barrios et al., 2020a).

Applications of MADM methods were not found in relocation problems. There are applications of MADM for industrial location (Farahani et al., 2010), but not for relocation problems. In an extensive literature review of AHP applications (Vaidya & Kumar, 2006), “allocations” is a major theme. However, no relocation application was referred. Their seven references on this theme were delimited to the problem of first location design or selection. An upper-to-date AHP review (Emrouznejad & Marra, 2017) did not apart “allocation” as an independent category. Likewise, several AHP applications on first location selection were included in the literature review.

A review of literature reviews (Zavadskas et al., 2014), undelimited to AHP, identified only one reference to the industrial location (Melo et al., 2009). Despite citing IR as “facility relocations”, that reference did not associate IR with MCDM methods. Another upper-to-date literature review (Mardani et al., 2015a) also identified a single reference on the relocation problem. However, this reference presents AHP applied to a relocation problem of agricultural technology (Cay & Yuan, 2013), not of industrial concern.

Therefore, there is a research gap in IR literature, since MADM and MCDM are well applied in several managerial decision problems, with emphasis on supply chain management (Ishizaka & Nemery, 2013; Khan et al., 2018; Mardani et al., 2015b; Ortiz-Barrios et al., 2020a, 2020b; Tramarico et al., 2015). Then, the following research questions arise:

  • Does MADM help managerial decision processes regarding IR?

  • What MADM method can be applied to support IR?

  • Would applications of different methods of MADM result in the same decision?

Applications of MADM methods in a real case of IR are presented in Sect. 5. This is another novelty of this paper: MADM applied to IR. As previously observed, mathematical models such as linear or non-linear regression and simulation have already been applied in IR studies. One main disadvantage of those models compared to MADM is the requirement for historical data. With MADM, knowledge from experts may be incorporated into the analysis.

For MADM application, a set of attributes was collected from the literature cited in this section (Chen et al., 2018; Jiang et al., 2018; Li & Wang, 2020; Liao & Chan, 2011; Liu et al., 2021b; Pappas et al., 2018; Tang et al., 2020; Zhang et al., 2017):

  • Attractive policy incentives.

  • Environmental aspects.

  • Labor shortage.

  • Need to cut production costs.

  • Physical proximity of high technology parks to industrial upgrading.

A five-component set of attributes may look too small for a strategic decision. However, size is not the most desirable characteristic of such a set. Mutually exclusiveness and collectively exhaustiveness (MECE) are the most desirable features for a set of attributes (Lee & Chen, 2018). In this sense, the analysis of benefits, opportunities, costs, and risks (BOCR Model) has been successfully applied with AHP (Wijnmalen, 2007). With the BOCR Model, one can assess all factors to decision-making: The certain (B and C), uncertain (O and R), positive (B and O), and negative (C and R) merits (Shih et al., 2014). Therefore, four attributes, exclusive and exhaustive could be better than more attributes not exclusive, or nor exhaustive.

3 Methodology

3.1 Generalities

Multiple criteria decision-making (MCDM) is a methodology that applies conflicting criteria in the solution of decision problems (Doumpos & Zopounidis, 2002). MCDM is divided into two branches (Farahani et al., 2010; Zavadskas et al., 2014): Multi-attribute decision-making (MADM) and multi-objective decision-making (MODM). Decision problems solved by MADM have a finite set of alternatives. In the extreme case, MADM deals with only two alternatives (Hwang & Yoon, 1981; Lombardi Netto, et al., 2020). MODM deals with larger sets of alternatives. In the extreme case, MODM deals with infinite alternatives. Therefore, this paper is on MADM, since a finite set of alternatives will be considered in an IR problem. This is due to the nature of IR problems. Usually, IR deals with a limited set of alternatives, for example, cities, regions, or countries.

This paper proposes the application of MADM to solve a decision problem of IR. Figure 1 presents the four-step research flow.

Fig. 1
figure 1

Research flow

For Step 1, IR may be defined as a “choice problem” (Ishizaka & Nemery, 2013; Roy, 1981). In this type of decision problem, the goal is to select the single best alternative, or a subset of alternatives, as restricted as possible. In IR, the relocation to a single place is desirable. Then, IR, as the choice problem, refers to decisions of “winner-take-all” (Lee et al., 2006).

Step 2 started with attributes identification in Sect. 2. These attributes were confirmed by the case presented in Sect. 4, where alternatives were also identified. As this research aimed the comparison of different MADM methods, these methods were applied with the same model, i.e., the same sets of attributes and alternatives.

AHP, MAUT, MAVT, and TOPSIS were the selected MADM methods, at first, by their leadership in MADM literature (Khan et al., 2018; Wallenius, et al., 2008). AHP, MAUT, MAVT, and TOPSIS were mainly selected because they are MADM methods developed for the solution of a choice problem (Ishizaka & Nemery, 2013). Besides some conceptual differences, all these methods result in decision vectors, \({\varvec{x}}\), where \({x}_{i}\) is the overall priority (AHP), joint utility (MAUT), aggregate value (MAVT), or closeness coefficient (TOPSIS) for Alternative \(i\). The best alternative is the one with the highest \({x}_{i}\).

Step 3 is presented in Sect. 5, with applications of AHP, MAUT, MAVT, and TOPSIS. Previously, in Sect. 4, the expert team is introduced. Four experts in the operations to be relocated provide input for MCDA. Results from methods application are presented in Sect. 5. Matter of fact, applications of different MADM methods may result in different decisions, even with the same data input (Liu et al., 2021a; Wang, 2015; Zamani-Sabzi et al., 2016; Zanakis et al., 1998; Zhang et al., 2015). Therefore, more than one MADM method will be applied and results compared, in the following step of this research.

Step 4 is presented in Sect. 5. The indication of the same best alternative from applications of two or more methods implies qualitatively similar results. Quantitative comparisons of results are also presented with two compatibility indices (Garuti and Saaty) and two correlation coefficients (Pearson and Spearman).

The remainder of this section presents concepts and details of MADM methods (Sects. 3.2 to 3.4), compatibility indices (Sect. 3.5), and correlation coefficients (Sect. 3.6).

3.2 Analytic hierarchy process

Developed by Saaty (1974, 1980), AHP is one the most applied MADM methods. There are AHP applications in diverse areas, like chemical engineering, computer science, ecology, energy sector, health sector, higher education sector, manufacturing, mathematical advances, supply chain management, and logistics (Emrouznejad & Marra, 2017; Tramarico et al., 2015).

The first step for AHP application is to model the decision problem in a hierarchy (Saaty & Rogers, 1976). Decision objective in the top, \(m\) attributes, usually named as criteria, in the middle of the hierarchy, and \(n\) alternatives in the bottom. Figure 2 presents an example of a hierarchy model for \(m=3\) and \(n=4\).

Fig. 2
figure 2

Three-level hierarchy for a decision problem with three attributes and four alternatives

The meaning of a hierarchical model is that the elements in a lower level need to be assessed regarding the elements in the higher level. Therefore, attributes shall be pairwise compared regarding decision objective, in this case, Industrial Relocation. Saaty Scale, more commonly referred to the Fundamental Scale of Absolute Numbers, as in Saaty (2013), is used in these pairwise comparisons. Saaty Scale is a linear 1–9 scale, with 1 for “equal importance”, 3 for “weak importance of one over another”, 5 for “strong importance”, 7 for “very strong or demonstrated importance”, and 9 for “absolute importance”. Intermediate values, as 2, 4, 6, and 8, and even rational numbers may be used, if needed (Saaty & Rogers, 1976).

Weights for the attributes, usually named as criteria priorities, are obtained normalizing the right eigenvector \({\varvec{w}}\) of the pairwise comparison matrix \({\varvec{A}}\), as in Eq. 1, where \({\lambda }_{\mathrm{max}}\) is its maximum eigenvalue.

$${\varvec{A}}{\varvec{w}}={\lambda }_{\mathrm{max}}{\varvec{w}}$$
(1)

Consistency checking is one of the great advantages of AHP against other MADM methods. A consistent pairwise matrix \({\varvec{A}}\) satisfies \({a}_{ij}={a}_{ik}{a}_{kj}\), resulting in \({\lambda }_{\mathrm{max}}=m\). Consistency index \(\mu \) is a measure of the consistency of a pairwise matrix, as in Eq. 2, where \(m\) is the number of attributes.

$$\mu =\frac{{\lambda }_{\mathrm{max}}-m}{m-1}$$
(2)

Consistency ratio \(CR\) is a better measure since it compares \(\mu \) with a random index \(RI\), computed by Oak Ridge Laboratory with more than 50,000 matrices (Saaty, 1980), as in Eq. 3.

$$CR=\frac{\mu }{RI}$$
(3)

Consistent matrices have \({\lambda }_{\mathrm{max}}=m\), then \(\mu =0\) and \(CR=0\). Inconsistent matrices have at least one comparison, and its reciprocal, \({a}_{ij}\ne {a}_{ik}{a}_{kj}\), resulting in \({\lambda }_{\mathrm{max}}>m\). It is desirable that \(CR\le 0.1\), then \({\varvec{A}}\) may be accepted, meaning “conformity with previous practice” or that decision-makers did not change their minds, when fulfilling a pairwise comparison matrix.

Alternatives shall be assessed regarding each attribute, resulting in their importance, likelihood, membership, or preference \({\varvec{x}}\). In AHP, these values of preferences are usually named as local priorities of alternatives (Saaty, 2013).

The same procedure may be used to obtain alternatives’ preferences regarding each attribute. This original procedure was named “relative measurement” (Saaty, 2013). Conversely, in “absolute measurement”, also named “ratings” (Saaty, 1986), alternatives are not pairwise compared to each other, instead of that, they are compared with standard levels (Salomon et al., 2016). Table 1 presents seven standard levels to assess alternatives regarding the attributes and proposal for their preference normalized in ideal synthesis (Salomon, 2016).

Table 1 Standard levels for relative measurement in AHP

One first limitation of AHP is the required independence of alternatives and attributes. Another limitation is on the total of alternatives and attributes. If.

\(m>9\), then attributes can be aggregated or grouped as sub-attributes which leads to one more level in the hierarchy. If \(n>9\), alternatives can be grouped, or absolute measurement may be applied.

Overall priorities of alternatives \({\varvec{y}}\) are obtained weighting local priorities of alternatives \({\varvec{x}}\) by weights of attributes \({\varvec{w}}\), as in Eq. 4.

$${\varvec{y}}={\varvec{w}}{\varvec{x}}$$
(4)

3.3 Multi-attribute utility and value theories

MAUT is indicated when decision-makers can construct utility functions (Ishizaka & Nemery, 2013; Oliveira & Dias, 2020). With MAUT a large number of quantitative and qualitative attributes can be evaluated (Chang, 2009). It is an MADM simpler method AHP since MAUT does not require pairwise comparisons (Doczy & AbdelRazig, 2017). Weights for attributes are directly assessed.

Preferences for alternatives come from utility functions \(U\), with codomain in interval \([\mathrm{0,1}]\). For each alternative, a marginal utility \({\varvec{u}}\) shall be determined regarding each attribute. For qualitative attributes, a five-point adapted Likert Scale, as in Table 2, can be adopted. In Table 2 a sixth point was added, with Null.

Table 2 Adapted Likert Scale for MAUT

This paper adopts linear utility functions, because they are simple to work, and because concepts as risks and uncertainty are not being considered. Then for instance, for three alternatives scored with 3, 4, and 5, from the Adapted Likert Scale, their marginal utilities are \(U\left(3\right)=0\), \(U\left(4\right)=0.5\), \(U\left(5\right)=1\), respectively.

MAVT is a variation of MAUT with value functions \(V\), instead of utility functions \(U\) (Oppio et al., 2018). \(U\) and \(V\) have the same codomain: interval \([\mathrm{0,1}]\). However, for each attribute, there will be always at least one alternative with \(U=0\). In MAVT, zero marginal values only result from null scores. Then, for Scores 3, 4, and 5, marginal values are \(V\left(3\right)=0.6\), \(V\left(4\right)=0.8\), and \(V\left(5\right)=1\), respectively.

3.4 Technique of order preference by similarity to the ideal solution

TOPSIS assesses the performances of alternatives through similarity with the ideal solution (Liern & Pérez-Gladish, 2020). As in MAUT and MAVT, weights for attributes are directly assessed, without pairwise comparisons. For qualitative attributes, Table 2 also can be used in TOPSIS. The last similarity with MAUT and MAVT is on the use of values for preferences of alternatives \({\varvec{x}}\). However, the decision matrix in TOPSIS has weighted values, as in Eq. 5, for \(i=1,\dots ,m\), and \(j=1,\dots ,n\).

$${x}_{ij}={w}_{j}{v}_{ij}$$
(5)

Besides the TOPSIS name referring to an ideal solution, this method also works with an anti-ideal solution, also referred to negative ideal solution (NIS). Positive ideal solution (PIS) \({{\varvec{a}}}^{+}\), and NIS \({{\varvec{a}}}^{-}\) can be obtained as in Eqs. 6 and 7, for \(j=1,\dots ,m\).

$${a}_{j}^{+}=\mathrm{max}({x}_{ij})$$
(6)
$${a}_{j}^{-}=\mathrm{min}({x}_{ij})$$
(7)

Then, Euclidean distances to NIS and PIS \({{\varvec{d}}}^{-}\) and \({{\varvec{d}}}^{+}\) must be obtained as in Eqs. 8 and 9, for \(i=1,\dots ,m\).

$${d}_{i}^{+}=\sqrt{\sum_{j=1}^{n}{({a}_{j}^{+}-{x}_{ij})}^{2}}$$
(8)
$${d}_{i}^{-}=\sqrt{\sum_{j=1}^{n}{({x}_{ij}-{a}_{j}^{-})}^{2}}$$
(9)

Finally, closeness coefficients \({c}_{i}\) are obtained as in Eq. 10, for \(i=1,\dots ,m\).

$${c}_{i}=\frac{{d}_{i}^{-}}{{d}_{i}^{-}+{d}_{i}^{+}}$$
(10)

When overall preferences for Alternative \(i\) are closer to PIS than NIS, then \({c}_{i}>0.5\).

3.5 Compatibility indices

Decision vectors are the main expected result with AHP, MAUT, MAVT, and TOPSIS applications. They have different names and concepts as “overall priority” for AHP, “joint utility” for MAUT, and “aggregate value” for MAVT, and “closeness coefficient” for TOPSIS. However, there are also some common features. Components ranging from 0 to 1 are the main mathematical similarity of AHP’s, MAUT’s, MAVT’s, and TOPSIS’s decision vectors.

As previously mentioned in Sect. 3, applications of different MADM methods may result in different decisions (Zanakis et al., 1998). On the other hand, not identical decision vectors may lead to the same decision. Mathematically, different vectors may be close or far. When two different vectors are close, they are compatible vectors (Saaty, 2013). Saaty Compatibility Index \(S\) between \({\varvec{x}}\) and \({\varvec{y}}\) is obtained as in Eq. 11, where \({a}_{ij}{=x}_{i}/{x}_{j}\), \({b}_{ij}={y}_{i}/{y}_{j}\), \({e}_{i}=1\), and \(\odot \) is the Hadamard Product operator, \({a}_{ij}\odot {b}_{ij}^{\mathrm{T}}= ={a}_{ij}{b}_{ji}\), for \(i=1,\dots ,m\) and \(j=1,\dots ,m\).

$$S=\frac{1}{{m}^{2}}{{\varvec{e}}}^{\mathrm{T}}{\varvec{A}}\odot {{\varvec{B}}}^{\mathrm{T}}{\varvec{e}}$$
(11)

The codomain of \(S({\varvec{x}},{\varvec{y}})\) is interval \([1,+\infty [\). Identical vectors \({\varvec{x}}={\varvec{y}}\) result in \(S=1\). Otherwise \(S>1\). For \(3\le m\le 5\), Saaty (2013) suggests that \(S\le 1.1\) indicates compatibility between \({\varvec{x}}\) and \({\varvec{y}}\).

Garuti Compatibility Index \(G\), between \({\varvec{x}}\) and \({\varvec{y}}\), is obtained as in Eq. 12.

$$G=\sum_{i=1}^{m}\left[\frac{\mathrm{min}\left({x}_{i},{y}_{i}\right)}{\mathrm{max}\left({x}_{i},{y}_{i}\right)}\frac{\left({x}_{i}+{y}_{i}\right)}{2}\right]$$
(12)

The codomain of \(G({\varvec{x}},{\varvec{y}})\) is interval \([\mathrm{0,1}]\). Identical vectors \({\varvec{x}}={\varvec{y}}\) also result in 1 for this compatibility index, or in \(G=1\). Otherwise \(G<1\). Ten-percent tolerance was also proposed for \(G\), as in Garuti (2019).

For \(n>5\), \(S\) is not indicated due to the sensibility to small components in \({\varvec{x}}\) and \({\varvec{y}}\). Garuti and Salomon (2012) presented two examples of validated AHP applications: distance of six cities to Philadelphia (Saaty, 1977) and electricity consumption by seven common household appliances (Whitaker, 2007). The example of distance estimation was also studied by Lipovetsky (2020), but not by Yoon (2020), which are two of the most recent works on compatibility indices. Lipovestky (2020) included indices \(G\) and \(S\), whereas Yoon (2020) limited to index \(S\).

For the distance estimation example, both indices \(G\) and \(S\) indicate compatible vectors. For the energy consumption example, \(S\approx 1.46\) and \(G\approx 0.92\). Therefore, indices \(G\) and \(S\) diverge on compatibility: they are compatible vectors, according to \(G\), but incompatible vectors, according to \(S\). AHP and real-distance vectors resulted in the same decision for the choice problem and the ranking problem. Garuti and Salomon (2012) concluded that \(G\) and \(S\) performed equally well for the example of distances to cities, but \(G\) performed better for the example of electricity consumption. Only two cases are too few to generalize one index superiority over another. This paper contributes to this literature with a new case.

3.6 Correlation coefficients

From Statistics, variance \(\mathrm{var}\left({\varvec{x}}\right)\) is a measure of one variable \({\varvec{x}}\) variation, obtained as in Eq. 13, where \(\overline{x}\) is the average, or arithmetic mean of \({x}_{i}\), for \(i=1,\dots ,m\).

$$\mathrm{var}\left({\varvec{x}}\right)=\frac{1}{m}\sum_{i=1}^{m}{({x}_{i}-\overline{x})}^{2}$$
(13)

Covariance \(\mathrm{cov }\left({\varvec{x}},{\varvec{y}}\right)\) is a measure of joint variation between two variables, \({\varvec{x}}\) and \({\varvec{y}}\) (Rice, 2007), obtained as in Eq. 14.

$$\mathrm{cov}({\varvec{x}},{\varvec{y}})=\frac{1}{m}\sum_{i=1}^{m}[({x}_{i}-\overline{x})({y}_{i}-\overline{y})]$$
(14)

Correlation is the degree of association between two variables (Asuero et al., 2006). Pearson correlation coefficient \(r\) between \({\varvec{x}}\) and \({\varvec{y}}\) is obtained as in Eq. 15.

$$r=\frac{\mathrm{cov }({\varvec{x}},{\varvec{y}})}{\sqrt{\mathrm{var}\left({\varvec{x}}\right)\mathrm{var}({\varvec{y}})}}$$
(15)

The codomain of \(r({\varvec{x}},{\varvec{y}})\) is interval \([-1, 1]\). Identical vectors \({\varvec{x}}={\varvec{y}}\) result in \(r=1\), referred to the perfect correlation between \({\varvec{x}}\) and \({\varvec{y}}\). Additionally, \(r=0\) means that there is no correlation between \({\varvec{x}}\) and \({\varvec{y}}\), and \(r<0\) results when there is a negative correlation between \({\varvec{x}}\) and \({\varvec{y}}\).

MADM deals with a finite number of alternatives. When the set of alternatives was little, for instance \(n=2, 3, 4\), there is a mathematical concerning with the use of \(r\). Correlation of too few observations may be unreliable. Therefore, that correlation should be checked by its statistical significance of the difference from zero, initially by calculating \({t}_{r}\), as in Eq. 16.

$${t}_{r}=\frac{r\sqrt{n-2}}{\sqrt{1-{r}^{2}}}$$
(16)

Then, it will be possible not only to indicate where such a correlation is higher or lower, but also it will be possible to test the quantitative reliability of the resulted value for \(r\). There are different statistical tests on the significance of correlation indices. This paper employs the Student’s \(t\)-distribution. In the two-tailed \(t\)-distribution, the significance level of \(\alpha =0.1\), results in \(t=2.92\). Therefore, for \({t}_{r}\ge 2.92\), \(r\) can be considered a reliable correlation index, with a probability at least equals to \(90\%\) (Venables & Ripley, 2002).

Spearman rank correlation coefficient \({r}_{s}\) between \({\varvec{x}}\) and \({\varvec{y}}\) is obtained as in Eq. 17, where \({p}_{i}=\mathrm{rank}({x}_{i}|{\varvec{x}})\), \({q}_{i}=\mathrm{rank}({y}_{i}|{\varvec{y}})\), for \(i=1,\dots ,m\).

$${r}_{s}=\frac{\mathrm{cov }({\varvec{p}},{\varvec{q}})}{\sqrt{\mathrm{var}\left({\varvec{p}}\right)\mathrm{var}({\varvec{q}})}}$$
(17)

4 Case study

A French industrial group has one of its major plants located in the municipality of Sao Paulo, the capital city of the state of Sao Paulo. This partially automated plant, with 400 employees, transforms 6000 tons of raw material into metal-mechanics products, annually. The group strategic plan identified the need for cost reductions of operations in Brazil. Therefore, the board of directors of the group considered relocating this plant.

The option of relocating the plant came from a sum of factors: The high cost of direct labor and respective fringe benefits; the high rental price of real estate for the plant site; understanding that the plant has reached the maximum operational efficiency, using all Lean Manufacturing concepts and tools. Another important factor was the decreasing competitiveness due to the launch of new products manufactured in China.

Table 3 presents the main attributes identified to this strategic decision.

Table 3 Attributes identified by the French group for industrial relocation

The five attributes identified in the literature review (Sect. 2) were combined in Attributes E, I, T, and W, with the “Need to cut production costs.” as a decision objective. E, I, T, and W is a four-element set of attributes. This set is mutually exclusive and collectively exhaustive since all elements of BOCR are considered. Environment is related to Risks, Infrastructure relates to Benefits, Transportation relates to Opportunities, and Workers relates to Costs.

Environment is a major issue nowadays. Brazil is at the center of several controversies. Forest fires and dam accidents have gained international repercussions in recent years. Fortunately, this is not the case of the studied company, which operates in the secondary sector of the economy. In other words, industries are not the greatest source of concern with Brazilian environmental laws. In that way, Infrastructure could be a more important criterion to guide an industrial relocation decision. However, as presented further in this section, every alternative considered for relocation has, at least, good infrastructure, in terms of electric power, gas and oil distribution, and also water and waste management.

Transportation is a major issue due to the Brazilian continental dimension. Ground transportation on paved highways is the primary modal for cargo and passengers. Rail and river transport is rarely used, as the air network is in full development. With more than 4000 airports and landing fields, Brazil has the second-largest network, only after the United States (Central Intelligence Agency, 2021).

Manpower for operating the plant only needs basic skills. They do not require a degree from a higher-level course. On the other hand, technical professionals, trained in the operation of machine tools, such as cutting and milling, are highly desirable.

Figure 3 presents a blank map of Brazil, with states borders. Four locations were identified from the company’s internal reports and surveys, referred by their airport codes: Recife (REC), Rezende (QRZ), Sao Jose dos Campos (SJK), and Sao Paulo (GRU). GRU and REC are two Brazilian top-ten populous cities. The other eight cities are included in red. The major Brazilian industrial districts are located in those cities.

Fig. 3
figure 3

Source: Wikimedia Commons, 2008)

Alternatives for industrial relocation (

The current location (GRU) is considered as an alternative. Still, this is a case of industrial relocation and not a case of industrial diversification. After all, even if the company chooses GRU, their products will be the same as the metal-mechanics industry.

Table 4 presents an assessment based on the company’s reports on the four alternatives, regarding the four attributes.

Table 4 Alternatives assessment from company’s reports

A team of experts in the plant operations was identified to provide data for the MADM application. The team of experts was composed of four persons:

  • Industrial Director, responsible for all activities in South America four plants, 35 years of experience in the company.

  • Senior Manufacturing Manager, responsible for two plants in Brazil, 25 years of experience.

  • Senior Supply Chain Manager, responsible for all logistic operations in South America, 20 years of experience.

  • Senior Engineering Manager, 18 years of experience in multinational companies.

This is the first MADM or MCDM application in a strategic decision by the French industrial group in Brazil. Therefore, the team of experts in plant operations was firstly introduced to MCDM general concepts. Then, details in AHP, MAUT, MAVT, and TOPSIS were presented as their applications evolved.

This is not group decision-making (GDM), as defined by Mu et al. (2020). Therefore, instead of each member of the team of experts providing a single, data were provided by consensus, in a single meeting. Then, aggregation was not an issue. On the other hand, as presented in Sect. 5, consistency indices of pairwise comparison matrices were quite satisfactory, under the 0.10 threshold.

5 Results

5.1 Application of AHP

Table 5 presents a pairwise comparisons matrix of attributes for industrial location. Weights were obtained with the normalized right eigenvector of the comparison matrix.

Table 5 Comparisons and weights of attributes for industrial relocation

The comparison matrix presented in Table 5 may be accepted, since its \({\lambda }_{\mathrm{max}}\approx 4.11\),

\(\mu \approx 0.040\), and \(CR=0.044.\) Table 6 presents local and overall priorities of alternatives for industrial relocation. Local priorities are obtained associating assessment from Table 4 with preferences from Table 1.

Table 6 Local and overall priorities of alternatives for industrial relocation

AHP application resulted in the highest overall priority for GRU. Then, the decision is not to relocate the plant from Sao Paulo (GRU). The second-best decision, relocate to Sao Jose dos Campos (SJK) is only 0.031 overall points behind GRU. Relocate to Recife (REC) or Rezende (QRZ) are also good decisions, with overall priorities higher than 0.8 and 0.7, respectively.

5.2 Applications of MAUT and MAVT

For MAUT, MAVT and TOPSIS applications weights of attributes should be the same as AHP’s priorities. After all, this paper aims to compare results from the applications of different MADM methods.

Table 7 presents scores of alternatives for industrial relocation, obtained associating assessment from Table 4 with scores from the adapted Likert Scale (Table 2).

Table 7 Scores of alternatives for industrial relocation

Table 8 presents marginal and joint utilities of alternatives for industrial relocation, obtained with linear utility functions of scores presented in Table 7.

Table 8 Marginal and joint utilities of alternatives for industrial relocation

MAUT application resulted in the highest joint utility for GRU. Therefore, the same decision from the AHP application: not to relocate the plant from Sao Paulo. Relocate to Sao Jose dos Campos (SJK) comes second, and relocation to Recife (REC) and Rezende (QRZ) had joint utilities lower than 0.5. This means that relocating to QRZ or REC might be not satisfying, considering the whole set of attributes.

Table 9 presentsmarginal and aggregate values of alternatives for industrial relocation. Marginal values \({v}_{ij}\) obtained normalizing scores \({s}_{ij}\) presented in Table 7, as \({v}_{ij}={s}_{ij}/{\mathrm{max}}_{i=1}^{m}{s}_{ij}\).

Table 9 Marginal and aggregate values of alternatives for industrial relocation

MAVT application results in the highest aggregate for GRU. Therefore, the same decision from AHP and MAUT applications: not to relocate the plant from GRU. Relocate to SJK comes second, relocate to QRZ in third, and relocate to QRZ in fourth and last. With all aggregate values higher than 0.5, relocation to all these cities may be valuable, considering the whole set of attributes.

5.3 Application of TOPSIS

Table 10 presents the decision matrix for TOPSIS, obtained with the decision matrix for MAVT with values weighted by the attribute’s weights.

Table 10 Decision matrix for TOPSIS

Table 11 presents Euclidean distances to NIS and PIS, respectively \({{\varvec{d}}}^{+}\) and \({{\varvec{d}}}^{-},\) and also presents closeness coefficients \({\varvec{c}}\) of alternatives for industrial relocation.

Table 11 Euclidean distances to NIS and PIS, and closeness coefficients of alternatives for industrial relocation

TOPSIS application resulted in the highest closeness coefficient for GRU. Therefore, the same decision from AHP and MAUT applications: not to relocate the plant from GRU. Relocate to SJK comes second, relocation to QRZ and REC had closeness coefficients lower than 0.5. This means that relocating to QRZ or REC might be not valuable, considering the whole set of attributes.

5.4 Qualitative comparison of MADM applications

Table 12 summarizes the results of AHP, MAUT, MAVT, and TOPSIS applications. GRU’s highest priority, utility, value, or closeness are highlighted in italics. QRZ’s and REC’s lower than 0.5 utility and closeness had bold because they are alternatives to be avoided according to MAUT and TOPSIS.

Table 12 Results of different MADM methods applications to industrial relocation

The same decision results from the four MADM methods: not to relocate the industrial plant from GRU. AHP, MAUT, and MAVT applications resulted in the same ranks: first GRU, second SJK, third QRZ, and fourth REC. With TOPSIS, REC was ranked third, and QRZ ranked fourth.

Summarizing all four methods led to the same decision. Almost the same ranks resulted from the four applications. AHP and MAVT were more benevolent with the last alternatives than MAUT and TOPSIS.

5.5 Quantitative comparison of MADM methods

Table 13 presents Saaty Compatibility Indices \(S\) between the results of AHP, MAUT, MAVT, and TOPSIS. Indices \(S\) between AHP and MAVT, and between MAUT and TOPSIS are highlighted in italics, because they are lower than the upper limit of 1.1, indicating compatibility between those pairs of methods.

Table 13 Saaty compatibility indices of different MADM methods applications to industrial relocation

Table 14 presents Garuti compatibility indices \(G\) between the results of AHP, MAUT, MAVT, and TOPSIS. Indices \(G\) between AHP and MAVT, and between MAUT and TOPSIS are also highlighted in italics, because they are greater than the lower limit of 0.9.

Table 14 Garuti compatibility indices of different MADM methods applications to industrial relocation

According to Tables 13 and 14, results from AHP and MAVT, and from MAUT and TOPSIS are compatible. Even resulting in the same decision of not relocating from Sao Paulo, results from AHP and MAUT, from MAUT and MAVT and from AHP and TOPSIS are indicated as not compatible by both indices \(G\) and \(S\).

Table 15 presents Pearson correlation coefficients \(r\) between the results of AHP, MAUT, MAVT, and TOPSIS. Coefficient \(r\) between AHP and TOPSIS is highlighted in italics because it is the only lower than 0.9.

Table 15 Pearson correlation coefficients of different MADM methods applications to industrial relocation

Table 16 presents \({t}_{r}\) values computed for \(r\) presented in Table 15.

Table 16 Testing results for Pearson correlation coefficients

As in Table 15, less correlated pair AHP–TOPSIS had the lower value for \({t}_{r}\), and the only \({t}_{r}<2.92\). Therefore, for AHP–TOPSIS correlation is the only not reliable, as presented in Sect. 3.5.

Table 17 presents Spearman rank correlation coefficients \({r}_{s}\) between the results of AHP, MAUT, MAVT, and TOPSIS. Coefficients \(r\) between TOPSIS and other methods are highlighted in italics because they are lower than 0.8.

Table 17 Spearman rank correlation coefficients of different MADM methods applications to industrial relocation

According to results from AHP, MAUT and MAVT are perfectly correlated, resulting in the same ranks. Results from TOPSIS are poorly correlated with the others, in this case.

5.6 Managerial implications

The applications of all four methods resulted in the same decision for the company: not to relocate from GRU. Moreover, this indication came from a multi-attribute model analyzed by different methods, i. e., considering diverse concepts such as Euclidean distance, priority, utility, and value. This unanimity reinforces and justifies the decision-making.

External and internal factors led the company’s top management to consider industrial relocation. Therefore, four major attributes were identified in the international literature for IR decisions: environmental, infrastructure, transportation, and workers. Company top managers assessed the attributes in a single meeting. With AHP their opinion was proved to have consistency, which was very welcomed by them. Therefore, managers are interested in applying AHP for other strategical decisions.

After collecting data from feasible locations, the decision of not relocating surprised managers, at a first glance. However, with the application of four different methods resulting the same, managers were convinced that was not the time for a strategic change as the industrial relocation.

After the application of MADM methods, but before the writing of this paper, the company decided to implement the decision. They announced to keep the plant located in Sao Paulo.

6 Conclusions

This paper presents a real-world case of industrial relocation (IR) solved by the applications of four MADM methods. The pre-COVID-19 case considered the relocation of an industrial plant in Sao Paulo to three other locations. However, considering attributes from IR literature, all four MADM methods result in not moving the plant from Sao Paulo. All methods resulted in Sao Jose dos Campos as the second location. Recife and Rezende resulted as good options for AHP and MAVT, but not for MAUT and TOPSIS.

The alternative of not relocating from Sao Paulo had the highest overall results as 0.914 overall priority with AHP, 0.820 joint utility with MAUT, 0.880 aggregate value with MAVT, and 0.723 closeness coefficient to the ideal solution with TOPSIS. On the other hand, alternatives to relocate to Recife or Rezende had joint utilities far lower than 0.5, indicating alternatives to be not considered, if the first alternative becomes unavailable.

Results were compared with compatibility indices and correlation coefficients. Results with AHP and MAVT were indicated as compatible and correlated. Results from TOPSIS had a lower compatibility correlation with results from other methods.

Hybridism is a trend in MADM and MCDM. This paper presents the use of compatibility indices and correlation coefficients, with concepts and indicators to promote discussion on results from different methods. Matter of fact, applications of different methods resulted in the same decision: Not relocate the plant. Ranks of the alternatives were also considered, even being not important result in this case. This consideration just emphasized the similarity of the results. Moreover, this analysis may be useful for future research.

Applications of new methods of MADM, as Best Worst Method (BWM), Complex Proportional Assessment (COPRAS), Potentially All Pairwise Rankings of all Possible Alternatives (PAPRIKA), or Visekriterijumska Optimizacija I Kompromisno Resenje (VIKOR) emerge as first sources for researches. Incorporating techniques as Delphi and Fuzzy may be useful if more alternatives or uncertainty were present in the decision problem. Additionally, it can be considered the inclusion of other attributes, dealing with political or social issues. Finally, the proposed methodology may be replicated in other industries.