Computing Economies of Scope Using Robust Partial Frontier Nonparametric Methods

This paper proposes a methodology to examine economies of scope using the recent order-α nonparametric method. It allows us to investigate economies of scope by comparing the efficient order-α frontiers of firms that produce two or more goods with the efficient order-α frontiers of firms that produce only one good. To accomplish this, and because the order-α frontiers are irregular, we suggest to linearize them by the DEA estimator. The proposed methodology uses partial frontier nonparametric methods that are more robust than the traditional full frontier methods. By using a sample of 67 Portuguese water utilities for the period 2002–2008 and, also, a simulated sample, we prove the usefulness of the approach adopted and show that if only the full frontier methods were used, they would lead to different results. We found evidence of economies of scope in the provision of water supply and wastewater services simultaneously by water utilities in Portugal.


Introduction
The existence of inadequate incentives in monopolistic markets jeopardizes the maximization of the organizations' efficiency and reduces the value for money [1].Although governments intervene many times to curb the market power and correct market failures, either through regulation or policies that promote competition, the presence of unexploited economies of scope and economies of scale hinders the public interest.
Economies of scope exist when the costs of producing two or more goods together by the same entity are lower than the costs of producing them separately by several entities [2].Thus, economies of scope are related to savings from the joint production of different goods [3], for example, due to the sharing of fixed costs [4].However, specialization (single product firms) may have advantages (diseconomies of scope), mainly in more complex products.
The results found in the literature are very different concerning the existence or non-existence of economies of scope (see, for instance, in the water sector [5]).Mainly two types of economies of scope have been investigated: vertical integration of various stages in the production chain, for example, from production to distribution; and the integration of different activities, such as electricity, water, wastewater, and gas supply services.The literature also adopts different types of methodologies to estimate this phenomenon [6].To investigate scope and scale economies, researchers have been using a wide range of performance evaluation methodologies, which usually are divided into two major groups: parametric and nonparametric methodologies.Although most studies in the literature apply parametric methodologies, the nonparametric ones have several advantages.They do not need to define, a priori, any functional form, neither do they consider too many assumptions to estimate the production (cost) technology.The traditional nonparametric methods (e.g., Data Envelopment Analysis-DEA) were adopted by Färe [7] and Färe et al. [8] to estimate the efficiency of joint production and, since then, several variants have been applied to evaluate economies of scope in many sectors, such as local public services [9], banking [10], hospitals [11], insurance [12], railways [13], electricity [14], and water [15].Most of the nonparametric methodologies adopted in the literature until now have been solely based on the full frontier methods (mainly DEA).However, due to their limitations, some recent studies employ bootstrapping methods [16] and other techniques [17] to mitigate the influence of some of their negative aspects (e.g., [9]).
Thus, to improve the use of nonparametric methodologies to investigate economies of scope, this paper proposes a methodology that uses robust partial frontier nonparametric methods.These recent nonparametric methods use only part of the sample to compute efficiency scores, thereby reducing the influence of outliers and extreme observations [18].They do not face the problem of the "curse of dimensionality" and have important statistical features and interesting econometric interpretations [19], which makes them more robust than the traditional full-frontier nonparametric methods, such as DEA [20] and the non-convex free-disposal hull (FDH) [21].Given all these advantages, the methodology proposed is an important and useful tool for estimating economies of scope.
A comparative method was used by applying the methodology proposed by [22], which is also based on partial frontier nonparametric methods.This methodology aims to analyze the influence of exogenous variables on the production process.To search for the existence of economies of scope through this methodology, the "type of utility" (grouped according to the number of services provided) was the exogenous variable [23] considered.In addition, a statistical test proposed by [19] was applied to assess whether the difference between the efficiency averages of the groups of multiproduct and single product firms is statistically significant or not.However, the results were not good enough and are not given here.We observed that the p-value formula leads to the acceptance of the null hypothesis (equality between the efficiency averages of the two groups of firms), even for samples formed by groups of firms with very different efficiency averages.
To demonstrate the superiority of the proposed methodology we applied it to two case studies: a real case and a simulated case.The real case study involves the search for the existence of economies of scope in the water sector in Portugal between drinking water and wastewater activities.The sample comprises 67 water utilities for the period 2002-2008, which serves about 63% of the Portuguese population.However, given that the real sample is not composed of all types of utilities necessary to search for economies of scope (it does not include utilities providing only wastewater services), we also applied the proposed methodology to a simulated case of known characteristics to prove its superiority.
The remaining of the paper is organized as follows.Section 2 explains the methodologies adopted.Section 3 describes the case studies, the data and the results obtained, and ends with the comparative analysis of the methodologies.Finally, Section 4 provides the main conclusions.

