Multicriteria Approach towards the Sustainable Selection of a Teahouse Location with Sensitivity Analysis

: The teahouse market has seen an expansion across various countries. In order to identify the most reasonable paths for development, the choice of location for the outlets needs to account for a number of conﬂicting criteria. Therefore, the multicriteria approach is required to effectively handle the location selection problem. In this paper, we develop a multicriteria framework for teahouse selection and apply it in the context of Lithuania. The indicator system is set up in order to capture the different aspects of the candidate locations. We also apply two multicriteria decision-making techniques (the evaluation based on distance from average solution (EDAS) method and the weighted aggregated sum product assessment with normalization (WASPAS-N) method) in order to ensure the robustness of the results. The weights of criteria were determined based on the expert survey. In addition, a Monte Carlo simulation was applied to check the sensitivity in changes of the criterion weights. The empirical application demonstrated validity of the proposed approach in choosing the optimal location of a teahouse.


Introduction
The hospitality sector is an important contributor to economic growth as it creates jobs, generates income, and facilitates international connections in multiple ways [1][2][3][4]. This sector is even more important for countries without fossil resource endowments. Lithuania may be a good case for the promotion of the hospitality sector. The food and beverage service is an important component of the hospitality sector. According to Statista [5], the turnover of the Lithuanian food and beverage service sector reached 512.9 million Euros as of 2014. The corresponding figure for 2013 was 474.7 million Euros. What is more, at the end of 2015, 980 enterprises for the manufacturing of beverages were in operation in Lithuania, with 20.9% of them being individual. The period of 2011-2015 saw an increase in the total number of enterprises by 16.3%, whereas the number of individual companies rose by 1.5% [6]. In order to ensure the further expansion of the Lithuanian hospitality industry, new forms of hospitality businesses need to be promoted. The teahouse could be a new niche for the hospitality sector. Statista [5] also reported the compound annual growth rate of 2.6% for the European tea market during the period of 2018-2023. Turning to Lithuania, the revenue generated in the tea subsector amounted to more than 31 million Euros in 2018 [7]. In addition, there has been a decline in the coffee market as the revenue generated there decreased by 4.6% according to Nielsen [8].
The teahouse can be considered as an appealing option for the hospitality sector development in Lithuania for the market for specialized teahouses in Vilnius has not been saturated yet. For instance, there are just nine teahouses in Vilnius as opposed to the 66 coffeehouses there, according to the results of a Google search. On the contrary, it is known that coffee market in Lithuania is dominated by chain cafes like Caffeine, Vero café, Presto, and so on. The case of Starbucks also showed that the coffee market in Lithuania is already saturated, as this brand decided not to come because of enormous competition [9]. The trend of the development of teahouses is indeed prevailing in Europe. People are becoming more concerned about the status of their health, and they are more prone to consider their food and drink rations. Some studies have suggested that tea is a healthier option and can cause a reduction in the probability of heart diseases as suggested by a survey by Gardner et al. [10]. What is more, it can interact with cancer according to Ek et al. [11]. Following the trend of healthier living is a great opportunity for entrepreneurs and investors. Therefore, the development of methodologies for the choice of a teahouse location can be of benefit to future investors and entrepreneurs.
The location selection is a multicriteria problem, which requires appropriate methodology. Multicriteria decision-making (MCDM) techniques can be applied to handle problems of location selection, taking multiple conflicting criteria into consideration. Location selection problems focused on infrastructure objects (airports, service facilities) have been addressed by means of MCDM techniques: Sennaroglu and Celebi [12] applied PROMETHEE (Preference Ranking Organization METHod for Enrichment of Evaluations) and VIKOR (Vlsekriterijumska Optimizacija I KOmpromisno Resenje) techniques for airport location selection; Nguyen Anh and Hoang Thi [13] utilized technique for order of preference by similarity to ideal solution (TOPSIS) for choosing the location of a distribution center; hybrid approaches involving analytic hierarchy process or analytic network process were applied by Chang et al. [14], Ertugrul and Karakaşoglu [15], Hu et al. [16], and Yıldız and Tüysüz [17]. The hospitality sector has also been analyzed in this regard by applying linguistic reasoning [18], fuzzy logics [19,20], non-compensatory measures [21], and hybrid MCDM techniques [22]. However, the problem of teahouse location selection has not received appropriate consideration in the literature. In order to promote the development of the hospitality sector in this regard, this paper proposes an MCDM-based framework.
The research aims to propose a framework based on MCDM methods for effectively handling the problem of location selection for a teahouse. The following research goals were implemented: (1) to analyze the main concepts underlying the MCDM methods; (2) to identify the specific methods to be adapted to the teahouse location selection problem; (3) to determine the criteria for solving the problem of the location selection; and (4) to apply the chosen methods when implementing the model for effective solution of the problem of the teahouse location selection. The criteria for location selection were determined following the literature review. The case of Vilnius (Lithuania) is considered in the empirical example with alternative locations chosen from the real estate database. The weights of criteria were set on the basis of the expert survey and Monte Carlo simulation.
The two MCDM techniques, namely the evaluation based on distance from average solution (EDAS) method and the weighted aggregated sum product assessment (WASPAS) method, were applied for the ranking of the alternatives. The use of these techniques allows the checking of the robustness of the results with regards to different normalization principles and utility functions. More specifically, the EDAS technique relies on the reference point approach, whereas the WASPAS technique applies value functions (additive and multiplicative utility functions). The proposed framework may help establish businesses in the tea market by taking into account the multiple alternatives and expert assessments. It is noteworthy that the proposed approach may be adapted to different contexts and applied in different locations. Methodologically, we revised the WASPAS technique by involving an additional stage of normalization when aggregating the additive and multiplicative utility functions.
The paper proceeds as follows. Section 2 comprises a review of the development and applications of the MCDM techniques. In this part, the theory and applications of MCDM for the location selection are discussed in detail. Section 3 focuses on an overview of the specific MCDM methods chosen for solving the location selection problem and their theoretical aspects. Section 4 presents the main results. Finally, Section 5 presents the concluding remarks.

