New families of bivariate copulas via unit weibull distortion

*Correspondence: aldhu1fa@cmich.edu Department of Statistics, Actuarial and Data Sciences, Central Michigan University, Mount Pleasant, MI 48859 USA Abstract This paper introduces a new family of bivariate copulas constructed using a unit Weibull distortion. Existing copulas play the role of the base or initial copulas that are transformed or distorted into a new family of copulas with additional parameters, allowing more flexibility and better fit to data. We present a general form for the new bivariate copula function and its conditional and density distributions. The tail behaviors are investigated and indicate the unit Weibull distortion may result in new copulas with upper tail dependence when the base copula has no upper tail dependence. The concordance ordering and Kendall’s tau are derived for the cases when the base copulas are Archimedean, such as the Clayton and Frank copulas. The Loss-ALEA data are analyzed to evaluate the performance of the proposed new families of copulas.


Introduction
Copulas serve numerous fields including insurance and finance. For example, (Frees and Valdez 1998) demonstrated their usefulness and explored practical applications such as estimation of joint life mortality and multi-decrement models. Nazemi and Elshorbagy (2012) implemented copula modeling to study the interdependence among hydrological data. The fitness of statistical models rests on its flexibility and more parameters may better accommodate various features in data. Construction of new families of copulas with better fitness has been of interest to researchers. In this paper, we provide a new distortion mechanism of copula construction and start by setting forth fundamentals and relevant literature below.
Let X and Y be continuous random variables with a joint distribution function H(x, y) = P(X ≤ x, Y ≤ y) and marginal cumulative distribution functions (cdf ) F(x) = P(X ≤ x) and G(y) = P(Y ≤ y). Sklar (1959) showed that there exists a unique copula C such that H(x, y) = C(F(x), G(y)) with joint probability density function (pdf ), h(x, y), given by h(x, y) = c (F(x), G(y)) f (x)g(y), where the copula pdf c (u, v) = ∂C (u, v)/∂u∂v, f (x) = dF(x)/dx = F (x), and g(y) = dG(y)/dy = G (y). Note the prime mark will be used to denote a derivative throughout the paper.
A bivariate copula can arise form a bivariate joint cdf. For example, the Gaussian copula is derived from the bivariate Gaussian distribution. Conversely, it can also be used to generate new bivariate probability distributions via (1); see Nelsen (2006) for summaries of methods of constructing copulas. Methods for constructing new bivariate joint distributions (Balakrishnan and Lai 2009) have also been adopted to create new copulas. For example, the framework of the Sarmanov-Lee distribution (Lee 1996) was utilized by Sharifonnasabi et al. (2018) and Cooray (2019) to construct new copulas. It is related to the bivariate FGM distribution introduced by Morgenstern (1956), given by h(x, y) = f (x)g(y) 1 + αh 1 (x)h 2 (y) , x, y ∈ R, where h 1 and h 2 are two functions satisfying certain conditions.
Distortion or transformation of existing copulas is another framework for forging new families of copulas. Valdez and Xiao (2011)proposed three kinds of distortion approaches: (1) distortion of the margins alone without altering the original copula structure; (2) simultaneous distortion of the margins and the copula structure; and (3) synchronized distortion of the copula and its margins. In this paper, we focus on the distortion of the third kind that acts on the copula and induces the copula defined in (2). A function T :[ 0, 1] →[ 0, 1] is said to be a distortion function if it is continuous and non-decreasing, not necessarily convex or concave, with T(0) = 0 and T(1) = 1. It is said to be admissible for a base or an initial copula C if the transformed copula C T (u, v) of the form is also a copula. Note that, as in Valdez and Xiao (2011), T is assumed to be strictly increasing such that T −1 exists and is continuous on [ 0, 1] . If the initial copula is Archimedean with generator φ, then C T is Archimedean with generator φ • T −1 ; see Di Bernardino and Rulliere (2013) and the right composition rule in Genest et al. (1998). A convex T is required for the admissibility; see Morillas (2005) or Theorem 3.3.3 in Nelsen (2006). Durante et al. (2010) showed T is admissible if T • exp : (−∞, 0) →[ 0, 1] is log-convex and suggested several distortion functions. The logconvex condition will be used to obtain the admissible parameter space for the proposed distortion. Samanthi and Sepanski (2019) constructed a new family of copulas via beta cdf distortion. The mixture of Max-infinitely divisible approach for constructing BB1-BB7 copulas in Joe (2015) is also a distortion method. Xie et al. (2019) presented a family of bivariate copulas by transforming an initial/base copula using two increasing functions. For more references, see Xie et al. (2019).
In this paper, we inaugurate a distributional distortion derived by a transformation of a Weibull random variable. This paper is organized as follows. In "Groundwork", we first lay some groundwork required for the derivation of properties of the new family of copulas. "The proposed unit weibull distortion" stages the unit Weibull (UW) distortion function and the admissibility conditions on the parameters. In "Unit-Weibull distorted copulas", the UW distorted copula distribution and its corresponding conditional and density distributions are formulated. Examples and possible limiting cases in parameters are presented. "Properties" investigates properties such as the tail dependence coefficients, tail orders, and concordance ordering. To assess its performance, the new UW-distorted copula model is applied to the Loss-ALAE data set in "Application" sections, followed by concluding remarks.

