scieee AI-readable full text Open interactive document viewer

Sharp convex bounds on the aggregate sums: An alternative proof

Yin, Chuancun,Zhu, Dan

Abstract

EconStor is a publication server for scholarly economic literature, provided as a non-commercial public service by the ZBW.

Full text

Yin, Chuancun; Zhu, Dan Article Sharp convex bounds on the aggregate sums: An alternative proof Risks Provided in Cooperation with: MDPI – Multidisciplinary Digital Publishing Institute, Basel Suggested Citation: Yin, Chuancun; Zhu, Dan (2016) : Sharp convex bounds on the aggregate sums: An alternative proof, Risks, ISSN 2227-9091, MDPI, Basel, Vol. 4, Iss. 4, pp. 1-8, https://doi.org/10.3390/risks4040034 This Version is available at: https://hdl.handle.net/10419/167899 Standard-Nutzungsbedingungen: Die Dokumente auf EconStor dürfen zu eigenen wissenschaftlichen Zwecken und zum Privatgebrauch gespeichert und kopiert werden. Sie dürfen die Dokumente nicht für öffentliche oder kommerzielle Zwecke vervielfältigen, öffentlich ausstellen, öffentlich zugänglich machen, vertreiben oder anderweitig nutzen. Sofern die Verfasser die Dokumente unter Open-Content-Lizenzen (insbesondere CC-Lizenzen) zur Verfügung gestellt haben sollten, gelten abweichend von diesen Nutzungsbedingungen die in der dort genannten Lizenz gewährten Nutzungsrechte. Terms of use: Documents in EconStor may be saved and copied for your personal and scholarly purposes. You are not to copy documents for public or commercial purposes, to exhibit the documents publicly, to make them publicly available on the internet, or to distribute or otherwise use the documents in public. If the documents have been made available under an Open Content Licence (especially Creative Commons Licences), you may exercise further usage rights as specified in the indicated licence. http://creativecommons.org/licenses/by/4.0/ risks Article Sharp Convex Bounds on the Aggregate Sums–An Alternative Proof Chuancun Yin * and Dan Zhu School of Statistics, Qufu Normal University, Qufu 273165, Shandong, China; [email protected] *Correspondence: [email protected]; Tel.: +86-537-4453221 Academic Editor: Mogens Steffensen Received: 8 May 2016; Accepted: 21 September 2016; Published: 29 September 2016 Abstract: It is well known that a random vector with given marginals is comonotonic if and only if it has the largest convex sum, and that a random vector with given marginals (under an additional condition) is mutually exclusive if and only if it has the minimal convex sum. This paper provides an alternative proof of these two results using the theories of distortion risk measure and expected utility. Keywords: comonotonicity; convex order; distortion risk measure; mutual exclusivity; stop-loss order 1. Introduction After years of efforts made by researchers, the study of sharp convex bounds on the sum of random variables (also known as aggregate sums) with given marginal distributions but unknown dependence structure has achieved many significant results. Mathematically, given an arbitrary Fréchet space R(F1 , . . . , Fn) of all random vectors having F1 , . . . , Fn as marginal distributions, the aim is to find two random vectors (Xm 1,...,Xm n)and (XM 1,...,XM n)belonging to R(F1,...,Fn), such that n ∑ i=1 Xm i≤cx n ∑ i=1 Xi≤cx n ∑ i=1 XM i for any (X1 , . .. , Xn)∈ R(F1 , . .. , Fn) , where ≤cx denotes the convex order. By definition, for a pair of random variables X and Y , we say that X is less than Y in the sense of convex order, denoted as X≤cx Y , if E f (X)≤E f (Y) for every convex function f , provided that expectations E f (X) and E f (Y) exist. In actuarial science, it is common to define convex order by using a stop-loss transform: X≤cx Y⇔EX =EY and X≤sl Y . Here X is said to precede Y in the stop-loss order sense, notation X≤sl Y, if and only if Xhas lower stop-loss premiums than Y: E(X−d)+≤E(Y−d)+,−∞<d<∞. A summary of other characterizations and properties of convex order can be found, for example, in [1,2]. Comonotonicity plays a crucial role in determining convex upper bound on aggregate sum. Let us recall the definition. For any X∈ R(F1,...,Fn),Xis said to be comonotonic if FX(x) = min 1≤k≤nFk(xk),∀x= (x1,x2,...,xn)∈Rn. Equivalently, X is comonotonic if and only if Xd = (F−1 1(U) , . .. , F−1 n(U)) , where U is a random variable uniformly distributed on the interval [ 0, 1 ] , denoted as U∼U[ 0, 1 ] . The concept of comonotonicity was introduced by Yaari [3] and Schmeidler [4]. For more details and other characterizations about the concept of comonotonicity and its applications in actuarial science and Risks 2016,4, 34; doi:10.3390/risks4040034 www.mdpi.com/journal/risks Risks 2016,4, 34 2 of 8 finance, we refer to the overview papers by Dhaene et al. [5,6] and more recently in [7]. Let S be the sum X1+ ... +Xn and Sc be the comonotonic sum Xc 1+ ... +Xc n , where (Xc 1 ,..., Xc n) is the comonotonic counterpart of X= (X1 ,..., Xn) . A well-known result between the sums S and Sc says that S≤cx Sc . Proofs of this fundamental result can be found in [8–12]. Müller [13] extended the result to higher dimensions as a special case of the concept of supermodular ordering. A simple geometric argument is given in [14], and Cheung [15] provided a new proof using the theory of majorization. The converse remains valid under the assumption that all marginal distribution functions are continuous and that the underlying probability space (Ω , F , P) is atomless. For more details, see [16]. This continuity assumption on the marginals was removed by Cheung [17]. A simple proof without the assumption that the underlying probability space (Ω , F , P) is atomless was given by Mao and Hu [18]. Some equivalent conditions on comonotonicity can be found in [19]. To summarize the above results, we arrive at the following theorem: Theorem 1. If (X∗ 1,...,X∗ n)∈ R(F1,...,Fn), then (X∗ 1,...,X∗ n)is comonotonic if, and only if X1+... +Xn≤cx X∗ 1+... +X∗ nfor all (X1,...,Xn)∈ R(F1,...,Fn). Now we focus on the lower convex bound of R(F1 , . .. , Fn) . When n= 2, the minimum sharp bound is obtained by the counter-monotonic scenario: F−1 1(U) + F−1 2(1−U)≤cx X1+X2for any (X1,X2)∈ R(F1,F2), where U∼U[ 0, 1 ] . Proofs for this assertion can be found in [9,19]. Moreover, Cheung and Lo [20] shows that the converse remains valid. However, the sharp lower convex bound for n≥ 3 is missing in general. Bernard et al. [21] gave an example showing that ∑n i=1Xi does not have a sharp lower convex bound. Sufficient conditions for the existence of a such sharp lower convex bound for some classes of distributions can be found in [19,21–24]. In another special case, when F1 , .. ., Fn are defined on [ 0, ∞) with ∑n i=1(1−Fi(0)) ≤1, the convex lower bound is obtained by the mutually exclusive scenario: X∗ 1+... +X∗ n≤cx X1+... +Xn for any (X1 , . .. , Xn)∈ R(F1 , . .. , Fn) , where (X∗ 1 , . .. , X∗ n)∈ R(F1 , . .. , Fn) and P(X∗ i> 0, X∗ j> 0 ) = 0 for all i6=j (see [25,26]). Mutual exclusivity can be considered as the strongest negative dependence structure in a multivariate setting. It was first studied in [26] when the marginals F1 , F2 , . .. , Fn are two-point distributions and in Dhaene and Denuit [25] in a more general setting. A revisited and further characterized treatment of mutual exclusivity can be found in [27]. A recent overview paper by Puccetti and Wang [19] introduced the concept of pairwise countermonotonicity, which is more general than that of mutual exclusivity. Moreover, several equivalent conditions on pairwise countermonotonicity have been provided (see Theorem 3.3 in [19]). Definition 1. (Definition 3.4 in [27]) Let X1 ,..., Xn be random variables with essential infima l1 ,..., ln and essential suprema u1,...,un, respectively. They are said to be (i) mutually exclusive from below if P(Xi>li,Xj>lj) = 0for all i 6=j; (ii) mutually exclusive from above if P(Xi<ui,Xj<uj) = 0for all i 6=j. The following theorem is concerned with mutually exclusive random variables and the minimal lower bound in convex order. Risks 2016,4, 34 3 of 8 Theorem 2. (Theorem 5.1 in [27]) Let X∗= (X∗ 1 , . . . , X∗ n) be a fixed random vector in R(F1 , . . . , Fn) ( n≥ 3) which satisfies ∑n i=1(1−Fi(li)) ≤1or ∑n i=1Fi(ui−)≤1. Then X∗is mutually exclusive if, and only if X∗ 1+... +X∗ n≤cx X1+... +Xn for all (X1,...,Xn)∈ R(F1,...,Fn). In this short note, we give alternative proofs of Theorems 1and 2in the next two sections. 2. Proof of Theorem 1 To prove Theorem 1, we need two useful lemmas. Here are some notations. Let FX be the cumulative distribution function of random variable X , and the decumulative distribution function is denoted by ¯ FX ; i.e., ¯ FX(x) = 1 −FX(x) = P(X>x) . A distortion function is defined as a non-decreasing function g:[ 0, 1 ]→[ 0, 1 ] such that g( 0 ) = 0 and g( 1 ) = 1. The distortion risk measure associated with distortion function gis defined by ρg[X] = Z+∞ 0g(¯ FX(x))dx +Z0 −∞[g(¯ FX(x)) −1]dx, for any random variable X , provided at least one of the two integrals above is finite. If X is a non-negative random variable, then ρgreduces to ρg[X] = Z+∞ 0g(¯ FX(x))dx. Obviously, a concave distortion function is continuous on ( 0, 1 ] and can only jump at 0. In view of Theorem 6 of Dhaene et al. [28], we know that for any concave distortion function g , one can rewrite ρg[X]as ρg[X] = Z[0,1]VaR1−q[X]dg(q), where VaRp[X]denotes the value-at-risk at level pof Xand is defined as VaRp[X] = inf{x∈R|FX(x)≥p},p∈(0, 1). The following theorem shows that stop-loss order can be characterized in terms of ordered concave distortion risk measures (see [29,30]). Here we provide a short proof. Lemma 1. For any random vector (X , Y) , we have that X≤sl Y if and only if their respective concave distortion risk measures are ordered: X≤sl Y⇔ρg[X]≤ρg[Y] for all concave distortion functions g . In particular, if E[X] = E[Y], then X ≤cx Y⇔ρg[X]≤ρg[Y]for all concave distortion functions g. Proof. For any concave distortion function g , it is differentiable at all but at most countably many points on [0, 1]. The distortion risk measure ρgcan be written as ρg[X] = Z1 0TVaRp[X]dµ(p), where µ(p):=Rp 0( 1 −α)dν(α) is a probability measure in which ν is defined by ν([ 0, p]) = g0( 1 −p) , TVaRpis the tail value-at-risk(also known as the expected shortfall), defined as TVaRp[X] = 1 1−pZ1 pVaRw[X]dw,p∈(0, 1), Risks 2016,4, 34 4 of 8 which is a distortion risk measure corresponding to the concave distortion function g(x) = min x 1−p, 1, 0 <p<1. The result follows as X≤sl Y⇔TVaRp[X]≤TVaRp[Y] for all p∈( 0, 1 ) (see Theorem 3.2 in [30]). The following subadditivity theorem can be found in [29], the bivariate case can be found in [31], see also [32]. Lemma 2. For any concave distortion function g and (X1,...,Xn)∈ R(F1,...,Fn), we have ρg[X1+... +Xn]≤ρg[X1] + ... +ρg[Xn]. Proof of Theorem 1. First we assume (X∗ 1 , . .. , X∗ n)∈ R(F1 , . .. , Fn) is comonotonic. For any concave distortion function gand (X1,...,Xn)∈ R(F1,...,Fn), by Lemma 2we have ρg[X1+... +Xn]≤ρg[X1] + ... +ρg[Xn]. (1) Comonotonicity of (X∗ 1,...,X∗ n)∈ R(F1,...,Fn)implies that (cf. Dhaene et al. [30]) ρg[X1] + ... +ρg[Xn] = ρg[X∗ 1+... +X∗ n]. (2) Therefore, combining (1) with (2), one has ρg[X1+... +Xn]≤ρg[X∗ 1+... +X∗ n], and the desired result follows from Lemma 1. To prove the other implication, we assume that (X∗ 1,...,X∗ n)∈ R(F1,...,Fn)and X1+... +Xn≤cx X∗ 1+... +X∗ nfor all (X1,...,Xn)∈ R(F1,...,Fn). From Lemma 1, we have that ρg[X1+... +Xn]≤ρg[X∗ 1+... +X∗ n], for any concave distortion function g. In particular, ρg[Xc 1+... +Xc n]≤ρg[X∗ 1+... +X∗ n], (3) where (Xc 1 ,..., Xc n) is the comonotonic counterpart of (X1 ,..., Xn) . On the other hand, by Lemma 2, we get ρg[X∗ 1+... +X∗ n]≤ρg[X∗ 1] + ... +ρg[X∗ n]. (4) Note that ρg[Xc 1+... +Xc n] = ρg[Xc 1] + ... +ρg[Xc n], (5) and ρg[Xc 1] + ... +ρg[Xc n] = ρg[X∗ 1] + ... +ρg[X∗ n]. (6) It follows from (3)–(6) that we have ρg[Xc 1+... +Xc n] = ρg[X∗ 1+... +X∗ n], (7) Risks 2016,4, 34 5 of 8 for any concave distortion function g . Therefore, Theorem 7 in [17] implies (X∗ 1 ,..., X∗ n) is comonotonic. This ends the proof of Theorem 1. 3. Proof of Theorem 2 To prove Theorem 2, we need two useful lemmas. Lemma 3gives a necessary and sufficient condition for the convex order of two random variables. Lemma 3. (Proposition 3.4.3 in [1]) Given two rvs X and Y, then the following statements are equivalent: (1) X≤cx Y. (2) E[v(X)] ≤E[v(Y)] for all convex functions v, such that the expectations exist. (3) E[v(X)] ≤E[v(Y)] for all functions v with v00 ≥0, such that the expectations exist. The following lemma, due to Cheung and Lo [33], will play a crucial role in the proof of Theorem 2. Lemma 4. (Theorem 3.1 in [33]) Let X1 , . . . , Xn be non-negative random variables, and f be a convex function such that E[f(∑n i=1Xi)] exists. (i) We have E"f n ∑ i=1 Xi!#≥ n ∑ i=1 E[f(Xi)] −(n−1)f(0); (ii) if f is strictly convex, then E"f n ∑ i=1 Xi!#= n ∑ i=1 E[f(Xi)] −(n−1)f(0) if, and only if X1,...,Xnare mutually exclusive random variables in the sense of Dhaene and Denuit [25]. Remark 1. We remark that the “if part" is still true when the function f is convex, but not necessarily strictly convex. Proof of Theorem 2. To prove Theorem 2, as in the proof to Lemma 3.6 in [27], there are three cases to consider. Recall that l1 ,..., ln are the essential infima of random variables X1 ,..., Xn , respectively (see Definition 1). Case 1. l1= ... =ln= 0. We assume (X∗ 1 ,..., X∗ n)∈ R(F1 ,..., Fn) is mutually exclusive. For any convex function uand (X1,...,Xn)∈ R(F1,...,Fn), by Lemma 4(i) we have E"u n ∑ i=1 Xi!#≥ n ∑ i=1 E[u(Xi)] −(n−1)u(0). (8) Thanks to Lemma 4(ii) and Remark 1, mutual exclusivity of (X∗ 1,...,X∗ n)implies that E"u n ∑ i=1 X∗ i!#= n ∑ i=1 E[u(X∗ i)] −(n−1)u(0). (9) Therefore, combining (8) with (9), and noting that E[u(X∗ 1)] + ... +E[u(X∗ n)] = E[u(X1)] + ... + E[u(Xn)], one has E"u n ∑ i=1 X∗ i!#≤E"u n ∑ i=1 Xi!#, from which and from Lemma 3, we deduce that X∗ 1+... +X∗ n≤cx X1+... +Xn Risks 2016,4, 34 6 of 8 for all (X1,...,Xn)∈ R(F1,...,Fn). To prove the other implication, we assume that (X∗ 1,...,X∗ n)∈ R(F1,...,Fn)and X∗ 1+... +X∗ n≤cx X1+... +Xnfor all (X1,...,Xn)∈ R(F1,...,Fn). From Lemma 3, we have that E[u(X∗ 1+... +X∗ n)] ≤E[u(X1+... +Xn)] for all convex functions u. In particular, E[u(X∗ 1+... +X∗ n)] ≤E[u(XM 1+... +XM n)] (10) where (XM 1 ,..., XM n) is the mutually exclusive counterpart of (X1 ,..., Xn) . On the other hand, by Lemma 4and Remark 1, we get E[u(X∗ 1+... +X∗ n)] ≥E[u(X∗ 1)] + ... +E[u(X∗ n))] −(n−1)u(0), (11) and E[u(XM 1+... +XM n)] = E[u(XM 1)] + ... +E[u(XM n))] −(n−1)u(0). (12) If (X∗ 1,...,X∗ n)is not mutually exclusive, then, for strict convex u, E[u(X∗ 1+... +X∗ n)] 6=E[u(X∗ 1)] + ... +E[u(X∗ n))] −(n−1)u(0). (13) Combining (10)–(12) with (13), we get E[u(X∗ 1)] + ... +E[u(X∗ n))] <E[u(XM 1)] + ... +E[u(XM n))], for any strict convex function u . This contradicts (X∗ 1 ,..., X∗ n) and (XM 1 ,..., XM n) having the same marginals. Thus, (X∗ 1,...,X∗ n)is mutually exclusive. Case 2. (X∗ 1 ,..., X∗ n) is mutually exclusive from below. For any (X1 ,..., Xn)∈ R(F1 ,..., Fn) , then Z:=Xi−li are non-negative random variables, and Z∗:=X∗ i−li are non-negative mutually exclusive random variables. Applying the result in Case 1, we obtain that (X∗ 1 ,..., X∗ n)is mutually exclusive ⇔ (X∗ 1−l1,...,X∗ n−ln)is mutually exclusive ⇔ n ∑ i=1 (X∗ i−li)≤cx n ∑ i=1 (Xi−li)⇔ n ∑ i=1 X∗ i− n ∑ i=1 li≤cx n ∑ i=1 Xi− n ∑ i=1 li ⇔ n ∑ i=1 X∗ i≤cx n ∑ i=1 Xi. Case 3. (X∗ 1 ,..., X∗ n) is mutually exclusive from above. For any (X1 ,..., Xn)∈ R(F1 ,..., Fn) , applying the result in Case 2, we have (X∗ 1 ,..., X∗ n) is mutually exclusive from above ⇔(−X∗ 1 ,..., −X∗ n) is mutually exclusive from below ⇔ − ∑n i=1X∗ i≤cx −∑n i=1Xi⇔∑n i=1X∗ i≤cx ∑n i=1Xi . The proof of Theorem 2is now complete. Acknowledgments: We are greatly indebted to Ruodu Wang and two anonymous referees for their insightful comments. The research was supported by the National Natural Science Foundation of China (No. 11171179, 11571198) and the Research Fund for the Doctoral Program of Higher Education of China (No. 20133705110002). Author Contributions: These authors contributed equally to this work. Conflicts of Interest: The authors declare no conflict of interest. Risks 2016,4, 34 7 of 8 References 1. Denuit, M.; Dhaene, J.; Goovaerts, M.; Kaas, R. Actuarial Theory for Dependent Risks: Measures, Orders and Models; John Wiley & Sons: Chichester, UK, 2005. 2. Shaked, M.; Shanthikumar, J.G. Stochastic Orders; Springer: New York, NY, USA, 2007. 3. Yaari, M.E. The dual theory of choice under risk. Econometrica 1987,55, 95–115. 4. Schmeidler, D. Integral representation without additivity. Proc. Am. Math. Soc. 1986,97, 255–261. 5. Dhaene, J.; Denuit, M.; Goovaerts, M.J.; Kaas, R.; Vyncke, D. The concept of comonotonicity in actuarial science and finance: Theory. Insur. Math. Econ. 2002,31, 3–33. 6. Dhaene, J.; Denuit, M.; Goovaerts, M.J.; Kaas, R.; Vyncke, D. The concept of comonotonicity in actuarial science and finance: Applications. Insur. Math. Econ. 2002,31, 133–161. 7. Deelstra, G.; Dhaene, J.; Vanmaele, M. An overview of comonotonicity and its applications in finance and insurance. In Advanced Mathematical Methods for Finance; Oksendal, B., Nunno, G., Eds.; Springer: Heidelberg, Germany, 2010; pp. 155–179. 8. Meilijson, I.; Nadas, A. Convex majorization with an application to the length of critical paths. J. Appl. Probab. 1979,16, 671–677. 9. Tchen, A.H. Inequalities for distributions with given marginals. Ann. Probab. 1980,8, 814–827. 10. Ru ¨ schendorf, L. Solution of a statistical optimization problem by rearrangement methods. Metrika 1983 ,30, 55–61. 11. Dhaene, J.; Goovaerts, M.J. Dependency of risks and stop-loss order. ASTIN Bull. 1996,26, 201–212. 12. Dhaene, J.; Goovaerts, M.J. On the dependency of risks in the individual life model. Insur. Math. Econ. 1997 , 19, 243–253. 13. Müller, A. Stop-loss order for portfolios of dependent risks. Insur. Math. Econ. 1997,21, 219–223. 14. Kaas, R.; Dhaene, J.; Vyncke, D.; Goovaerts, M.; Denuit, M. A simple geometric proof that comonotonic risks have the convex-largest sum. ASTIN Bull. 2002,32, 71–80. 15. Cheung, K.C. Comonotonic convex upper bound and majorization. Insur. Math. Econ. 2010,47, 154–158. 16. Cheung, K.C. Characterization of comonotonicity using convex order. Insur. Math. Econ. 2008,43, 403–406. 17. Cheung, K.C. Characterizing a comonotonic random vector by the distribution of the sum of its components. Insur. Math. Econ. 2010,47, 130–136. 18. Mao, T.T.; Hu, T.Z. A new proof of Cheung’s characterization of comonotonicity. Insur. Math. Econ. 2011 ,48, 214–216. 19. Puccetti, G.; Ruodu Wang, R. Extremal dependence concepts. Stat. Sci. 2015,30, 485–517. 20. Cheung, K.C.; Lo, A. Characterizations of counter-monotonicity and upper comonotonicity by (tail) convex order. Insur. Math. Econ. 2013,53, 334–342. 21. Bernard, C.; Jiang, X.; Wang, R. Risk aggregation with dependence uncertainty. Insur. Math. Econ. 2014 ,54, 93–108. 22. Wang, B.; Wang, R. The complete mixability and convex minimization problems for monotone marginal distributions. J. Multivar. Anal. 2011,102, 1344–1360. 23. Wang, B.; Wang, R. Joint mixability. Math. Op. Res. 2016,41, 808–826. 24. Jakobsons, E.; Han, X.; Wang, R. General convex order on risk aggregation. Scand. Actuar. J. 2016 ,8, 713–740. 25. Dhaene, J.; Denuit, M. The safest dependence structure among risks. Insur. Math. Econ. 1999,25, 11–21. 26. Hu, T.Z.; Wang, Z.Q. On the dependence of risks and the stop-loss premiums. Insur. Math. Econ. 1999 ,24, 323–332. 27. Cheung, K.C.; Lo, A. Characterizing mutual exclusivity as the strongest negative multivariate dependence structure. Insur. Math. Econ. 2014,55, 180–190. 28. Dhaene, J.; Kukush, A.; Linders, D.; Tang, Q. Remarks on quantiles and distortion risk measures. Eur. Actuar. J. 2012,2, 319–328. 29. Dhaene, J.; Wang, S.; Young, V.; Goovaerts, M. Comonotonicity and maximal stop-loss premiums. Bull. Swiss Assoc. Actuar. 2000,2, 99–113. 30. Dhaene, J.; Vanduffel, S.; Tang, Q.; Goovaerts, M.; Kaas, R.; Vyncke, D. Risk measures and comonotonicity: A review. Stoch. Models 2006,22, 573–606. 31. Denneberg, D. Non-Additive Measure and Integral; Kluwer Academic Publilshers: Dordrecht, The Netherlands, 1994. Risks 2016,4, 34 8 of 8 32. Wang, S.; Dhaene, J. Comonotonicity, correlation order and premium principles. Insur. Math. Econ. 1998 ,22, 235–242. 33. Cheung, K.C.; Lo, A. General lower bounds on convex functionals of aggregate sums. Insur. Math. Econ. 2013,53, 884–896. © 2016 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).