Next Article in Journal
Different Ground Subsidence Contributions Revealed by Integrated Discussion of Sentinel-1 Datasets, Well Discharge, Stratigraphical and Geomorphological Data: The Case of the Gioia Tauro Coastal Plain (Southern Italy)
Next Article in Special Issue
Experiencing a Food Production Site Using Wearable Devices: The Indirect Impact of Immersion and Presence in VR Tours
Previous Article in Journal
Nature Interpretation as an Environmental Educational Approach in Visitor Management; The Application Dilemma for Different Target Groups at Masai Mara National Reserve, Kenya
Previous Article in Special Issue
Exploring the Role of Norms and Habit in Explaining Pro-Environmental Behavior Intentions in Situations of Use Robots and AI Agents as Providers in Tourism Sector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bootstrapped DEA and Clustering Analysis of Eco-Efficiency in China’s Hotel Industry

1
Department of International Business, National Taiwan University, Taipei 10617, Taiwan
2
College of Business Administration, Fujian Business University, Fuzhou 350016, China
3
Department of International Business Administration, Chinese Culture University, Taipei 11114, Taiwan
4
New Huadu Business School, Minjiang University, Fuzhou 350108, China
5
School of Economics and Management, Guangdong University of Petrochemical Technology, Maoming 525000, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(5), 2925; https://doi.org/10.3390/su14052925
Submission received: 30 January 2022 / Revised: 25 February 2022 / Accepted: 1 March 2022 / Published: 2 March 2022

Abstract

:
As one of the world’s largest and fastest growing industries, tourism is facing the challenge of balancing growth and eco-environmental protection. Taking tourism CO2 emissions as undesirable outputs, this research employs the bootstrapping data envelopment analysis (DEA) approach to measure the eco-efficiency of China’s hotel industry. Using a dataset consisting of 31 provinces in the period 2016–2019, the bootstrapping-based test validates that the technology exhibits variable returns to scale. The partitioning around medoids (PAM) algorithm, based on the bootstrap samples of eco-efficiency, clusters China’s hotel industry into two groups: Cluster 1 with Shandong as the representative medoid consists of half of the superior coastal provinces and half of the competitive inland provinces, while Cluster 2 is less efficient with Jiangsu as the representative medoid. Therefore, it is suggested that the China government conduct a survey of only Shandong and Jiangsu to approximately capture the key characteristics of the domestic hotel industry’s eco-efficiency in order to formulate appropriate sustainable development policies. Lastly, biased upward eco-efficiencies may provide incorrect information and misguide managerial and/or policy implications.

1. Introduction

The tourism sector is one of the largest and fastest growing industries in the world. According to the Report on World Tourism Economy Trends [1], global tourism revenues were only USD 2 billion in 1950, but hit USD 5.8 trillion in 2019, which is equivalent to 6.7% of global gross domestic product (GDP). In addition to tourism revenues, global total tourist arrivals (including domestic and inbound tourist arrivals) were 25 million in 1950, but reached 12.31 billion in 2019 for an increase of 4.6% over 2018. As their demand highly correlates with the number of tourists seeking to stay overnight, hotels are one of the fastest-growing segments of the tourism sector. Furthermore, the hotel industry provides similar products and services that are visible and easily imitated by competitors; thus, it is highly competitive. Therefore, the measurement of hotel efficiency is a vital theme in tourism research.
Rapid economic development is often accompanied by considerable environmental costs such as environmental pollution. Since the 1970s, the international community has called on companies to actively practice social responsibility and promote sustainable development; for example, the United Nations recommends that companies disclose their environmental, social, and governance (ESG) practices by 2030 [2]. Sustainable development should meet the needs of the present without sacrificing the ability of future generations to meet their needs and also strive to achieve balanced development between human beings and nature [3]. The tourism sector contributed approximately 5% of all man-made CO2 emissions in 2005, with accommodation accounting for 21% of overall emissions, second only to transportation at 75% [4]. Thus, sustainability is clearly desirable in this important global industry.
The hotel industry can impact the environment through waste disposal such as food waste [5,6], energy consumption, and water use. For instance, the International Hotel Environment Initiative [7] estimates that each guest generates an average of 0.8–1.2 kg of waste per day, which doubles on check-out day [8,9]. In addition, more than one-third of this waste is food waste, of which about three-quarters is still edible [10]. Given the increasing concern about ecological environment (eco-environment) issues, tourism is facing the challenge of balancing growth and ecological protection, also known as the challenge of sustainable tourism [11,12].
Since its reform and opening up in 1979, China has experienced rapid development, becoming the world’s largest emerging economy and now its second largest economy. Its tourism industry has also experienced rapid growth [13]. The Golden Week holiday, which first started in 1999, has further accelerated the development of China’s tourism industry [14]. With its beautiful scenery, unique culture, and rich cultural relics, China has become the world’s third largest inbound tourism destination since 1990 [15] and it was the second largest tourism market in 2018 [16]. However, shortening the sector’s development time has come with high environmental costs. In order to comply with the global trend of ecological environmental protection and solve its deteriorating domestic environmental problems, China has actively promoted sustainable development in recent years. Striking a balance between tourism growth and eco-environmental protection is clearly a major challenge for China.
The concept of eco-efficiency (ecological efficiency) is to create more goods and services while using fewer resources and reducing their impact on the natural environment [17,18]. It focuses on the practices of efficient use of resources for economic and environmental progress [19] and provides quantitative guidance for sustainable tourism growth [20,21]. In recent years, eco-efficiency has attracted more and more attention from tourism researchers [20,22], but many studies have used a partial factor framework in their assessment, ignoring the substitution relationships between factor inputs and leading to biased estimates [23]. Data envelopment analysis (DEA), as a non-parametric and total-factor analysis approach, is essentially a linear programming model that evaluates the efficiencies of decision-making units (DMUs) by calculating the best multiplier for inputs and outputs. It has been widely used in tourism eco-efficiency assessments [21,24,25].
Initially proposed by Charnes et al. [26] and called the CCR model, DEA assumes that production technology exhibits constant returns to scale (CRS). Banker et al. [27] relaxed the CRS assumption to account for variable returns to scale (VRS) and called it the BCC model. If the production technology is CRS, then both the CCR model and the BBC model will generate consistent estimators, but the latter is less efficient than the former. However, when technology exhibits VRS, then the latter still remains consistent, whereas the former becomes inconsistent [28]. Traditional DEA approaches employ linear programming techniques for solutions and thus do not directly offer clear statistical inferences.
To correctly measure the eco-efficiency of the hotel industry, this research extends the bootstrapping procedure initially proposed by Simar and Wilson [28], to statistically test the property of returns to scale for China’s hotel industry by considering undesirable environmental emissions. After recognizing the appropriate returns to scale, we take a similar bootstrapping approach to obtain bias-corrected eco-efficiencies and then analyze China’s hotels by location. Lastly, the bootstrapped DEA samples can help us to cluster hotels that are more similar to each other in the same group compared to those in other groups.

2. Literature Review