The Preliminaries for MCDM and Its Application for the Location Selection Problems
This section comprises the two main parts. The first subsection is devoted to a more general discussion on the nature of the MCDM techniques and the main theoretical principles to be followed. The second subsection discusses the possibilities of the application of the MCDM techniques for solving the location selection problem, which is a focal point of this research.

MCDM Methods
In this subsection, we briefly present the major strands in MCDM. MCDM problems have several common features. First, MCDM considers multiple objectives represented by the corresponding criteria, and the decisionmakers are allowed to establish conflicting relationships among these criteria. Second, the criteria to be considered are expressed in different dimensions (units of measurement), which does not allow for a trivial aggregation of the decision matrix. Third, the MCDM problems consider multiple alternatives that represent socio-economic decisions at different levels. The presence of the aforementioned properties implies that one has to handle multidimensional problems. In this case, decision making can rely on different types of the MCDM methods [23]. The MCDM methods are primarily aimed at quantifying the underlying utilities and ranking the available alternatives. However, different MCDM methods can result in different outcomes (i.e., the same set of alternatives might be attributed with different rankings). This can be attributed to different mathematical operations employed by the considered methods. Therefore, the issue of identifying the most suitable MCDM method for a particular case still exists [24].
The MCDM methods can be further divided into multi-objective decision-making (MODM) methods and multi-attribute decision-making (MADM) methods [25]. The MODM methods deal with continuous optimization problems (i.e., the best alternative is chosen from infinite set of alternatives). However, integer programming might be applied to a discrete set of alternatives. Conversely, the MADM methods can only handle discrete sets of alternatives. However, the literature usually refers to the MADM methods as the MCDM methods (we follow this strand in the other sections of this paper too).
The MADM methods can be grouped into those based on the value function, reference point, and outranking [26]. The techniques based on the value function include simple additive weighting [27], weighted product method [28], multi-attribute utility theory (MAUT) [29], and the weighted aggregated sum product assessment (WASPAS) [30]. These techniques aggregate the normalized values of the decision matrix by taking the weighting associated with the criteria into consideration. The analytic hierarchy process (AHP) by Saaty [31] can also be considered as a value function-based approach [32]. The technique for order of preference by similarity to ideal solution (TOPSIS) [33] and VIKOR [34] can be given as examples of the reference point-based approaches. In this case, the distances among the alternatives and ideal solutions (reference points) are measured. TOPSIS relies on Euclidean distance, whereas VIKOR involves Manhattan and Chebyshev distances. Certain generalizations are possible. For instance, the Multi-Objective Optimization by Ratio Analysis (MOORA) technique [35] involves both the reference point approach and value function approach. The evaluation based on distance from average solution (EDAS) approach [36] considers the average solution as the reference point. Finally, the outranking approaches include methods based on preference relationships (e.g., ELECTRE (ELimination Et Choix Traduisant la REalité) [37] or PROMETHEE [38]).
As regards the continuous optimization and the MODM methods, Marler and Arora [39] presented a survey on MODM techniques. The different utility functions can be employed to aggregate the objective functions in MODM. Therefore, one can arrive at a weighted-sum model or a weighted-product model, among other formulations. Lexicographic models define the ordering of the objectives without quantitative considerations. The objective functions can also be defined in a number of ways. For instance, Romero et al. [40] linked the goal programming, compromise programming, and reference point formulation. Sawik [41][42][43] applied weighted-sum and lexicographic approaches in connection with mathematical programming for decision analysis.
For more detailed reviews on the use of different MCDM techniques one can consult papers by Kumar et al. [25], Diaz-Balteiro et al. [44], Zavadskas, Turskis [45], Song et al. [46,47], and Marler and Arora [39]. As we focus on the location selection for a teahouse, the set of alternatives is finite and determined by the supply of the premises. Therefore, we apply the MADM approach, and the term MCDM will refer to the MCDM methods in the following sections. Among the possible techniques for MADM, there is a need for choosing the methods that allow the aggregating of the decision information by applying different principles (i.e., normalization, utility functions) in order to ensure the robustness of the analysis.