Groundwork
In this section, we describe notation, definitions and some known results that would be applied to distorted copulas; see Joe (2015) for more details. From (1), a copula contains the dependence structure between two random variables and links a bivariate distribution function to its univariate marginal cdf 's. It has the following properties: is a continuous, strictly decreasing and convex function such that φ(1) = 0 and φ(0) = ∞.
Tail dependence coefficients are measures of extremal dependence that quantify the dependence in the lower-left-quadrant tail or upper-right-quadrant tail of a bivariate distribution. Let U and V be two unit uniform random variables with a joint copula cdf The lower tail dependence coefficient, λ L , is defined as the limit value of the conditional probability of U ≤ u given V ≤ u as u → 0 + and can be calculated as lim u→0 + C(u, u)/u. The upper tail dependence coefficient, λ U , is defined as the limit value of the conditional probability of U > u given V > u as u → 1 − . It can be simplified as lim Let T be an admissible distortion function, then the induced copula is of the form displayed in (2). Since T is a distortion function and by L'Hopital's rule, the lower tail dependence coefficient for a T-distortion induced copula is given by where t(u) = dT(u)/du, if the lower tail dependence coefficient λ L of the initial copula C and the limit of t (C(u, u))/t(u) at 0 + exist. Since lim u→1 − T(u) = 1, with the substitution of v = T −1 (u) and by L'Hopital's rule, the upper tail dependence coefficient is given by if the upper tail dependence coefficient λ U of the initial copula C and the limit of t (C(u, u))/t(u) at 1 − exit.
Let f 1 and f 2 be two functions. If lim u→u 0 f 1 (u)/f 2 (u) = 1, we denote it by f 1 (u) ∼ f 2 (u) as u → u 0 . A positive function f 1 defined on (0, ∞) is regularly varying at 0 with index γ , in which case we write f 1 ∈ R(γ ), if for some real number γ it satisfies lim x→0 f 1 (sx)/f 1 (x) = s γ for all s > 0.
A function f 1 is said to be slowly varying if γ = 0. Karamata's Characterization Theorem (Bingham et al. 1989) states that every regularly varying function f 1 with index γ is of the form f 1 (x) = x γ (x), where is a slowly varying function. Buldygin et al. (2006) derived that if f 1 (x) is regularly varying at 0 (or ∞) with an order of γ ∈ R, then For a bivariate copula C, if C(u, u) ∼ u κ L (u) as u → 0 + , where (u) is slowly varying at 0 + , then κ L is referred to as the lower tail order of the copula C. Let C(u, v) be the survival copula. The upper tail order is defined as κ U if C(u, u) ∼ u κ U * (u) as u → 0 + for some slowly varying function * (u). When κ L = 2 and (u) → q as u → 0 + , for some positive q, the variables are near independent in the lower tail. If 1 < κ L < 2, the variables are positively associated and have intermediate tail dependence. The case κ L = 1 corresponds to the usual tail dependence coefficient λ L ∈ (0, 1) with lim u→0 + C(u, u)/u = lim u→0 + (u). Similar conclusions are made for the upper tail and κ U ; see Hua and Joe (2013) for more details.

The proposed unit weibull distortion
By the definition of a distortion function, a continuous cdf with domain I is a distortion function. In this section, we define the unit Weibull cdf and examine its admissibility.
Let W be a non-negative continuous random variable with cdf G(.) and pdf g(·). Define Z = exp(−W ). Then, the cdf of Z is given by which is the survival function of W evaluated at − log z. It is related to the expres-sionḠ(− log z), whereḠ is a survival function, suggested by Durante et al. (2010). If W is a Weibull random variable, we name Z as a unit Weibull (UW) random variable with a support of the unit interval I. Let G be the Weibull cdf given by G(w) = 1 − exp (−bw a ) , a, b > 0, w ≥ 0, then the cdf of Z is given by The UW quantile function and UW pdf are given by, respectively, To find the admissibility of the distributional distortion function T, we employ the following proposition shown in Durante et al. (2010).
is log-convex, then the function C T in (2) is a copula.
The following corollary specifies constraints on the parameter values in the UW distortion to ensure the admissibility.
Corollary 1 Let T(u) and T −1 (u) for u ∈[ 0, 1] be the UW-distortion and quantile functions in (7) and (8), respectively. Then, the function C T in (2) is a copula if 0 < a ≤ 1 and b > 0. respectively. Since the second derivative l is non-negative if 0 < a ≤ 1 and b > 0, the conclusion follows from Proposition 1.