The tourism industry consists of several different sub-sectors such as hotels, transportation, travel agencies, and so on. The market to attract tourists is highly competitive, and so the measurement of tourism efficiency is an important theme in the tourism literature. Over the past few decades, many studies have employed DEA to evaluate tourism efficiency [29,30,31,32]. Since the hotel industry is a vital sector in tourism, there is quite an extensive amount of literature on hotel efficiency [33,34,35,36,37]. An early study by Banker and Morey [29] measured the technical and scale efficiencies of 60 restaurants in a fast food chain, while Johns et al. [33] were the first to compare hotels using DEA. Since then, many scholars have applied DEA models to evaluate the efficiency of the hotel industry in various regions, such as the United States [38], central Europe [39], European countries [40], China [32,41], Portugal [30,42], Ecuador [43], Oman [44], Taiwan [45,46], Tunisia [31], and more.
Sustainable development must not only fulfill the needs of the present without sacrificing the ability of future generations, but must also aim to accomplish balanced development between human beings and nature [3]. In view of this, the concept of eco-efficiency offers a quantitative direction for enhancing intensive development while achieving sustainability. The simplest way to measure eco-efficiency is by taking the ratio between economic value and environmental load [47]; however, it is generally insufficient to measure eco-efficiency in a single factor framework [23]. Some studies have used composite environmental loads or adverse ecological impacts to calculate eco-efficiency indices from two or more environmental aspects [48,49,50], but Kuosmanen and Kortelainen [51] argued that the eco-environment itself cannot produce any output and must be combined with inputs such as labor and capital to produce economic output. Therefore, Reinhard et al. [52] proposed that the total-factor framework, DEA is suitable to overcome the above limitations. Since then, many scholars have chosen DEA models to assess eco-efficiency [53,54,55,56,57,58].
As tourism continues to grow globally and its adverse impact on the environment through energy consumption and pollutant emissions becomes more prominent [59], the debate on how to achieve sustainable tourism growth has become a major focus of the tourism economy literature [60]. Gossling et al. [20] first introduced the concept of eco-efficiency into the field of sustainable tourism. From then onwards, tourism eco-efficiency, which seeks to offer more tourism products and services while reducing resource consumption and pollutant emissions, has offered quantitative guidance for the sustainable growth of tourism [24]. Many research studies have assessed and analyzed tourism eco-efficiency [60,61,62] and eco-productivity [24,25], but to our knowledge no studies on tourism and/or the hotel industry have used a bootstrapped DEA approach for clustering studies of eco-efficiency. Hence, extending the research on eco-efficiency in tourism via an appropriate clustering study is of great significance.

3. Methodology

Conventional DEA approaches rely on linear programming techniques for solutions and therefore are unable to directly show any statistical inference. Thus, this study extends Simar and Wilson’s [28] bootstrapping procedure to perform a hypothesis test of returns to scale by considering undesirable by-products, obtains bias-corrected eco-efficiencies, and then tests for comparisons of two means. This section provides the data to be used herein and the input–output variables.

3.1. The Bootstrap-Based Test of Returns to Scale

Suppose that there are H decision making units (DMUs). Each DMU employs n inputs v = ( v 1 , , v n ) + n to produce m outputs u = ( u 1 , , u m ) + m . The CCR model proposed by Charnes et al. [26], uses the smallest free disposal convex set covering all the data to estimate the production possibility set as follows:
Ψ = { ( v ,   u ) | v V λ ,   u U λ ,   λ 0 }
where V = ( v 1 , , v H ) , U = ( u 1 , , u H ) , λ is an (N × 1) vector of intensity variables, and 0 is an (N × 1) vector of zeros.
Many industries are characterized by multiple outputs, and though they may produce desirable outputs, they also generate some undesirable or unwelcome by-products. For example, the production of paper is desirable, but both biochemical oxygen and sulfur oxides are undesirable in paper and pulp mills [63]; performing loans are desirable while non-performing loans are undesirable in the banking industry [64,65]. Undesirable outputs do not show the characteristics of strong disposability, as they are only weakly disposable, and hence cannot be handled by the CCR model.
Färe and Grosskoph [66] subsequently extended the CCR model by incorporating q undesirable outputs b = ( b 1 , , b q ) + q as follows:
Ψ c r s = { ( v ,   b ,   u ) |   v V λ ,     b = B λ ,     u U λ ,     λ 0 }
where B = ( b 1 , , b H ) . The superscript “crs” indicates that production technology exhibits constant returns to scale (CRS), which is an assumption implicit in the CCR model. The inequality constraints for desirable output vector u and input vector v reflect the property of strong disposability, while the equality constraint for undesirable output vector b reveals the characteristics of weak disposability. Given desirable and undesirable outputs, the input-oriented eco-efficiency of the CCR model for DMU h, β ^ h c r s , is defined as:
β ^ h c r s = m i n   β ,   λ { β |   ( β v h , b h , u h ) Ψ c r s }
The CRS assumption requires that all DMUs work at optimal scale, while some DMUs may actually operate below or above optimal scale. Hence, inefficiency can arise from inappropriate scaling as well as the unsuitable mixing of inputs and outputs. The BCC model, suggested by Banker et al. [27], permits production technology to have variable returns to scale (VRS). Mathematically, we just need to add the convexity constraint 1 λ = 1 in Equation (2), where 1 = ( 1 , , 1 ) is an (N × 1) vector of ones. The corresponding input-oriented eco-efficiency is:
β ^ h v r s = m i n   β ,   λ { β |   ( β v h , b h , u h ) Ψ v r s }
where
Ψ v r s = { ( v ,   b ,   u ) |   v V λ ,     b = B λ ,     u U λ ,     1 λ = 1 ,     λ 0 }
The superscript “vrs” indicates that the production technology is variable returns to scale (VRS).
If the production technology exhibits CRS globally, then both β ^ h c r s and β ^ h v r s are consistent estimators of eco-efficiency, but the former is more efficient than the latter. However, if the production technology displays VRS in some locations, then β ^ h c r s becomes inconsistent, while β ^ h v r s remains consistent [28]. Nevertheless, conventional CCR and BCC models rely on linear programming techniques to measure eco-efficiencies and are therefore unable to directly perform any statistical inference to examine the property of returns to scale. This study extends the bootstrapping procedure of Simar and Wilson [28,67] by considering undesirable environmental emissions to generate bootstrap samples, thereby testing the returns to scale. We summarize the bootstrapping procedure as follows.
(i)
Use the appropriate DEA model to calculate ρ ^ h (= 1 / β ^ h ) (h = 1, 2, …, H).
(ii)
Randomly draw a sample of size H with replacement from the following set:
{ ( 2 ρ ^ 1 ) , , ( 2 ρ ^ H ) , ρ ^ 1 , ,   ρ ^ H }
to obtain { ρ 1 * , ,   ρ H * } .
(iii)
Calculate ρ ˜ h * = ρ ¯ * + ( τ h * ρ ¯ * ) / 1 + ( γ / σ ^ * ) 2 , h = 1 , , H , where τ h * = ρ h * + γ ε h , γ = 1.06 m i n { σ ^ ,   R ^ / 1.34 } H 1 / 5 , σ ^ and R ^ are respectively the standard deviation and the interquartile range of { ρ ^ 1 , ,   ρ ^ H } , ε h draws randomly from a standard normal distribution, and ρ ¯ * and σ ^ * are the mean and standard deviation of { ρ 1 * , ,   ρ H * } , respectively.
(iv)
Calculate the bootstrap estimate β ^ h * (h = 1, 2,…, H ) by the appropriate DEA model with technology (V*, B, U), where V* = ( v 1 * , , v H * ) , v h * = ρ ˜ h * * v h / ρ ^ h , and ρ ˜ h * * = 2 ρ ˜ h * if ρ ˜ h *   <   1 and ρ ˜ h * * = ρ ˜ h * otherwise.
(v)
Re-do the above steps (ii)~(iv) L times to get bootstrap estimates { β ^ h l * } h = 1 H for l = 1 , , L .
Because the assumption of CRS is more restrictive than that of VRS, we specify the null hypothesis H0 and alternative hypothesis Ha as follows.
H0: the technology set is global CRS versus Ha: the technology set is VRS.
The appropriate test statistic, suggested by Simar and Wilson [67], is:
T = H 1 h = 1 H β ^ h c r s / β ^ h v r s
where β ^ h c r s / β ^ h v r s provides an estimate of the distance between the CRS and VRS estimated frontiers in the input direction at the point where ( v h , b h , u h ) is projected onto the estimated frontiers. By construction, β ^ h c r s β ^ h v r s , means rejecting the null hypothesis H0 if the test statistic T is too small.
To perform the hypothesis test empirically, we simulate L bootstrap samples of size H under H0 (the production technology is global CRS). In other words, we evaluate eco-efficiencies in step (i) by the CCR model, and for each bootstrap sample l we compute the bootstrap estimate T l * to obtain the bootstrap sample { T l * } l = 1 L . Since we reject the null hypothesis H0 if the test statistic T is small enough, the p-value is approximated by the proportion of bootstrap samples with values T l * less than t o (the observed value of T under H0 and the observed dataset):
p - value l = 1 L I ( T l * t o ) / L ,
where I(α) = 1 if α is true and 0 otherwise.