Location Selection and MCDM
There has been a body of literature on the application of MCDM techniques in the domain of location selection [48,49]. The MADM was applied by implementing different techniques for discrete sets of alternatives (locations). The fuzzy logic and preference relations were involved into the choice of the distribution center location by Chen [48] in order to account for uncertainties in the decision information. The triangular fuzzy numbers were applied to map the linguistic variables to the fuzzy values. The fuzzy numbers were used to represent the ratings, weights of criteria, and the resulting utility values. The fuzzy preference matrix was established in order to compare the alternatives via the step-wise approach. Cagri Tolga et al. [49] tackled the retail location selection problem by proposing a fuzzy MCDM procedure relying on the analytic network process. The application of the fuzzy MCDM in the hospitality sector can be found in the study by Chou et al. [19]. The latter study offered a fuzzy MCDM procedure for choosing the hotel location considering such variables as geographical properties, traffic conditions, hotel requirements and performance in terms of operation management. Tzeng et al. [50] applied AHP and VIKOR for restaurant location selection. Ziemba et al. [51] and Ziemba [52] proposed a fuzzy PROMETHEE-based technique for selection of a wind farm.
The problem of location selections has also been addressed by means of the MODM. The integer programming model was also set up by Jovanovic [53] for optimization of the electricity grid. Specifically, the decision variables included both the location and capacity of the transformer. The parameters of the model included investment costs (per annum), voltage fluctuations, load of feeders, and substations. The data envelopment analysis was combined with the integer programming model [54] in order to facilitate the identification of the most promising location for investment and implementation of a project. Indeed, data envelopment analysis defined the Pareto efficient frontier based on multiple criteria, whereas the integer programming problem can identify the most desirable combination of the decisions with regards to the efficiency frontier. Another application of the integer programming problem can be found in the study by Bhaumik [55]. The latter study aimed to identify the nodes of a distribution network to be shut down given the locations of the distributors and retailers. Tzeng and Chen [56] operationalized the choice of number and locations of fire stations within an international airport via the fuzzy multi-objective model. In order to tackle the complexity surrounding the fuzzy problem, a genetic algorithm was applied and juxtaposed to the enumeration method. Arora and Arora [57] and Serrano-Hernandez et al. [58] applied mathematical programming for plant location problems in different contexts. Ho et al. [59] combined AHP and multi-choice goal programming for solving the problem of restaurant location choice.
There has also been a body of literature focusing on the outlet (e.g., store or restaurant) location problems. The main principles underlying the issue of the retail location problem have been discussed by Nelson [60]. The rents attributed to stores were analyzed by Applebaum and Cohen [61] with regards to the store sites. The use of the multiplicative competitive interactive approach was advocated by Jain and Mahajan [62] when analyzing the competitive environment in the retail activities. The approach for measurement of the service quality (SERVQUAL) was applied by Hing and Lee [63] in the context of restaurant operation. The quality of the hospitality management services was assessed by Olsen and Roper [64]. The latter study also looked into the strategic management issues amidst the effects of globalization, multinationals, and corporate strategies. An approach integrating fuzzy numbers and the analytic hierarchy process was proposed by Kuo et al. [65] in order to select the most suitable location for a convenience store. The optimal location of the new franchises among the existing ones has been tackled by means of multi-objective integer programming [66]. Such criteria as customer flow and the presence of competing outlets were considered in their model. Therefore, the location selection problem is topical across different sectors and levels of aggregation. The problems of the location selection rely on different theoretical preliminaries and research methods. The use of the MCDM techniques is needed to effectively handle the multiple effects of the feasible solutions. In the context of the present study, the choice of the teahouse location can be formulated as an MCDM problem.

