New Classes of Distortion Risk Measures and Their Estimation

: In this paper, we present a new method to construct new classes of distortion functions. A distortion function maps the unit interval to the unit interval and has the characteristics of a cumulative distribution function. The method is based on the transformation of an existing non-negative random variable whose distribution function, named the generating distribution, may contain more than one parameter. The coherency of the resulting risk measures is ensured by restricting the parameter space on which the distortion function is concave. We studied cases when the generating distributions are exponentiated exponential and Gompertz distributions. Closed-form expressions for risk measures were derived for uniform, exponential, and Lomax losses. Numerical and graphical results are presented to examine the effects of the parameter values on the risk measures. We then propose a simple plug-in estimate of risk measures and conduct simulation studies to compare and demonstrate the performance of the proposed estimates. The plug-in estimates appear to perform slightly better than the well-known L-estimates, but also suffer from biases when applied to heavy-tailed losses.


Introduction
Risk measures or premium calculation principles are statistical tools for the calculations of the insurance price corresponding to a risk in the actuarial literature; see Fischer et al. (2018) for a discussion on recent risk measures with application to credit risk.They have also been developed to assess the reserve capital required to cover unexpected losses and ensure financial stability.For example, the value at risk (VaR) at level p is a quantile-based risk measure such that VaR p (X) = inf[x|F(x) ≥ p] = F −1 (p), 0 ≤ p ≤ 1, where X is a non-negative loss or risk random variable with a cumulative distribution function (cdf) F(X) = P(X ≤ x), x > 0. The conditional tail expectation (CTE) is another risk measure defined as the average value of losses beyond the VaR p value, i.e., CTE p = E[X|X ≥ VaR p (X)].Another popular premium principle is the Gini shortfall; see Furman et al. (2017) and Eugene et al. (2021).It is given by GS p = CTE p + λTGini p , where λ is the loading parameter and TGini p = 2 1 p F −1 (u)(2u − 1 − p)du/(1 − p) 2 .It is said to be more comprehensive as it combines the average serenity and the variability of a loss distribution tail.Artzner et al. (1997Artzner et al. ( , 1999) ) addressed the desired behaviors and characteristics of risk measures.More specifically, here, we are concerned about the coherency.Let ρ(X) be a risk measure associated with a loss random variable X.A risk measure is coherent if it satisfies the following four axioms: (i) monotonicity: if Y ≤ X ⇒ ρ(Y) ≤ ρ(X); (ii) subadditivity: ρ(X + Y) ≤ ρ(X) + ρ(Y); (iii) positive homogeneity: for any c > 0, ρ(cX) = cρ(X); (iv) translation invariance: for any c > 0, ρ(X + c) = ρ(X) + c.Subadditivity specifies that a diversified portfolio reduces the overall risk profile.It is well known that the VaR does not satisfy the subadditivity axiom and, hence, is not coherent; see Denuit et al. (2005) for Risks 2023, 11, 194.https://doi.org/10.3390/risks11110194https://www.mdpi.com/journal/risksexamples.In contrast, the CTE developed by Rockafellar and Uryasev (2002) is a coherent risk measure.
A distortion risk measure introduced by Denneberg (1994) calculates the insurance premium by transforming or distorting the accumulative distribution function of the loss variable.A distortion function is a non-decreasing function mapping the unit interval [0,1] to the unit interval such that g(0) = 0 and g(1) = 1.The risk-adjusted distortion risk measure denoted by ρ or ρ(X) for a continuous loss X is given by where the survival function S(x) = 1 − F(x), f (x) = dF(x)/dx, and g (•) is the first derivative of g(•).Based on (1), the distortion risk measure ρ can be interpreted as the mean of a random variable Y with cdf 1 − g(S(x)).It reflects the notion in Yaari (1987) that probabilities are liable to be in the decision-maker's perception (see Pflug 2009) and that the distortion risk can be seen as the expected utility with a utility function of xg (S(x)) with respect to the loss distribution.According to (2), the distortion risk measure is also a spectral risk measure, which is a quantile-based risk measure that takes the form of 1 0 F −1 (t)ψ(t)dt, where ψ(t) ≥ 0, 1 0 ψ(t)dt = 1, and ψ (t) ≥ 0. The function ψ represents the user's risk attitude; see Acerbi (2002) and Dowd et al. (2008).For the applications of distortion risk measures, see Sereda et al. (2010) and Bihary et al. (2020).Wang (1995Wang ( , 2000) ) proposed two classes of distortion operators or transformation for pricing financial and insurance risk: power distortion or the proportional hazards transform defined by g(w) = w a , 0 ≤ a ≤ 1, and the Wang transform defined by g(w) = Φ Φ −1 (w) + λ , where Φ(•) is the standard normal cdf and λ is a scalar.The power distortion brings about the dual-power distortion function g(w) = 1 − (1 − w) b , b ≥ 1.It is shown that, if the distortion function is concave, then the resulting risk measure is coherent.Wirch and Hardy (1999) considered the beta distribution distortion given by where Note that the beta distortion is concave when 0 ≤ a ≤ 1 and b ≥ 1.The distortion utilizes the beta probability density function (pdf) as its generating function.Setting a = 1, the beta distortion yields the dual-power distortion.It is the power distortion when b = 1.Applying the framework in (3), Samanthi and Sepanski (2018) proposed new classes of distortion functions by replacing the beta pdf with other pdf's.Yin and Zhu (2018) presented three methods for constructing new classes of distortion functions: compositing two distortions, convex linear mixing of distortions, and employing copula cdf.Minasyan (2020) studied two new classes of financial risk measures defined by the power function of the VaR and CTE.Minasyan (2021) introduced the concept of variance distortion, i.e., distorting the variance instead of the mean, by using popular existing distortion functions that yield the VaR and CTE.More references can be found in Minasyan (2021).
Let X 1 , X 2 , • • • , X n be independent and identically distributed (iid) loss random variables with the pdf, cdf, and survival function given by f , F, and S, respectively.Grounded by (2), Jones and Zitikis (2003) and Jones and Zitikis (2007) proposed the following empirical L-estimator, which is a linear combination of order statistics with weights being the score function g(1 − t) : The estimator is nonparametric, intuitive, and simple to implement.Due to the scarcity of data in the tails for a heavy-tailed loss, bias-correction methods have been proposed.For example, Brahimi et al. (2012) introduced a new reduced-biasestimator for heavy-tailed losses; see the references therein for more bias-correction methods.
The organization of the paper is as follows.In Section 2, we stage the method employed to produce new distortion functions and demonstrate the method by employing exponentiated exponential and Gompertz distributions as the generating functions.In Section 3, closed-form expressions for distortion risk measures are derived for uniform, exponential, Lomax, and Weibull loss distributions.Numerical and graphical comparisons are also included.In Section 4, we propose a plug-in, empirical-distribution-based estimator of the risk measures.The estimator does not require parametric assumptions on the loss distribution.Simulations were conducted to compare and demonstrate the performance of the proposed estimator.Concluding remarks are made in Section 5.

Proposed Methods and Distortions
We first advance the method, similar to those in Aldhufairi and Sepanski (2020) and Aldhufairi et al. (2020), for constructing a distortion.The main idea stems from the fact that a cdf with a domain of the unit interval by definition is also a distortion function.
Let Y be a non-negative continuous random variable with cdf G. Consider the following two transformations of the random variable Y: (5) The transformed random variables W and W 1 both have a support of I = [0, 1], and their respective cdfs are given by for w ∈ I.Both H(•) and K(•) are continuous and non-decreasing with H(0) = K(0) = 0 and H(1) = K(1) = 1.Let G (y) = dG(y)/dy be the pdf of Y.Note that K(w) in (7) can also be derived from the framework: motivated by the framework in (3); see Samanthi and Sepanski (2018).Since W or W 1 is continuous, its quantile function or inverse cdf is non-decreasing and may also potentially serve as a distortion for the purpose of constructing new risk measures.While the cdfs and inverse cdfs of W and W 1 are distortions, they must be concave to produce coherent risk measures.To ensure the concavity, we restricted the parameter spaces where the second derivative of the functions is non-positive.
We next showcase four new admissible distortions deriving from the cases when Y has an exponentiated exponential distribution and a Gompertz distribution.A distortion function is said to be admissible if it yields coherent risk measures, i.e., if it is concave.The parameter space over which the distortion is concave is presented for each of the four distortions.