3.2. Correcting Eco-Efficiency Bias

The eco-efficiency obtained by DEA, by construction, is biased upward [68]. Although it is consistent, the bias will disappear at a slow rate [69]. Let θ ^ be an estimator of θ. The bias is defined as:
B I A S ( θ ^ ) = E ( θ ^ ) θ
where E ( θ ^ ) is the mathematical expectation of θ ^ . Note that if θ ^ is an unbiased estimator of θ, then B I A S ( θ ^ ) = 0 since E ( θ ^ ) = θ .
To correct the eco-efficiency bias empirically, we first apply the bootstrap procedure, introduced in Section 3.1, for L replications to obtain the bootstrap sample { β ^ h l * } l = 1 L for h = 1 , , H . The bias-corrected estimator of β ^ h c is then:
β ^ h c = β ^ h B I A S ^ ( β ^ h )
where B I A S ^ ( β ^ h ) = β ^ ¯ h * β ^ h , and β ^ ¯ h * is the average value of the bootstrap sample { β ^ h l * } l = 1 L . Nevertheless, Efron and Tibshirani [70] indicated that correcting for the bias also causes additional noise and suggested that the bias-correction should not be used unless:
|   B I A S ^ ( β ^ h )   | > s t d ( β ^ h * ) / 4
where s t d ( β ^ h * ) = ( L 1 l = 1 L ( β ^ h l * ) 2 ( β ^ ¯ h * ) 2 ) 0.5 .

3.3. The Bootstrap Test for Comparison of Two Means

Suppose that we want to investigate whether the mean eco-efficiencies of group 1 and group 2 are equal. Null hypothesis H0 and alternative hypothesis H1 are thus specified as:
H 0 :    μ 1 = μ 2    v e r s u s     H 1 :    μ 1 μ 2 ,
where μ 1 and μ 2 are the mean eco-efficiencies of group 1 and group 2, respectively. The corresponding test statistic is:
μ ^ = h group   1 β ^ h / H 1 h group   2 β ^ h / H 2
If we have the sampling distribution of μ ^ , then it is straightforward to construct the 100(1−2α)% confidence interval for μ ( = μ 1 μ 2 ). Nevertheless, we generally do not know the finite sampling distribution of μ ^ . The bootstrap sample { μ ^ l * } l = 1 L provides an empirical distribution of ( μ ^ * μ ^ ) and allows us to locate the values a α * and b 1 α * such that:
Pr ( a α * μ ^ * μ ^ b 1 α * ) = ( 1 2 α )
Because ( μ ^ μ ) approximates ( μ ^ * μ ^ ) , we use the bootstrap approximation:
Pr ( a α * μ ^ μ b 1 α * )   ( 1 2 α ) ,
so as to obtain the corresponding 100(1−2α)% confidence interval for μ ( = μ 1 μ 2 ) as:
[ ( μ ^ b 1 α * ) ,     ( μ ^ a α * ) ]
If the above interval does not include zero, then we reject H0 at the 200α% level of significance.

3.4. Data and Input–Output Variables

The dataset, obtained from the Yearbook of China Tourism Statistics and the China Statistical Yearbook on the Environment, consists of 31 provinces over the period 2016–2019. Since we have 4 years of data, there are 124 observations. In addition, all nominal variables are deflated by the GDP deflator with 2017 as the base year.
The basic concept of measuring eco-efficiency is to maximize the outputs by minimizing inputs and environmental damage. Based on previous studies, we consider three desirable outputs, revenues from guest rooms, revenues from food and beverage (F&B), and other revenues [45,71,72,73,74], and three input variables, employees, number of guest rooms, and fixed assets [71,72,74,75].
Operations in the hotel industry generate negative impacts on the environment (undesirable by-products), because it is an intensive energy user [76] and responsible for significant amounts of food waste [5,6], and so on. The carbon footprint method, which measures carbon dioxide (CO2) emissions, is an important approach for assessing the environmental impact of tourism [77]. Following the suggestion of Gossling [78] and Qiu et al. [79], this study takes the bottom-up approach to measure CO2 emissions by using the following formula:
E i = 365 × δ × B E D i × l i
where E i is the accommodation carbon footprint in province i, B E D i is total number of available beds in province i, l i is the annual average room occupancy rate in province i, and δ is CO2 emission per bed per night and set to be 2.458 g. Table 1 reports the summary statistics of the inputs and outputs used herein.
The output and input variables in the DEA model should satisfy the property of isotonicity—that is, increased inputs cannot reduce outputs. Table 2 presents the Pearson correlation coefficients of the input and output variables. All values are significantly positive at the 0.1% level, indicating that our selected input and output variables do indeed meet the property of isotonicity.

4. Results and Discussion

This section discusses the test for returns to scale, the eco-efficiencies of China’s hotel industry, coastal versus inland hotels, and the cluster analysis.

4.1. Bootstrap-Based Test of Returns to Scale

In any study of production efficiency, a vital question concerns whether the returns to scale of the underlying technology are increasing, constant, or decreasing. If the production technology displays CRS globally, then both β ^ c r s and β ^ v r s are consistent, while the latter is less efficient than the former. This implies that we should use the CCR model to estimate eco-efficiency. On the other hand, if the production possibility set exhibits VRS at some locations, then β ^ v r s still remains consistent, but β ^ c r s turns out to be inconsistent [28]. This situation suggests that the BCC model is more suitable for measuring eco-efficiency than the CCR model.
Conventional DEA models typically rely on linear programming techniques for solutions and thus cannot generally provide a formal statistical test of returns to scale. In addition, even though the DEA estimator under mild assumptions is consistent, it converges slowly and is known as the curse of dimensionality [80]. Löthgren and Tambour [81] implemented a native bootstrap method to perform statistical analyses, but unfortunately the native bootstrap method yields inconsistent estimates [82]. This study extends the bootstrap approach proposed by Simar and Wilson [28] to implement a statistical test of returns to scale by considering undesirable CO2 emissions.
This study uses R software for empirical analysis. Table 3 shows that based on 2000 replications for the bootstrapped procedure, the observed test statistic of 0.9722 is less than the critical value of 0.9773 for a 1% level test (p-value = 0.0005). Therefore, we reject the null hypothesis in which the production technology is globally CRS. Hence, we apply the BCC model to estimate the eco-efficiency and analyze relevant properties of China’s hotel industry.

4.2. General Analysis

Eco-efficiency, by construction, is biased and asymptotically consistent at a slow rate. Following the suggestion of Simar and Wilson [67], this study applied the bootstrap procedure for 2000 replications to obtain the bootstrap sample and to correct the bias. However, Efron and Tibshirani [70] suggested that the bias-correction should not be used unless the condition of Equation (10) is satisfied. Since bootstrapping results do fulfill the condition of Equation (10) for each DMU (the mean values of |   B I A S ^ ( β ^ h )   | and s t d ( β ^ h * ) / 4 are 0.0514 and 0.0081, respectively), we shall use Equation (9) to correct the bias of eco-efficiencies.
Table 4 shows the descriptive statistics of the BCC eco-efficiency β ^ and the bias-corrected eco-efficiency β ^ c . The average value of β ^ , obtained from Equation (4), is 0.9105, implying on average that China’s hotel industry could decrease the number of employees by about 9%. Nevertheless, since eco-efficiencies are biased upward, the mean value of the bias-corrected eco-efficiency β ^ c at 0.8591 indicates that the hotel industry on average could potentially reduce its employees by 14% instead of 9%. In addition, 35 DMUs (28.23%) are BCC-efficient ( β ^ = 1), or operating on the production frontier, while the maximum value of bias-corrected eco-efficiency β ^ c is 0.9609, suggesting that there is no efficient DMU in our sample after correcting for bias.