Framework for Multicriteria Choice of the Teahouse Location
Deciding which location is the best one for expanding a business is a complicated process, and it includes multiple conflicting factors [67]. From the techniques available for solving MCDM problems, one needs to choose the suitable ones for addressing particular problems. In this study, EDAS and WASPAS-N methods are to be used for solving the location problem. In this section, the two MCDM techniques will be presented in a more detailed manner. The main properties of the underlying utility functions will be discussed, which allow one to understand the major differences between these two techniques.
This study focuses on the case of teahouse location selection. This problem is a multicriteria one, as multiple conflicting objectives arise when comparing different locations. Specifically, the economic, social, and environmental criteria can be regarded as those determining the attractiveness of a certain location. The criteria defining the objectives of the business are also expressed in different dimensions. Finally, the importance of different objectives might vary. All these circumstances imply the need for multi-criteria analysis.
In this section, the framework for the multicriteria choice of a teahouse location is presented. First, the two MCDM methods are presented. Second, the approach for sensitivity analysis involving simulation of the random vectors of the criteria weights is discussed. This kind of sensitivity analysis enables one to look into the stability of the vector of ranks as the vector of weights is changed. Note that we also provide the weights of criteria based on an expert survey. Experts were given a questionnaire based on a Likert scale (scores from 1 to 5).