Nonparametric Methods
This section may be divided by subheadings.It should provide a concise and precise description of the experimental results along with their interpretation, as well as the experimental conclusions that can be drawn.Although nonparametric methodologies are scarcely used in the literature to compute economies of scope, the current study adopted them because of their great advantages.As stressed earlier, they do not require any a priori assumption for the functional form of the frontier of the production (cost) technology.Furthermore, they allow for the use of multiple inputs and outputs, do not require information relative to input or output prices, identify best practices and are easy to explain to different stakeholders (see, for example, Fried et al. [24] for an overview of the benchmarking techniques).Nevertheless, nonparametric methodologies also have some drawbacks, especially the full frontier methods, such as DEA and FDH estimators.They are particularly sensitive to extreme data and outliers and suffer from the problem of the "curse of dimensionality", thus requiring the use of a large number of observations and a small quantity of inputs and outputs to avoid inaccurate results [19].Additionally, they are deterministic methods with a non-statistical nature, which means that all the deviations from the efficient frontier are seen as inefficiency and, therefore, they do not consider the existence of errors.
However, with the appearance of the probabilistic formulation of the production process proposed by [25] a new group of nonparametric methods gained prominence: the partial frontier nonparametric methods.These nonparametric methods, which include the order-α and order-m methods, use only part of the observations in the sample (α (probability), in the case of order-α and m observations, in the case of the order-m method) to determine efficiency scores.Therefore, they are less sensitive to extreme data and outliers and more robust than the traditional nonparametric full frontier methods.In addition, they do not suffer from the problem of the "curse of dimensionality", have important statistical features and an interesting econometric interpretation and allow for the easy inclusion of exogenous variables in efficiency computation [19].
According to Cazals et al. [25], the production process can be defined by the joint distribution function of inputs (X) and outputs (Y) and the efficiencies can be obtained from a conditional distribution function resulting from the decomposition of that joint distribution function.Later on, Daraio and Simar (2005) proposed an alternative probabilistic formulation to the production process suggesting the following probability function, which represents the probability of an observation operating at the level (x,y) being dominated: which can be decomposed, in a context of input orientation, into a conditional distribution function of X, F X|Y px|yq, and into a survivor function of Y, S Y pyq: This study adopted an input-orientation because the observations under study (water utilities) aim at rationalizing the quantity of inputs for a given level of outputs.Following Daraio and Simar [19], the order-α input efficiency score, for example, for an observation operating at the level (x,y) can be defined in terms of these probabilities as: which can be estimated by replacing the unknown conditional distribution function of X by its empirical version: FX|Y,n pθx|yq " where I(k) is the indicator function that takes the value of I pkq " 1 if k is true or I(k) = 0, otherwise.This conditional distribution function of X can be computed from the following expression ( The corresponding order-α frontier is defined as the input level not exceeded by p1 ´αq ˆ100% of the observations among the population producing at least a level y of outputs.The order-α method was introduced by [26] for the univariate case and later extended to the multivariate setting by Daouia and Simar (2007).
Thus, the order-α efficiency estimates are done by previously fixing the quantity (probability 1 ´α) of observations that are above the frontier.An observation operating at the level px, yq with an order-α input efficiency score θ α px, yq " 1 is efficient at the level α ˆ100%, being dominated by the observations producing more outputs than y with a probability 1-α.Additionally, an observation with an order-α input efficiency score θ α px, yq ă 1 has to reduce its inputs to the level θ α px, yqx to reach the input efficient frontier of level α ˆ100%.In contrast, if θ α px, yq ą 1 indicates that the observation px, yq can increase its inputs by factor θ α px, yq to reach the same frontier, being considered as super-efficient with respect to the order-α frontier.The order-α efficiency converges to the FDH efficiency score when α Ñ 1 [19].Figure 1 provides an example of a partial frontier and of a full frontier method.Observations A and B are efficient relatively to the full frontier and partial frontier, respectively.Observation C is only inefficient for the full frontier and observation D is inefficient for both frontiers.In this paper, we apply the order-α method because it is more appealing.Compared with the order-m method, the order-α method is more intuitive with respect to the amount (probability 1-α) of observations that are above the frontier and, among all the nonparametric methods, this is the least sensitive to extreme data [27].

Proposed Methodology
As mentioned above, economies of scope are related to savings originating from joint production of goods or services.This type of economies is related to the concept of subadditivity of the cost function and, therefore, to the concept of natural monopoly.A natural monopoly exists whenever a single firm produces a good or service at a lower cost than two or more firms that produce the same good or service separately, that is, whenever the cost function is subadditive [2].A cost function C pyq is said to be subadditive at the output level y if for any and all quantities of outputs y i such that ř k i"1 y i " y, if ř k i"1 Cpy i q ą Cpyq.The economies of scope imply a restricted form of subadditivity which occurs when there is orthogonality between the output vectors of the specialized firm, that is, if y i ¨yj " 0, i ‰ j.
According to [28], economies of scope can also be investigated from an efficiency point of view.That is, if we observe that, for example, the joint production of two or more goods or services is more efficient than their separate production, we can deduce that there are economies of scope.This is evaluated by comparing the joint efficient production frontier with the separate efficient production frontiers.If the efficient frontier of the joint production dominates the efficient frontier of separate production there are economies of scope, otherwise there are diseconomies of scope.Following the ideas of Morita [28], Cummins et al. [12], and Growitsch and Wetzel [13], we propose a methodology to test for economies of scope using the robust order-α frontier nonparametric method instead of the traditional (and less robust) full frontier methods.This proposed methodology intends to take advantage of all the benefits (described above) of robust partial frontier nonparametric methods and, more importantly, also attenuate the classic problem of the comparison of groups with different sample sizes (which induce different sizes for the bias) which occurs in the application of the full frontier nonparametric methods.This last target is achieved given the probabilistic nature of these robust partial frontier methods.
We considered the provision of two services (for example, water services (Ws) and wastewater services (Ww)) and the existence of three types of firms (n Ws firms that only provide the water service, n Ww firms that only provide wastewater services and n WsWw firms that provide both services simultaneously).To obtain the efficient frontier of separate production (S Ws`Ww ) we created n Ws ˆnWw fictitious firms (Ws+Ww) (for more details about how the fictitious firms are created, see Morita [28]) combining the real Ws and (mostly fictitious) Ww firms through the sum of their corresponding inputs and outputs.
For each firm with joint production (WsWw multiproduct firms), we determined the efficiency relative to the frontier of its own group (S WsWw ) (θ) and the efficiency relative to the frontier of the fictitious Ws+Ww firms (S Ws`Ww ) (θ j ) and compared them with the following ratio: This ratio allows us to assess which frontier dominates the other and, consequently, to conclude about economies or diseconomies of scope.
Suppose, for example, that the frontiers of the multiproduct WsWw firms and the fictitious Ws+Ww firms are the ones represented in Figure 2. Considering an input orientation, the multiproduct firm M P , for example, has an efficiency relative to the frontier of the fictitious Ws+Ww firms equal to θ MP,j " 0D1{0P1 and an efficiency relative to the frontier of the WsWw multiproduct firms equal to θ MP " 0C1{0P1, to which corresponds the ratio (8): Given that 0D1 ą 0C1, the multiproduct firm MP will have a ratio ηMP ą 1 indicating that its own frontier (S Ws +S Ww ) is located above the separate production frontier (S Ws S Ww ) and, therefore, its own frontier dominates the separate production frontier, meaning that this firm displays economies of scope.
For the multiproduct firm MQ, the ratio is: As 0F1 ă 0E1, the ratio to the multiproduct firm MQ is ηMQ ă 1, which means that its own frontier (S Ws S Ww ) is dominated by the separate production frontier (S Ws +S Ww ), so that the multiproduct firm MQ presents diseconomies of scope.Thus, multiproduct firms with a ratio η ą 1 have economies of scope and multiproduct firms with a ratio η ă 1 exhibit diseconomies of scope.If one assumes an output orientation, economies of scope can be evaluated in the same manner as above, relating the efficiency of the multiproduct firm frontier (θ) to the efficiency of the separate production frontiers (θ j ) in an output-output space: η " θ j {θ.In the same way as before, the multiproduct firms with a ratio η ą 1 exhibit economies of scope and multiproduct firms with a ratio η ă 1 have diseconomies of scope.See, for example, the M P multiproduct firm of Figure 3: it has an efficiency relative to the frontier of the separate production θ MP,j " 0P 2 {0D 2 and an efficiency relative to the frontier of the multiproduct firms θ MP " 0P 2 {0C 2 and, thus, a ratio ηMP " θ j {θ " 0C 2 {0D 2 ą 1, meaning that the M P firm displays economies of scope.Moreover, the multiproduct firm M Q has the ratio ηMQ " θ j {θ " `0Q 2 {0F 2 ˘{ `0Q 2 {0E 2 ˘" 0E 2 {0F 2 ă 1, meaning that this firm presents diseconomies of scope.Note that if one considers constant returns to scale (CRS), the output-oriented efficiencies coincide with the input-oriented efficiencies and, therefore, economies of scope may be valued also from an output-output space for an input orientation.To evaluate the consistency of economies of scope, we estimated the previous ratios for each multiproduct firm by using robust order-α efficiencies for different values of α.We estimated the various order-α frontiers of the groups of multiproduct firms and fictitious separate firms and compared them.However, since these frontiers are irregular (Figure 1), and to avoid situations as the one illustrated in Figure 4 within the red circle, where it is hard to assess the dominance of a frontier over the other, we applied the DEA method to linearize order-α frontiers of the two groups of firms (Figure 5) (instead of linearizing order-α frontiers as in Carvalho and Marques [29]; in which the convexity of the production set is not imposed).This procedure was performed by excluding the super-efficient firms (that is, the firms that are above the order-α frontier) in the corresponding order-α frontiers and then estimate the DEA frontier of firms only located below the order-α frontier.Thus, after the two linearized order-α efficient frontiers of the joint and separate production for several values of α have been estimated, we computed the DEA efficiencies θ of multiproduct firms relative to their own linearized frontier (joint production frontier-S WsWw ) and their DEA efficiencies θ j relatively to the linearized frontiers of the separate production (S Ws`Ww ) (considering as the reference set the group of separate firms located below the corresponding order-α frontier) and afterwards compared by the ratio η.For example, the efficiency θ j relative to the separate production frontier (S Ws`Ww ) for the multiproduct firm MP was estimated from the following linear programming problem: min θ j s.t.
where x il and y rl are inputs and outputs of n separate firms (located below the order-α frontier), m is the number of inputs, s the number of outputs, and X iP and Y rP are the inputs and outputs of WsWw multiproduct firm MP, respectively.In the case of considering variable returns to scale (VRS) the additional restriction ř n l"1 λ l " 1 is imposed.As highlighted by Chavas and Kim [30] economies of scope are also influenced by economies of scale.Thus, to evaluate the possible influence of the economies of scale, we applied the DEA method adopting a CRS technology and then the DEA method adopting a VRS technology.To overcome the problem of infeasibility in efficiency calculation when considering VRS [31], we used the methodology proposed by [32], which adopts both input and output orientation super-efficiency models to fully characterize the firms' super-efficiency.
The methodology here suggested is robust and very useful, since in most real situations the samples are contaminated by extreme observations and outliers, sometimes leading to biased conclusions.By using this methodology, we can assess if the conclusions are affected by extreme observations or outliers and check if the results are consistent by varying the values of α and evaluate the existence of economies of scope individually.

Conditional and Unconditional Efficiency Scores
To compare with our method, we applied the methodology proposed by Daraio and Simar [22], which is also based on the order-α method.This methodology aims to analyze the influence of exogenous variables on the production process and consists in a linearized nonparametric regression of the ratios between the conditional efficiencies order-α and the unconditional efficiencies order-α, Q z α " θα,n px, y|zq{ θα,n px, yq, on the exogenous variable (usually referred to as Z).The influence of exogenous variables on the production process is interpreted according to the evolution of this linearized nonparametric regression all over the values of Z.In an input-orientation context (as in the current study) when the linearized nonparametric regression shows a growing tendency as the values of Z increase, it means that the exogenous variable penalizes the efficiency.Conversely, if the regression presents a decreasing slope, it means that the exogenous variable is favorable to efficiency.
In order to obtain the conditional efficiencies, it is enough to constrain the production process (defined by the expression (1)) to a given value of the exogenous variable Z = z [22]: Allowing us, subsequently, to obtain an estimator for the conditional order-α input efficiency: where FX|Y,Z,n pθ x|y, zq " ř n i"1 I pX i ď θx, Y i ě yq K ppz ´Zi q{hq{ ř n i"1 I pY i ě yqK ppz ´Zi q{hq, I pkq is the indicator function that takes the value of I pkq " 1 if k is true or I pkq " 0 otherwise, K p¨q is the kernel function for continuous variables, and h the appropriate bandwidth.This empirical version of the conditional distribution function of X is equivalent to [19]: where M y , χ i , and χ y pjq are described previously.Additionally, Z y rjs represents the observation Z i corresponding to the order statistic χ y pjq , Q y,z " ř n i"1 I pY i ě yqK ppz ´Zi q {hq, and So the nonparametric estimator of the conditional order-α input efficiency measure given that Z = z is: Obtaining conditional efficiencies involves the estimation of the non-standard conditional distribution function FX|Y,Z,n , which requires the use of linearizing techniques for the exogenous variables.Such linearizing techniques still require the choice of a kernel function and the determination of a bandwidth.In this research, we used the following kernel function proposed by [33] for the environmental variables (which are "unordered" categorical variables) instead of K ppz ´Zi q {hq (which is used for continuous variables): To obtain the optimal bandwidths, we used the likelihood cross validation based on a k-nearest neighbor (k-NN) method.This method consists in evaluating the leave-one-out kernel density estimate of Z ( f p´iq k pZ i q) for a set of values of k and find the value of k which maximizes the function: where f p´iq k pZ i q " 1{ rλ pn ´1qs ř n j"1,j‰i K `Zj , Z i , λ ˘, where λ is the local bandwidth and chosen so that there will be k points Z j verifying ˇˇZ j ´Zi ˇˇă λ [22].
Thus, to find economies of scope, the "type of utility" was considered as an environmental variable and its influence on the production process was examined by analyzing the linearized nonparametric regression of the ratio Q z α " θα,n px, y|zq{ θα,n px, yq on the environmental variable (Z).

Economies of Scope and the Portuguese Water Sector
In Portugal, the drinking water supply and wastewater collection and treatment are under the responsibility of the municipalities.They can provide these activities directly or indirectly by municipal companies (which can have a private partner as minority shareholder) or by the private sector through concession contracts.There are utilities that provide only water, and others providing both water and wastewater.Some utilities also include urban waste and some also provide other activities, such as transportation, gardening, and swimming pools.The predominant model is the water and wastewater utility.However, for several years it has been discussed what the optimal market structure is, thus acknowledging the importance of economies of scale and scope.In 1993 the water sector was reformed and was unbundled (separating the production from the distribution in water and collection from treatment in wastewater) spending millions of euros to specialize the utilities.In 2001 the sector was reformed again and, from then on, the orientation was to amalgamate the utilities and not to specialize them.Again, millions of euros were spent.Recently the last government has announced the reform (again) of the water sector market structure, aiming at its optimization (see Marques [34] for a detailed analysis of the Portuguese water sector).
Economies of scope are quite important for utilities, in general, and for water utilities, in particular.In addition to the similar features of network industries (water, wastewater, electricity, gas, etc.) and the fact that they are under the responsibility of the same public body, i.e., the municipality, there are additional potential savings concerning the final customers of the water utilities (metering, billing, collection, etc.).These are some of the reasons why the model of multi-utilities is so popular in certain countries (e.g., Germany or Switzerland).In Portugal, as in other countries, water utilities are strictly regulated and the regulator defines targets for their performance.The targets should not be the same for a water utility or for a water and wastewater utility.Therefore, the existence of (dis)economies of scope should be accounted for due to their significant economic consequences in the tariff setting.These are some of the motivations to investigate the presence of economies of scope in the water utilities and particularly in the Portuguese water sector.
So, from the proposed methodology we investigated the presence of economies of scope among water utilities (Ws) and water and wastewater utilities (WsWw) in Portugal.The sample is made up of 67 water utilities of the retail segment for the period between 2002 and 2008, covering almost 63% of the Portuguese population.Altogether the sample is composed of 39 WsWw utilities, 10 Ws utilities and 18 water utilities that provide water supply, wastewater, and other services, such as solid waste collection.From the latter utilities, and from the WsWw utilities, were built Ww fictitious utilities, drawing from these utilities the inputs and outputs corresponding to the wastewater activity.This simplification was done in order to apply the proposed methodology, which requires the existence of all types of utilities in the market study: utilities that provide only one of the outputs and utilities that provide all outputs.Note that over the period studied some water utilities have changed concerning the services provided.Figure 6 shows the distribution of observations according to the types of services provided during the period studied.In this case study we selected three "consensual" input variables [35,36]: labor cost, other operational costs, and capital costs (in thousands of euros); and the following outputs: volume of delivered water, volume of treated and collected effluents in wastewater service (in thousands of m 3 ) and the number of water customers and of wastewater customers.The output "volume of delivered water" gathers the volume of delivered water in the retail and wholesale segments.The outputs "volume of treated effluents" and "volume of collected effluents" are also aggregated into a single output.Further analysis showed that the costs associated with the wastewater treatment are usually approximately two times higher than the collection costs, which led us to create the following composite output: "1/3 volume of collected effluent +2/3 volume of treated effluent".Tables 1 and 2 describe some statistics of these input and outputs variables for the Ws and WsWw utilities, respectively.All of the data were obtained from the annual activity and account reports published by the water utilities and from annual reports of the Portuguese regulator.In Portugal there is only one utility responsible for the wastewater service (Ww utilities).Therefore, we had to assume additional fictitious Ww utilities that were constructed as follows: from each water utility of the sector that provides the wastewater service along with other services, we removed the inputs and outputs corresponding to the wastewater service.Later, it was assumed that these inputs and outputs are the variables that characterize the activity of these fictitious Ww utilities.This "drawing" of inputs and outputs was performed in a careful manner in order not to lead to biased results, especially the costs (inputs) which in Portugal are not displayed per activity undertaken by water utilities.That is, this "drawing" of inputs and outputs was carried out bearing in mind that if there is an underestimation of the costs associated with the activity of wastewater it can lead to diseconomies of scope (benefiting the Ws+Ww utilities in terms of efficiency compared with the WwWs utilities) or overestimation can lead to economies of scope (benefiting the WsWw utilities in terms of efficiency).So, this estimation of costs associated with the wastewater activity was performed taking into account the extent of wastewater activity in each utility represented by the amount of treated and collected effluents and by the revenue generated by the wastewater activity, and also bearing in mind a study [37] which showed that in Portugal, on average, in the WsWw utilities, 36% of the total costs are allocated to the activity of wastewater and the remaining 64% allocated to the water supply.This study applies the nonparametric method (a shared input DEA model) that allows to estimate for each utility the share of the total costs allocated to each service or output produced and to evaluate the respective partial cost efficiencies associated with each service (see, for example, [38,39]).Table 3 presents some statistics of inputs and outputs for these fictitious Ww utilities.

Simulated Case Study
Since the sample of the real case study of the Portuguese water sector does not have the ideal characteristics to search for economies of scope through the proposed methodology, because it does not include all kind of utilities, we also incorporated in this paper a simulated case study with known characteristics.
In this simulated sample it was assumed that in a particular sector there are utilities A (50 utilities) providing the product A, utilities B (40 utilities) providing the product B, and utilities AB (100 utilities) providing both products.All utilities use two inputs (X 1 and X 2 ) (for example, labor costs and capital costs) to produce the outputs A and B and the different groups of utilities (A, B, and AB) face different technologies, but all of them governed by a Cobb-Douglas function, with the following expressions, respectively: For utilities A and B it was assumed that all inputs are randomly drawn from a uniform distribution between [0.1, 2] while for the utilities AB it was assumed that all inputs are randomly drawn from a uniform distribution between [0.25,5].u is a small perturbation which is drawn from an exponential distribution with mean 1/50 for all groups of utilities.The coefficients and the exponents in the above Cobb-Douglas functions were assigned so that small size utilities AB would present economies of scope and large size utilities AB would present diseconomies of scope (as it is usually observed in the real world), and in such a way that diseconomies of scope could occur in this case study.For diseconomies of scope to occur it was assumed that, on average, the integrated utilities AB use more inputs and produce less outputs than the disintegrated utilities (A+B) (Table 4) or, in other words, on average, the integrated utilities AB are less efficient than the disintegrated utilities A+B.
Both in this case study, and in the previous one, an input orientation was adopted, since the aim of the utilities was to decrease the inputs consumed for a given level of outputs produced.

Proposed Methodology
For the case study of the Portuguese water sector the results show that the frontier of the WsWw utilities dominates mostly the frontier of the Ws+Ww utilities for several α values, since the average of η (as is its median) is always greater than one (Table 5).This means that WsWw utilities with economies of scope predominate (considering the order-α quantile frontier (α = 0.999) of the separate utilities and the order-α quantile frontiers (α = 0.99, α = 0.995 and α = 0.999) of the group of joint production firms).For example, for the case where the order-α quantile frontier of the separate utilities (Ws + Ww utilities) is α = 0.999 and the order-α quantile frontier of the group of joint production firms (WsWw utilities) is α = 0.995, on average, WsWw utilities must reduce their inputs by about 2.8% (1-0.972) to reach the frontier of separate utilities (Ws + Ww utilities); that is, to become as efficient as separate utilities (Ws+Ww utilities), while reaching their own frontier, they will have to reduce their inputs by about 17.0% (1-0.830).Conversely, if the order α-quantile frontier of the group of joint production firms (WsWw utilities) move to an α = 0,999, on average, WsWw utilities need to further reduce their inputs by about 22.8% (1-0.772) to become efficient in relation to their own group of utilities (WsWw ).This means that WsWw utilities are, on average, more efficient than separate utilities (Ws+Ww utilities), which means that there are economies of scope in providing these two types of services for these values of α.
Table 5. Statistics of order-α input efficiencies scores and ratios η for WsWw utilities obtained in the analysis of economies of scope for several α values and for CRS technology (for the order-α (α = 0.999) quantile frontier of the Ws+Ww utilities).As we can see in Table 5, as the value of α increases, tending to 1, that is, as we move towards the full frontier (FDH), the average of ratios η increases, being always higher than unity, thereby increasing the percentage of utilities (Table 5) and observations (Figure 7) with economies of scope.This sensitivity analysis, making the values of α change, proves the consistency of the results.An analysis using only full frontier methods leads to the conclusion that in the Portuguese water sector most of the WsWw utilities (76% corresponding to 34 WsWw utilities) present economies of scope, and the remainder (24% which corresponds to 11 WsWw utilities) have diseconomies of scope (Table 5).However, an analysis using the latest and robust partial frontier nonparametric methods shows that when we neglect the extreme utilities (which constitute the full FDH frontier), the amount of WsWw utilities with economies of scope decreases (for example, to around 64% when considering α = 0.995).
A sensitivity analysis for the Ws+Ww frontier was carried out, but as in this case study the sample of observations Ws+Ww is extremely compact, for α values less than 0.999, a huge amount of observations were also above the frontier, leading to a situation whose relevance does not justify being studied.
The simulated case study proves, as expected, the predominance of diseconomies of scope, since the ratios η prevail for various order-α partial frontiers (Table 6 and Figure 8).As in the previous case study, results for α ă 0.987 partial frontiers of utilities AB and for α ă 0.999 partial frontiers of utilities A+B, are not given because in these cases a large part of the utilities are above the corresponding frontiers and, for this reason, they have no interest.From the proposed method it is demonstrated that as the frontier of utilities AB moves towards the full frontier ( α Ñ 1 ) will dominate the frontier of utilities A+B and, consequently, the number of utilities with economies of scope increases.According to the results in a CRS situation and when 9% of utilities AB are above the frontier of the joint production (α " 0.989), 80% utilities AB show diseconomies of scope (Table 6).However, when the frontier of multioutput utilities AB moves towards the full frontier (for the α " 0.99), leaving only 1% of the utilities AB above the frontier, the number of utilities AB with diseconomies of scope drops to about 45% (Table 6).An analysis considering a VRS situation proves that, in fact, the diseconomies of scope dominate as it gets closer to the full frontier 6 and Figure 8).This means that if partial frontiers had not been applied in a CRS situation and, if only full frontiers were applied (such as DEA), we would have concluded the opposite, the predominance of economies of scope, which would not correspond to reality.After an analysis of the sample it was found that 9% of the utilities AB (which are above the frontier α " 0.989) are actually extreme observations which consume a very small amount of inputs and also produce a low level of output; therefore, considering them leads to the opposite conclusions.

Conditional and Unconditional Efficiency Scores
To test for economies of scope among Ws and WsWw utilities in the Portuguese water sector, we carried out two analyses.First, we analyzed the influence of the exogenous variable "type of utility" on the production process considering the previous two groups of utilities, the group of multiproduct utilities (WsWw) to which corresponds the categorical variable with value 1, and a group of fictitious utilities (separate firms) (Ws+Ww) to which corresponds the categorical variable with value 2. In a second analysis, we examined the influence of the exogenous variable ("type" of utility) just for the real utilities (Ws and WsWw utilities).We considered a categorical variable with values: 1 for Ws utilities and 2 for WsWw utilities.To perform this last analysis it was necessary to build the following two compound outputs which gather outputs from each of the activities provided: "Delivered water volume (retail and wholesale) + 1/3 volume of collected wastewater + 2/3 volume of treated wastewater" (these ratios were assumed since in Portugal the wastewater treatment activity is roughly a weight cost twice the activity of wastewater collection) and "Water customers + Wastewater customers".
According to the results presented in Figures 9 and 10 the simultaneous provision of water and wastewater services is favorable to performance (at the 99% confidence level), displaying the existence of economies of scope between these activities.When studying the influence of the exogenous variable "type of utility" on the production process of the WsWw multiproduct utilities (1) and the group of fictitious utilities (separate utilities) (Ws+Ww) (2), the existence of economies of scope is clearer.The presence of economies of scope becomes more pronounced (steeper nonparametric regressions) as we move away from a unit value for α.This was also observed in the second analysis when the influence of the exogenous variable on the production process of the real utilities, Ws utilities (1) and WsWw utilities (2) were investigated.However, this second analysis leads to opposite conclusions for higher values of α or, in other words, when it coincides with the full FDH frontier.This may be directly related to the fact that two different types of utilities are compared in this second analysis, particularly a group of multiproduct utilities (WsWw) with a group of utilities that provide only one service (water supply) and that can, therefore, be responsible for skewing the results.This means that this analysis may not be the most appropriate one to find economies of scope, because it compares two completely different types of utilities, some utilities are characterized by an output while others are characterized by two outputs.Regarding the simulated case study, the application of this methodology has also concluded that utilities A and B are more efficient than utilities AB (Figure 11) and, therefore, diseconomies of scope dominate.However, this finding is extracted from the frontiers α ď 0.989, as it is not possible to draw any conclusions from the frontiers near of full frontier (α -1) about the existence or absence of economies of scope, since the nonparametric regression is horizontal.This fact highlights the great usefulness of partial frontiers in the search for economies of scope.The use of this methodology to look for economies of scope needs to be done with care.See, for example, the above analysis in which economies of scope are investigated between utilities A and AB.In the computation of order-α efficiencies, the utilities A can be considered less efficient than utilities AB because they have zero output from output B. To overcome this drawback, we investigated economies of scope using this methodology among the groups of utilities AB and A+B.However, the results do not allow us to prove clearly the existence of scope diseconomies.Only for α " 0.999 case (frontier near the full frontier) is it possible to conclude that there are scope diseconomies but it is still difficult to prove the existence of diseconomies, given the low slope of the non-parametric regression (Figure 12).

Discussion
The results obtained are very encouraging demonstrating the potential of the proposed methodology compared with the other methodologies usually applied in the literature of nonparametric methods.Although the second methodology (which allows the analysis of the influence of the exogenous variables on the production process) is not the most appropriate one to investigate economies of scope, it has also been used in the literature to explore such economies [40].We argue that the methodology proposed is the most useful since it enjoys all the advantages of the recent partial frontier nonparametric methods.It is a robust tool to investigate economies of scope because it allows us to assess the influence of extreme observations and outliers on the results and does not face the problem of the "curse of dimensionality", especially when α is farthest from the unit.In addition, our methodology is closer to the definition of economies of scope because it compares the production function of firms that produce two or more goods with the production function of firms that produce only one good.Moreover, not only does it provide an overview of the whole sample but it also enables the evaluation of each observation, and allows us to further assess the robustness of the results obtained by varying the value of α.The results obtained in the simulated case study proves the usefulness and potentiality of the methodology.These results show that if only full-frontier non-parametric methods (such as DEA) were used, we would have concluded quite the opposite.That is, we would have concluded for the predominance of economies of scope when, in fact, diseconomies of scope prevail.The results obtained by the proposed methodology demonstrate that this is caused by 9% of extreme observations, which use a very small amount of inputs and also produce a low level of outputs.
Regarding the case study of the Portuguese water sector, the results obtained by the methodology proposed show the existence of economies of scope between water and wastewater activities in Portugal.These results are in line with most studies in the literature, which find economies of scope between water and wastewater services (e.g., [41][42][43][44]), especially when dealing with small companies.The fact that the water utilities sector in Portugal is characterized by many small water utilities which are vertically unbundled (not integrated) in the retail segment might be one of the main justifications for the existence of economies of scope.The management model of the water utilities in the retail segment in Portugal, contrary to most countries around the world, follows a model which is vertically unbundled (not integrated); that is, in most cases there are different actors supplying the "wholesale" segment (which includes intake, treatment, and transportation in water supply, and transportation and wastewater treatment in the wastewater activity) and the "retail" segment (which includes water distribution in water supply and wastewater collection in the wastewater activity).This last fact may also be one more justification for the existence of economies of scope.Thus, the joint provision of water supply and wastewater collection services may provide a cost advantage, especially in the sharing and joint utilization of resources (inputs), such as human resources (management, administrative labor, and so on), machinery, and information technology, the sharing of technical knowledge, and better coordination of management and maintenance of both networks (which are sometimes interconnected), as well as savings in fixed costs of production, such as purchasing and marketing [45].

Conclusions
This paper investigates the use of nonparametric methodologies to estimate the existence of economies of scope, proposing a methodology based on partial frontier methods.These methods are more robust than the full frontier methods, mainly because the partial frontier methods are less sensitive to extreme data and outliers and do not face the problem of the "curse of dimensionality".
The results of this research show that our methodology to investigate economies of scope is more robust than the ones that use full frontier methods.Firstly, it is closest to the true concept of economies of scope and, secondly, it uses the most robust nonparametric methods, also allowing us to assess the robustness of the results obtained by varying the value of α.In addition, this methodology allows us to evaluate the existence of economies of scope individually for each observation.
Concerning the Portuguese water utilities case study, we found the existence of economies of scope among water utilities and water and wastewater utilities, proving that there could be cost reductions by fostering synergies in providing these two services simultaneously in Portugal.For example, for the order-α quantile frontier of the separate utilities (α = 0.999) and for the order-α quantile frontier of the group of joint production utilities (α = 0.995), on average, WsWw utilities must reduce their inputs in only about 2.8% to become as efficient as separate utilities, while to reach their own frontier they will have to reduce their inputs in about 17.0%.
Therefore, this paper provided evidence of the existence of economies of scope in the water sector in Portugal.Among other reasons, this can be very helpful for decision-makers to reform the water sector since the market structure in Portugal is considered very inefficient.In a period when it is very important, and even compulsory, to recover all the costs by the water utilities, it is socially and politically key to mitigate the impact of tariff increases and market structure reform can be a good contributor for this.
Thus, the policy implications of this investigation refer to the need of Portuguese policy-makers to promote the merger between the wholesale and the retail utilities, the provision, whenever possible, of water and wastewater together and encourage small utilities, for example, with less than 50,000 inhabitants, to provide other services besides water supply and wastewater collection and treatment, such as solid waste, street cleaning, transportation, and gardening, among others.

Figure 1 .
Figure 1.Example of a partial frontier (order-α) (blue) and of a full DEA frontier (orange).

Figure 3 .
Figure 3. Separate production frontiers (S Ws`Ww ) and joint production frontiers (S WsWw ) in the output-output space.

Figure 4 .
Figure 4. Example of order-α joint production and separate production frontiers.

Figure 5 .
Figure 5. Linearized frontiers of joint and separate production.

Figure 6 .
Figure 6.Distribution of observations according to the types of services provided for the period 2002-2008.

Figure 7 .
Figure 7. Percentage of observations with economies and diseconomies of scope, considering CRS (a) and VRS (b) technologies.

Figure 8 .
Figure 8. Ratios η as a function of the output YA for various partial frontiers of utilities AB, for CRS (a) and VRS (b) technologies.

Table 1 .
Statistical features of the model variables relative to water utilities (Ws).

Table 2 .
Statistical features of the model variables relative to water and wastewater utilities (WsWw).

Table 3 .
Statistical features of the variables concerning the activity of collection and treated wastewater.

Table 4 .
Statistical features of the variables used in the simulated case study.

Table 6 .
Average of ratios η and percentage of utilities AB with economies and diseconomies of scope for various quantile frontiers of the utilities AB and for the quantile frontier (α = 0.999) of the utilities A+B.