An international comparison of productivity change in the textile and clothing industry: a bootstrapped Malmquist index approach

Firms in the textile and clothing industry operate in competitive international markets characterized by the liberalized trade after the removal of multi-fiber agreement quotas in 2005, and have to address rapid changes in consumer preferences and production technology. Hence, improving competitiveness is crucial for firm survival. Competitiveness of the sector often depends on its firms meeting their production potential. This paper analyzes productivity changes in the textile and clothing industry worldwide during the period 1995–2004. A bootstrapped Malmquist approach is used to identify the respective contributions of technical change, technical efficiency change, and scale efficiency change. Moreover, differences in productivity changes across different groups of firms are statistically assessed. Our results show a relatively small overall productivity increase for both textile and clothing firms due to positive technical change, despite declines in technical and scale efficiency. Furthermore, our results indicate that productivity and its components differ for textile firms and clothing firms, for firms in countries that benefited and did not benefit from the quotas’ elimination, and for firms in different regions.


Introduction
The textile and clothing industry is facing great challenges due to its highly competitive nature and the ongoing globalization and internationalization of the textile and clothing markets. For over 30 years, trade relations in the textile and clothing sector were dominated by the multi-fiber agreement (MFA). The agreement played a crucial role in protecting producers from the developed world against competition from developing countries by allowing countries to impose quotas on imports of different categories of textile and clothing products. During the Uruguay Round, it was agreed to remove the quotas progressively over a 10-year period until 2005. Additionally, firms in the textile and clothing industry face rapid changes in production technology and consumer preferences toward more variety, higher quality, and more fashionable products (Owen 2001;Stengg 2001). Over the last two decades, the sector restructured and modernized its production, resulting in a large number of firm exits in the industry and a reduction of the industry's workforce by more than 40 % in Europe and approximately 30 % in the USA. In addition, the Asian textile and clothing industry was restructuring and making considerable adjustments to reduce employment, return to profitability, decrease the number of firms, and increase efficiency (Yeung and Mok 2004). Although the trends of decreasing employment and number of firms are expected to continue in the international textile and clothing industry, the industry remains a relatively important sector of the economies of many developed and developing countries. In 2004, the US textile and clothing industry accounted for 8 % of the nation's total value added in manufacturing; in Japan, China, and Italy, the industry accounted for more than 9, 11 and over 11 % of total manufacturing value added, respectively (World Bank 2005).
Pressures related to technological change and internationalization have affected the textile and clothing industry in recent decades, which has considerably altered the environment in which firms operate. The globalization and liberalization of markets increases competition, which could force firms to operate more efficiently in order to survive. In addition, the expanded set of customer tastes could imply that firms need to modify their technologies to produce new goods. While these changes might result in technical progress, they do not necessarily increase technical efficiency or productivity, especially when firms improve the quality of output to meet market demand for variety and fashion. In addition, the new environment of the textile and clothing industry can privilege some of the sub-sectors while disadvantaging others; this disparity could result in differences in efficiency or technical change between the textile and clothing industries. The manufacturing of clothing is largely low-tech and labor-intensive, while the textile industry is more technology-intensive and requires less unskilled labor (Mittelhauser 1997;Datta and Christoffersen 2005). However, the clothing sector also has a segment of high-quality fashion characterized by the use of modern technology. Furthermore, the textile industry is less flexible in terms of adjusting to consumer tastes than clothing (Nordås 2004). Finally, at the firm level, the contributions of efficiency, scale, and technical change depend on the firm's regional location due to the region-specific characteristics of this industry.
Although productivity growth is essential for survival of textile and clothing firms that operate in international markets and in a dynamic environment, it is unclear what causes productivity growth. Existing studies analyzing the evolution of productivity of textile and clothing companies have used both parametric (e.g., Stochastic production frontier) and nonparametric (e.g., Data Envelopment Analysis) approaches. Previous studies have focused on the analysis of the textile and clothing industry in individual countries. Datta and Christoffersen (2005) focused on the USA and found productivity growth of 2.1 % that was mainly driven by technical change. Wadud (2007) analyzed Australian textile and clothing firms over the period 1995-1998 and found both years with productivity decline and years with productivity increases. Goaïed and Ayed-Mouelhi (2000) and Ayed-Mouelhi and Goaïed (2003) focused on Tunisia; they found both technical regress and technical progress, depending on the method applied. Several studies focused on countries in Asia, such as Leung (1998), Färe et al. (2001), andChen (2003), who found high productivity growth rates of 7, 6, and 9 % for Singapore, Taiwan, and China, respectively. In all cases, productivity growth was mainly driven by technical change. Lower productivity growth rates were found by Mahadevan (2002) for Malaysia (0.3 %) and by Joshi and Singh (2010) for India (1.7 %). Margono and Sharma (2006) and Kong et al. (1999) reported productivity declines in Indonesia and China, respectively. Goncharuk (2007) found a large productivity decline of 7.29 % for Ukraine. Taymaz and Saatçi (1997) found rapid technical progress (6 %) for Turkish textile and clothing firms.
Previous studies analyzing productivity change in the textile and clothing industry have generally focused on a specific country and do not assess the differences across countries/regions. An international comparison of productivity change within this sector has not been made thus far. A systematic study that compares the evolution of productivity over time for textile and clothing firms across several countries/regions worldwide is fundamentally important for assessing the competitiveness of this sector in the globalized market.
Given this background, the objective of this paper is to analyze productivity changes of firms in the textile and clothing industry. The empirical application focuses on panel data over the period 1995-2004 of larger stock exchange listed firms in the textile and clothing industry in Asia, Europe, North America, Middle and South America, and the Rest of the World. The time period considered in this study covers the period of increasing trade liberalization in the textile and clothing industry when MFA quotas were gradually removed. We use the Malmquist index and follow Ray and Desli (1997) to decompose productivity growth to identify the respective contributions of technical efficiency change, scale efficiency change, and technical change. Unlike previous studies of the textile and clothing sector, our research undertakes a worldwide analysis of the textile and clothing sector and accounts specifically for scale efficiency change. Furthermore, our study applies the bootstrap method for computing the Malmquist index to obtain statistical inference. The bootstrap approach used in this study extends that of Simar and Wilson (1999) by implementing the Ray and Desli (1997) decomposition of productivity growth.
The remainder of this paper is organized as follows. Section 2 elaborates on the bootstrapped Malmquist approach and its decomposition. This is followed by the presentation of the data (Sect. 3) and discussion of results (Sect. 4). The final section offers concluding comments.