Indicator System
The identification of the best location is a critical step in cost/benefit analysis framework for any coffee/teahouse, restaurant or other food and drink service outlet [68]. Nevertheless, there exists a variety of various different goals underlying the problem of location selection. A teahouse not only provides a selection of different types of tea or snacks but is also important in the sense of communication and social interaction in general. The success of a teahouse is led by an appropriate location. It goes without saying, that the level of the convenience of service and, hence, the number of potential customers is critically related to the choice of the place of operation. Therefore, selecting the location of an outlet is of high significance when deliberating on buying or renting suitable premises for food or drink businesses. In most cases, the issue of selecting location can rest on the past experience or the opinion from experts [69].
To mitigate the shortcomings of intuition judgements, or even to control costs and indicate perks of alternative locations, it is handy to adopt a scientific method. The location selection of a teahouse is quite diverse from common selection decisions. This is because of the complexity of the issue; the problem deals with customer uniqueness, competitors, or even costly decoration. It is known that the tea market in Lithuania is not as saturated as the coffee market; however, there are still many issues that must be faced. Indeed, in order to select an eligible site for a teahouse business, both quantitative and qualitative criteria are required to be considered simultaneously [70]. For these reasons, this study presents a quantitative assessment tool to assist tea market players in choosing the most promising location for acquiring a property for a teahouse in terms of both quantitative and qualitative aspects.
Location selection is usually initiated once the need for additional capacity is realized [71]. This triggers the process of searching for the optimal solution with regard to the new premises. The factors that determine the overall utility of a certain alternative (i.e., location) generally fall in the following groups: market available, transportation network, labor supply, surrounding site, availability of raw materials and services, access to utilities, the presence of regulations, and the community environment. Nevertheless, these categories are not fully suitable for solving the teahouse location problem. Thus, it is clear that the issue to be tackled in this study is a bit more specific.
The inflow of customers and the prospective growth of the business is impacted by such factors as access to the public transport (buses, trolleybuses) and the availability of parking spaces for private transportation near the teahouse, landscaping, distance from public facilities, and the size of the property and commercial area. This research, thus, seeks to develop a MCDM framework for ranking of the alternative locations for a teahouse. The framework includes the following: decision criteria, alternative locations, and weighting of the criteria. The empirical example considered is the teahouse location problem for a teahouse in Vilnius, Lithuania. The literature review was carried out in order to identify the criteria for the empirical research. The resulting set of indicators is outlined in Table 1. Table 1. Criteria for the choice of a teahouse location.

Criteria Sources
1. Rent cost [50] 2. Property area [59,72] 3. Distance to scenery [73] 4. Public transportation [50] 5. Pedestrian flow [50,65] 6. Parking capacity [50,65] 7. Number of competitors [50,65] 8. Number of crimes in the surrounding area [19,74] 9. Distance from public facilities [19,50] 10. Outdoor advertisement [50,65,73] 11. Distance from garbage containers [50] It is well known that the cost of property and its area are the main criteria for location selection; nevertheless, the decision also depends on many other criteria. In this paper, the cost of the property is expressed as the amount of money that must be spent to rent the space per month, and the property area is expressed as the inner space of the available area for rent (m 2 ). Focusing on the customers, teahouse customers arrive on foot or by public or private transport. Therefore, these options should be accounted for during the analysis. The pedestrian flow is included as a criterion in order to identify the possible flow of clients that arrive on foot. The availability of parking spaces is also included a criterion in order to identify the possibilities for arrival of customers relying on private transport. In this case, the capacity of the parking space within 500 m around the candidate location is considered. The access to public transport is also included in the model in order to capture the possible customer flow relying on the use of public transport. The access to public transport is quantified by considering the frequency of arrival of public transport in vicinity of the candidate location. Distance to scenery is yet another factor that can increase the attractiveness of the candidate place, even though the current customer flows are subdued. Therefore, this criterion is also included in the analysis. Similarly, the number of crimes in the area is also likely to shape the flows of customers. The distance from the public facilities is included into the model in order to account for possible customer inflow.
The economic performance of the teahouse is impacted by the competitive environment. Therefore, the number of competitors enters into the MCDM model. It is measured as the number of similar teahouses within the radius of 500 m [62]. The possibilities for outdoor advertising depend on the regulations prevailing in the area. Indeed, this might also imply changes in the competitiveness of the outlet. We include the area of possible outdoor advertisement as the corresponding indicator. The disposing of the waste generated during the process of operation is related to both economic and environmental dimensions (and, eventually, the social one) of the teahouse business. Therefore, the criteria for distance to the garbage containers can be considered as yet another measure for the sustainability of the operation of a teahouse.