4.3. Analysis of Coastal and Inland Hotel Industries

Ever since the launch of economic reforms in 1978, China has advocated a progressive reform policy and shifted from a planned economy to a market economy. A series of preferential policies have shown obvious regional tendencies. The economic development of China was first concentrated in the Pearl River Delta region during the 1980s, then the Yangtze River Delta in the 1990s, and the BoHai Rim since the end of the 20th century. Hence, economic development in China is mainly concentrated in the coastal provinces (including Beijing, Tianjin, Liaoning, Hebei, Shandong, Jiangsu, Shanghai, Zhejiang, Fujian, Guangdong, Guangxi, and Hainan). As such, we might expect that the hotel industry in coastal provinces (hereafter, coastal hotels) outperform those in the inland provinces (hereafter, inland hotels). Figure 1 exhibits the provinces of the coastal and inland regions.
Table 5 shows the summary of the eco-efficiencies for the coastal and inland hotel industries. The average BCC eco-efficiency β ^ of coastal hotels is 0.9293, which is greater than that of inland hotels at 0.8986, while the mean bias-corrected eco-efficiency β ^ c of coastal hotels is 0.8676, which is only a little higher than that of inland hotels at 0.8537. These results imply that China’s hotel industry could potentially reduce the number of employees by approximately 14%, regardless of location. Moreover, 35 BCC-efficient DMUs operate on the production frontier: 20 are coastal DMUs and 15 are inland DMUs.
Although the average eco-efficiencies of the hotel industry reveal discrepancies among different locations, we need additional analysis to examine whether they are significantly different. Some research studies have used analysis of variance (ANOVA) to investigate whether the mean eco-efficiencies among groups are equal [65,83,84]. However, ANOVA assumes that all populations are normal and have equal variances. Some studies have applied the non-parametric test [85], which depends on the properties of large samples [86]. Daraio and Simar [69] pointed out that even though BCC eco-efficiency estimators are asymptotically consistent, they converge very slowly.
We employed the bootstrap procedure for 2000 replicates to acquire bootstrap samples and to construct confidence intervals for the mean eco-efficiency differences between coastal and inland hotels. Table 6 indicates that the 90% and 95% confidence intervals do contain zero values. Hence, we conclude that the mean eco-efficiencies between coastal and inland hotels do not reveal significant differences at the 10% level (p-value = 0.24). Note that the p-values of the non-parametric test (whether mean eco-efficiencies between coastal and inland hotels are equal) for BCC eco-efficiency is 0.0248—that is, coastal hotels statistically outperform inland hotels at the 5% level of significance. This suggests that biased upward eco-efficiencies may deliver incorrect information and misguide managerial and/or policy implications.

4.4. Cluster Analysis

We have so far investigated the eco-efficiency of the hotel industry by coastal and inland regions, but this is unlikely to fully characterize China’s hotel industry. Cluster analysis is a multivariate method that aims to classify similar objects (such as regions, provinces, companies, and so on) in the same group based on a set of variables. It involves a large number of techniques, methods, and algorithms that can be applied in various fields, including banking, marketing, tourism, etc. Several studies have used it to complement DEA studies [87,88,89].
The k-medoids or partitioning around medoids (PAM) algorithm is a clustering approach used to map a distance matrix into k clusters. It is like the k-means algorithm, but is based on the search for k representative objects, or medoids, among the objects in the dataset [90]. It is also more robust to noise and outliers as compared to k-means since it minimizes a sum of general pairwise dissimilarities instead of a sum of squared Euclidean distances. Furthermore, the medoids are robust representations of the cluster centers, whereas the cluster center in k-means is not necessarily one of the objects (it is the average between the objects in the cluster).
The bootstrap samples of eco-efficiency can help us to construct a distance matrix based on the probability that the order of two objects may be reversed [89,91]. Let β ¯ h be the mean eco-efficiency of the hth object over the sample periods. If the jth object completely outperforms the hth object, then Pr ( β ¯ h β ¯ j > 0 ) = 0. Contrarily, if the hth object completely outperforms the jth object, then Pr ( β ¯ h β ¯ j > 0 ) = 1. Both cases reveal the largest dissimilarity between the hth object and the jth object, that is, both objects never switch in rank over all of the bootstrap samples.
In contrast to those two extreme cases, Pr ( β ¯ h β ¯ j > 0 ) = 0.5 discloses the smallest dissimilarity between the hth object and the jth object. In other words, both objects have an equal likelihood of outperforming each other. Hence, we measure the distance between the hth object and the jth object by:
D p ( h , j ) = | Pr ( β ¯ h β ¯ j > 0 ) 0.5 | × 200
The values of D p ( h , j ) are between 0 and 100. Higher values of D p ( h , j ) imply a larger dissimilarity between the hth object and the jth object. Given the bootstrap samples of the eco-efficiency, we use the distance metric D p ( h , j ) to compute the (symmetric) dissimilarity matrix. The k-medoids clustering minimizes the average dissimilarity of objects to their closest representative object (medoid).
We applied the “cluster” package in R software to map a distance matrix into k clusters. Figure 2a,b are bivariate cluster plots for k = 2 and k = 3, respectively. All 31 objects are standardized and represented by points in the plot, using the first two principal components. An ellipse is drawn around each cluster. Figure 2b shows that there are considerable observations in the overlapping area of Cluster 1 and Cluster 2. Hence, we chose k = 2 for the cluster analysis of China’s hotel industry.
Table 7 indicates that Cluster 1 has higher mean values for both eco-efficiencies β ^ and β ^ c (0.9738 and 0.9018) than Cluster 2 (0.8583 and 0.8239). In addition, Cluster 1 contains all the efficient provinces of the BCC model, while the corresponding maximum eco-efficiency in Cluster 2 is 0.9307. In Cluster 1, the coastal and inland provinces are evenly divided. Hence, we can characterize the samples of Cluster 1 as half being superior coastal provinces and half of those competitive provinces being located in the inland region, while Cluster 2 is less efficient than Cluster 1 and mainly covers inland provinces.
One useful clustering application is to characterize the clusters by means of typical or representative objects. The k-medoids algorithm is based on the search for k representative medoids with a dataset whose average dissimilarity to all the objects in the cluster is minimal. A medoid is not only the representative object of a cluster and the most centrally located point in that cluster, but it is also actually an observation and can be used economically or conveniently to represent that cluster. In other words, these medoids provide a characterization of the clusters and also offer a chance for investigation by using only a small set of k medoids instead of a large set to start off with. Table 7 shows that the representative medoids for Cluster 1 and Cluster 2 are Shandong and Jiangsu, respectively. Hence, we are able to conveniently analyze these two provinces so as to properly capture the main features of the eco-efficiency of China’s hotel industry. Figure 3 displays the provinces in Cluster 1 and Cluster 2.
The empirical results show that the mean eco-efficiencies do not reveal significant differences between coastal and inland hotels, which may suggest that governments can simply plan or execute national policies to promote the sustainable development of the hotel industry, regardless of regional and/or economic development. However, cluster analysis indicates that the eco-efficiency of China’s hotel industry presents considerable differences, and is divided into at least two clusters. Therefore, the central government should formulate or implement eco-environmental policies according to the degree of eco-efficiency in order to effectively promote the sustainable development of its domestic hotel industry.