Methodology
In this section, we present the input-oriented Malmquist index and its decomposition. The input orientation in the Malmquist index is appropriate for the textile and clothing industry as this sector is embedded in a competitive market, where outputs are defined by demand and limited only by the utilization of the available resources. The input-oriented Malmquist index is based on the input distance function introduced by Shephard (1970). The production technology of a firm producing output y t and using input x t is represented by the production set and the input sets: The input distance function for period t is defined as follows: where λ is the value of the input distance function. Using t +1 instead of t in (3), one can obtain D t+1 (x t+1 , y t+1 ). In addition to these two distance functions, the computation of Malmquist requires the formulation of two mixed period distance functions. One concerns the firm at time t + 1 in relation to the technology at t The second refers to the firm at time t in relation to the technology at t + 1 D t+1 (x t , y t ), which can be obtained by simply swapping the superscript on D with those of x and y in (4). Assuming constant returns to scale (CRS, indicated with the subscript C), Färe et al. (1992) define an input-oriented Malmquist index between time period t and t + 1 as 1 : 1 It is worth pointing out that the geometric mean version of the adjacent Malmquist productivity index, which we apply in this paper, does not satisfy the circularity test. The seriousness of the failure to satisfy this test has not been resolved in the literature (Fried et al. 2008). A weaker version of time-neutrality of technical change is necessary and sufficient for this property to hold. For details, see Pastor and Lovell (2007).
Values of the Malmquist input-oriented index that are less than 1 indicate improvements in productivity in period t +1 compared to period t, whereas values greater than 1 indicate decreases in productivity in period t + 1 compared to period t. Although one well-known pitfall of the Malmquist index is that it only measures a local technical change and does not have the total factor productivity (TFP) interpretation (see O'Donnell (2010) for more details), it is widely applied in empirical research for measuring productivity growth in a variety of industrial sectors. The literature proposes different decompositions of the Malmquist index. The most commonly used decompositions are the formulations proposed by Färe et al. (1992) and Färe et al. (1994). Färe et al. (1992) decomposes the Malmquist index to identify two sources of productivity change: technical efficiency change and technical change. Technical efficiency change reflects the movement of a firm toward the CRS frontier (technical efficiency increase) or away from the CRS frontier (technical efficiency decrease). Technical change measures the shift of the CRS frontier itself from one period to the next. Färe et al. (1994) further decompose the efficiency change index into two indexes: one that measures the change in efficiency with respect to the variable returns to scale (VRS) technology and one that assesses the scale efficiency change, which is the comparison of the distance function gauged relative to the CRS with that relative to the VRS technology. However, the CRS assumption in the technical change component in this formulation is problematic: the estimation of the technology would be statistically inconsistent if the underlying true technology does not exhibit CRS. Hence, it is useful to estimate the changes in technology by the changes in the VRS frontier, as proposed by Ray and Desli (1997). Whereas the Ray and Desli formulation is widely supported in the literature 2 , its empirical application is limited (examples of its usage include Kumar and Basu 2008;Guan and Chen 2010;and Mahlberg and Url 2010). Ray and Desli (1997) break down the Malmquist index into three sources of productivity change: Note that the subscript V refers to the VRS technology. Scale efficiency change is the product of scale efficiency as in Färe et al. (1994) and the ratio of technical change and technical change itself under CRS and VRS technologies, respectively. The technical change component is the geometric mean of two ratios of VRS distance functions reflecting period t and period t + 1 technologies. The first ratio measures technical change as the ratio computed using period t +1 quantities; the second ratio is computed using period t quantities. The technical change component represents the productivity change that can be attributed to a movement in the VRS frontier between two time periods. If any of these components is smaller than 1, then a firm has improved its technical efficiency, scale efficiency, or technology; components equal to 1 suggest no change; and components larger than 1 suggest a decrease in technical efficiency, scale efficiency, or technology.
It is worth pointing out that one could also consider applying the formulations and decompositions of Malmquist index that control for capital fixity. Such an approach could be appropriate for the technology-intensive textile industry that is characterized by a large share of fixed inputs that cannot be adjusted rapidly. However, our database concerns firms in regions across the world, and the textile industry is more capital intensive in e.g., North America and Europe, but less capital intensive in Asia. Therefore, the capital fixity assumption could be more suitable for textile firms in North America and Europe and less for Asia. In order to be consistent across the entire dataset of textile firms and due to the limitation of data as it is not possible to do the separate estimations for every region, the issue of correcting for input fixity is left for future work.
We calculate the component distance functions of the Malmquist index using the linear programming technique known as Data Envelopment Analysis (DEA) (Charnes et al. 1978;Banker et al. 1984). The Malmquist computation requires solving eight linear programs for each decision making unit i (DMU i ): four single-period (under VRS and CRS) and four mixed period distance functions (under VRS and CRS). Assume k DMUs, each of which uses N inputs (x t ni ) to produce M outputs (y t mi ). The distance function D t C (x t , y t ) is calculated by solving the following linear program: In this program, z i reflects the activity levels associated with inputs and outputs of DMU i used to construct the linear segments of the CRS frontier. Note that D t+1 C (x t+1 , y t+1 ) is obtained using the analogous program. The mixed period distance function D t C (x t+1 , y t+1 ), which projects quantities in period t + 1 on the production frontier from period t, is given by: Again, the mixed period distance function D t+1 C (x t , y t ), which projects quantities in period t on the production frontier in period t + 1, is obtained analogously. Adding the restriction k i=1 z i = 1 to the linear program in (7) or (8) gives the distance function under the VRS technology.
To perform statistical inference for the Malmquist index and its components, bootstrap techniques can be used. Bootstrapping is defined as a repeated simulation of the data-generating process achieved by resampling from the data and applying the original estimator to each simulated sample; resulting estimates imitate the original unknown sampling distribution of the estimators of interest. We follow Simar and Wilson's (Simar andWilson 1999, 2000) algorithm adapted for the Ray and Desli (1997) decomposition. Although the bootstrap has been applied to the Malmquist index before (e.g., Gitto and Mancuso 2012), it has so far been rarely applied to the decomposition of the Malmquist index proposed by Ray and Desli (1997). The algorithm is summarized in the following steps: 1) From the original data set, compute the distance functions composing the solving the linear programming model (7) and its variations for t +1, mixed period and VRS distance functions.

2) Use a smoothed bootstrap to generate a pseudo-dataset
. . , k by drawing with replacement from the joint density of original distance functions. This density is estimated using a bivariate kernel density and univariate reflection method adapted to the bivariate case; this approach accounts for the possible temporal correlation arising from the time series dimension of the panel data. When the joint density is already estimated, drawing a pseudo-sample involves solving the additional linear programs. For observations for which one of the linear programs gives infeasible solutions, the bootstrap values are deleted and the procedure of drawing with replacement is repeated 3 . 3) Construct the bootstrap estimates of distance functions that compose the Malmquist index and its decomposition indices for each DMU i = 1, . . . , k by applying the original estimators to the pseudo-dataset derived from step 2. Use these estimates to compute bootstrap estimates of the Malmquist index and its components for