EDAS Method
The evaluation based on distance from average solution (EDAS) method was proposed by Keshavarz Ghorabaee et al. [36]. The EDAS technique can be identified from the other methods in that it relies on the measurement of the distances from the average solution rather than the ideal ones. In this regard, the decision space is less impacted by the outlying observations. The measures of the positive distance from the average (PDA) and negative distance from the average (NDA) are applied in order to define the relative performance of the alternatives considered with respect to multiple criteria. Indeed, the alternatives showing higher (resp. lower) values of PDA (resp. NDA) are more preferred. Keshavarz Ghorabaee et al. [75] extended the EDAS technique into the fuzzy environment.
Let there be n alternatives and m criteria. The EDAS method then proceeds in the following way: Step 1. The decision matrix is constructed. For instance, the rows represent the alternatives and the columns represent the criteria as follows: where X ij denotes the value defining the rating of the alternative i with respect to the criterion j with i = 1, 2, . . . , n and j = 1, 2, . . . , m, where n and m are the numbers of alternatives and criteria, respectively.
Step 2. For each criterion, calculate the average solution: where each component of AV is obtained as follows: Step 3. The relative positions against the average solution are obtained in the form of PDA and NDA.
The following matrices store the results: The elements of the matrices defined by Equations (4) and (5) depend on the type of the criteria. If the j-th criterion is benefit one, then the following is true: If the j-th criterion is cost one, then the following is true.
where the positive and negative distances for alternative i with regards to criterion j are denoted by PDA ij and NDA ij , respectively.
Step 4. The positive and negative distances obtained for each alternative are aggregated by applying the weighted sum. Therefore, SP i and NP i denote the two weighted sums: where w j is the weight associated with criterion j.
Step 5. The aggregates obtained in Step 4 are normalized with respect to the maximum values in the following manner: Step 6. Each of the alternatives is attributed with the utility score, which is obtained as the average of the two normalized aggregates: where 0 ≤ AS i ≤ 1.
Step 7. The alternatives are ranked in descending order of the utility scores from Step 6.

WASPAS-N Method
The weighted aggregated sum product assessment (WASPAS) technique is based on two types of aggregation, namely the weighted sum and weighted product model [30,76]. Indeed, these two concepts are related to the arithmetic and geometric means, respectively. The combination of the two approaches is appealing in the context of the measurement of sustainability, as the weighted arithmetic mean allows for complete substitution among the arguments (criteria in the case of MCDM), whereas the substitution is limited in the case of the weighted geometric mean. The combination of the two approaches, therefore, offers an intermediate case. The WASPAS technique has been applied for energy system management [77,78], third-party logistics providers [79] and indoor construction [80], among others. Jahan [81] extended the WASPAS technique for target-based normalization.
The use of the arithmetic and geometric means involves certain methodological issues. Specifically, the geometric mean is not greater than the arithmetic mean. Therefore, we modify the WASPAS technique and propose normalizing the two means with respect to the maximum values. We term the proposed technique WASPAS-N.
Step 1. The MCDM problem is formulated as a decision matrix as follows: where the number of alternatives is represented by n and the number of criteria is given by m; x ij denotes the rating of alternative i against criterion j.
Step 2. The decision matrix X is normalized by applying the linear normalization. This procedure depends on the type of criteria.
For benefit criteria, (16) and for cost criteria, where x ij is the normalized value of x ij .
Step 3. The WASPAS relies on the arithmetic and geometric weighted means. The relative additive utility of the i-th alternative is obtained as the weighted arithmetic mean: where the importance of the j-th criterion is represented by w j . The relative multiplicative utility of the i-th alternative is obtained as the weighted geometric mean: Step 4. We then offer a modified calculation procedure to account for different levels of the arithmetic and geometric means by calculating the overall relative utility as follows: Step 5. The alternatives are ordered in descending order of Q i .