Distortions via Exponentiated Exponential Distribution
Let Y be an exponentiated exponential random variable.Its cdf is G(y) = (1 − e −αy ) θ , where y > 0, θ > 0, α > 0. We refer to cdf G as the generating distribution.Then, the cdfs of the transformed variables in (6) are given by The respective inverse cdfs of H e (w) and K e (w) are In this case, H e and K −1 e are of the same functional form with different parameterizations, so are K e and H −1 e .The distortion H e (•) can also be formulated with the Kumaraswamy pdf instead of the beta pdf in (3), which was studied in Samanthi and Sepanski (2018).The distortion H e will be cited as the Kumaraswamy distortion below.
Since the support of transformed variables in (5) is the unit interval, we will refer to K e as the unit-exponentiated exponential (UEE) distortion below.Note that the power and dual-power distortions are special cases of the beta distortion in (3), the Kumaraswamy and UEE distortions.
The notations of the parameters in (8) or K e are changed in the following lemma so that, throughout the paper, θ consistently represents a parameter with a constraint of greater than or equal to one.
Proof.The first and second derivatives of K e (w) with respect to w are given by . In this case, if α ≤ 1, K e is non-positive for all w in I.If θ < 1, there is no α > 0 value such that K o (w) ≤ 0 for all w in I.

Distortions via Gompertz Distribution
Let Y be a Gompertz random variable with cdf G(y) = 1 − e −θ(e αx −1) , y > 0, α > 0, θ > 0. The resulting distortions H o (w) and K o (w) from ( 6) and ( 7) and their inverse functions are The first and second derivatives of H o (w) are, respectively, There exists no admissible parameter space for α and θ on which H o (w) is non-positive for all w ∈ I. Similarly, one can show that the same conclusion holds for K −1 o .The cdfs of the transformed variables defined in (5) are distortions, but may not be admissible, e.g., H o .When G is the Weibull cdf, the resulting distortions are not admissible, though not shown here.
In what follows, the distortion K o will be called the unit Gompertz (UG) distortion and H −1 o the unit Gompertz quantile (UGQ) distortion.
Proof.The respective first and second derivatives of K o (w) are given by Proof.The respective first and second derivatives of In summary, applying the exponential transformation to negative random variables with the exponentiated exponential and Gompertz distributions, we obtain four new admissible distortions.Table 1 summarizes the four lemmas in this section.
Table 1.Proposed distortions and admissible parameter spaces.

Distortion Function Form Admissible Parameter Space
Kumaraswamy Figure 1 shows how the two new distortions UG and UGQ behave at varying θ and α values.Assuming that a risk-neutral agent would not distort the survival distribution, Belles-Sampera et al. (2016) used the area under a distortion function as a measure of global risk attitude.The area under a concave curve on [0, 1] is always greater than half, and a larger area indicates a higher level of global risk-tolerant attitude.
For the UG distortion, when α = 0.5, Graph (a) indicates that a higher level of global risk-tolerant attitude and risk aversion is associated with a larger θ value.Setting θ = 6 (Graph (b)), a larger α value corresponds to a larger level of global risk-tolerant attitude.When α = 10 and θ = 6, the area under the curve is close to 1, which reflects a very conservative global risk-tolerant attitude.Similarly, for the UGQ distortion, the choice of a larger θ at a fixed α value or a larger α at a fixed θ value reflects a higher level of global risk-tolerant attitude.
From (2), the derivative of the distortion g (w), 0 < w ≤ 1, i.e., the slope of the tangent line to a distortion curve at w, is the assigned weight to the loss of S −1 (w).For example, when α = 0.5, at a small extreme survival value w, the slopes of the tangent lines to the curves increase as θ increases.That is, a greater weight is assigned to a large extreme loss as θ increases, indicating a higher level of risk aversion.One, therefore, expects to obtain a larger distortion risk measure for a larger θ; see Tables 2 and 3 below

Examples of Distortion Risk Measures
In this section, we calculate the distortion risk measures defined in (1) for uniform, exponential, Lomax, and Weibull losses using the distortions in Table 1.The derivations of the Kumaraswamy and UEE risk measures for the exponential and Lomax losses can be found in Samanthi and Sepanski (2018).

