Characterizations of the Logistic and Related Distributions

Chin-Yuan Hu and Gwo Dong Lin

arXiv:1803.06620v1 [math.ST] 18 Mar 2018

National Changhua University of Education and Academia Sinica

Abstract. It is known that few characterization results of the logistic distribution were
available before, although it is similar in shape to the normal one whose characteristic properties have been well investigated. Fortunately, in the last decade, several authors have made
great progress in this topic. Some interesting characterization results of the logistic distribution have been developed recently. In this paper, we further provide some new results by
the distributional equalities in terms of order statistics of the underlying distribution and
the random exponential shifts. The characterization of the closely related Pareto type II
distribution is also investigated.

AMS subject classiﬁcations: Primary 62E10, 62G30, 60E10.
Key words and phrases: Characterization, order statistics, stochastic order, logistic distribution, exponential distribution, Pareto type II distribution.
Short title: The logistic and related distributions
Postal addresses: Chin-Yuan Hu, Department of Business Education, National Changhua
University of Education, Changhua 50058, Taiwan. (E-mail: buhuua@gmail.com)
Gwo Dong Lin, Institute of Statistical Science, Academia Sinica, Taipei 11529, Taiwan.
(E-mail: gdlin@stat.sinica.edu.tw)

1

1. Introduction
The logistic distribution is similar to a normal distribution in shape (Mudholkar and
George 1978) and has an explicit closed form, so it has some advantages in practical applications. As remarked by Kotz (1974), few characterizations of the logistic distribution were
available before, but recently, some interesting results have been developed. In this paper,
we will further provide some more new results by properties of order statistics.
We ﬁrst introduce some notations. Let X obey the distribution F , denoted by X ∼ F.
Let {Xj }nj=1 be a random sample of size n from distribution F and denote the corresponding
order statistics by X1,n ≤ X2,n ≤ · · · ≤ Xn,n . The distribution function of Xk,n is denoted
by Fk,n . It is known that Fk,n is the composition of Bk,n−k+1 and F (see, e.g., Hwang and
Lin 1984), where Bα,β is the beta distribution with parameters α, β > 0, namely,
  Z F (x)
n
Fk,n (x) = Bk,n−k+1(F (x)) = k
tk−1 (1 − t)n−k dt, x ∈ R ≡ (−∞, ∞), (1)
k
0
Z
Γ(α + β) u α−1
Bα,β (u) =
t (1 − t)β−1 dt, u ∈ [0, 1].
(2)
Γ(α)Γ(β) 0
On the other hand, for Y ∼ G, we say that X is less than or equal to Y in the usual
stochastic order, denoted by X ≤st Y, if F ≤ G, where F (x) = 1 − F (x) = Pr(X > x).
Let us start with an interesting simple example. Clearly, for general distribution F , we
2

have X1,2 ≤st X because F1,2 = 1−F ≥ F by (1). One possibility to adjust this ”inequality”
is to choose a nonnegative random variable Z, independent of X and Xj ’s, such that
d

X = X1,2 + Z

(3)

or
d

X − Z = X1,2 ,

(4)

d

where = means equality in distribution. One might think that the solutions of the distributional equations (3) and (4) are the same, but this is not true in general, because
the characteristic function of Z is not equal to the reciprocal of that of −Z, namely,
E[eitZ ] 6= (E[e−itZ ])−1 , t ∈ R, in general. For example, if Z has the standard exponential
distribution E, then the solution of (3) is a logistic distribution F (x) = 1/[1+e−(x−µ) ], x ∈ R,
2

where µ ∈ R is a constant, while the solution of (4) is a negative (or reversed) exponential
distribution F (x) = e(x−µ)/2 , x ≤ µ, where µ ∈ R is a constant (see, e.g., George and Mudholkar 1982, Lin and Hu 2008, and Ahsanullah et al. 2011, and note that the smoothness
conditions on F therein are redundant due to Lemmas 1–3 below).
Throughout the paper, let U and ξ obey the uniform distribution U on [0, 1] and the
standard exponential distribution E, respectively. Moreover, let {Uj }nj=1 and {Uj′ }nj=1 be two
random samples of size n from U, and let {ξj }nj=1 and {ξj′ }nj=1 be two random samples of size
n from E. All the above random variables X, U, ξ, Xj , Uj , Uj′ , ξj and ξj′ , j = 1, 2, . . . , n, are
assumed to be independent from now on.
Mimicking the above characterization approaches (3) and (4), several authors have considered the general stochastic inequality Xk,n ≤st Xk+1,n and solved the distributional
d

d

equations (a) Xk+1,n = Xk,n + aξ and (b) Xk,n = Xk+1,n − bξ, or, more generally, (c)
d

Xk,n + aξ1 = Xk+1,n − bξ2 , where a and b are nonnegative constants. In particular, the
equality
Xk,n +

1
1
d
ξ1 = Xk+1,n − ξ2
n−k
k