4) Repeat steps 2 and 3 B times to obtain a set of estimates:
There are two complementary ways of performing the statistical inference on computed indices: estimating bias-corrected indices or estimating confidence intervals. In this paper, we compute the confidence intervals because the bias-corrected estimator may be less reliable when it has a higher standard deviation than the original estimator (Simar and Wilson 1999). The confidence intervals for the Malmquist index are calculated by sorting the valuesM . . , B and deleting α 2 × 100 -percent of the elements at either end of this sorted array. Then, we set −b * α and −â * α â * α ≤b * α equal to the endpoints of the resulting array. As a result, the estimated (1 − α) -percent confidence interval is formulated as follows: The same procedure applies to the estimation of the confidence intervals for technical efficiency change, technical change, and scale efficiency change. It should be noted though that confidence intervals are calculated only for observations for which the Malmquist index is well defined, i.e., observations for which all distance functions have feasible outcomes. Therefore, when an infeasibility of mixed period distance functions occurs, the confidence intervals for such observations are not calculated. The interpretation of the confidence intervals is as follows: the Malmquist index for the ith DMU is significantly different from 1 at the 1, 5, or 10 % level if the 99, 95, and 90 % interval in (9) does not include 1. In this paper, the Malmquist index and its decomposition are estimated using FEAR 1.1 package 4 using B = 2000 bootstrap replications.