5. Conclusions

Global climate change, which is mainly characterized by climate warming, has increasingly become a serious challenge for mankind. Tourism is one of the world’s largest and fastest growing industries, and its eco-environmental issues have increasingly attracted widespread attention from governments and the international community. Hence, the tourism industry is facing the challenge of balancing growth and eco-environmental protection. Taking tourism CO2 emissions as undesired outputs, this study utilized the bootstrapping procedure to measure the eco-efficiency of China’s hotel industry. After identifying appropriate returns to scale by using the bootstrapping-based test, we applied a similar bootstrapping approach to correct for biases in eco-efficiency and to perform cluster analysis.
The dataset consists of 31 provinces over the period 2016–2019, thus covering 124 observations. The test for returns to scale does support the research hypothesis that production technology exhibits VRS. Furthermore, biased upward eco-efficiencies may provide incorrect information and misguide managerial and/or policy implications. The k-medoids clustering, based on the bootstrap samples of eco-efficiency, partitions China’s hotel industry into two clusters: Cluster 1 contains half of the better coastal provinces and half of those competitive provinces in the inland region, whereas Cluster 2 is less efficient than Cluster 1, covering about two-thirds of the inland provinces. In addition, the representative medoids of Cluster 1 and Cluster 2 are Shandong and Jiangsu, respectively.
One limitation of this study is the use of the carbon footprint method to measure CO2 emissions as an undesirable by-product of the hotel industry. The operational impacts of this industry on the ecological environment consist of diverse components such as energy consumption, water use, and waste disposal. They may have distinct impacts on the environment, but a single index assumes that their effects are homogeneous. Future studies could construct a multi-indicator framework to measure the undesirable by-products of the hotel industry so as to comprehensively evaluate its eco-efficiency. In addition, future research could also apply conditional non-parametric frontier methods to analyze how external environmental variables that are neither inputs, nor desirable outputs, nor undesirable by-products affect tourism eco-efficiency.

Author Contributions