Unit-Weibull distorted copulas
We below present the copula distribution, pdf, and conditional pdf, and derive limiting cases for the proposed new family of UW distorted copulas. Applying (7) and (8), the copula C T in (2) is of the following general form where 0 < a ≤ 1 and b > 0. When a = 1 and b = 1, then T(u) = u, i.e., the initial copula is not distorted. The initial copula is a member of the proposed family of copulas. When a = 1, the UW distortion is the power distortion.
If the initial copula C is Archimedean with a strict generator function of φ, then C T is Archimedean with generator given by Example 1 UW-Clayton copula. Consider the Clayton copula expressed as copula has the following expression and its generator is given by copula derived by Joe (2015) is a special case when b = 1.
Example 2 UW-Gumbel copula. Consider the Gumbel copula expressed as C (u, v; The UW-Gumbel copula has the following expression Note the parameters θ and a cannot be identified separately. Reparameterizing by setting θ/a = δ, we see that the UW distortion of the Gumbel copula returns the Gumbel copula and does not yield a new family of copulas. This example prompts us to consider the UW distortion of extreme-value bivariate copulas such that C u 1/m , v 1/m m = C(u, v), for integers m ≥ 1, and are of the form where A(·) is convex and satisfies certain constraints; see Gudendorf and Segers (2010).

Aldhufairi and Sepanski Journal of Statistical Distributions and Applications
(2020) 7:8 Page 6 of 20 The parameter b originated from the UW distortion disappears.
Example 3 UW-independence copula. Consider the independence copula expressed as C(u, v) = uv. The UW-independence copula has the following expression Therefore distorting the independence copula results in the Gumbel copula. That is, the proposed UW distortion gives another genesis of the Gumbel copula.
The UW-Frank copula has the following expression

Conditional distribution and copula density
The conditional density C(u|v) plays a key role in simulating bivariate data linked by a copula C since the conditional distribution P(X ≤ x|Y = y) = ∂C (F(x), G(y)) /∂v and C(u|v) = ∂C(u, v)/∂v. A general algorithm to generate draws from a bivariate copula C using the conditional distribution approach (see Nelsen (2006)) is described as follows. (i) Generate two independent uniform random values (u 1 , v) and (ii) solve C (u 2 |u 1 ) − v = 0 for u 2 . The desired pair is (u 1 , u 2 ) . We will also obtain the copula density function needed for computing the maximum likelihood estimates of parameters. (10) can be rewritten as The conditional cdf and copula pdf can be respectively derived by The derivatives of H with respect to x and x with respect to u are The conditional cdf is therefore given by u .

Aldhufairi and Sepanski Journal of Statistical Distributions and Applications
(2020) 7:8 Page 7 of 20 The derivation of the copula density c T (u, v) in (12) requires the calculation of ∂ 2 H/∂y∂x, which is lengthy and tedious therefore not presented here. Next, we consider the case when the initial copula is Archimedean with generator φ. In this case, let x = φ T −1 (u) and y = φ T −1 (v) . Note that x and y are decreasing and map [ 0, 1] to [ 0, ∞] such that φ(1) = 0 and φ(0) = ∞. Then, whereH is a bivariate survival function with univariate margins in (9). The conditional cdf and pdf of the UW distorted copula can be obtained from (12) and the following: where Example 5 UW-Clayton Copula. The Clayton copula, see Example 1, is Archimedean in the following components into (13) and (14) to obtain the conditional distribution and density of the UW-Clayton copula: (13) and (14) to derive the conditional distribution and density of the UW-Frank copula, the following expressions will be required: φ φ −1 (x + y) = θ 2 e x+y e −θ − 1 −1 1 + e x+y e −θ − 1 −1 .

Limiting cases
When a = 1, the UW distortion function becomes T(u) = u b , the power distortion, and the UW distortion results in copulas of the form C T (u, v) = C u 1/b , v 1/b b . Proposition 2 below is not applicable to the case when the initial copula is an extreme-value copula, for the power distortion doesn't produce a new family of copulas in this case; see Example 2.
Proposition 2 Let C T be the unit-Weibull distorted copula in (10), where 0 < a ≤ 1 and b > 0. Then, C T approaches the independence copula when b → ∞ and a → 1.
Proof Let r = 1/b, x = e −(−r log u) 1/a , y = e −(−r log v) 1/a , and A r = C(x, y). The derivative A r = dA r /dr is given by Hopital's Rule and chain rule, the limit of the exponent term in (10) have that x → 1, y → 1, and A r → 1. When a → 1, lim A r = log(uv) since C 2|1 and C 1|2 are conditional distributions. Therefore, lim b→∞ C T (u, v) = exp(log(uv)) = uv.
Proposition 2 provides the limit of the UW copulas when b → ∞ and a → 1 without specifying the initial copula. In the following, we find the limiting copulas in the parameter θ originated from the initial copula for families of UW-Clayton and UW-Frank copulas.
Example 7 Consider UW-Clayton copula in Example 1. By the same arguments for the limit of the Clayton copula in Joe (2015), Therefore, the UW-Clayton copula of the form T C T −1 (u), T −1 (v) , by Example 3, The limit of UW-Clayton copulas as θ → ∞ is the Gumbel copula. When b = 1, the UW-Clayton begets the BB3 copula. Therefore, UW-Clayton copulas approach the comonotonicity copula when θ → ∞ or a → 0 + , and the Gumbel family when θ → 0 + .
As θ → 0 + , by (16) and the facts that e −θ ∼ 1 − θ, (1 + s/n) n ∼ e s as n → ∞, using the Archimedean representation of the UW-Frank copula, the limit of the UW-Frank copula is which is the Gumbel copula.

Properties
We obtain the tail dependence coefficients and tail orders for the UW-distorted copulas, and study the tail concordance in the parameters for the UW-Clayton copulas. (7).

Consider the distortion T(u)
When b = 1, T turns into the form suggested by Durante et al. (2010), and when b = 1 and a = 1, T is the identity distortion, i.e., no distortion is applied. Definitions of tail orders can be found in "Groundwork" section. We note here the joint survival function (u, v) and Talyor's series approximation of (1 + u) a ∼ 1 + au, log(1 − u) ∼ −u, e u ∼ 1 + u, as u → 0.
Therefore, we have that Below we assume that the lower tail coefficient λ L = 0 when κ L > 1 and the upper tail coefficient λ U = 0 when κ U > 1 for the initial copula. Let the subscript T notation denote the properties of the UW distorted copulas, e.g., λ T,U is the lower tail coefficient for the UW distorted copulas.

by L'Hopital's rule and lim
We obtain below tail orders of C T . By (5) and (18) and since T −1 (u) = exp − −b −1 log u 1/a and lim u→0 + log (u)/ log u = 0 (Bingham et al. 1989), Note that for s > 0, sT(v))/v, and lim v→0 + q(v) = 1. One can then show the exponential term in (22) is slowly varying by definition. Note that if a = 1 and b = 1, then (22) returns the assumption that C(u, u) ∼ u κ L (u) as u → 0 + . Using the approximations in (18) and (19), Therefore, for the upper tail order, by (18), indicating that UW distorted copulas have an upper tail order of 1 except when a = b = 1.
The following table summarizes tail orders and dependence coefficients for the family of UW-distorted copulas with the initial copulas being BB1, Clayton, Frank, and Gaussian copulas, where θ, δ, and ρ are the parameters in the initial copulas with formulas shown in Joe (2015).
Density contour plots with standard normal margins for various combinations of (a, b) are shown in Fig. 1. The parameter θ is chosen so that the initial copula has Kendall's tau of 2/7 or -2/7. As indicated by Proposition 3 or Table 1, the family of UW-distorted copulas not only preserves the tail dependence of the initial copula but also can accommodate upper tail dependence. Unlike the initial Frank copula with a = b = 1, the resulting UW-Frank copulas are asymmetric. The graphs also reflect the results in Proposition 3, the upper tail dependence becomes stronger as a decreases.