Sensitivity Analysis
It is known that a number of factors can affect the outcomes of MCDM, and one of the major factors is the weights of criteria. Due to the necessity to ascertain if the alterations in weight vector induce alterations in the resulting orderings of the alternatives, one can apply multiple weighting schemes. What is more, the stability of the analysis may be impacted by the number of criteria and alternatives. Because of the difficulties of modeling and reporting these relationships, the robustness of the MCDM-based ranking is tested by exploiting the Monte Carlo procedure. Here, the Monte Carlo technique is utilized due to the need to know whether modifications in criteria weights can impact the consistency of the results. In this paper, we assume that criteria weights are derived from variables randomly drawn from a uniform distribution: where b = 1, 2, . . . , B keeps track of the Monte Carlo replications, B denotes the total number of replications carried out, and i.i.d. means identically and independently distributed. Accordingly, the Monte Carlo procedure generates B vectors of values used to derive the weights based on the uniform distribution. The following formula shows how the elements within the draws are normalized: The resulting values w b j then serve as the weights of criteria in replication b. In this paper, we use B = 1000.

Results
This section presents results of the analysis. We include the decision matrix in this section as it required much additional research when combining different sources of information. We then proceed to the ranking of the alternative locations by applying the two MCDM techniques. The ranking based on Monte Carlo simulation is contrasted to the one based on the expert survey.

Construction of the Decision Matrix
The criteria and the associated directions of optimization (minimization or maximization) for location selection of a teahouse are displayed in Table 2. We chose Vilnius, the capital city of Lithuania, as the case study, as this is the most appealing option for investment given the development trends of the economy and population flows. Six alternative locations were identified for consideration when locating a teahouse: Location 1-space for rent on Aušros Vartų Street, Location 2-space for rent in Vilniaus Street, Location 3-space for rent in Dominikonų Street, Location 4-space for rent in Vingrių Street, Location 5-space for rent in Bazilijonų Street, and Location 6-space for rent in Didžioji Street.
The decision matrix was then established by assigning the criterion values for each alternative ( Table 3). The data on rent cost, property area, and parking capacity criteria were collected from the official website for the advertisement of premises to let in Vilnius (Aruodas.lt). The data on pedestrian flow were collected from JCDecaux pedestrian flow research based on expert evaluations. The data on public transportation were collected by counting the frequency of buses and trolleybuses stopping by the alternatives per day (Trafi.lt). The data on distances to sceneries and public facilities were collected by measuring the distances on the virtual Vilnius map (maps.vilnius.lt). The data on the distances to garbage containers were collected by measuring the distances on the virtual map of the garbage containers in Vilnius (gis.v-planas.com). The data on the number of crimes were collected by counting the number of crimes on the virtual map of crime activity in Vilnius (ird.lt). The data on outdoor advertisement possibilities were collected by the virtue of the interactive map of Vilnius (maps.vilnius.lt).   Note: the weights of criteria were provided by the expert survey.
After constructing the matrix, the weights of criteria were calculated and further analysis was carried out both by using the ratings provided by the experts (the questionnaire was based on the Likert scale) and by applying the Monte Carlo simulation procedure to the WASPAS method and EDAS method. In this research, EDAS and WASPAS methods were applied by using Monte Carlo simulation and expert survey weights. Using the previously presented procedures for implementation of the WASPAS and EDAS methods, the optimal solution for the case study was determined, and the results based on each of the methods were compared.