Conceptualization, Y.L. and A.-C.L.; methodology, Y.L.; software, Y.L.; validation, Y.L., Y.Z. (Yueru Zhang) and Y.Z. (Yiting Zhan); formal analysis, Y.L., A.-C.L., Y.-Y.Y. and W.-C.L.; investigation, Y.Z. (Yueru Zhang) and Y.Z. (Yiting Zhan); resources, W.-C.L.; data curation, W.-C.L.; writing—original draft preparation, Y.L., A.-C.L., Y.-Y.Y. and W.-C.L.; writing—review and editing, Y.L., A.-C.L. and Y.-Y.Y.; visualization, Y.Z. (Yueru Zhang) and Y.Z. (Yiting Zhan); supervision, Y.L. and Y.-Y.Y.; project administration, Y.-Y.Y. and W.-C.L.; funding acquisition, W.-C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Report on World Tourism Economy Trends. 2020. Available online: https://cn.wtcf.org.cn/special/2020/0110/4988.html (accessed on 26 January 2022).
  2. Sustainable Stock Exchanges. Sustainable Stock Exchanges Initiative: Model Guidance on Reporting ESG Information to Investors. 2015. Available online: www.sseinitiative.org/wp-content/uploads/2015/09/SSE-Model-Guidance-on-Reporting-ESG.pdf (accessed on 29 February 2020).
  3. World Commission on Environment and Development. Our Common Future; Oxford University Press: Oxford, UK, 1987. [Google Scholar]
  4. World Tourism Organization; International Transport Forum. Transport-Related CO2 Emissions of the Tourism Sector—Modelling Results; UNWTO: Madrid, Spain, 2019. [Google Scholar] [CrossRef]
  5. Malefors, C.; Callewaert, P.; Hansson, P.A.; Hartikainen, H.; Pietilainen, O.; Strid, I.; Strotmann, C.; Eriksson, M. Towards a baseline for Food-Waste quantification in the hospitality Sector-Quantities and data processing criteria. Sustainability 2019, 11, 3541. [Google Scholar] [CrossRef] [Green Version]
  6. Amicarelli, V.; Aluculusei, A.C.; Lagioia, G.; Pamfilie, R.; Bux, C. How to manage and minimize food waste in the hotel industry? An exploratory research. Int. J. Cult. Tour. Hosp. Res. 2021. [Google Scholar] [CrossRef]
  7. International Hotel Environmental Initiative. Hotels Case: Community Action and Responsibility for the Environment; IHEI: London, UK, 2002. [Google Scholar]
  8. Pirani, S.I.; Arafat, H. Solid waste management in the hospitality industry: A review. J. Environ. Manag. 2014, 146, 320–336. [Google Scholar] [CrossRef] [PubMed]
  9. Abdulredha, M.; Al Khaddar, R.; Jordan, D.; Kot, P.; Abdulridha, A.; Hashim, K. Estimating solid waste generation by hospitality industry during major festivals: A quantification model based on multiple regression. Waste Manag. 2018, 77, 388–400. [Google Scholar] [CrossRef] [PubMed]
  10. Williams, P.; Leach, B.; Christensen, K.; Armstron, G.D.; Hawkins, R.P.; Lane, A.J.G.; Scholes, P. The Composition of Waste Disposed of by the UK Hospitality Industry; Report RES093-001; Waste & Resources Action Programme (WRAP): Banbury, UK, 2011. [Google Scholar]
  11. Hall, C.M. Framing behavioural approaches to understanding and governing sustainable tourism consumption: Beyond neoliberalism, “nudging” and “green growth”? J. Sustain. Tour. 2013, 21, 1091–1109. [Google Scholar] [CrossRef]
  12. Bella, G. Estimating the tourism induced environmental Kuznets curve in France. J. Sustain. Tour. 2018, 26, 2043–2052. [Google Scholar] [CrossRef]
  13. Yang, Y.; Liu, Z.-H.; Qi, Q. Domestic tourism demand of urban and rural residents in China: Does relative income matter? Tour. Manag. 2014, 40, 193–202. [Google Scholar] [CrossRef] [Green Version]
  14. York, Q.Y.; Zhang, H.Q. The determinants of the 1999 and 2007 chinese golden holiday system: A content analysis of official documentation. Tour. Manag. 2010, 31, 881–890. [Google Scholar] [CrossRef]
  15. Deng, T.; Li, X.; Ma, M. Evaluating impact of air pollution on China’s inbound tourism industry: A spatial econometric approach. Asia Pac. J. Tour. Res. 2017, 2, 771–780. [Google Scholar] [CrossRef]
  16. World Travel and Tourism Council. Economic Impact Reports. 2019. Available online: https://www.wttc.org/economic-impact/country-analysis/ (accessed on 14 October 2020).
  17. Verfaillie, H.A.; Bidwell, R. Measuring Eco-Efficiency. In A Guide to Reporting Company Performance; World Business Council for Sustainable Development: Geneva, Switzerland, 2000. [Google Scholar]
  18. Amicarelli, V.; Rana, L.R.; Lombardi, M.; Bux, C. Material flow analysis and sustainability of the Italian meat industry. J. Clean. Prod. 2021, 299, 126902. [Google Scholar] [CrossRef]
  19. United Nations. Eco-Efficiency Indicators: Measuring Resource-Use Efficiency and the Impact of Economic Activities on the Environment; ST/ESCAP/2561; United Nations Publication: Bangkok, Thailand, 2009. [Google Scholar]
  20. Gossling, S.; Peeters, P.; Ceron, J.P.; Dubois, G.; Patterson, T.; Richardson, R.B. The eco-efficiency of tourism. Ecol. Econ. 2005, 54, 417–434. [Google Scholar] [CrossRef]
  21. Zha, J.; Yuan, W.; Dai, J.; Tan, T.; He, L. Eco-efficiency, eco-productivity and tourism growth in China: A non-convex metafrontier DEA-based decomposition model. J. Sustain. Tour. 2020, 28, 663–685. [Google Scholar] [CrossRef]
  22. Hadjikakou, M.; Miller, G.; Chenoweth, J.; Druckman, A.; Zoumides, C. A comprehensive framework for comparing water use intensity across different tourist types. J. Sustain. Tour. 2015, 23, 1445–1467. [Google Scholar] [CrossRef] [Green Version]
  23. Huang, J.; Yang, X.; Cheng, G.; Wang, S. A comprehensive eco-efficiency model and dynamics of regional eco-efficiency in China. J. Clean. Prod. 2014, 67, 228–238. [Google Scholar] [CrossRef]
  24. Liu, J.; Zhang, J.; Fu, Z. Tourism eco-efficiency of Chinese coastal cities-analysis based on the DEA-Tobit model. Ocean. Coast Manag. 2017, 148, 164–170. [Google Scholar] [CrossRef]
  25. Peng, H.; Zhang, J.; Lu, L.; Tang, G.; Yan, B.; Xiao, X.; Han, Y. Eco-efficiency and its determinants at a tourism destination: A case study of Huangshan National Park, China. Tour. Manag. 2017, 60, 201–211. [Google Scholar] [CrossRef]
  26. Charnes, A.; Cooper, W.W.; Rhodes, E. Measuring the Efficiency of Decision-making Units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  27. Banker, R.D.; Charnres, A.; Cooper, W.W. Some Models for Estimation of Technical and Scale Inefficiencies in Data Envelopment Analysis. Manag. Sci. 1984, 30, 1078–1092. [Google Scholar] [CrossRef] [Green Version]
  28. Simar, L.; Wilson, P.W. Non-Parametric Tests of Returns to Scale. Eur. J. Oper. Res. 2002, 139, 115–132. [Google Scholar] [CrossRef]
  29. Banker, R.D.; Morey, R.C. Efficiency analysis for exogenously fixed inputs and outputs. Oper. Res. 1986, 34, 513–521. [Google Scholar] [CrossRef]
  30. Barros, C.P.; Mascarenhas, M.J. Technical and allocative efficiency in a chain of small hotels. Int. J. Hosp. Manag. 2005, 24, 415–436. [Google Scholar] [CrossRef]
  31. Hathroubi, S.; Peypoch, N.; Robinot, E. Technical efficiency and environmental management: The Tunisian case. J. Hosp. Tour. Manag. 2014, 21, 27–33. [Google Scholar] [CrossRef]
  32. Xia, B.; Dong, S.; Zhao, M.; Li, Z.; Li, F.; Li, Y.; Cheng, H. Analysis of economic efficiency and eco-efficiency of Chinese star hotels based on SBM model. Earth Environ. Sci. 2018, 190, 12066. [Google Scholar] [CrossRef] [Green Version]
  33. Johns, N.; Howcroft, B.; Drake, L. The use of data envelopment analysis to monitor hotel productivity. Prog. Tour. Hosp. Res. 1997, 3, 119–127. [Google Scholar] [CrossRef]
  34. Assaf, A.; Barros, C.P. Performance analysis of the Gulf hotel industry: A Malmquist index with bias correction. Int. J. Hosp. Manag. 2011, 30, 819–826. [Google Scholar] [CrossRef]
  35. Assaf, A.; Cvelbar, L.K.; Pahor, M. Performance drivers in the casino industry: Evidence from Slovenia. Int. J. Hosp. Manag. 2012, 32, 149–154. [Google Scholar] [CrossRef]
  36. Mendieta-Peñalver, L.F.; Perles-Ribes, J.F.; Ramón-Rodríguez, A.B.; Such-Devesa, M.J. Is hotel efficiency necessary for tourism destination competitiveness? An integrated approach. Tour. Econ. 2018, 24, 3–26. [Google Scholar] [CrossRef]
  37. Lado-Sestayo, R.; Fernández-Castro, Á.S. The impact of tourist destination on hotel efficiency: A data envelopment analysis approach. Eur. J. Oper. Res. 2019, 272, 674–686. [Google Scholar] [CrossRef]
  38. Anderson, R.I.; Fok, R.; Scott, J. Hotel Industry Efficiency: An Advanced Linear Programming Examination. Am. Bus. Rev. 2000, 18, 40–48. [Google Scholar]
  39. Ilić, I.; Petrevska, I. Using DEA method for determining tourism efficiency of Serbia and the surrounding countries. Hotel. Tour. Manag. 2018, 6, 73–80. [Google Scholar] [CrossRef]
  40. Soysal-Kurt, H. Measuring Tourism Efficiency of European Countries by Using Data Envelopment Analysis. Eur. Sci. J. 2017, 13, 1857–7881. [Google Scholar] [CrossRef] [Green Version]
  41. Sun, Y.; Hou, G.; Huang, Z.; Zhong, Y. Spatial-Temporal Differences and Influencing Factors of Tourism Eco-Efficiency in China’s Three Major Urban Agglomerations Based on the Super-EBM Model. Sustainability 2020, 12, 4156. [Google Scholar] [CrossRef]
  42. Barros, C.P.; Alves, F.P. Productivity in tourism industry. Int. Adv. Econ. Res. 2004, 10, 215–225. [Google Scholar] [CrossRef]
  43. Higuerey, A.; Viñan-Merecí, C.; Malo-Montoya, Z.; Martínez-Fernández, V. Data Envelopment Analysis (DEA) for Measuring the Efficiency of the Hotel Industry in Ecuador. Sustainability 2020, 12, 1590. [Google Scholar] [CrossRef] [Green Version]
  44. Oukil, A.; Channouf, N.; Al-Zaidi, A. Performance evaluation of the hotel industry in an emerging tourism destination: The case of Oman. J. Hosp. Tour. Manag. 2016, 29, 60–68. [Google Scholar] [CrossRef]
  45. Lee, Y.L.; Kuo, S.H.; Jiang, M.Y.; Li, Y. Evaluating the performances of Taiwan’s international tourist hotels: Applying the directional distance function and meta-frontier approach. Sustainability 2019, 11, 5773. [Google Scholar] [CrossRef] [Green Version]
  46. Yin, P.; Chu, J.; Wu, J.; Ding, J.; Yang, M.; Wang, Y. A DEA-based two-stage network approach for hotel performance analysis: An internal cooperation perspective. Omega 2020, 93, 102035. [Google Scholar] [CrossRef]
  47. Zhang, B.; Bi, J.; Fan, Z.; Yuan, Z.; Ge, J. Eco-efficiency analysis of industrial system in China: A data envelopment analysis approach. Ecol. Econ. 2008, 68, 306–316. [Google Scholar] [CrossRef]
  48. Mickwitz, P.; Melanen, M.; Rosenstrom, U.; Seppala, J. Regional eco-efficiency indicators—A participatory approach. J. Clean. Prod. 2006, 14, 1603–1611. [Google Scholar] [CrossRef]
  49. Ji, D. Evaluation on China’s regional eco-efficiency: Based on ecological footprint methodology. Contemp. Econ. Manag. 2013, 35, 57–62. [Google Scholar]
  50. Wang, J. A comprehensive evaluation of China’s provincial ecological pressure and ecological efficiency: Based on ecological footprint method. J. Contemp. Financ. Econ. 2016, 8, 3–15. [Google Scholar]
  51. Kuosmanen, T.; Kortelainen, M. Measuring eco-efficiency of production with data envelopment analysis. J. Ind. Ecol 2005, 9, 59–72. [Google Scholar] [CrossRef]
  52. Reinhard, S.; Lovell, C.A.K.; Thijssen, G.J. Environmental efficiency with multiple environmentally detrimental variables; estimated with SFA and DEA. Eur. J. Oper. Res. 2000, 121, 287–303. [Google Scholar] [CrossRef]
  53. Mahlberg, B.; Luptacik, M. Eco-efficiency and eco-productivity change over time in a multisectoral economic system. Eur. J. Oper. Res. 2014, 234, 885–897. [Google Scholar] [CrossRef] [Green Version]
  54. Lampe, H.W.; Hilgers, D. Trajectories of efficiency measurement: A bibliometric analysis of DEA and SFA. Eur. J. Oper. Res. 2015, 240, 1–21. [Google Scholar] [CrossRef]
  55. Woo, C.; Chung, Y.; Chun, D.; Seo, H.; Hong, S. The static and dynamic environmental efficiency of renewable energy: A Malmquist index analysis of OECD countries. Renew. Sust. Energy Rev. 2015, 47, 367–376. [Google Scholar] [CrossRef]
  56. Broadstock, D.C.; Li, J.; Zhang, D. Efficiency snakes and energy ladders: A (meta-) frontier demand analysis of electricity consumption efficiency in Chinese households. Energy Policy 2016, 91, 383–396. [Google Scholar] [CrossRef] [Green Version]
  57. Yu, Y.; Huang, J.; Zhang, N. Industrial eco-efficiency, regional disparity, and spatial convergence of China’s regions. J. Clean. Prod. 2018, 204, 872–887. [Google Scholar] [CrossRef]
  58. Chen, L.; Huang, Y.; Li, M.J.; Wang, Y.M. Meta-frontier analysis using cross-efficiency method for performance evaluation. Eur. J. Oper. Res. 2020, 280, 219–229. [Google Scholar] [CrossRef]
  59. Sofield, T.; Li, S. Tourism governance and sustainable national development in China: A macro-level synthesis. J. Sustain. Tour. 2011, 19, 501–534. [Google Scholar] [CrossRef]
  60. Han, J. Carrying capacity of low carbon tourism environment in coastal areas from the perspective of ecological efficiency. J. Coast. Res. 2018, 83, 199–203. [Google Scholar] [CrossRef]
  61. Yang, G.H.; Li, P.; Zheng, B.A.; Zhang, Y.Q. GHG emission-based eco-efficiency study on tourism itinerary products in Shangri-La, Yunnan Province, China. Curr. Issues Tour. 2008, 11, 604–622. [Google Scholar] [CrossRef]
  62. Rungsrisawat, S. Image of tourism destinations-principles of macromarketing and environmental efficiency: The case of Thailand. Int. J. Ecol. Dev. 2017, 32, 133–140. [Google Scholar]
  63. Färe, R.; Grosskopf, S.; Lovell, C.K.; Yaisawarng, S. Derivation of Shadow Prices for Undesirable Outputs: A Distance Function Approach. Rev. Econ. Statis. 1993, 75, 374–380. [Google Scholar] [CrossRef]
  64. Chang, C.-C. The Nonparametric Risk-Adjusted Efficiency Measurement: An Application to Taiwan’s Major Rural Financial Intermediaries. Am. J. Agric. Econ. 1999, 81, 902–913. [Google Scholar] [CrossRef]
  65. Li, Y.; Hu, J.L.; Liu, H.W. Non-performing Loans and Bank Efficiencies: An Application of the Input Distance Function Approach. J. Stat. Manag. Syst. 2009, 12, 435–450. [Google Scholar] [CrossRef]
  66. Färe, R.; Grosskopf, S. New Directions: Efficiency and Productivity; Springer: New York, NY, USA, 2004. [Google Scholar]
  67. Simar, L.; Wilson, P.W. Sensitivity Analysis of Efficiency Scores: How to Bootstrap in Nonparametric Frontier Models. Manag. Sci. 1998, 44, 49–61. [Google Scholar] [CrossRef] [Green Version]
  68. Banker, R.D. Maximum Likelihood, Consistency and Data Envelopment Analysis: A Statistical Foundation. Manag. Sci. 1993, 39, 1265–1273. [Google Scholar] [CrossRef]
  69. Daraio, C.; Simar, L. Advanced Robust and Nonparametric Methods in Efficiency Analysis: Methodology and Applications; Springer: New York, NY, USA, 2007. [Google Scholar]
  70. Efron, B.; Tibshirani, R.J. An Introduction to the Bootstrap; Chapman and Hall: New York, NY, USA, 1993. [Google Scholar]
  71. Wang, F.C.; Hung, W.T.; Shang, J.K. Measuring Pure Managerial Efficiency of International Tourist Hotels in Taiwan. Serv. Ind. J. 2006, 26, 59–71. [Google Scholar] [CrossRef]
  72. Shyu, J.; Hung, S.C. The True Managerial Efficiency of International Tourist Hotels in Taiwan: Three–stage Data Envelopment Analysis. Serv. Ind. J. 2012, 32, 1991–2004. [Google Scholar] [CrossRef]
  73. Chen, C.M.; Chang, K.L.; Langelett, G. How Demand Uncertainty and Market Concentration Affect Long–term Price Instability. Int. J. Hosp. Manag. 2014, 37, 146–149. [Google Scholar] [CrossRef]
  74. Guo, X.; Wei, W.; Li, Y.; Wang, L.Y. A Study of Different Types of Air Pollutants on the Efficiency of China’s Hotel Industry. Int. J. Environ. Res 2019, 16, 4319. [Google Scholar] [CrossRef] [Green Version]
  75. Shieh, H.S.; Hu, J.L.; Gao, L.Y. Tourist Preferences and Cost Efficiency of International Tourist Hotels in Taiwan. Int. J. Mark. Stud. 2014, 6, 1918–7203. [Google Scholar] [CrossRef] [Green Version]
  76. World Tourism Organization; United Nations Environment Programme. Climate Change and Tourism: Responding to Global challenges; World Tourism Organization: Madrid, Spain, 2008. [Google Scholar]
  77. Dong, X.; Zhang, J.; Zhang, J.; Cheng, S. A Critical review on several issues of regional tourism-related carbon emissions or its carbon footprint. Acta Ecol. Sin. 2016, 36, 554–568. [Google Scholar]
  78. Gossling, S. Global environmental consequences of tourism. Glob. Environ. Chang. 2002, 12, 283–302. [Google Scholar] [CrossRef]
  79. Qiu, X.; Fang, Y.; Yang, X.; Zhu, F. Tourism Eco-Efficiency Measurement, Characteristics, and Its Influence Factors in China. Sustainability 2017, 9, 1634. [Google Scholar] [CrossRef] [Green Version]
  80. Kneip, A.; Simar, L.; Wilson, P.W. Asymptotics for DEA Estimators in Non-Parametric Frontier Models; Discussion Paper No. 0317; Institut de Statistique, Universite’ Catholique de Louvain: Louvain-la-Neuve, Belgium, 2003. [Google Scholar]
  81. Löthgren, M.; Tambour, M. Bootstrapping the DEA-based Malmquist Productivity Index. Appl. Econ. 1999, 31, 417–425. [Google Scholar] [CrossRef]
  82. Simar, L.; Wilson, P.W. Statistical inference in nonparametric frontier models: The state of the art. J. Product. Anal. 2000, 13, 49–78. [Google Scholar] [CrossRef]
  83. Aly, H.Y.; Grabowski, R.; Pasurka, C.; Rangan, N. Technical, Scale, and Allocative Efficiencies in U.S. Banking: An Empirical Investigation. Rev. Econ. Stat. 1990, 72, 211–218. [Google Scholar] [CrossRef]
  84. Li, Y. DEA Efficiency Measurement with Undesirable Outputs: An Application to Taiwan’s Commercial Banks. Int. J. Serv. Technol. Manag. 2005, 6, 544–555. [Google Scholar] [CrossRef]
  85. Li, Y.; Chen, Y.K.; Chien, F.S.; Lee, W.C.; Hsu, Y.C. Study of Optimal Capital Adequacy Ratios. J. Product. Anal. 2016, 45, 261–274. [Google Scholar] [CrossRef]
  86. Wackerly, D.D.; Mendenhall, W., III; Scheaffer, R.L. Mathematical Statistics with Applications, 7th ed.; Duxbury Press: Pacific Grove, CA, USA, 2008. [Google Scholar]
  87. Amin, G.R.; Emrouznejad, A.; Rezaei, S. Some Clarifications on the DEA Clustering Approach. Eur. J. Oper. Res. 2011, 215, 498–501. [Google Scholar] [CrossRef]
  88. Dai, X.; Kuosmanen, T. Best-Practice Benchmarking Using Clustering Methods: Application to Energy Regulation. Omega 2014, 42, 179–188. [Google Scholar] [CrossRef]
  89. Li, Y. Analyzing Efficiencies of City Commercial Banks in China: An Application of the Bootstrapped DEA Approach. Pac-Basin Financ. J. 2020, 62, 101372. [Google Scholar] [CrossRef]
  90. Kaufman, L.; Rousseeuw, P.J. Finding Groups in Data: An Introduction to Cluster Analysis; Wiley: New York, NY, USA, 1990. [Google Scholar]
  91. Hirschberg, J.G.; Lye, J.N. Clustering in a Data Envelopment Analysis Using Bootstrapped Efficiency Scores. In Department of Economics—Working Papers Series 800; The University of Melbourne: Melbourne, Australia, 2001. [Google Scholar]