Uniform Loss
When the loss variable is uniform over the interval [0, 2b], then the respective cdf and survival function are F(x) = x/2b and S(x) = 1 − x/2b for 0 ≤ x ≤ 2b.In this case, we use the formula of ρ = g S(x) dx in (1) to calculate the distortion risk measures.
The Kumaraswamy risk measure for the uniform loss is given by, with substitution When θ = 1, the Kumaraswamy risk measure is the PH or power risk measure of 2b/(1 + α) for 0 < α ≤ 1; see also Wang (1995).
The UEE risk measure, with substitution s = For the UG distortion, with substitution s = (x/2b) −α , the UG risk measure is given by 2b where E(a, z) = ∞ 1 e −zs /s a ds is the generalized exponential integral.The UGQ distortion risk measure is, with substitution The distortion risk measures derived for a uniform loss on [0, 2b] in this section are summarized in Table 2 below.The distortion risk measures by power and dual-power distortions in Wang (1995) are special cases of the Kumaraswany risk measures.They are presented not only as a strategy for double-checking the calculations, but also due to their popularity.

Distortion
Parameter Space Risk Measures

Exponential Loss
When the loss variable X has an exponential distribution with mean b, the respective cdf, survival function, and their inverses are Below, the formulas in (2) will be used and the derivative g of each distortion in Table 1 will be required.
The Kumaraswamy risk measure for an exponential loss can also be found in Samanthi and Sepanski (2018).Briefly, from (2), it is given by, with substitution where the integral is the expected value of the logarithm of a beta variable with parameters 1 and θ and the digamma function ; see https://en.wikipedia.org/wiki/Digamma_function, (accessed on 6 October 2021).
The UG risk measure is given by, with substitution where E is the generalized exponential integral.It is shown in Section 3.4 that this UG exponential risk measure as a special case of UG Weibull risk measures being finite.
The UGQ risk measure using the formula ρ = g(S(x))dx is When α ≥ 1, the risk measure is not finite.When 0 Table 3 summarizes the distortion risk measures for an exponential loss with mean b.It is obvious that a larger mean of b results in a larger distortion risk measure.As the distortion parameter θ increases, the resulting distortion risk measure increases.As indicated in Figure 1, a larger θ represents a higher level of risk aversion.The effects of the distortion parameter α do not have the same pattern across the distortions considered in Table 3.For example, for 0 ≤ α < 1, as α increases, the power and Kumaraswamy risk measures for exponential losses decrease while the UGQ exponential risk measure increases.

Distortion
Parameter Space Risks Measure For an exponential loss with a mean of 50, the 3D Figure 2 displays the UEE and UG distorted risk measures involving the sum of an infinite series, numerically computed using R, with the parameter values α and θ in their respective admissible spaces.The plots for the UEE and UG risk measures allow us to gain better insight into the effect of the parameters on the risk measures.A larger θ value appears to yield a larger UEE or UG exponential risk measure, which is not at all transparent by examining the formulas.The effect of increasing α seems less dramatic when alaha is large for the UG distortion.