also characterizes the logistic distribution. (See AlZaid and Ahsanulla 2003, Wesolowski and
Ahsanulla 2004, and Ahsanulla et al. 2010, 2012 for equations arising from Xk,n ≤st Xk+r,n
with 1 ≤ r ≤ n − k.) Besides, the distributional equations arising from (a) X ≤st Xn,n ,
(b) X1,n ≤st Xn,n and (c) Xm,m ≤st Xn,n , where m < n, were investigated by Zykov and
Nevzorov (2011), Ananjevskii and Nevzorov (2016) as well as Berred and Nevzorov (2013),
respectively.
In this paper we will solve the distributional equations arising from stochastic inequalities:
(i) Xk,n ≤st Xk,n−1 , (ii) Xk,n−1 ≤st Xk+1,n , (iii) Xk,n ≤st Xk,k , (iv) X1,k ≤st Xn−k+1,n , (v)
Xk,n ≤st Xk,n−m, and (vi) Xm−k,n−k ≤st Xm,n .
To do this, some useful lemmas are given in the next section. The main characterization
results are stated and proved in Sections 3 and 4. For simplicity, we ﬁrst deal with the closely
related Pareto type II distribution in Section 3, and then the logistic distribution in Section
4. Finally, we pose an open problem in Section 5.
2. Lemmas
3

We need some lemmas in the sequel. Lemma 1(i) was given without proof in Lukacs (1970,
p. 38), but has been ignored in the literature. We provide a proof here for completeness.
Lemma 1. (i) Let Y and Z be two independent random variables. If Y has an absolutely
continuous distribution, then so does Y + Z, regardless of the distribution of Z.
(ii) If, in addition to the assumptions in (i), both Y and Z are positive random variables,
then the product Y Z has an absolutely continuous distribution.
Proof. Let F , G and H be the distributions of Y + Z, Y and Z, respectively. Then
Z ∞
F (x) =
G(x − z)dH(z), x ∈ R.
−∞

Since G is absolutely continuous, we have that for each ε > 0, there exists a δ > 0 such
P
Pn
that nj=1 [G(bj ) − G(aj )] < ε if
j=1 (bj − aj ) < δ, where aj < bj ≤ aj+1 < bj+1 , j =
1, 2, . . . , n − 1. This in turn implies that for the above {(aj , bj )}nj=1 ,
n
X

[F (bj ) − F (aj )] =

Z

∞

n
X

[G(bj − z) − G(aj − z)]dH(z) <

−∞ j=1

j=1

Z

∞

ε dH(z) = ε.
−∞

Hence, part (i) is proved. To prove part (ii), we recall ﬁrst that both the logarithmic and
exponential functions are absolutely continuous, and that the composition preserves the
property of absolute continuity. Then consider log(Y Z) = log Y + log Z and use part (i)
to conclude that log(Y Z) has an absolutely continuous distribution, and hence so does Y Z.
This completes the proof.
It is known that the inverse function of an absolutely continuous function with positive
derivative almost everywhere is not necessarily absolutely continuous. However, we have the
following useful result.
Lemma 2. Let F be an absolutely continuous distribution on [0, 1] and F ′ (x) = f (x) > 0
on (0, 1). Then the inverse function of F is itself an absolutely continuous distribution.
Proof. By the assumptions, F is a strictly increasing and continuous function from [0, 1]
to [0, 1], so is its inverse function F −1 . This implies that F −1 is a continuous distribution on [0, 1]. Moreover,

d
F −1 (t)
dt

= 1/f (F −1(t)) is a positive measurable function on (0, 1)

(see, e.g., Shorack and Wellner 1986, pp. 8-9). By changing variables x = F −1 (t), we have
R 1 d −1
R1
[ F (t)]dt = 0 [1/f (x)] · f (x)dx = 1. Therefore, the distribution function F −1 has no
0 dt

singular part and is absolutely continuous. The proof is complete.
4

Lemma 3. If the distribution Fk,n of order statistic Xk,n is absolutely continuous, then so
is the underlying distribution F of X.
Proof. Recall (see (1)) that Fk,n (x) = Bk,n−k+1(F (x)), x ∈ R, where Bk,n−k+1, deﬁned in
(2), is the beta distribution Bα,β with parameters α = k and β = n − k + 1, and has a
−1
positive continuous density function on (0, 1). Therefore, F = Bk,n−k+1
◦ Fk,n is absolutely

continuous by Lemma 2. The proof is complete.
Lemma 4. Let 1 ≤k 
< n. Then we have the following identities:

n
n−1
n−k
k n−k
F F
(i) Fk,n − Fk+1,n =
F kF
, (ii) Fk,n − Fk,n−1 =
, and
k
k−1

n−k
n−1
(iii) Fk,n−1 − Fk+1,n =
F kF
.
k
Proof. For parts (i) and (ii), see, e.g., David and Shu (1978) as well as David and Nagaraja
(2003, p. 23), while part (iii) follows from parts (i) and (ii) and the identity:
  
 

n
n−1
n−1
=
+
.
k
k−1
k
Denote the left and the right extremities of F by ℓF and rF , respectively. It is known that
if the absolutely continuous F satisﬁes the functional equation F ′ (x) = F (x)(1 − F (x)), x ∈
(ℓF , rF ), then F is a logistic distribution. In the next lemma we extend this result.
Lemma 5. Let r, θ > 0, a ∈ [0, 1], and let F be an absolutely continuous distribution
function satisfying xa F ′ (x) = θF (x)(1 − F r (x)) for x ∈ (ℓF , rF ).
(i) If a = 1, then ℓF = 0, rF = +∞ and
1/r

λxrθ
, 0 ≤ x < ∞,
F (x) =
1 + λxrθ
where λ is a positive constant.
(ii) If a ∈ [0, 1), then ℓF = −∞, rF = +∞ and
F (x) =

rθ 1−a
λ exp( 1−a
x )
rθ 1−a
x )
1 + λ exp( 1−a

!1/r

, −∞ < x < ∞,

where λ is a positive constant.
Proof. Deﬁne the increasing function G(x) = (1 − F r (x))−1 − 1 from (ℓF , rF ) onto (0, ∞).
Then G′ (x) = rF r−1 (x)F ′ (x)(1 − F r (x))−2 , and hence
xa G′ (x) = rθG(x), x ∈ (ℓF , rF ).
5

(a) If a = 1, solving the above equation leads to G(x) = λxrθ , x ∈ (ℓF , rF ), for some constant
λ > 0. On the other hand, we have, by the deﬁnition of G, that

1/r 
1/r
G(x)
λxrθ
F (x) =
=
, x ∈ (ℓF , rF ),
1 + G(x)
1 + λxrθ
and hence, ℓF = 0 and rF = +∞, because F is a distribution function. This proves part (i).
rθ 1−a
x ), x ∈ (ℓF , rF ), for some constant
(b) If a ∈ [0, 1), we have instead G(x) = λ exp( 1−a

λ > 0. The required result then follows from both the deﬁnition of the function G and the
fact that F is a distribution function. The proof is complete.
Some equalities (in distribution) of the next lemma are essentially due to Nevzorov (2001,
Lecture 3), but we provide here an alternative and possibly simpler proof.
Lemma 6. Let ξk,n (Uk,n , resp.) be the k-th smallest order statistic of a random sample
of size n from the standard exponential distribution E (the uniform distribution U, resp.).
Then the following statements are true.
n−k+1
(i) The Laplace transform of ξk,n is Lξk,n (s) = n−k+1+s
·
d P
1
ξk−j+1, where 1 ≤ k ≤ n.
(ii) ξk,n = kj=1 n−j+1

n−k+2
n−k+2+s

n
· · · n+s
, s ≥ 0.

d

′
(iii) ξm,n = ξm−k,n−k + ξk,n
, where 1 ≤ k < m ≤ n.
d Q
1/(n−j+1)
(iv) Un−k+1,n = kj=1 Uj
, where 1 ≤ k ≤ n.
d

′
(v) Uk,n = Uk,m−1 · Um,n
, where 1 ≤ k < m ≤ n.

 
n R E(x) k−1
Proof. Recall that the distribution of ξk,n is Fξk,n (x) = k
t (1 − t)n−k dt, x ≥ 0,
k 0
where E(x) = Fξ (x) = 1 − e−x , x ≥ 0. Then the Laplace transform of ξk,n is
 Z ∞
Z ∞
n
−sξk,n
−sx
Lξk,n (s) = E[e
]=
e dFξk,n (x) = k
e−(n−k+1+s)x (1 − e−x )k−1 dx, s ≥ 0.
k
0
0
By integration by parts, it follows from the above integral that
 
Z ∞
k−1
n
Lξk,n (s) = k
e−(n−k+2+s)x (1 − e−x )k−2 dx = · · ·
k n−k+1+s 0
 
Z ∞
k−2
1
k−1
n
·
···
·
e−(n+s)x dx
= k
k n−k+1+s n−k+2+s
n−1+s 0
n−k+2
n−1
n
n−k+1
·
···
·
, s ≥ 0.
=
n−k+1+s n−k+2+s
n−1+s n+s
This proves part (i), which in turn implies parts (ii) and (iii) by the fact that E[e−s(ξ/k) ] =
d

k/(k + s), s ≥ 0. Part (iv) follows from part (ii) because Uj = exp(−ξk−j+1 ) and the order
6

d

d

′
statistic Un−k+1,n = exp(−ξk,n ). To prove part (v), we have Un−m+1,n = Un−m+1,n−k · Un−k+1,n

by using part (iii), and then reset k = n − m + 1. The proof is complete.
Lemma 7. Let Y obey the Pareto type II (or log-logistic) distribution G(y) = y/(1+y), y ≥
0. Let {Yj }nj=1 , independent of U and {Uj }nj=1 , be a random sample of size n from G. Then
we have the following equalities in distribution:
d

d

(i) 1/Y = Y and in general, 1/Yk,n = Yn−k+1,n , where 1 ≤ k ≤ n.
d

(ii) Yk,n−1 = Yk,n /U 1/(n−k) , where 1 ≤ k ≤ n − 1.
d

(iii) Yk,n−m = Yk,n/Un−k−m+1,n−k , where 1 ≤ k ≤ n − m.
d

(iv) Yk,n−1 = Yk+1,n · U 1/k , where 1 ≤ k ≤ n − 1.
d

(v) Ym−k,n−k = Ym,n · Um−k,m−1 , where 2 ≤ k + 1 ≤ m ≤ n.
Proof. It is easy to check part
 To prove the remaining parts, recall that the distri (i).
n R G(y) k−1
bution of Yk,n is Gk,n (y) = k
t (1 − t)n−k dt, y ≥ 0. Then we have H(y) ≡
k 0
R1
Pr(Yk,n /U 1/(n−k) ≤ y) = 0 Gk,n (yu1/(n−k))du. By changing variables,
Z 1
Z 1
n−k
H(y) =
Gk,n (yt)dt
= (n − k)
Gk,n (yt)tn−k−1dt, y ≥ 0.
0

0

Therefore, Gk,n−1 = H iﬀ, by diﬀerentiation,
n Z 1 
n+1

1
1
=
tn−1 dt, y ≥ 0,
n
1+y
1
+
yt
0
which is, however, a special case of the well-known identity

β1
β1 +β2
Z 
1
Γ(β1 + β2 ) 1
1
=
tβ1 −1 (1 − t)β2 −1 dt, y ≥ 0
1+y
Γ(β1 )Γ(β2 ) 0 1 + yt
(see Gradshteyn and Ryzhik 2014, p. 314). This proves part (ii). Part (iii) follows from part
(ii) by iteration and Lemma 6(iv), while part (iv) follows from either parts (i) and (ii) (letting
d

d

k∗ = n − k) or Lemma 8(iii) below because Y = exp(X) and U = exp(−ξ). Finally, we prove
part (v) by using part (iv), iteration and Lemma 6(iv) again. The proof is complete.
Lemma 8. Let X obey the standard logistic distribution F (x) = 1/[1 + exp(−x)], x ∈ R.
Then we have the following equalities in distribution:
d

(i) Xk,n−1 = Xk,n +

1
ξ,
n−k

where 1 ≤ k ≤ n − 1.

d

(ii) Xk,n−m = Xk,n + ξm,n−k , where 1 ≤ k ≤ n − m.
d

(iii) Xk,n−1 = Xk+1,n − k1 ξ, where 1 ≤ k ≤ n − 1.
7

d

(iv) Xm−k,n−k = Xm,n − ξk,m−1 , where 2 ≤ k + 1 ≤ m ≤ n.
d

d

Proof. The results follow from Lemma 7 by noting that (a) X = log Y, (b) ξ = − log U and
d

(c) − log Uk,n = ξn−k+1,n for all 1 ≤ k ≤ n. The proof is complete.
For the proof of the next lemma, see Lin and Hu (2008, Lemma 5).
Lemma 9. Let f and g be two functions real analytic and strictly monotone in [0, ∞).
Assume that for each n ≥ 1, the n-th derivatives f (n) and g (n) are strictly monotone in some
interval [0, δn ). Let {xn }∞
n=1 be a sequence of positive real numbers converging to zero. If
f (xn ) = g(xn ), n = 1, 2, . . . , then f = g.
3. Characterizations of the Pareto type II distribution
We start with the Pareto type II distribution which is easier to handle, and recall that
the uniform order statistic Uk,n ∼ Bk,n−k+1.
Theorem 1. Let Y ∼ G be a positive random variable and let 1 ≤ k ≤ n − 1 be ﬁxed
integers. Let Y1 , Y2 , . . . , Yn be n independent copies of Y, and let U, independent of {Yi }ni=1 ,
be a random variable with uniform distribution on [0, 1]. Then the distributional equality
d

Yk,n−1 = Yk,n /U 1/(n−k)

(5)

holds iﬀ G is a Pareto distribution G(y) = λy/(1 + λy), y ≥ 0, where λ > 0 is a constant.
d

Proof. The suﬃciency part follows from Lemma 7(ii) because (λY )ℓ,m = λYℓ,m for λ > 0 and
for all 1 ≤ ℓ ≤ m. To prove the necessity part, we note, by Lemma 1(ii), that the distribution
Gk,n−1 of Yk,n−1 is absolutely continuous, and so is G by Lemmas 2 and 3. Rewrite (5) as
Z 1
Gk,n−1(y) =
Gk,n (yu1/(n−k) )du, y ≥ 0.
0

By changing variables t = yu1/(n−k), we have
Gk,n−1(y) = (n − k)y

−(n−k)

Z

y

Gk,n (t)tn−k−1 dt, y > 0.
0

Taking diﬀerentiation leads to
yG′k,n−1(y) = (n − k)[Gk,n (y) − Gk,n−1(y)], y > 0.
With the help of (1) and Lemma 4(ii), (6) is equivalent to
yG′(y) = G(y)[1 − G(y)], y ∈ (ℓG , rG ).
8

(6)

Finally, Lemma 5(i) with r = θ = 1 completes the proof.
Corollary 1.

Under the same assumptions of Theorem 1, the distributional equality

d

Yk,n−1 = Yk,n/U 1/α holds for some α > 0, iﬀ G is a Pareto distribution with G(y) =
1/[1 + λy α/(n−k) ], y ≥ 0, where λ is a positive constant.
α/(n−k)

Proof. Consider Yi′ = Yi

α/(n−k)

′
for i = 1, 2, . . . . Then Yk,n
= Yk,n

, and apply Theorem

d

′
′
1 to the case: Yk,n−1
= Yk,n
/U 1/(n−k) .

Corollary 2. Under the same assumptions of Theorem 1, the distributional equality
d

Yk,n−1 = Yk+1,n · U 1/k
holds iﬀ G is a Pareto distribution G(y) = λy/(1 + λy), y ≥ 0, where λ > 0 is a constant.
Proof. The suﬃciency part is a consequence of Lemma 7(iv). To prove the necessity
d

∗
part, denote Y ∗ = 1/Y. Then Yℓ,m
= (1/Y )ℓ,m = 1/Ym−ℓ+1,m for all 1 ≤ ℓ ≤ m. By
d

d

∗
assumptions, we have the equality 1/Yk,n−1 = 1/Yk+1,n · 1/U 1/k , or, equivalently, Yn−k,n−1
=
∗
Yn−k,n
· 1/U 1/(n−(n−k)) . It follows from Theorem 1 (letting k∗ = n − k) that Y ∗ has a Pareto

distribution G∗ (y) = λ∗ y/(1 + λ∗ y), y ≥ 0, for some constant λ∗ > 0. This in turn implies
that Y has the Pareto distribution G(y) = λy/(1 + λy), y ≥ 0, where λ = 1/λ∗. We can
d

prove the last claim directly, or by using Lemma 7(i), because λ∗ Y ∗ = 1/(λ∗ Y ∗ ) = Y /λ∗
having the standard Pareto type II distribution. The proof is complete.
Theorem 2. Let Y ∼ G be a positive random variable and let 1 ≤ k ≤ n − 1 be ﬁxed
integers. Let Y1 , Y2 , . . . , Yn be n independent copies of Y, and let B, independent of {Yi }ni=1 ,
be a random variable having beta distribution Bα,β with parameters α = 1 and β = n − k,
that is, FB (u) = 1 − (1 − u)n−k , u ∈ [0, 1]. Assume further that limy→0+ G(y)/y = λ > 0.
Then the distributional equality
d

Yk,k = Yk,n /B

(7)

holds iﬀ G is the Pareto distribution G(y) = λy/(1 + λy), y ≥ 0.
d

Proof. The suﬃciency part follows from Lemma 7(iii) with n − m = k and the fact B =
U1,n−k . To prove the necessity part, we note ﬁrst that G is absolutely continuous by Lemmas
1–3, and then rewrite (7) as the functional equation:
Z 1 Z G(yu)  
n k−1
k
G (y) =
k
t (1 − t)n−k dtdFB (u), y ≥ 0.
k
0
0
9

(8)

Now, it suﬃces to prove that the solution of equation (7) is unique under the smoothness
condition on the distribution. Namely, if the absolutely continuous distribution F on (0, ∞)
satisﬁes limy→0+ F (y)/y = λ > 0 and
k

F (y) =

Z

1
0

Z

F (yu)

0

 
n k−1
t (1 − t)n−k dtdFB (u), y ≥ 0,
k
k

then we will prove that F = G. From (8) and (9) it follows that
Z 1
1
k
k
|F (y) − G (y)| ≤
|F k (yu) − Gk (yu)|dFB (u), y ≥ 0,
E[B k ] 0
 
n
k
where E[B ] = 1/
. Deﬁne the bounded function
k
g(y) =

(9)

(10)

F k (y) − Gk (y)
, y > 0, and g(0) = lim+ g(y) = 0,
y→0
yk

and the increasing function
h(y) = sup g(t), y > 0, and h(0) = lim+ h(y) = 0.
y→0

0≤t≤y

By (10), we see that
g(y) ≤

Z

1

g(uy)dH(u), y ≥ 0,

(11)

0

where H(u) = (1/E[B k ])

Ru
0

tk dFB (t), u ∈ [0, 1]. Now, by (11) and the deﬁnition of the

increasing function h, we have
Z 1
Z
h(y) ≤
h(uy)dH(u) ≤ h(y)
0

1

dH(u) = h(y), y ≥ 0.

0

This in turn implies that h is a constant function and hence h(y) = 0, y ≥ 0, because
h(0) = h(0+ ) = 0. Consequently, g(y) = 0, y ≥ 0, and F = G. The proof is complete.
The next result is the counterpart of Theorem 2 for the minimum order statistics.
Corollary 3. Under the same setting in Theorem 2 with the condition on G replaced by
limy→0+ G(1/y)/y = 1/λ > 0 (equivalently, limy→+∞ yG(y) = 1/λ > 0), the distributional
d

equality Y1,k = Yn−k+1,n · B holds iﬀ G is the Pareto distribution G(y) = λy/(1 + λy), y ≥ 0.
d

Proof. Use Lemma 7(v), Theorem 2 and the fact (1/Y )ℓ,m = 1/Ym−ℓ+1,m for all 1 ≤ ℓ ≤ m.
We now further extend Theorem 2 under some stronger smoothness conditions.
10

Theorem 3. Let Y ∼ G be a positive random variable and let n, m, k be three ﬁxed positive
integers with 1 ≤ k ≤ n − m. Let Y1 , Y2 , . . . , Yn be n independent copies of Y, and let B1 ,
independent of {Yi }ni=1 , be a random variable having beta distribution Bα,β with parameters
α = n − m − k + 1 and β = m. Assume further that the distribution function G satisﬁes the
following conditions:
(i) G is real analytic and strictly increasing in [0, ∞) and for each i ≥ 1, its i-th derivative
G(i) is strictly monotone in some interval [0, δi ).
(ii) limy→0+ [Gk (y) − (λy)k ]/(λy)k+1 = −k for some positive constant λ.
Then the distributional equality
d

Yk,n−m = Yk,n /B1

(12)

holds iﬀ G is the Pareto distribution G(y) = λy/(1 + λy), y ≥ 0.
d

Proof. The suﬃciency part follows from Lemma 7(iii) and the fact B1 = Un−m−k+1,n−k . To
prove the necessity part, we note ﬁrst that G is absolutely continuous as before, and then
we rewrite (12) as the functional equation:
Z 1
Gk,n−m(y) =
Gk,n (yu)dFB1 (u), y ≥ 0.

(13)

0

Now, it suﬃces to prove that the solution of equation (12) is unique under the smoothness
condition on the distribution. Namely, if the absolutely continuous distribution F on (0, ∞)
satisﬁes the above conditions (i) and (ii) and
Z 1
Fk,n−m(y) =
Fk,n (yu)dFB1 (u), y ≥ 0,

(14)

0

then we will prove that F = G.
Deﬁne the increasing function H(y) = max{F (y), G(y)}, y ≥ 0. From (1) it follows that
for any a > 0 and 0 ≤ y ≤ a,

 Z F (y)
n−m
|Fk,n−m (y) − Gk,n−m (y)| = k
tk−1 (1 − t)n−m−k dt
k
G(y)


1
n−m
≥ k
(1 − H(y))n−m−k F k (y) − Gk (y)
k
k


n−m
≥
(1 − H(a))n−m−k F k (y) − Gk (y) .
k
11

(15)

On the other hand, we have
 
n
F k (y) − Gk (y) , y ≥ 0.
|Fk,n (y) − Gk,n (y)| ≤
k
Combing (12)–(16) leads to


n−m
(1 − H(a))n−m−k F k (y) − Gk (y) ≤ |Fk,n−m(y) − Gk,n−m(y)|
k
 Z 1
Z 1
n
≤
|Fk,n (yu) − Gk,n (yu)| dFB1 (u) ≤
F k (yu) − Gk (yu) dFB1 (u).
k
0
0

(16)

(17)

Deﬁne the bounded increasing function
h(y) = sup
0<t≤y

F k (t) − Gk (t)
, y > 0,
tk+1

and h(0) = lim+ h(y) = 0.
y→0

Then from the inequality (17) it follows that for any a > 0,


 Z 1
n−m
n
n−m−k
(1 − H(a))
h(y) ≤
h(yu)uk+1dFB1 (u)
k
k
0
 
 
Z 1
n
n
h(y)
≤
h(y)E[B1k+1], 0 ≤ y ≤ a.
uk+1dFB1 (u) =
k
k
0

  
n−m
n
Recall that E[B1k ] =
/
. Then rewrite the inequality (18) as follows:
k
k
E[B1k ](1 − H(a))n−m−k h(y) ≤ h(y)E[B1k+1], 0 ≤ y ≤ a, a > 0.

(18)

(19)

We claim that there exists a y0 > 0 such that h(y0 ) = 0. Otherwise, we have, by (19),
E[B1k ](1 − H(a))n−m−k ≤ E[B1k+1 ], ∀ a > 0,
which in turn implies, by letting a → 0+ , that E[B1k ] ≤ E[B1k+1 ], a contradiction. Therefore,
h(y0 ) = 0 for some y0 > 0 and hence F (y) = G(y) for y ∈ [0, y0]. By Lemma 9 and the
assumptions on F and G, we conclude that F = G. The proof is complete.
Corollary 4. Let Y ∼ G be a positive random variable and let n, m, k be three ﬁxed positive
integers with k + 1 ≤ m ≤ n. Let Y1 , Y2 , . . . , Yn be n independent copies of Y, and let B2 ,
independent of {Yi }ni=1 , be a random variable having beta distribution Bα,β with parameters
α = m − k and β = k. Assume further that the distribution function G∗ of 1/Y satisﬁes the
following conditions:
12

(i) G∗ is real analytic and strictly increasing in [0, ∞) and for each i ≥ 1, its i-th derivative
(i)

G∗ is strictly monotone in some interval [0, δi ).
(ii) limy→0+ [Gk∗∗ (y) − (y/λ)k∗ ]/(y/λ)k∗ +1 = −k∗ for some positive constant λ, where k∗ =
n − m + 1.
d

Then the distributional equality Ym−k,n−k = Ym,n · B2 holds iﬀ G is the Pareto distribution
G(y) = λy/(1 + λy), y ≥ 0.
d

Proof. Use Lemma 7(v), Theorem 3 and the fact (1/Y )ℓ,m = 1/Ym−ℓ+1,m for all 1 ≤ ℓ ≤ m.
In summary, for a positive random variable Y ∼ G, we have the following characteristic
properties of the Pareto distribution G(y) = λy/(1 + λy), y ≥ 0, where λ is a positive
constant (compare with Lemma 7).
d

1. Yk,n−1 = Yk,n/U 1/(n−k) .
d

2. Yk,n−1 = Yk+1,n · U 1/k .
d

d

3. Yk,k = Yk,n /B (equivalently, Yk,k = Yk,n /U1,n−k ).
d

d

4. Y1,k = Yn−k+1,n · B (equivalently, Y1,k = Yn−k+1,n · U1,n−k ).
d

d

5. Yk,n−m = Yk,n /B1 (equivalently, Yk,n−m = Yk,n /Un−m−k+1,n−k ).
d

d

6. Ym−k,n−k = Ym,n · B2 (equivalently, Ym−k,n−k = Ym,n · Um−k,m−1 ).
Here, the random variables U ∼ U, B ∼ B1,n−k , B1 ∼ Bn−m−k+1,m , B2 ∼ Bm−k,k , and on
the RHS of each equality, the two random variables are independent.
4. Characterizations of the logistic distribution
We are now ready to provide characterization results of the logistic distribution.
Theorem 4. Let X ∼ F and let 1 ≤ k ≤ n − 1 be ﬁxed integers. Then the distributional
equality
d

Xk,n−1 = Xk,n +

1
ξ
n−k

holds iﬀ F is a logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R, where µ ∈ R is a constant.
Proof. Take Yi = exp(Xi ), U = exp(−ξ) and λ = e−µ . Then the result follows immediately
from Theorem 1.
Corollary 5. Let X ∼ F, α > 0 and let 1 ≤ k ≤ n − 1 be ﬁxed integers. Then the
d

distributional equality Xk,n−1 = Xk,n + α1 ξ holds iﬀ F is a logistic distribution F (x) =
1/{1 + e−[α/(n−k)](x−µ) }, x ∈ R, where µ ∈ R is a constant.
13

Corollary 6. Let X ∼ F, α > 0 and let 1 ≤ k ≤ n − 1 be ﬁxed integers. Then the
distributional equality
d

Xk,n−1 = Xk+1,n −

1
ξ
α

(20)

holds iﬀ F is a logistic distribution F (x) = 1/[1 + e−(α/k)(x−µ) ], x ∈ R, where µ ∈ R is a
constant.
d

Proof. Use Corollary 5 and the fact that Xℓ,m = −(−X)m−ℓ+1,m for all 1 ≤ ℓ ≤ m.
d

The counterpart of (20), namely, Xk+1,n = Xk,n−1 + α1 ξ, and the two-sided case: Xk,n +
d

aξ1 = Xk,n−1 + bξ2 (see Corollary 5), where α, a, b > 0, were investigated by Wesolowski and
Ahsanullah (2004). All the solutions of these two equations are exponential distributions.
The next result improves and extends Theorem 6 of Lin and Hu (2008) by an approach
diﬀerent from the previous method of intensively monotone operator (Kakosyan et al. 1984).
Theorem 5.

Let 1 ≤ k ≤ n − 1 be ﬁxed integers.

Assume that X ∼ F satisﬁes

limx→−∞ F (x)/ex = e−µ for some constant µ ∈ R. Then the distributional equality
d

Xk,k = Xk,n + ξn−k,n−k
holds iﬀ F is the logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R.
Proof. The suﬃciency part follows from Lemma 8(ii), while the necessity part is a consequence of Theorem 2.
The next result is the counterpart of Theorem 5 for the minimum order statistics.
Corollary 7. Let 1 ≤ k ≤ n − 1 be ﬁxed integers.

Assume that X ∼ F satisﬁes

limx→−∞ F (−x)/ex = eµ for some constant µ ∈ R. Then the distributional equality
d

X1,k = Xn−k+1,n − ξn−k,n−k
holds iﬀ F is the logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R.
d

Proof. Use Lemma 8(iv), Theorem 5 and the fact Xℓ,m = −(−X)m−ℓ+1,m for all 1 ≤ ℓ ≤ m.
Using Theorem 3 and Lemma 8(ii), we further extend Theorem 5 to the following.
Theorem 6. Let n, m, k be three ﬁxed positive integers with 1 ≤ k ≤ n − m and let X ∼ F
satisfy limx→−∞ [e−k(x−µ) F k (x) − 1]/ex−µ = −k for some constant µ ∈ R. Assume further
that the distribution function G of exp(X1 ) is real analytic and strictly increasing in [0, ∞)
14

and that for each i ≥ 1, its i-th derivative G(i) is strictly monotone in some interval [0, δi ).
Then the distributional equality
d

Xk,n−m = Xk,n + ξm,n−k
holds iﬀ F is the logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R.
As before, Theorem 6 and Lemma 8(iv) together lead to the following.
Corollary 8. Let n, m, k be three ﬁxed positive integers with k + 1 ≤ m ≤ n and let
X ∼ F satisfy limx→−∞ [e−k∗ (x+µ) (F (−x))k∗ − 1]/ex+µ = −k∗ for some constant µ ∈ R, where
k∗ = n−m+1. Assume further that the distribution function G∗ of exp(−X1 ) is real analytic
(i)

and strictly increasing in [0, ∞) and that for each i ≥ 1, its i-th derivative G∗ is strictly
monotone in some interval [0, δi ). Then the distributional equality
d

Xm−k,n−k = Xm,n − ξk,m−1

(21)

holds iﬀ F is the logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R.
In summary, for a random variable X ∼ F, we have the following characteristic properties
of the logistic distribution F (x) = 1/[1 + e−(x−µ) ], x ∈ R (compare with Lemma 8). Here,
µ ∈ R, ξ ∼ E, and on the RHS of each equality, the two random variables are independent.
d

1. Xk,n−1 = Xk,n +

1
ξ.
n−k

d

2. Xk,n−1 = Xk+1,n − k1 ξ.
d

3. Xk,k = Xk,n + ξn−k,n−k .
d

4. X1,k = Xn−k+1,n − ξn−k,n−k .
d

5. Xk,n−m = Xk,n + ξm,n−k .
d

6. Xm−k,n−k = Xm,n − ξk,m−1 .
5. Open problem
Finally, we would like to pose the following problem in which part (i) is the counterpart
of (21) for exponential distribution, and in part (ii), the ﬁrst two cases, m = k + 1, k + 2,
have been solved by AlZaid and Ahsanullah (2003) and Ahsanullah et al. (2010).
Problem. Let X ∼ F and let 1 ≤ k < m ≤ n be ﬁxed integers. Then solve the general
d

d

distributional equations: (i) Xm,n = Xm−k,n−k + ξk,n and (ii) Xm,n = Xk,n + ξm−k,n−k .
15

Acknowledgments. The authors would like to thank the Editor-in-Chief and the Referee
for helpful comments and constructive suggestions which improve the presentation of the
paper.
References
Ahsanullah, M., Berred, A. and Nevzorov, V.B. (2011). On characterizations of the exponential distributions. J. Appl. Statist. Sci., 19, 37–43.
Ahsanullah, M., Nevzorov, V.B. and Yanev, G.P. (2010). Characterizations of distributions
via order statistics with random exponential shifts. J. Appl. Statist. Sci., 18, 297–305.
Ahsanullah, M., Yanev, G.P. and Onica, C. (2012). Characterizations of logistic distribution through order statistics with independent exponential shifts. Economic Quality
Control, 27, 85–96.
AlZaid, A.A. and Ahsanullah, M. (2003). A characterization of the Gumbel distribution
based on record values. Commun. Statist. - Theory and Methods, 32, 2101–2108.
Ananjevskii, S.M. and Nevzorov, V.B. (2016). On families of distributions characterized
by certain properties of ordered random variables. Vestnik St. Petersburg University:
Mathematics, 49, 197–203.
Berred, A. and Nevzorov, V.B. (2013). Characterizations of families of distributions, which
include the logistic one, by properties of order statistics. Journal of Mathematical
Sciences, 188, 673–676.
David, H.A. and Nagaraja, H.N. (2003). Order Statistics, 3rd edn. Wiley, New Jersey.
David, H.A. and Shu, V.S. (1978). Robustness of location estimators in presence of an
outlier. In: Contributions to Survey Sampling and Applied Statistics, ed. H.A. David,
pp. 235–250. Academic Press, New York.
George, E.O. and Mudholkar, G.S. (1982). On the logistic and exponential laws. Sankhyā,
Ser. A, 44, 291–293.
16

Gradshteyn, I.S. and Ryzhik, I.M. (2014). Table of Integrals, Series, and Products, 8th edn.
Elsevier, New York.
Hwang, J.S. and Lin, G.D. (1984). Characterizations of distributions by linear combinations
of moments of order statistics. Bulletin of the Institute of Mathematics, Academia
Sinica, 12, 179–202.
Kakosyan, A.V., Klebanov, L.B. and Melamed, J.A. (1984). Characterization of Distributions by the Method of Intensively Monotone Operators. Springer, New York.
Kotz, S. (1974). Characterizations of statistical distributions: a supplement to recent
surveys. Int. Statist. Rev., 42, 39–65.
Lin, G.D. and Hu, C.-Y. (2008). On characterizations of the logistic distribution. J. Statist.
Plann. Infer., 138, 1147–1156.
Lukacs, E. (1970). Characteristic Functions, 2nd edn. Hafner Pub. Co., New York.
Mudholkar, G.S. and George, E.O. (1978). A remark on the shape of the logistic distribution. Biometrika, 65, 667–668.
Nevzorov, V.B. (2001). Records: Mathematical Theory. Translations of Mathematical
Monographs, Vol. 194, Amer. Math. Soc., Rhode Island.
Shorack, G.R. and Wellner, J.A. (1986). Empirical Processes with Applications to Statistics.
Wiley, New York.
Wesolowski, J. and Ahsanullah, M. (2004). Switching order statistics through random
power contractions. Aust. N.Z. J. Statist., 46, 297–303.
Zykov, V.O. and Nevzorov, V.B. (2011). Some characterizations of families of distributions,
including logistic or exponential ones, by properties of order statistics. Journal of
Mathematical Sciences, 176, 203–206.

17