Figure 1. Provinces of the coastal and inland regions.
Figure 1. Provinces of the coastal and inland regions.
Sustainability 14 02925 g001
Figure 2. Bivariate cluster plots of partitioning around medoids. (a) Bivariate Cluster Plot of PAM for k = 2; (b) Bivariate Cluster Plot of PAM for k = 3.
Figure 2. Bivariate cluster plots of partitioning around medoids. (a) Bivariate Cluster Plot of PAM for k = 2; (b) Bivariate Cluster Plot of PAM for k = 3.
Sustainability 14 02925 g002
Figure 3. Provinces in Cluster 1 and Cluster 2.
Figure 3. Provinces in Cluster 1 and Cluster 2.
Sustainability 14 02925 g003
Table 1. Descriptive statistics of inputs and outputs.
Table 1. Descriptive statistics of inputs and outputs.
VariablesMeanMedianStd. Dev.MinMax
Inputs
 Number of employees (thousand)58.39148.89447.1325.077245.588
 Number of guest rooms (thousand)43.98542.19728.4632.565164.458
 Fixed assets (RMB billion)15.64511.25513.4781.50161.481
Undesirable Output
 CO2 emissions (thousand tons)3.7233.4712.5050.17613.877
Desirable Outputs
 Revenues from guest rooms (RMB billion)6.9584.8226.7940.40031.840
 Revenues from F&B (RMB billion)4.6943.2154.3270.19718.681
 Other revenues (RMB billion)1.7050.9462.2770.0699.882