Concordance ordering
A copula is said to be positively ordered (u, v; β) whenever the parameters α ≤ β for all u, v ∈ I; see Nelsen (2006) for more details.
Proposition 4 If the initial copula is positively or negatively ordered by its parameter, then the unit-Weibull distortion preserves the concordance order in the parameter of the initial copula.
Proof If the initial copula C is positively ordered, then, for θ 1 ≤ θ 2 , That is, the family of T-distortion induced copulas is also positively (negatively) ordered by the parameter θ originated from the initial copula if the initial copula is positively (negatively) ordered by the parameter.
To examine the concordance ordering in the parameters a and b introduced by the UW distortion for UW-Clayton copula, we present the following corollary and lemma; see Schweizer and Sklar (1983) or Nelsen (2006).
Corollary 2 Let C 1 and C 2 be Archimedean copulas with generators φ 1 and φ 2 , respectively. Then C 1 ≺ C 2 holds if one of the following conditions is satisfied (i) φ 1 • φ −1 2 Table 1 Examples of tail orders and dependence coefficients

Example 9
Consider the family of the UW-Clayton copulas in Example 1. Below we show the family of UW-Clayton copulas is negatively ordered by the parameter b but not a. It is Archimedean, see Example 5, with generator and inverse generator given by respectively. We wish to use Corollary 2 to show the concordance order in the parameter and b. Define h b (u) to be The first derivative and second derivatives of h b are give by That is, the family of the UW-Clayton copulas is negatively ordered by the parameter b. Define h a (u) to be The first derivative and second derivatives of h a are give by h a (u) = a 2 θ 1−a 2 /a 1 a 1 e θ 1−a 2 /a 1 [log(1+θu)] a 2 /a 1 log(1 + θu) The UW-Clayton is negatively ordered by the parameter a if (23) is nonpositive for a 2 < a 1 and all u ∈ I. As we will see from Fig. 2, the UW-Clayton is not ordered by the parameter a for all θ and b values.

Measures of concordance
In this section, we explore two widely known scale-invariant measures of concordance or association: Spearman's rho and Kendall's tau. If X and Y are continuous random variables with copula C, then the Spearman's rho and Kendall's tau can be expressed as, respectively, For an Archimedean copula (Genest and MacKay 1986), Kendall's tau is also given by  (u, v)/∂u and C 1|2 (u, v) = ∂C(u, v)/∂v. Numerical integration methods can be employed to compute the concordance measures. If the inital copula is Archimedean with generator φ(u), from (11) and (24), (u) = φ T −1 (u) and (u) = φ T −1 (u) /t T −1 (u) . The Kendall's tau for a UW distorted copula is given by Proposition 5 Let X and Y be random variables with copula C T , a UW distorted copula of the form in (10), where the initial copula is Archimedean with generator φ(·). Then, the Kendall's tau between X and Y can be expressed as Proof: From (9) and (25) Example 10 Kendall's tau of the UW-Clayton copula. When the initial copula is the Clayton copula, Proposition 5 gives With (26), one can readily write R programs to compute Kendall's tau values at various parameter values. Figure 2 plots values of Kendall's tau for various parameter values. The plot for b = 8 indicates when θ is large, e.g., θ = 30, the family of UW-Clayton copulas is not ordered by the parameter a as the resulting tau values are not monotone in a; see also Example 9.
Example 11 Kendall's tau of the UW-Frank copula. From Example 6 and Proposition 5, with substitution v = θe −u , the Kendall's tau of a UW-Frank copula is given by for θ = 0, 0 < a ≤ 1, and b > 0. We compute Kendall's tau coefficients at various parameter values using the formula in (27) and produce Fig. 3. While not mathematically shown due to its tediousness, the plot for θ = 10 illustrates that the family of UW-Frank copulas is not ordered by the parameter a as the resulting tau values are not monotone in a.