Data
This study uses the database of publicly traded firms in the international textile and clothing industry. This database combines information collected from three sources: COMPUSTAT, DATASTREAM, and OSIRIS. Although three different sources are used, the consistency of the database is guaranteed, because we use publicly available accounting data that are highly harmonized across countries and the information is audited for most firms. Linking these databases is undertaken using the international codes that identify listed firms such as ISIN (International Securities Identifying Number) and SEDOL (Stock Exchange Daily Official List). In addition, firms are classified as either 'textile firms,' or 'clothing firms.' COMPUSTAT and OSIRIS already divide the companies into these categories: COMPUSTAT uses a classification similar to the SIC code 5 , while OSIRIS is based on the NACE Rev. 2 code 6 . In DATASTREAM, all firms fall into one general category of "Clothing and accessories," following the ICB classification 7 . The classification by either textile or clothing industry is obtained using the extended business description provided by DATASTREAM or the firms' webpages. We extracted the firm's balance sheet and profit and loss account from the databases.
Initially, 2,588 firms are extracted from three databases, which are matched such that overlapping observations are removed. This results in a sample of 1,693 firms. After deleting observations with missing information on variables that are necessary to compute Malmquist index, the sample further reduces to 1,064 firms (5,337 observations). Finally, after removing the outliers 8 , the final sample consists of 1,016 textile and clothing firms that operated at least two consecutive years during the period from 1995 to 2004. The panel is unbalanced and it sums up to 4,985 observations. The firms in the sample are generally larger firms that are listed on the stock exchange and whose main activity is textile and clothing manufacturing. Our sample is almost equally distributed between textile and clothing firms: 2,584 observations are from the clothing industry and 2,401 are from the textile industry. The data come from 39 countries in Asia, Europe, North America, Middle and South America, and the Rest of the World (Australia, Oceania, and Africa). Table 1 provides details of the geographic location and the activity profile (textile or clothing) of the companies included in the sample, along with the information on the number of outliers removed. As it can be appreciated from the table the outliers account for 6.6 % of the initial sample. 5 SIC (Standard Industrial Classification) is a number that identifies a company's primary business. In COMPUSTAT number 2290 represents "Apparel and Accessories," while number 2500 refers to "Textiles." 6 NACE Rev. 2 is the European Union classification of economic activities. According to this code, number 13 represents "Manufacture of textiles," while number 14 refers to "Manufacture of wearing apparel." 7 ICB (Industry Classification Benchmark) is a classification standard for trading and investment decisions. The category of "Clothing and accessories" (3763) falls under Consumer Goods (3000)-Personal and Household goods (3700)-Personal goods (3760). 8 It is common knowledge that the DEA is very sensitive to the presence of outliers in the sample (Timmer 1971) and one might be concerned that firms display extreme behavior in the input-output proportions that they use. For the purpose of the paper we want to reduce the number of infeasibilities of the Malmquist index, therefore, we choose the outlier detection method that results in the smallest number of infeasibilities, but may have less proven qualities as an outlier detection method compared to other methods. In this way, we eliminate the firms which may or may not be outliers. Initially, the outliers are determined using the method of Wilson (1993). However, it proves insufficient for us, since after removing outliers we could not get feasible solutions for the Malmquist index for approximately 50 % of sample observations. Therefore, we switch to another method which is frequently applied in the literature in which the outliers are determined using ratios of output to input: an observation is defined as an outlier if the ratio of output over any of the three inputs is outside the interval of the median plus and minus three times the interquartile range. There are also alternative methods of outliers' detection; these include, for example, the methods proposed by Simar (2003) or Ondrich and Ruggiero (2002).  (1) 10 (1) 26 (2) Estonia 2 (0) 6 (1) 8 (1) Finland 30 (0) 0 (0) 30 (0) France 103 (4) 122 (6) 225 (10) Germany 87 (1) 154 (4) 241 (5) Greece 59 (19) 15 (3) 74 (22) Hong Kong 41 (3) 111 (12) 152 (15) Hungary 0 (0) 6 (0) 6 (0) India 17 (2) 12 (0) 29 (2) Indonesia 84 (4) 44 (2) 128 (6) Italy 55 (5) 93 (2) 148 (7) Japan 421 (21)   Following Piesse and Thirtle (2000), Wu et al. (2007) and Kapelko and Oude Lansink (2014) 9 , we consider one output and three inputs: costs of goods sold, tangible fixed assets, and number of employees. Output is measured as total sales and other operating revenues. Costs of goods sold involve all costs allocated to production. Tangible fixed assets are measured as the total amount of property, the plant and equipment in monetary terms, and net of depreciation. The number of employees is the annual average of full-time workers. All variables in monetary units (output, costs of goods sold, and tangible fixed assets) are initially expressed in local currencies in current prices 10 . To make these variables comparable across countries, they are divided by the purchasing power parity (PPP) of the local currency to the US dollar. To obtain constant 1995 prices, producer price indices (PPI) are used 11 . Both PPP and PPI are obtained from the World Bank, the International Monetary Fund and the Organization for Economic Co-operation and Development. Table 2 provides the descriptive statistics of the variables for the 1995-2004 time period separately for the textile and clothing companies.
The Malmquist index is calculated for each firm between two consecutive years of the 1995-2004 time period. Moreover, we perform the computations of the Malmquist index separately for the groups of textile and clothing firms because the production technology differs between these industries. However, we are not able to do the separate estimations for every region analyzed because we do not have enough firms in our database for some regions and some periods of time; the data are particularly limiting for following regions: the Rest of the World, and Middle and South America. Table 3 summarizes the geometric means of the Malmquist index and its decomposition separately for the two samples of textile and clothing firms. It should be noted that some of the mixed period linear programs used to compute the Malmquist index may not have a feasible solution. One possible solution to this problem is to omit the observations with infeasible solutions in the computation of averages, which is the strategy we have followed. In line with the recommendation of Briec and Kerstens (2009), we report the infeasibilities that occurred in the empirical application in Table  3. Out of 4985 observations, only 21 observations have an infeasible solution for one of the mixed period linear programs; this represents 0.4 % of the entire sample. Furthermore, Table 3 reports the results of non-parametric test of stochastic dominance (Kolmogorov-Smirnov (K-S) test statistics) for the differences between textile and clothing industries, which was developed specifically for efficiency analysis by Delgado et al. (2002) 12 . Table 3 suggests that in the period 1995-2004, the textile and clothing sectors have experienced small average annual productivity growth of 0.3 and 1 %, respectively. The productivity increase is largely driven by average annual technological progress of 1.5 and 2 % for the textile and clothing sectors, respectively. The positive contribution of technological progress to productivity growth is partially offset by a deterioration of technical efficiency and scale efficiency. Technical efficiency change and scale efficiency change are slightly more negative in the textile sector (−0.7 and −0.5 %, respectively) than in the clothing sector (−0.6 and −0.3 %, respectively).