Note: All nominal variables are deflated by the GDP deflator with 2017 as the base year.
Table 2. Pearson correlation coefficients between input and output variables.
Table 2. Pearson correlation coefficients between input and output variables.
OutputsRevenues from
Guest Rooms
Revenues from
F&B
Other
Revenues
CO2
Emissions
Inputs
Number of
employees
0.92650.95370.82710.8864
(<0.001)(<0.001)(<0.001)(<0.001)
Fixed assets0.84350.85920.79060.9911
(<0.001)(<0.001)(<0.001)(<0.001)
Number of
guest rooms
0.94110.88970.91490.8967
(<0.001)(<0.001)(<0.001)(<0.001)
Notes: The values in parentheses are p-values. All correlation coefficients are significant at the 0.1% level.
Table 3. Bootstrap-based test of returns to scale.
Table 3. Bootstrap-based test of returns to scale.
t o p-ValueCritical Values
α = 0.01α = 0.05α = 0.10
H0: CRS
Ha: VRS
0.97220.00050.97730.98020.9814
Note: We use 2000 replications for the bootstrapped procedure.
Table 4. Descriptive statistics of eco-efficiencies.
Table 4. Descriptive statistics of eco-efficiencies.
MeanMedianStd. Dev.MinMax
BCC Eco-Efficiency β ^ 0.91050.91340.08130.70851
Bias-corrected Eco-Efficiency β ^ c 0.85910.87630.06350.67760.9609
Notes: (i) We use 2000 replications for the bootstrapped procedure. (ii) Bootstrapping results do fulfill the condition of Equation (10) for each DMU, and thus we correct the bias of eco-efficiencies.
Table 5. Summary of eco-efficiencies for the coastal and inland hotel industries.
Table 5. Summary of eco-efficiencies for the coastal and inland hotel industries.
MeanMedianStd. Dev.MinMaxNo. of Provinces
(DMUs)
BCC Eco-Efficiency β ^
Overall0.91050.91340.08130.7085131 (124)
Coastal Hotels0.92930.95420.07980.7583112 (48)
Inland Hotels0.89860.89680.08040.7085119 (76)
Bias-corrected Eco-Efficiency β ^ c
Overall0.85910.87630.06350.67760.960931 (124)
Coastal Hotels0.86760.88310.05960.73220.958812 (48)
Inland Hotels0.85370.86520.06570.67760.960919 (76)
Table 6. The bootstrap test for coastal versus inland hotels.
Table 6. The bootstrap test for coastal versus inland hotels.
Mean Eco-EfficiencyThe 100(1 − 2α)% Confidence Intervalsp-Value
β ^ β ^ c α = 0.025α = 0.05
Coastal Hotels0.92930.8676[−0.0085, 0.0343][−0.0050, 0.0324]0.24
Inland Hotels0.89860.8537
Notes: (1) We use 2000 replications for the bootstrapped procedure. (2) The p-values of the non-parametric test for the BCC eco-efficiency is 0.0248.
Table 7. Results of partitioning around medoids for K = 2.
Table 7. Results of partitioning around medoids for K = 2.
Cluster 1Cluster 2
β ^ β ^ c β ^ β ^ c
Mean0.97380.90180.85830.8239
Median0.97070.90160.88030.8377
Std. Dev.0.02270.01920.05130.0482
Min0.93740.86980.77750.7407
Max10.92790.93070.9137
Representative
Province
ShandongJiangsu
0.96640.89440.88990.8514
Number of Provinces1417
  Coastal Provinces75
  Inland Provinces712
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Liu, A.-C.; Yu, Y.-Y.; Zhang, Y.; Zhan, Y.; Lin, W.-C. Bootstrapped DEA and Clustering Analysis of Eco-Efficiency in China’s Hotel Industry. Sustainability 2022, 14, 2925. https://doi.org/10.3390/su14052925

AMA Style

Li Y, Liu A-C, Yu Y-Y, Zhang Y, Zhan Y, Lin W-C. Bootstrapped DEA and Clustering Analysis of Eco-Efficiency in China’s Hotel Industry. Sustainability. 2022; 14(5):2925. https://doi.org/10.3390/su14052925

Chicago/Turabian Style

Li, Yang, An-Chi Liu, Yi-Ying Yu, Yueru Zhang, Yiting Zhan, and Wen-Cheng Lin. 2022. "Bootstrapped DEA and Clustering Analysis of Eco-Efficiency in China’s Hotel Industry" Sustainability 14, no. 5: 2925. https://doi.org/10.3390/su14052925

APA Style

Li, Y., Liu, A. -C., Yu, Y. -Y., Zhang, Y., Zhan, Y., & Lin, W. -C. (2022). Bootstrapped DEA and Clustering Analysis of Eco-Efficiency in China’s Hotel Industry. Sustainability, 14(5), 2925. https://doi.org/10.3390/su14052925

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop