In various research areas related to decision making, problems and their solutions frequently rely on certain functions being monotonic. In the case of non-monotonic functions, one would then wish to quantify their lack of monotonicity. In this paper we develop a method designed specifically for this task, including quantification of the lack of positivity, negativity, or sign-constancy in signed measures. We note relevant applications in Insurance, Finance, and Economics, and discuss some of them in detail.
Non-monotonic functionssigned measuresHahn and Jordan decompositionsweighted premiumrisk measuregain–loss ratio28E0526A4862P0597M30Introduction
In various research areas such as Economics, Finance, Insurance, and Statistics, problems and their solutions frequently rely on certain functions being monotonic, as well as on determining the degree of their monotonicity, or lack of it. For example, the notion of profit seekers in Behavioural Economics and Finance is based on increasing utility functions, which can have varying shapes and thus characterize subclasses of profit seekers (e.g., [9]). In Reliability Engineering and Risk Assessment (e.g., [13]), a number of notions such as hazard-rate and likelihood-ratio orderings rely on monotonicity of the ratios of certain functions. The presence of insurance deductibles and policy limits often change the pattern of monotonicity of insurance losses (e.g., [3]). In the literature on Statistical Inference, the so-called monotone-likelihood-ratio family plays an important role.
Due to these and a myriad of other reasons, researchers quite often restrict themselves to function classes with pre-specified monotonicity properties. But one may not be comfortable with this element of subjectivity and would therefore prefer to rely on data-implied shapes when making decisions. To illustrate the point, we recall, for example, the work of Bebbington et al. [1] who specifically set out to determine whether mortality continues to increase, or starts to decelerate, after a certain species-related late-life age. This is known in the gerontological literature as the late-life mortality deceleration phenomenon. Naturally, we refrain from elaborating on this complex topic and refer for details and further references to the aforementioned paper.
Monotonicity may indeed be necessary for certain results or properties to hold, but there are also many instances when monotonicity is just a sufficient condition. In such cases, a natural question arises: can we depart from monotonicity and still have valid results? Furthermore, in some cases, monotonicity may not even be expected to hold, though perhaps be desirable, and so developing techniques for quantifying the lack of monotonicity becomes of interest. Several results in this direction have recently been proposed in the literature (e.g., [5, 21, 15, 16]; and references therein). All in all, these are some of the problems, solutions to which ask for indices that could be used to assess monotonicity, or lack of it. In the following sections we introduce and discuss several such indices, each designed to reveal different monotonicity aspects.
We have organized the rest of the paper as follows. In Section 2 we present a specific example, driven by insurance and econometric considerations, which not only motivates present research but also highlights the main idea adopted in this paper. In Section 3 we introduce and discuss indices of lack of increase, decrease, and monotonicity. We also suggest a convenient numerical procedure for a quick calculation of the indices at any desirable precision. In Section 4 we introduce orderings of functions according to the values of their indices and argue in favour of using normalized versions of the indices. In Section 5 we introduce a stricter notion of ordering, and in Section 6 we develop a theory for quantifying the lack of positivity (or negativity) in signed measures. Section 7 concludes the paper.
Motivating example
Insurance losses are non-negative random variables X≥0, and the expected value E[X] is called the net premium. All practically-sound premium calculation principles (pcp’s), which are functionals π assigning non-negative finite or infinite values to X’s, are such that π[X]≥E[X]. The latter property is called non-negative loading of π.
As an example, consider the following ‘dual’ version of the weighted pcp ([7, 8]; and references therein):
πw[X]=E[F−1(U)w(U)]E[w(U)],
where U is the random variable uniform on [0,1], F−1 is the quantile function of X defined by the equation F−1(p)=inf{x:F(x)≥p}, and w:[0,1]→[0,∞] is an appropriately chosen weight function, whose illustrative examples will be provided in a moment. We of course assume that the two expectations in the definition of πw are well-defined and finite, and the denominator is not zero.
When dealing with insurance losses, researchers usually choose non-decreasing weight functions in order to have non-negative loading of πw. For example, w(t)=1{t>p}p\}$]]> for any parameter p∈(0,1) is non-decreasing, and it turns πw into the average value at risk, also known as tail conditional expectation. Another example is w(t)=ν(1−t)ν−1 with parameter ν>00$]]>. If ν∈(0,1], then w is non-decreasing, and πw becomes the proportional hazards pcp [18, 19]. If ν≥1, then w is non-increasing, and πw becomes the (absolute) S-Gini index of economic equality (e.g., [23]; and references therein).
Note that the pcp πw can be rewritten as the weighted integral
πw=∫01F−1(t)w∗(t)dt,
with the normalized weight function w∗(t)=w(t)/∫01w(u)du, which is a probability density function, because w(t)≥0 for all t∈[0,1] and ∫01w(u)du∈(0,∞). This representation of πw connects our present research with the dual utility theory ([20, 17]; and references therein) that has arisen as a prominent counterpart to the classical utility theory of von Neuman and Morgenstern [14]. It is also important to mention that in Insurance, integral (2) plays a very prominent role and is known as the distortion risk measure (e.g., [6]; and references therein).
In general, the function w and its normalized version w∗ may not be monotonic, as it depends on the shape of the probability density function of the random variable W∈[0,1] in the following reformulation of πw:
πw[X]=E[F−1(W)].
In the econometric language, E[F−1(W)] means the average income (assuming that X stands for ‘income’) possessed by individuals whose positions on the society’s income-percentile scale are modelled by W, which is, naturally, a random variable.
Hence, we are interested when E[F−1(W)] is at least E[W] (insurance perspective) or at most E[W] (income inequality perspective). In view of equations (1) and (2), our task reduces to verifying whether or not
cov[F−1(U),w(U)]≥0.
If the function w happens to be non-decreasing (as in insurance), then bound (4) holds (e.g., [11]). For example, with parameter λ>00$]]>, the weight function w(s)=sλ leads to the size-biased pcp, w(s)=eλs to the Esscher pcp, w(s)=1−e−λs to the Kamps pcp, and the already noted function w(s)=1{s>λ}\lambda \}$]]> leads to the average value at risk.
As already noted, broader contexts than that of classical insurance suggest various shapes of probability distortions and thus lead to functions w that are not necessarily monotonic. Indeed, in view of representation (3), the average income of individuals depends on the distribution of W, whose shape is governed by societal opportunities that the individuals are exposed to. A natural question arises:
What shapes of w could ensure property (4)?
To answer this question in an illuminating way, let w be absolutely continuous and such that w(0)=0, which are sound assumptions from the practical point of view. Denote the density of w by w′, and let 1{u>t}t\}}$]]> be the indicator. We have
cov[F−1(U),w(U)]=cov[F−1(U),∫0Uw′(t)dt]=∫01cov[F−1(U),1{U>t}]w′(t)dt.t\}}\big]{w^{\prime }}(t)\hspace{0.1667em}\mathrm{d}t.\end{array}\]]]>
The function v(t)=cov[F−1(U),1{U>t}]t\}}]$]]> is non-negative (e.g., [11]), and so is also the integral θ=∫01v(t)dt, which makes v(t)/θ a probability density function. Denote the corresponding distribution function by V. We have the equations
cov[F−1(U),w(U)]=θ∫01w′dV=θ∫01w′∘V−1dλ,
where λ denotes the Lebesgue measure. Hence, property (4) is equivalent to the bound ∫01gdλ≥0 with g=w′∘V−1, and it can of course be rewritten as the inequality
∫01g+dλ≥∫01g−dλ,
where g+=max{g,0} and g−=max{−g,0}. Bound (5) says that, in average with respect to the Lebesgue measure, the positive part g+ must be larger than the negative part g−.
Those familiar with asset pricing will immediately see how inequality (5), especially when reformulated as
∫01g+dλ∫01g−dλ≥1
is connected to the gain–loss ratio [2], as well as to the Omega ratio [10]. Furthermore, we can reformulate bound (6) as
∫01g+dλ∫01|g|dλ≥12,
with the ratio on the left-hand side being equal to 0 when the function g is non-positive, and equal to 1 when g is non-negative. Hence, bound (7) says that the ratio must be at least 1/2, which means that, in average, the function g must be more positive than negative.
The above interpretations have shaped our considerations in the present paper, and have led toward the construction of monotonicity indices that we introduce and discuss next.
Assessing lack of monotonicity in functions
We are interested in assessing monotonicity of a function g0 on an interval [a,b]. Since shifting the function up or down, left or right, does not distort its monotonicity, we therefore ‘standardize’ g0 into
g(x)=g0(x+a)−g0(a)
defined on the interval [0,y] with y=b−a. Note that g satisfies the boundary condition g(0)=0. We assume that g0 and thus g are absolutely continuous.
Let Fy denote the set of all absolutely continuous functions f on the interval [0,y] such that f(0)=0. Denote the total variation of f∈Fy on the interval [0,y] by ‖f‖y, that is, ‖f‖y=∫0y|f′|dλ. Furthermore, let Fy+ denote the set of all f∈Fy that are non-decreasing. For any g∈Fy, we define its index of lack of increase (LOI) as the distance between g and the set Fy+, that is,
LOIy(g)=inf{∫0y|g′−f′|dλ:f∈Fy+}.
Obviously, if g is non-decreasing, then LOIy(g)=0, and the larger the value of LOIy(g), the farther the function g is from being non-decreasing on the interval [0,y]. For the function g0, which is where g originates from, the LOI of g0 on the interval [a,b] is
LOI[a,b](g0):=LOIy(g).
(Throughout the paper we occasionally use ‘:=,’ when a need arises to emphasize that certain equations are by definition.) Determining LOIy(g) using definition (9) is not straightforward. To facilitate the task, we next give an integral representation of the index.
The infimum in definition (9) is attained at a functionf1∈Fy+such thatf1′=(g′)+, and thusLOIy(g)=∫0y(g′)−dλ.When g originates fromg0via equation (8), we haveLOI[a,b](g0)=∫ab(g0′)−dλ.
Theorem 1, though easy to prove directly, follows immediately from a more general result, Theorem 2 of Section 6, and we thus do not provide any more details. The index of lack of decrease (LOD) is defined analogously. Namely, in the computationally convenient form, it is given by the equation
LODy(g)=∫0y(g′)+dλ.
When g originates from g0 via equation (8), we have
LOD[a,b](g0)=∫ab(g0′)+dλ.
In turn, the index of lack of monotonicity (LOM) of g is given by the equation
LOMy(g)=2min{LOIy(g),LODy(g)},
and when g originates from g0 via equation (8), we have
LOM[a,b](g0)=2min{LOI[a,b](g0),LOD[a,b](g0)}.
The reason for doubling the minimum on the right-hand sides of definitions (12) and (13) will become clear from properties below.
The index LOIy is translation invariant, that is, LOIy(g+α)=LOIy(g) for every constant α∈R. The index LODy is also translation invariant.
The index LOIy is positively homogeneous, that is, LOIy(βg)=βLOIy(g) for every constant β≥0. The index LODy is also positively homogeneous.
LOIy(βg)=(−β)LODy(g) for every negative constant β<0, and thus in particular, LOIy(−g)=LODy(g).
LOIy(g)+LODy(g)=‖g‖y. Consequently, we have the following two observations:
LOIy(g) and LODy(g) do not exceed ‖g‖y, with both indices achieving the upper bound. Namely, LOIy(g)=‖g‖y whenever (g′)+≡0, and LODy(g)=‖g‖y whenever (g′)−≡0.
min{LOIy(g),LODy(g)}≤‖g‖y/2 and thus LOMy(g)≤‖g‖y, which justifies the use of the factor 2 in definitions (12) and (13).
An illustrative example follows.
Consider the functions g0(z)=sin(z) and h0(z)=cos(z) on [−π/2,π]. Neither of them is monotonic on the interval, but a visual inspection of their graphs suggests that sine is closer to being increasing than cosine, which can of course be viewed as a subjective statement. To substantiate it, we employ the above introduced indices. First, by lifting and shifting, we turn sine into g(x)=1−cos(x) and cosine into h(x)=sin(x). Since g′(x)=sin(x) and h′(x)=cos(x), we have ∫03π/2(g′)−dλ=1and∫03π/2(g′)+dλ=2,∫03π/2(h′)−dλ=2and∫03π/2(h′)+dλ=1. Consequently,
LOI[−π/2,π](g0)=LOI3π/2(g)=1,LOI[−π/2,π](h0)=LOI3π/2(h)=2,LOD[−π/2,π](g0)=LOD3π/2(g)=2,LOD[−π/2,π](h0)=LOD3π/2(h)=1,LOM[−π/2,π](g0)=LOM3π/2(g)=2,LOM[−π/2,π](h0)=LOM3π/2(h)=2.
Hence, for example, sine is at the distance 1 from the set of all non-decreasing functions on the noted interval, whereas cosine is at the distance 2 from the same set. Note also that the total variations ∫03π/2|g′|dλ and ∫03π/2|h′|dλ of the two functions are the same, equal to 3. This concludes Example 1.
Example 1 is based on functions for which the four integrals in equations (14) and (15) are easy to calculate, but functions arising in applications are frequently quite unwieldy. For this, we need a numerical procedure for calculating the integral ∫0yH(f′)dλ for various transformations H, such as H(x)=x− and H(x)=x+. A convenient way is as follows:
∫0yH(f′)dλ≈∑n=1NyNH(Ny{f(nNy)−f(n−1Ny)})=∑n=1Nb−aNH(Nb−a{f0(a+nN(b−a))−f0(a+n−1N(b−a))}),
where f0 is the underlying function (on the interval [a,b]) and f is the shifted-and-lifted function (on the interval [0,y]) defined by the equation f(x)=f0(x+a)−f0(a) for all x∈[0,y] with y=b−a.
The use of the integral ∫0yH(f′)dλ in Note 1 hints at the possibility of distances other than L1-norms when defining indices. One may indeed wish to de-emphasize small values of f′ and to emphasize its large values when defining indices. In a simple way, this can be achieved by taking the p-th power of (f′)−, (f′)+ and f′ for any p≥1. This argument leads us to the Lp-type counterpart of minimization problem (9) which can be solved along the same path as that used in the proof of Theorem 1. It is remarkable that the minimizing function does not depend on the choice of the metric and remains equal to (f′)+. Consequently, for example, the following Lp-index of lack of increase arises:
LOIp,y(f)=(∫0y((f′)−)pdλ)1/p.
Other than Lp-norms can also be successfully explored, and this is important in applications, where no p-th power may adequately (de-)emphasize parts of f′. The phenomenon has prominently manifested, for example, in Econometrics (e.g., [4, 22]; and references therein). In such cases, more complexly shaped functions are typically used, including those H:[0,∞)→[0,∞) with H(0)=0 that give rise to the class of Birnbaum–Orlicz (BO) spaces.
Monotonicity comparisons and normalized indices
The index values in Example 1 suggest that on the noted interval, sine is more increasing than cosine, because sine is closer to the set Fy+ than cosine is. In general, given two functions g1 and h1 on the interval [0,y], we can say that the function g1 is more non-decreasing than h1 on the interval whenever LOIy(g1)≤LOIy(h1). Likewise, we can say that the function g2 is more non-increasing than h2 on the interval [0,y] whenever LODy(g2)≤LODy(h2). But an issue arises with these definitions because for a given pair of functions g and h, the property LOIy(g)≤LOIy(h) may not be equivalent to LODy(g)≥LODy(h). Though it may look strange at first sight, this non-reflexivity is natural because the total variations of the functions g and h on the interval [0,y] may not be equal, and in such cases, comparing non-monotonicities of g and h is not meaningful.
If, however, the total variations of g and h are equal on the interval [0,y], then LOIy(g)≤LOIy(h) if and only if LODy(g)≥LODy(h). This suggests that in order to achieve this ‘if and only if’ property in general, we need to normalize the indices, which gives rise to the following definitions
LOIy∗(g)=∫0y(g′)−dλ∫0y|g′|dλandLODy∗(g)=∫0y(g′)+dλ∫0y|g′|dλ,
of the normalized indices of lack of increase and decrease, respectively. Obviously,
LOIy∗(g)+LODy∗(g)=1,
and thus LOIy∗(g)≤LOIy∗(h) if and only if LODy∗(g)≥LODy∗(h). Furthermore, the normalized index of lack of monotonicity is LOMy∗(g):=2min{LOIy∗(g),LODy∗(g)}, which we rewrite as
LOMy∗(g)=1−|g(y)|∫0y|g′|dλ
using the identity min{u,v}=(u+v−|u−v|)/2 that holds for all real numbers u and v. In terms of the original function g0 on the interval [a,b], we have
LOI[a,b]∗(g0)=∫ab(g0′)−dλ∫ab|g0′|dλ,LOD[a,b]∗(g0)=∫ab(g0′)+dλ∫ab|g0′|dλ,
and
LOM[a,b]∗(g0)=1−|g0(b)−g0(a)|∫ab|g0′|dλ.
Obviously, LOI[a,b]∗(g0)=LOIy∗(g), LOD[a,b]∗(g0)=LODy∗(g), and LOM[a,b]∗(g0)=LOMy∗(g). To illustrate the normalized indices, we continue Example 1.
Recall that we are dealing with the functions g0(z)=sin(z) and h0(z)=cos(z) on the interval [−π/2,π]. We transform them into the functions g(x)=1−cos(x) and h(x)=sin(x) on the interval [0,3π/2]. From equations (14) and (15), we see that the total variations ∫03π/2|g′|dλ and ∫03π/2|h′|dλ are equal to 3, and so we have the equations:
LOI[−π/2,π]∗(g0)=LOI3π/2∗(g)=13,LOI[−π/2,π]∗(h0)=LOI3π/2∗(h)=23,LOD[−π/2,π]∗(g0)=LOD3π/2∗(g)=23,LOD[−π/2,π]∗(h0)=LOD3π/2∗(h)=13,LOM[−π/2,π]∗(g0)=LOM3π/2∗(g)=23,LOM[−π/2,π]∗(h0)=LOM3π/2∗(h)=23.
The numerical procedure of Note 1 can easily be employed to calculate these normalized indices. This concludes Example 2.
Next are properties of the normalized indices.
The three indices LOIy∗(g), LODy∗(g), and LOMy∗(g) are normalized, that is, take values in the unit interval [0,1], with the following special cases:
LOIy∗(g)=0 if and only if (g′)−≡0, that is, when g is non-decreasing everywhere on [0,y], and LOIy∗(g)=1 if and only if (g′)+≡0, that is, when g is non-increasing everywhere on [0,y].
LODy∗(g)=0 if and only if (g′)+≡0, that is, when g is non-increasing everywhere on [0,y], and LODy∗(g)=1 if and only if (g′)−≡0, that is, when g is non-decreasing everywhere on [0,y].
LOMy∗(g)=0 if and only if g is either non-decreasing everywhere on [0,y] or non-increasing everywhere on [0,y], and LOMy∗(g)=1 if and only if LOIy∗(g)=LODy∗(g) (recall equation (16)).
The three indices LOIy∗(g), LODy∗(g), and LOMy∗(g) are translation invariant, that is, LOIy∗(g+α)=LOIy∗(g) for every constant α∈R, and analogously for LODy∗(g) and LOMy∗(g).
The three indices LOIy∗(g), LODy∗(g), and LOMy∗(g) are positive-scale invariant, that is, LOIy∗(βg)=LOIy∗(g) for every positive constant β>00$]]>, and analogously for LODy∗(g) and LOMy∗(g).
Moreover, LOMy∗(g) is negative-scale invariant, that is, LOMy∗(βg)=LOMy∗(g) for every negative constant β<0, and thus, in general, LOMy∗(βg)=LOMy∗(g) for every real constant β≠0 (recall equation (17)).
LODy∗(−g)=LOIy∗(g), and thus LODy∗(βg)=LOIy∗(g) for every negative constant β<0.
We next use the above indices to introduce three new orderings:
The function g is more non-decreasing than h on the interval [0,y], denoted by g≥I,yh, if and only if LOIy∗(g)≤LOIy∗(h).
The function g is more non-increasing than h on the interval [0,y], denoted by g≥D,yh, if and only if LODy∗(g)≤LODy∗(h).
The function g is more monotonic than h on the interval [0,y], denoted by g≥M,yh, if and only if LOMy∗(g)≤LOMy∗(h).
We see that on the interval [0,y], the function g is more non-decreasing than h if and only if the function g is less non-increasing than h. In other words, g≥I,yh is equivalent to g≤D,yh, which we have achieved by introducing the normalized indices.
Stricter notion of comparison: a note
One of the fundamental notions of ordering random variables is that of first-order stochastic dominance (e.g., [6, 12, 13]). Similarly to this notation, our earlier introduced orderings can be strengthened by first noting that the integral ∫0y(g′)−dλ is equal to ∫0∞Sy−(z∣g)dz, where the function
Sy−(z∣g)=λ{x∈[0,y]:(g′)−(x)>z}z\big\}\]]]>
counts the ‘time’ that the function (g′)− spends above the threshold z during the ‘time’ period [0,y]. Likewise, we define the ‘plus’ version Sy+(z∣g). This lead us to the following definitions of ordering functions according to their monotonicity.
g is more strictly non-decreasing than h on the interval [0,y], denoted by g≥SI,yh, whenever
Sy−(z∣g)∫0y|g′|dλ≤Sy−(z∣h)∫0y|h′|dλfor allz>0.0.\]]]>
g is more strictly non-increasing than h on the interval [0,y], denoted by g≥SD,yh, whenever
Sy+(z∣g)∫0y|g′|dλ≤Sy+(z∣h)∫0y|h′|dλfor allz>0.0.\]]]>
Obviously, if g≥SI,yh, then g≥I,yh, and if g≥SD,yh, then g≥D,yh.
Assessing lack of positivity in signed measures
We now take a path in the direction of general Measure Theory. Namely, let Ω be a set equipped with a sigma-algebra, and let M denote the set of all (signed) measures ν defined on the sigma-algebra. Furthermore, let M+⊂M be the subset of all positive measures. Given a signed-measure ν∈M, we define its index of lack of positivity (LOP) by the equation
LOP(ν)=inf{‖ν−μ‖:μ∈M+},
where ‖·‖ denotes the total variation. Specifically, with (Ω−,Ω+) denoting a Hahn decomposition of Ω, let (ν−,ν+) be the Jordan decomposition of ν. Note that ν− and ν+ are elements of M+. The variation of ν is |ν|=ν−+ν+, and its total variation is ‖ν‖=|ν|(Ω). The following theorem provides an actionable, and crucial for our considerations, reformulation of the index LOP(ν).
The infimum in the definition ofLOP(ν)is attained on the unique element ofM+, which is the measureν+, and thus theLOPindex can be written asLOP(ν)=‖ν−‖.
Since ν+∈M+, we have the bound LOP(ν)≤‖ν−ν+‖, which can be rewritten as LOP(ν)≤‖ν−‖. To prove the opposite bound LOP(ν)≥‖ν−‖, we proceed as follows. For any μ∈M+, we have the bound
‖ν−μ‖≥|ν−μ|(Ω−).
Since ν=ν+−ν− and ν+(A)=0 for every A⊂Ω−, the right-hand side of bound (19) is equal to |ν−+μ|(Ω−), which is not smaller than ν−(Ω−). The latter is, by definition, equal to ‖ν−‖. This establishes the bound LOP(ν)≥‖ν−‖ and completes the proof of equation (18). We still need to show that μ=ν+ is the only measure μ∈M+ such that the equation
‖ν−μ‖=‖ν−‖
holds. Note that ‖ν−‖=ν−(Ω−). Since
‖ν−μ‖≥|ν−+μ|(Ω−)+|ν+−μ|(Ω+),
in order to have equation (20), the right-hand side of inequality (21) must be equal to ν−(Ω−). This can happen only when μ(Ω−)=0 and |ν+−μ|(Ω+)=0, with the former equation implying that the support of μ must be Ω+, and the latter equation implying that μ must be equal to ν+ on Ω+. Hence, μ=ν+. This finishes the proof of Theorem 2. □
Similarly to the LOP index, the index of lack of negativity (LON) of ν∈M is given by the equation
LON(ν)=‖ν+‖,
and the corresponding index of lack of sign (LOS) is
LOS(ν)=2min{‖ν−‖,‖ν+‖}.
The reason for doubling the minimum on the right-hand side of equation (23) is the same as in the more specialized cases discussed earlier. Namely, due to the equation ‖ν−‖+‖ν+‖=‖ν‖, we have that min{‖ν−‖,‖ν+‖} does not exceed ‖ν‖/2. Hence, for LOS to be always between 0 and ‖ν‖, just like LOP and LON are, we need to double the minimum. Finally, we introduce the normalized indices
LOP∗(ν)=‖ν−‖‖ν‖,LON∗(ν)=‖ν+‖‖ν‖,
and
LOS∗(ν)=2min{LOP∗(ν),LON∗(ν)},
whose values are always in the unit interval [0,1].
Conclusion
In this paper, we have introduced indices that, in a natural way, quantify the lack of increase, decrease, and monotonicity of functions, as well as the lack of positivity, negativity, and sign-constancy in signed measures. In addition to being of theoretical interest, this research topic also has practical implications, and for the latter reason, we have also introduced a simple and convenient numerical procedure for calculating the indices without resorting to frequently unwieldy closed-form expressions. The indices satisfy a number of natural properties, and they also facilitate the ranking of functions according to their lack of monotonicity. Relevant applications in Insurance, Finance, and Economics have been pointed out, and some of them discussed in greater detail.
Acknowledgments
We are indebted to Professor Yu. Mishura and two anonymous reviewers for insightful comments and constructive criticism, which guided our work on the revision. We gratefully acknowledge the grants “From Data to Integrated Risk Management and Smart Living: Mathematical Modelling, Statistical Inference, and Decision Making” awarded by the Natural Sciences and Engineering Research Council of Canada to the second author (RZ), and “A New Method for Educational Assessment: Measuring Association via LOC index” awarded by the national research organization Mathematics of Information Technology and Complex Systems, Canada, in partnership with Hefei Gemei Culture and Education Technology Co. Ltd, Hefei, China, to Jiang Wu and RZ.
ReferencesBebbington, M., Lai, C.D., Zitikis, R.: Modelling deceleration in senescent mortality. 18, 18–37 (2011). MR2770696. doi:10.1080/08898480.2011.540173Bernardo, A.E., Ledoit, O.: Gain, loss, and asset pricing. 108, 144–172 (2000). doi:10.1086/262114Brazauskas, V., Jones, B.L., Zitikis, R.: Trends in disguise. 9, 58–71 (2015). doi:10.1017/S1748499514000232Chakravarty, S.R.: Extended Gini indices of inequality. 29, 147–156 (1988). MR0954119. doi:10.2307/2526814Davydov, Y., Zitikis, R.: An index of monotonicity and its estimation: a step beyond econometric applications of the Gini index. 63 (special issue in memory of Corrado Gini), 351–372 (2005). MR2276056Denuit, M., Dhaene, J., Goovaerts, M., Kaas, R.: . Wiley, New York (2005)Furman, E., Zitikis, R.: Weighted premium calculation principles. 42, 459–465 (2008). MR2392102. doi:10.1016/j.insmatheco.2007.10.006Furman, E., Zitikis, R.: Weighted pricing functionals with application to insurance: an overview. 13, 483–496 (2009). MR2595160. doi:10.1080/10920277.2009.10597570Gillen, B., Markowitz, H.M.: A taxonomy of utility functions. In: Aronson, J.R., Parmet, H.L., Thornton, R.J. (eds.) . Springer, New York (2009)Keating, C., Shadwick, W.F.: In: , The Finance Development Centre Limited. London, UK (2002)Lehmann, E.L.: Some concepts of dependence.. 37, 1137–1153 (1966). MR202228. doi:10.1214/1177699260.1966.10597570Levy, H.: . Springer, New York (2006). MR2239375Li, H., Li, X.: Stochastic Orders in Reliability and Risk. In: . Springer, New York (2013). MR3157470. doi:10.1007/978-1-4614-6892-9von Neuman, J., Morgenstern, O.: . Princeton University Press, Princeton, NJ (1944). MR0011937Qoyyimi, D.T., Zitikis, R.: Measuring the lack of monotonicity in functions. 39, 107–117 (2014). MR3307986Qoyyimi, D.T., Zitikis, R.: Measuring association via lack of co-monotonicity: the LOC index and a problem of educational assessment. 3, 83–97 (2015). MR3418658. doi:10.1515/demo-2015-0006Quiggin, J.: . Kluwer, Dordrecht (1993).Wang, S.S.: Insurance pricing and increased limits ratemaking by proportional hazards transforms. 17, 43–54 (1995). MR1363642. doi:10.1016/0167-6687(95)91054-PWang, S.S.: Premium calculation by transforming the layer premium density. 26, 71–92 (1996). doi:10.2143/AST.26.1.563234Yaari, M.E.: The dual theory of choice under risk. 55, 95–115 (1987). MR0875518. doi:10.2307/1911158Yitzhaki, S., Schechtman, E.: . Springer, New York (2013). MR3012052. doi:10.1007/978-1-4614-4720-7Zitikis, R.: Asymptotic estimation of the E-Gini index. 19, 587–601 (2003). MR1997934. doi:10.1017/S0266466603194042Zitikis, R., Gastwirth, J.L.: Asymptotic distribution of the S-Gini index. 44, 439–446 (2002). MR1934733. doi:10.1111/1467-842X.00245