Results
Productivity growth occurred in almost all years; the only exceptions were 2002/2003 and 2003/2004 for the textile sector and 1997/1998 for the clothing sector. Technological progress is also found in almost all years for the textile and clothing sectors. The finding of technological regress in some years suggests that technological change has eliminated some productive options that were previously available for the textile and clothing firms. For example, the increasing liberalization of trade due to the removal of MFA quotas might have forced the companies to undertake additional investments to adapt to the international standards and absorb additional costs with-  (11) 2574 (10) ***, **, * Significant differences between T and C at the 1, 5 and 10 % level, respectively out a direct productivity increasing impact. The decrease of scale efficiency of textile and clothing firms in almost all years indicates that scale of operation has, on average, become less beneficial. Overall, the growth of the gap between efficient and inefficient firms and the technological advances suggest that much of the decrease in technical and scale efficiency can be attributed to the failure of textile and clothing firms to adapt to the technological improvements made by some of their competitors; that is, a few textile and clothing firms are innovators and shift the frontier, while most of the other firms fail to adapt to the technological improvements and fall behind. This result could also explain the large number of exits from the textile and clothing industry in the past decades.
To appropriately compare the differences in indices between textile and clothing firms, we need to look at the results of test of stochastic dominance in Table 3. These findings indicate that for all periods except for 1995/1996, the null hypothesis of equality of distributions of textile versus clothing firms with regard to the Malmquist index is rejected at all conventional significance levels. Furthermore, there are no significant differences in efficiency change between sectors for three periods of time: 1996/1997, 2000/2001, and 2002/2003. Additionally, there are no significant differences in scale efficiency change across textile and clothing firms: for four periods of time (1995/1996, 1996/1997, 1998/1999, and 2000/2001), the null hypothesis of equal distributions is not rejected. This result reflects the fact that in the textile and clothing industry, scale efficiency change does not depend on the activity profile. On the other hand, technical change differs across subsectors, as evidenced by the fact that the null hypothesis of equality of distributions is rejected for all periods. As indicated in Table 3, relative to the textile industry, clothing makes more technological improvements in almost all periods of time. These improvements might explain the higher productivity growth observed in this sector for the 1995-2004 period as a whole. All these findings reflect the textile and clothing firms' ability to expand their production possibilities through innovation, while the overall productivity growth demonstrates that the improvement in textile and clothing companies' technology offsets the increases in technical and scale inefficiencies. Figure 1 illustrates the differences in the distributions of the Malmquist index and efficiency change, technical change, and scale efficiency change between textile and clothing firms. Because of space limitations, the figure displays only some time periods.
The graphs show that the distributions of the Malmquist index and technical change of clothing firms are always below the analogous distributions of textile companies, demonstrating the lower levels of these indices for clothing firms. The graphs suggest that clothing firms made more improvements in productivity and technology than textile companies.
To interpret the changes in the Malmquist index and its components in terms of their statistical significance, we need to analyze the confidence intervals computed by the bootstrap procedure. Table 6 in the Appendix reports the number of firms with significant and insignificant (at the 1 % critical level) growth or decline of the Malmquist index and its components separately for the textile and clothing sectors. As it was outlined before, for these firms for which some of the Malmquist components are infeasible to compute (the infeasibility of mixed period distance functions), the bootstrap's confidence intervals are not estimated as well; therefore, the total numbers reported in Table 6 are equal to the sample size after removing infeasibilities as reported in Table 3. The findings in this table reveal that for the majority of both textile and clothing firms, the changes in productivity, efficiency, technology, and scale are significantly different from 1 at the 1 % significance level. In addition, for a relatively large number of firms, stagnation in productivity and efficiency is found; this result is especially prevalent in the clothing sector. The results in the table show that the number of textile and clothing firms experiencing a significant productivity growth is almost equal to the number of firms experiencing productivity decline. Hence, it is clear that a large number of individual firms have experienced decline, even though Table 3 suggests that productivity has grown on average. The results for technical change in this table also display a larger number of firms experiencing technological progress than regress. Textile and clothing firms experiencing technical efficiency and scale efficiency decline outnumber those experiencing improvement.
4.1 Productivity change and country position regarding MFA import quotas' removal Next, we want to identify the impact of quotas' removal on productivity growth and its components. We compare two groups of firms: the firms in countries which benefited from the removal of quotas and these that did not benefit from the removal, based on the value of country exports before and after the removal of quotas. In particular, the calcu- K-S test 0.114*** 0.067 0.066 0.202*** ***, **, * Significant differences between firms that benefited and that did not benefit from the quotas' removal at the 1, 5 and 10 % level, respectively lation compares the mean value of sample country world exports for textile/clothing for the pre-2005 period to that observed in 2005, as published by the World Trade Organization 13 . Furthermore, these analyses are undertaken by separating the dataset into three periods that indicate the progressive removal of quotas as established by the World Trade Organization. The first period ("First stage") ranges from 1995/1996 till 1996/1997, reflecting the integration of products representing not less than 16 % of country 1990 total imports of textile and clothing. The second period ("Second stage") ranges from 1997/1998 till 2001/2002 when not less than further 17 % was integrated.  Table 4 presents the results for Malmquist productivity growth and its components in the first, second, and third phase of removal of quotas, together with K-S test statistics for the differences between firms in countries that benefited and that did not benefit from the elimination of quotas. A major finding from this analysis is that firms in countries that benefited from the removal of quotas and these that did not benefit respond differently to opening up to trade. Firms that benefited from the quotas' removal improve their performance over time with regard to productivity change, efficiency change, and technical change. It suggests that the increase in export opportunities of firms results in productivity, efficiency, and technology gains. On the contrary, the results reveal the worsening of performance across firms in countries that did not benefit from the quotas' removal after barriers to trade were drastically reduced. In particular, for these firms productivity growth is becoming smaller along the stages of relaxing of quota restrictions. Initial efficiency growth, when still considerable quota restrictions are presented in the market, turns into the efficiency decline in two last stages of opening up to trade.
Although still these firms experience positive technical change, this increase in technical change is becoming smaller when barriers to trade are completely removed. For both types of firms, scale efficiency change remains similar along the stages of quotas' removal. Therefore, the change in trade policy results in the improvement of performance for firms in countries that increase the exports of textile and clothing products and causes the worsening of performance for firms in countries that see the decrease in their value of exports. However, the encountered differences between the productivity growth and its components between firms in countries that benefited and did not benefit from the quotas' removal are not always statistically significant: results of Kolmogorov-Smirnov two-sided test show that this is the case for all stages of quotas' removal with regard to productivity change, for two initial stages for efficiency and technical change and for the last stage for scale efficiency change component. Table 5 compares the differences in geometric means of the Malmquist index and its components across firms with headquarters 14 in Asia, Europe, North America, Middle and South America, and the Rest of the World. To test the statistical significance of the differences across regions, a Kruskal-Wallis (K-W) nonparametric rank test is employed 15 . Table 5 also reports the growth of R&D intensity in each region between 1995 and 2004. R&D intensity reflects the R&D expenditures of the country as a percentage of GDP 16 . Growth in R&D intensity in each region is calculated based on the level of R&D intensity of countries in the sample as obtained from the World Bank. Table 5 shows important differences in productivity and its components among regions. Although Table 3 shows that the sample of all textile firms, on average, experienced productivity growth in the period 1995-2004, the results in Table 5 show that the growth only occurred in Asia and North America. Textile firms in Europe, Middle and South America and particularly those in the Rest of the World, experienced a productivity decrease. In all regions, efficiency change made a negative contribution and technical change made a positive contribution to productivity growth. Scale efficiency change made a negative contribution in all regions except the Rest of the World. The results for clothing firms in different regions are in line with the findings for the entire sample of firms, i.e., clothing firms on average experienced a productivity increase in all regions despite the negative contributions of efficiency change and scale efficiency change in most regions. Furthermore, technical progress occurred in all regions for both textile and clothing firms, but is particularly large for textile firms in Asia, and 14 Rather than analyzing productivity growth based on the locations of firms' headquarters, it would be interesting to perform this analysis based on the location where most production takes place because it is well known that textile and clothing companies in developed countries are moving the production to less developed countries to benefit from cost advantages and to increase efficiency. Unfortunately, data to do such an analysis are not available for us. 15 The test of stochastic dominance cannot be applied for this purpose as it is designed to compare only two groups of results. The application of the Kruskal-Wallis test to analyze frontier shifts is discussed in Sueyoshi and Aoki 2001. 16 Unfortunately, data on R&D expenses at firm's level are not available in the datasets that are used in this study.  Table 5 show that the Malmquist index, technical change and efficiency change significantly differ between regions for textile firms, whereas scale efficiency change does not significantly differ. These results for textile firms show that technical change, technical efficiency change and improvement of productivity depend on the geographical location of the firm. Firms in Asia and North America display a better performance than firms in Europe, Middle and South America, and the Rest of the World. For clothing firms, there are significant differences between regions only for scale efficiency change 17 .

Productivity change in regions
To graphically reflect the differences in indices across regions, Fig. 2 shows the distributions of the Malmquist index, efficiency change, technical change and scale efficiency change for some periods of time. Due to space limitations, we only present the graphs of Asia and Europe for textile firms. Figure 2 shows that the Malmquist index and technical change of textile firms differ between these two regions; distributions for 17 Although statistically significant differences are found in productivity growth and its decomposition between regions, the reader should keep in mind that within textile and clothing sectors firms are still heterogeneous. A statistically significant difference between regions might reflect differences in the composition of the subgroups. Asia lie below those for Europe, while distributions for efficiency and scale efficiency change do not markedly differ. Table 7 in the Appendix reports the bootstrap findings grouped by regional location. This table shows that the majority of textile and clothing firms in all analyzed regions have experienced significant (at the 1 % critical level) changes in the Malmquist index and its components. It is also worthwhile to note that a large number of textile and clothing firms from Asia, Europe and North America stagnated in terms of efficiency; additionally, a large number of North American firms experienced stagnation in terms of overall productivity.

Conclusions and future research
This paper analyzed productivity change and its decomposition into technical change, technical efficiency change and scale efficiency change of firms in the textile and clothing industry in the context of a gradual removal of MFA quotas on textile and clothing trade. In the face of the new liberalized market, improving productivity is an important condition for competitiveness of the textile and clothing industry. The empirical application used panel data for larger stock exchange listed textile and clothing firms from Asia, Europe, North America, Middle and South America, and the Rest of the World. In the 1995-2004 time period, productivity in the textile and clothing industries increased by an average of 0.3 and 1 % per year, respectively. The productivity increase was mainly driven by technical progress, which was partially offset by a decrease in technical efficiency and scale efficiency in both industries. A deeper analysis of productivity growth using the bootstrap method showed that productivity has significantly increased for more than half of the textile and clothing firms. This paper also finds that firms in countries that benefited from the quotas' removal improved their performance, while firms in countries that did not benefit from the quotas' removal worsened their performance in the period of quotas' elimination. Productivity growth and its components also differ depending on a firm's geographic region. Textile firms in North America and Asia experienced productivity growth on average in the period under investigation, whereas productivity declined for firms in Europe, Middle and South America, and the Rest of the World. Clothing firms in all regions experienced productivity growth on average. Another interesting finding is that the firms in regions that achieve the fastest technical change (Asia, Middle and South America, and the Rest of the World), are the ones that have the highest growth in R&D expenses.
The patterns of productivity, efficiency and technical change observed in this study reflect the rapid adjustments in the production environment in the past decade in the textile and clothing industry. In response to increased competition and new opportunities due to the removal of barriers to trade, the innovative firms have adapted quickly, while others have responded slowly and fallen behind. Of course, we are also cautious in interpreting our results: other forces may have contributed to these changes. To corroborate the findings, future research will need to test whether productivity change is correlated with, for example, the economic conditions of the region/country or the size of firms. Furthermore, in future work, we would need to treat the firms from the Australia and Oceania region separately from those located in Africa. This advance would require more data on companies from these regions. Furthermore, if we are able to collect more data in some regions, future research would allow us to compute the Malmquist index and its components for each region separately; this approach would allow for region-specific production technologies. Future research could also analyze productivity growth and its components of post-quota textile and clothing period to compare it with the current results during the trade reform. With regard to the methodology used, because one of the shortcomings of Malmquist index is that it only measures the local technical change, a promising avenue for future study would be to apply an index that has a TFP interpretation, such as the Hicks-Moorsteen index (O'Donnell 2010(O'Donnell , 2012a. Another methodological extension for future research, if the data allows for this, would be to model capital as a quasi-fixed input in the technology, e.g., following the proposals of Ouellette and Vierstraete (2004) or De Borger and Kerstens (2000).
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

Appendix
See the Tables 6 and 7  Table 6 Summary of bootstrap results for Malmquist index and its decomposition for textile and clothing subsectors, 1995-2004 (significant differences from 1 at 1 % level)