Lomax Loss
The cdf, survival function, and their inverses, with parameters a > 0 and b > 0, of a Lomax loss variable X are given by, for x > 0 and 0 < u ≤ 1, respectively.The Lomax loss has a mean of b/(a − 1).The calculations of the Kumaraswamy and UEE distortion risk measures for a Lomax loss are shown in Samanthi and Sepanski (2018).By (2), (11), and with substitution s = u −α , the UG distortion risk measure for a Lomax loss is where E is the generalized exponential integral, since (1 The UGQ distortion risk measure for a Lomax loss is given by Note that, by L'Hopital's rule, for 0 Therefore, the risk measure in ( 12) is not well defined.
The mean of a Lomax loss increases as a decreases and b increases.Intuitively, one expects the parameters a and b to have the same effect on the distortion risk measures, which can be verified from the formulas except for the UEE distortion in Table 4. Varying the α and θ values, we plot the UEE and UG risk measures for a Lomax loss with a mean of 50 with (a, b) = (12.61,580.40) in Figure 3.The figure shows that the UEE Lomax risk measure increases as θ increases and α decreases and that the UG Lomax risk measure increases as θ or α increases.When the risk aversion indicator θ increases, all distorted risk measures increase.

Distortion
Parameter Space Risk Measures

Weibull Loss
When the loss variable X follows a Weibull distribution, its respective cdf and survival function are The hazard rate function of the Weibull distribution with parameters b and c is (c/b)x (c−1) .It is a decreasing function when c < 1, marking a heavy-tailed distribution, and increasing when c > 1, marking a light-tailed distribution.The characteristics of its hazard function make the Weibull distribution an adequate model for a variety of applications, such as weather forecasting, insurance modeling, and financial risk analysis; see Frees (2018).
For the Weibull loss, the Kumaraswamy risk measure is defined to be since Bernoulli's inequality states that (1 + t) r ≥ 1 + rt for t ≥ −1 and r ≥ 1, and (1 + t) r ≤ 1 + rt for t ≥ −1 and 0 ≤ r ≤ 1.With a substitution of x c , the integral ∞ 0 e −αx c /b dx = Γ(1/c)(b/α) 1/c /c for c > 0 and is finite.Thus, the Kumaraswamy distortion risk measure for a Weibull loss is finite.By the binomial expansion of (1 When α = 1, we obtain the dual-power distortion risk measure given by When c = 1, the Weibull loss reduces to the exponential loss.Using the Newton series for the digamma function, the risk measure in ( 13) is equal to which coincides with the result for the exponential loss in Table 3 ) θ for the admissible parameter space of 0 < α ≤ 1, θ ≥ 1, and by applying Bernoulli's inequality, the UEE risk measure for the Weibull loss, we obtain that The UG Weibull risk measure has no closed form and needs to be computed numerically.Below, we show that the UG Weibull risk measure is finite.It is defined to be 14) since e t ≥ 1 + t or 1 − e t ≤ −t by the Taylor series.The integrand in A 1 is bounded on the unit interval; therefore, A 1 in ( 14) is finite.Let z = 1 − e −1/b .With the substitution t = (1 − e −x c /b ) ∈ (1 − e −1/b , 1) and x = [−b ln(1 − t)] 1/c for 1 < x < ∞ and using integration by parts, A 2 is given by The limit term in ( 15) can be shown to be 0 by repeatedly applying L'Hopital's rule and limit laws.By the Cauchy-Schwarz inequality for integrals, we obtain that Combining ( 14), ( 15), ( 16), and ( 17), we conclude that the UG Weibull risk measure is finite.
From Table 5, the power distortion churns out a risk measure that is θ 1/c times the mean loss.For a heavy-tailed Weibull loss with c < 1, the distorted risk measure would increase by a higher magnitude than the case with c > 1.By Table 5 and Figure 4, the distortion risk measure increases as θ increases or α decreases for the UEE and Kumaraswany distortions.

Distortion
Parameter Space Risk Measures

Numerical Analyses and Estimation
In this section, we numerically compute distortion risk measures for uniform, exponential, Lomax, and Weibull losses.Note that the parameter θ in Table 1 assumes values of at least 1.We then propose a non-parametric, plug-in estimator for the distortion risk measure and conducted simulations to compare the proposed estimator with the empirical L-estimator in (4).

Numerical Results for Distortion Risk Measures
The loss distributions considered are uniform on the interval (0, 100), exponential with mean 50, Lomax with parameters (12.61, 580.40), and Weibull with parameters (5, c = 0.50) and (412.20,c = 1.50).The loss distribution parameters were selected so that all loss distributions had a mean of 50.We fixed the mean value of the loss distributions and inspected how distortion risk measures vary with the tailedness, measured by the kurtosis, of a loss distribution.A normal distribution has a kurtosis of 3 and an exponential distribution of 6.The Lomax loss parameters were chosen to have a kurtosis of 12.For the Weibull loss, we used c = 0.5 and c = 1.5, corresponding to heavy-tailed and light-tailed distributions with kurtoses of 84.72 and 1.39, respectively.Note that the UGQ distorted Lomax risk measures are not well-defined; see Table 4.The VaR and CTE values at various percentile levels for each loss distribution are reported as benchmark comparisons and can provide insights into the right tails of the loss random variables under study.They are supplied in Table 6.The risk premium, defined as the difference between the premium, e.g., the CTE, and the expected value, is larger for a loss with a heavier tailedness.Table 6.VaR and CTE values at levels of 0.25, 0.5, 0.75, 0.95, and 0.99 for uniform (0, 100), exponential (0.02), Lomax (12.61, 580.40),Weibull (5, 0.5), and Weibull (412.20,1.5) losses.Table 7 tabulates the beta, Kumaraswamy, and UEE risk measures for the losses examined in Section 3 and Table 8 the UG and UGQ risk measures for some admissible parameter values.The dash notation in Table 8 connotes that the resulting risk measure is not finite for UGQ Lomax risk or not well-defined due to the lack of admissibility.The beta and Kumaraswamy distortions both have the power and dual-power distortions as their special cases when θ = 1 or α = 1; therefore, they are expected to have the same risk measures in theses cases.While this conclusion seems trivial, we computed them to possibly detect programming errors.