Application
In this section, the Loss-ALAE insurance data set is analyzed to evaluate the performance of the proposed UW distorted copula models. It is easily accessible from the R copula package. The loss variable is the general liability claims and the allocated loss adjustment expenses (ALAE) is attributable to the settlement of individual claims (e.g., lawyer's fees, claims investigation expenses). The summary statistics, including the standard deviation (SD), the first quartile (Q 1 ) and the third quartile (Q 3 ) are reported in Table 2.
To visualize the relationship, scatter plots in Fig. 4 are constructed on the real dollar scale and on the log scale. There seems to be an upper tail dependence in the data.
Using the notation in (1), the log-likelihood function for the data (x i , y i ) n i=1 is given by where α 1 and α 2 are the parameters in the marginal distributions, and θ, a and b are parameters in the copula function. Rather than a full maximum likelihood estimation, one of the more attractive estimation methods is the two-stage maximum likelihood estimation, also know as inference function for margins (IFM); see Joe (1997). The IFM first obtains estimates, α 1 and α 2 , of parameters in the marginals by maximizing n i=1 log f (x i ; α 1 ) and n i=1 log g (y i ; α 2 ) , and then computes estimates of the parameters θ, a, and b by maximizing To determine the appropriate marginals for the ALAE and loss variables, (Frees and Valdez 1998) and (Frees 2018) overlaid the fitted Pareto cdf and the empirical cdf and found the two curves are reasonably close to each other for both variables. Since we ignore the mild censoring in the loss variable, we reassess the fit using P-P plots. The P-P plots in Fig. 5 indicate that the Pareto margins fit the data well. Therefore, we model both of the marginals using the Pareto distribution with the cdf of the from 1 − [1 + (x/α i1 )] −α i2 , i = 1, 2, where α i1 is the scale parameter and α i2 is the shape parameter. From the loss marginal, α 1 = ( α 11 , α 12 ) = (16228.15, 1.238), and for ALAE, α 2 = ( α 21 , α 22 ) = (15133.6, 2.223). Plugging in α 1 and α 2 in (28), we then ascertain estimates in Table 3 by maximizing (28). Genest et al. (1998) ignored the censoring in the loss data, so did we in this paper. The pseudo log-likelihood estimation maximizes (28) using nonparametric, empirical distribution estimates of F and G. The results are of little differences and therefore not reported here.
Based on the scatter plots, we select survival Clayton (SClayton) (a 180 degrees rotation of the Clatyon copula) and Gumbel copulas, in addition to Frank and Gaussian copulas. The two-parameter survival BB1 (SBB1) is selected by the SelectCopula() in R Vine package as the best-fit bivariate copula model. The R optim() function is employed to compute the IFM estimates.
In addition to standard errors (SE) of the estimators, we also report the resulting loglikelihoods (IFML) in (28) and AIC for comparison purposes. The parameter θ (and δ for the survival BB1 copula) is the one originated from the initial copula, a and b are parameters injected by the UW distortion. The estimate of the correlation parameter ρ in th Gaussian copula is denoted by θ in Table 3. The estimate of θ resulting from fitting the initial copula is used as a starting point in the optim() for fitting UW-distorted models and the Beta(B)-distorted copulas (Samanthi and Sepanski 2019). The dash symbol (-) is used for models with a lower or an upper tail dependence coefficient of 0. Note that the sample Kendall's tau is 0.315. The Cramer-von Mises goodness-of-fit test statistic (CvM) for copulas in Genest et al. (2009) measures the sum of square deviations between the empirical cdf and an estimated copula cdf. Larger CvM values are less desirable. The bootstrap approach detailed in Genest et al. (2009) was used to calculate p-values. Note the sample size is 1500 and we used 1000 bootstrap replications. The results, e.g., the CvM for the UW-SBB1 copula is 0.296 with a p-value of 0.127, are tabulated in Table 4. Apart from the Frank, Gaussian, and SClayton copulas, the copulas listed in Table 4 provide an adequate fit to the data.  The distortion induced copulas, as expected, outperform the initial copulas in terms of log-likelihoods. There is a sizable improvement in the log-likelihoods of the UW-Frank, UW-Gaussian, Beta-Frank, and Beta-Gaussian copulas over the Frank and Gaussian copulas. It may due to the fact that the distorted copulas can accommodate the upper tail dependence in the data. As indicated in Table 4, distortions can improve the goodnessof-fit in terms of the CvM statistic. While more parameters are expected to yield better log-likelihood results, the AIC that penalized for having more parameters indicates that the two-parameter survival BB1 is the winner among the chosen copulas for fitting this particular data set. However, the estimated Kendall's tau calculated based on the estimated survival BB1 model seems to deviate from the sample Kendall's tau more than the ones based on the estimated UW-SBB1 and B-SBB1 models. The upper tail coefficient estimates from models with better performance in terms of IFML seem to suggest the upper tail dependence exists in the data. The standard errors of estimates from the UW-distorted copulas are smaller than those from Beta distorted copulas, although they perform comparably in terms of the AIC and log-likelihood.

Concluding remarks
This paper constructs a new family of copulas by employing the unit Weibull distributional distortion function. With an additional two parameters in the unit Weibull distortion, the new family of copulas allows for more modeling flexibility and versatility. Note also that the initial copula is a special case of the UW-distorted copulas and