Ranking the Alternatives
The WASPAS and EDAS techniques were applied without assuming any specific weight vectors during the Monte Carlo procedure. The results are given in Table 4. As it was stated before, the uniform distribution was applied in order to generate the weights. The rows in Table 4 represent the ranks assigned for different alternatives during the Monte Carlo simulation, whereas the columns represent different locations considered in the MCDM problem. The most preferable location is attributed with rank 1. Therefore, higher values within the cells indicate higher probabilities for assigning a certain rank for a particular location. By considering the highest frequencies observed within each column in Table 4, one can establish the most probable ranking in case the weights are perturbed. Obviously, the rankings present in Table 4 induce no conflicting situations among the alternatives as each column has unique most probable value (rank) independently on the MCDM technique applied. The results indicate that stability of the most preferable alternative is rather high (74%) for both MCDM techniques. In both cases, Location 1 is identified as the most preferable alternative. The probability of ranking Location 1 as the second most preferable alternative is 0.23 for the WASPAS technique and 0.21 for the EDAS technique. In the case of the WASPAS technique, the ranking of Locations 3 and 4 is the least certain as the corresponding stabilities are below 50%). In the case of the EDAS method, the ranking of Locations 1, 4, and 6 turned out to be the least certain (the stability was below 50%). In particular, the ranking of Location 4 under the EDAS technique suggested a scattered pattern. These results indicate that the use of single MCDM technique and single weight vector might render the ranking unreliable.
Besides the Monte Carlo simulation, we also opted for the expert survey. A total of five experts holding executive positions in the real estate or hospitality sector with work experience of not less than 10 years was surveyed. They were asked to rate the importance of the criteria considered on a five-point Likert scale. The resulting average scores were normalized with respect to their sum in order to retrieve the resulting weight vector (Table 3 presents the resulting weights). The rankings based on the Monte Carlo procedure and expert survey are summarized in Table 5.
As one can note from Table 5, both the Monte Carlo simulation and the expert survey resulted in the same optimal solution: Location 2 was identified as the most preferable location for a teahouse. Indeed, Location 2 is located on Vilnius Street. The pedestrian flow per day is the highest there. What is more, the property area and public transportation frequency showed good performance of this particular alternative. We can see the difference between the Monte Carlo simulation for WASPAS and Monte Carlo simulation for EDAS with regards to ranks 4-6. The results based on the weights provided by the expert survey are invariant to the choice of the MCDM technique. In any case, Location 2 remained the most preferable one.

Conclusions
This study presented the case of MCDM application to teahouse location selection in Vilnius, Lithuania. The results of the research are important as the tea market is not as saturated yet as compared with the substitute markets (e.g., coffee market) in Lithuania. From the theoretical viewpoint, the present study revises the WASPAS technique and proposes WASPAS-N involving normalization of the additive and multiplicative utility functions.
The empirical results of the present study are of importance for prospective investors, entrepreneurs, and enterprisers. More specifically, the set of criteria that determine the success of the possible locations for a teahouse has been established. The proposed framework allows imputing the data available and comprehensively comparing the possible alternatives for the future business plans.
The proposed research framework involved the two MCDM methods (i.e., WASPAS and EDAS). The application of the two methods allowed us to increase the robustness of the analysis. The same objective was further tackled by implementing the sensitivity analysis based on the Monte Carlo simulation.
Several possible locations in Vilnius were considered in the empirical analysis. By applying the developed framework and collecting the relevant data from multiple information sources, the most promising alternatives have been identified. It turned out that the city center is where the most appealing alternative for a teahouse is located (Location 2, Vilnius Str.). Specifically, the latter location features the highest pedestrian flow and public transportation frequency when compared with the other alternatives. Also, the largest area is available in this instance.
The results indicate that the optimal decision does not depend on the technique applied (i.e., WASPAS or EDAS). What is more, the same results are obtained if weights provided by the experts are applied rather than assuming random weights during the Monte Carlo simulation. Therefore, the results are robust to both random changes in the weights and switching from random weighting to expert survey.
As regards further research, fuzzy logics could be involved in the analysis to allow for the use of imprecise data. Such an approach would allow the capturing of more complex information patterns when solving economic problems. On the other hand, different MCDM techniques can be applied for further analysis, especially those based on mathematical programming (weighted-sum approach, reference point method, or lexicographic approach). In this paper, we considered eleven criteria for choosing the optimal location for a teahouse. This allows the considering of multiple facets of the location selection problem. However, some criteria might provide little additional information (i.e., redundant criteria might appear). In the future research, the sensitivity analysis could aim at checking the impacts of the reduction of the number of criteria.