Level
Table 7.The beta, Kumaraswamy, UEE risk measures for uniform (0, 100), exponential (0.02), Lomax (12.61, 580.40),Weibull (5, 0.5), and Weibull (412.20,1.5) losses.The distortion risk measures in Tables 7 and 8 appear to increase as θ increases, as shown in the previous section.As α increases, the resulting risk measures appear to decrease for the beta, Kumaraswamy, and UEE distortions and increase for the UG and UGQ distortions if well-defined.The loss distributions with a larger kurtosis or tailedness are more sensitive to the distortion in the sense that the distortion will yield a larger risk measure.The UGQ distortion measures seem to be more sensitive than the other distortions to the loss kurtosis.As we can see for the exponential and Weibull (c = 0.5) losses, they can be much larger than the VaR and CTE.Both the VaR and CTE values of the four loss distributions fall in between the ranges of the distortion measures produced by the beta, Kumaraswamy, UEE, and UG distortions presented, which may be seen as an indication that those distortions have practicability.

Estimation of Distortion Risk Measures
Instead of (2), we considered the equation in (1) such that ρ = ∞ 0 xg S(x) f (x)dx = E Xg S(X) , the expected value of the function Xg S(X) .An intuitive estimator of an expected value is the sample mean.Let F(x) = ∑ n j=1 I(X j ≤ x)/n denote the empirical cdf.Then, the survival function S(X) can be estimated by The proposed plug-in estimator is given by No parametric assumptions were made about the underlying loss distribution.
Below, we employed simulations for the time being to preliminarily investigate the behaviors of the proposed estimator and compare it with the empirical L-estimator in (4) for the exponential loss and the Lomax losses with parameters of 5 and 200, both with a mean of 50.The exponential distribution had a kurtosis of 6 and the Lomax distribution a kurtosis of 70.8.Five-hundred simulations were run for each of sample sizes n = 50, 100, 500, and 1000 for various distortions.For each distortion, we ran them for three pairs of admissible α and θ values.The true theoretical distortion risk measures, the sample mean, the bias defined as the difference between the sample mean and the theoretical risk measure, and the standard deviations of estimates from the 500 simulations are displayed in Tables 9 and 10.
As one would expect, the bias and standard deviation decreased as the sample size increased in Tables 9 and 10.As seen in the previous section, as θ increased, the resulting distortion risk measure increased, so were the biases and standard deviations of the Lestimator and plug-in estimators.A distortion of a loss distribution with a larger kurtosis would result in a larger risk measure and also a larger bias and standard deviation in the two estimators.When the sample size was 500 and 1000, in general, the proposed plug-in estimates appeared to perform slightly better than the L-estimates in terms of the bias and standard deviation.Based on the results from the sample size of 1000, both the L-estimator and plug-in estimators seemed to be consistent for the distortions with parameter values of α = 1 and θ = 3.Though the proposed plug-in estimator had smaller biases in the majority of the cases, both estimators performed poorly in estimating the UGQ risk measures for the exponential loss.

Figure 1 .
Figure 1.UG distortion cures are displayed in (a) and (b), and UGQ in (c) and (d) for varying θ or α.

Figure 2 .
Figure 2. The 3D graphs of UEE and UG risk measures for an exponential loss with mean 50.

Figure 3 .
Figure 3.The 3D graphs of UEE and UG risk measures for a Lomax loss with mean 50.

Table 4 .
Distortion risk measures for a Lomax loss with parameters (a, b).

Figure 4 .
Figure 4.The 3D graphs of the Kumaraswamy, UEE, and UG risk measures of Weibull loss with mean 50.

Table 3 .
Distortion risk measures for an exponential loss with mean b.

Table 5 .
Distortion risk measures for a Weibull loss with parameters b and c.

Table 9 .
Sample means, biases and standard deviations in parentheses of L-estimates and plug-in estimates from 500 simulations for exponential loss with mean 50.

Table 10 .
Sample means, biases and standard deviations parentheses of L-estimates and plug-in estimates from 500 simulations for Lomax loss with parameters (5, 200).