Projecting Statistical Functionals
Projecting Statistical Functionals
, Springer
Tomasz Rychlik
Institute of Mathematics
Polish Academy of Sciences
Chopina 12
87100 Toruil
Poland
trychlik@[Link]
9 8 7 654 3 2 1
Tomasz Rychlik
January, 2001
Contents
Preface v
List of Tables ix
2 Basic Notions 11
2.1 Elements of Hilbert Space Theory 11
2.2 Statistical Linear Functionals . . . 16
2.3 Restricted Families of Distributions. 25
3 Quantiles 33
3.1 General and Symmetric Distributions . . . . . . . . . . . . 33
3.2 Distributions with Monotone Density and Failure Rate . .. 36
3.3 Distributions with Monotone Density and Failure Rate on
the Average . . . . . . . . . . . . . 44
3.4 Symmetric Unimodal Distributions 50
3.5 Open Problems . . . . . . . . . . . 54
1.1 Introduction
This work presents a method of using projections of functions onto convex
cones in Hilbert spaces for determining sharp bounds on values of statistical
functionals over general and restricted families of distributions, expressed in
terms of moment parameters of the distributions. The method is based on
representing the statistical functionals and families of distributions as fixed
elements and convex cones, respectively, in a common real Hilbert space.
Then the norm of the projection of the element onto the cone provides the
optimal bound. The distribution for which the bound is attained is derived
by a simple transformation of the projection.
The advantage of the projection method lies in providing definite an-
swers for numerous simply stated but nontrivial problems of theoretical
and practical importance. The method enables us to optimally evaluate
various random objects in terms of the first two moments which are the
most classical parameters describing the population. Natural restrictions on
the random structure of observations are allowed. Although some results
presented here were proven earlier by means of specific tools, our unified
approach provides simpler proofs and indicates mutual relations among
various evaluations. On the other hand, for numerous problems presented
in the book, we do not see alternative ways of solving them without refer-
ence to projections. The idea of projections has been effectively exploited
in various aspects of statistical inference, for instance in the least squares
and minimum distance estimation. The majority of applications are based
2 1. Introduction and Notation
where
/;:n(X) = n(; =:)xi- 1 (1- x)n- j
is the density function of the jth order statistic from the LLd. standard uni-
form sample. This may be interpreted as the inner product of the centered
quantile function F-l_P.F with /;:n -1 in the Hilbert space £2([0, 1), dx) of
the square integrable functions on the unit interval. Applying the Schwarz
inequality to (1.1), and noting that IIF-l - [Link] = UF, we obtain
(1.2)
The equality is attained here if the arguments of the inner product are
proportional; that is,
F- 1 (x) - P.F = a[/;:n(x) -1] (1.3)
for some a ~ O. All the possible centered quantile functions of distributions
with a finite variance form the convex cone of nondecreasing functions
integrating to 0 in £2([0, 1),dx). The cone is further denoted by Co.
rr j = n, then the right-hand side of (1.3) actually increases, and inte-
grates to O. Therefore bound (1.2) is then tight, and (1.3) enables us to
determine a (power) distribution that attains the bound (see Gumbel [36],
and Hartley and David [38]). Otherwise /;:n -1 is not nondecreasing, and
(1.2) cannot be sharp.
The most natural idea is to replace /;:n - 1 by the closest element of
the family of centered quantile functions Co. We call it the projection of
/;:n - 1 onto CO and denote it by pO(/j:n - 1). The idea does work here
and in many other problems considered in the monograph. Namely,
[EFXj:n - P.F]/UF $ B°(j,n) = IlpO(fj:n -1)11 (1.4)
for all F with u} < 00, and we get the equality in (1.4) if
F-l(X) - P.F = P°(fj:n -1)(x) E CO,
(1.5)
UF IIP°(fj:n - 1)11
1.1 Introduction 3
1.2 Notation
We tried to set a unified notation for the whole work. A list of symbols is
presented below. Only notation used locally is not included there. We gen-
erally preserved standard symbols. However, some other notation may look
strange at the first glance. For instance, numerous convex cones of nmc-
tions were considered in the text. They were denoted by C with various
superscripts and subscripts whose meaning is explained in the list below
(see also Section 2.3). The projection operator onto a given cone is written
as P with the same upper and lower indices. Moreover, the indices appear
in the notation of sharp bounds on functionals that are obtained by means
of projections on specified convex cones. For instance, B°(j, n) calculated
in Example 1, denotes the sharp mean-variance bound on the expectation
of the jth order statistic of an independent identically distributed sample
of size n with arbitrary marginal distribution that has a finite variance. The
bound is obtained by means of projection pO of a properly chosen functional
(see Section 2.2) onto the convex cone Co of quantile functions of all distri-
butions with finite variance, centered about the respective mean. Convex
cones of quantile functions for families of distributions satisfying various
restrictions and respective projection operators need more sophisticated
notation. Nevertheless, we tried to introduce the symbols in a coherent
manner and believe that the reader shall get used to them.
List of symbols
(1l, (. , .)) (real) Hilbert space with inner product (.,.)
11·11 = (., .)1/2 - norm in (1l, (. ,.))
L2([a, d), w(x) dx) Hilbert space of functions 9 : [a, d) I-t !R
t
satisfying g2(X)W(X) dx < 00 with inner
a d
product (g, h) = fa g(x)h(x)w(x) dx
for a positive weight nmction w(x)
C¥.,(3, ... scalars, parameters of functions
C¥.*,(3., ••. optimal scalars, parameters of projections
g,h, .. . functions, elements of Hilbert spaces
G,H, .. . antiderivatives of g, h, .. . , respectively
gH(x) = g(H(x)) -composition of functions H andg
H(x) greatest convex minorant of H(x)
h(x) (right) derivative of H(x)
l(x) constant function equal to 1
lA(X) indicator of A (= 1 on A and 0 elsewhere)
x+ = max{x,O} - positive part of a number
lxJ = max{k :5 x: kinteger}-floorofa number
6 l. Introduction and Notation
fi.W{P) = Ewmin{X,p}
v?v{P) = Ew[min{X,p}]2
u?v{P) = Varwmin{X,p}
11w{a,p) = Ew{X - a) 1 (ti,d) (X)
1?~{a,p) = Varw{X - a) 1 (ti,d) (X)
f/w{P) = 11w{P, P) = Ew{X - P)+
v?v{P) = Ew[{X _ P)+]2
J~{P) = 1?~{P,P) = Varw{X - P)+
'iJw{P) = Ew{P-X)+
t?~{P) = Varw{P - X)+
XI, ... ,Xn. ... independent identically distributed (LLd.)
random variables with distribution function F
YI, ... ,Yn, ... possibly dependent identically distributed
random variables with common marginal F
X j :n jth order statistic of Xl, .. . , X n, 1 :5 j :5 n
Y;:n jth order statistic of YI , .. . , Yn, 1 :5 j :5 n
Ej=l CjXj :n L-statistic of independent sample
Ej=l Cj Y;:n L-statistic of dependent sample
C = (CI, ... ,en) vector of coefficients of L-statistic
Ln nth occurrence time of (first) record (increase
in sequence of sample maxima Xj:j, j ~ 1)
Rn = XLn - nth value of (first) record,
L~) nth occurrence time of (kth) record (increase
in sequence of kth greatest order statistics
X j +1-k:j, j ~ k), k ~ 1
~) = XL~.)+1-k:L~.) - nth value of kth record
/;:n{x) = n(j=Dxj - I {I- x)n- j l[O,lj (x) - density
function of jth order statistic of LLd. stan-
dard uniform sample of size n, 1 :5 j :5 n,
expectation functional for Xj:n
Fj:n{x) = E:=j (:)xk{1 - x)n-k, 0 :5 x :5 1,
- distribution function of /;:n{x)
Gj:n{X) distribution function of jth order statistic
of dependent identically distributed sample
of size n
Gj:n{x) = (nx + 1- j)l[(j-I)/n,lj{x)/{n + 1- j) - sto-
chastically largest distribution function
of jth order statistic of dependent sample
of size n with standard uniform marginal
gj:n{X) = nl[(j-l)/n,l) {x)/{n + 1 - j) - density func-
tion of Gj:n, expectation functional for Y;:n
8 1. Introduction and Notation
c·b(:::S.)W = {g E Cw
' : 9 is convex (concave) on r"w,dw)}
lfO
- family of compositions of symmetric W (x)
with centered quantile functions of C' for
symmetric P t. (js)W
In particular:
This also shows that (2.1) becomes the equality iff 9 = 0, h = 0, or 9 = a:h
for some a: > o. We use (2.1) for verifying that the function
Ilhll = (h,h)1/2
h I---t
defines a norm in 1i. H (1i, II ·ID is complete, then
(1i, (. ,.)) is called the
Hilbert space. The Riesz representation theorem asserts that every linear
continuous functional defined on a Hilbert space can be written as
Th(g) = (g, h), 9 E 1i,
12 2. Basic Notions
for some h E 1l. By (2.1) again, 11Th II = Ilhll. One can see that the normal-
ized nonzero functional Th(g)/lIgll, 0 :F 9 E 1l, attains its maximum IIhll
at 9 = ah with a > O.
In numerous statistical problems, it is important to maximize a linear
normalized functional over a convex cone in a Hilbert space. We say that
C C 1l is a convex cone if f, 9 E C implies that af + f3g E C for arbitrary
a, f3 ~ o. If h E C, then the solution of our restricted maximization problem
coincides with that of the general one. Otherwise we show that h should
be replaced by its projection Ph onto C, that is, the element of C that is
least distant from h. This can be deduced from the following theorem (cf.
Balakrishnan [9, Section 1.4]).
Theorem 1 If h is an arbitrary element of a real Hilbert space 1l and C is
a closed convex cone in 1l, then there exists a unique projection Ph of h
onto C that is characterized by two relations
Setting 9 = aPh for some positive a and using (2.3), the equality holds
in (2.5). If Ph = 0 then, due to (2.2), Th is nonpositive on C, and clearly
Th(Ph) = O. Due to the fundamental significance of Theorem 1 for our
further considerations, we recall its proof here.
Ilgn; h _ gm 2- hl1
2
+ Ilgn; h + gm 2- hl1 2
For arbitrary f > 0, the right-hand side of (2.6) is not greater than D2 + f
if nand m are large enough. Since (gn + gm)/2 E C, the latter term of the
left-hand side is not less than D2. Therefore we have
the former being concluded from the triangle inequality, imply that actually
IIPh-hll =D.
Now we verify that (2.4) completely characterizes Ph. For arbitrary gEe
we define a function Dg on the unit interval [0,1] as
Dg(O:) = IIh - (1 - o:)Ph - o:gll2
= IIh - Phll 2 + 20:(h - Ph, Ph - g) + 0:2 II Ph _ g1l2.
We have
Now we use the fact that C is a convex cone. Plugging 9 = o:Ph into (2.4),
we obtain
(0: - 1)(Ph, h - Ph) ::; 0, 0: ~ 0,
which yields (Ph, h - Ph) = O. This proves (2.3), and in combination
with (2.4) gives (2.2) •
]1/2
fa d
g(x)h(x)w(x) dx::;
[
fa d
g2(x)w(x) dx fa d
h 2(x)w(x) dx (2.10)
14 2. Basic Notions
EXAMPLE 2. Let
C+ = {h E L2([a,d),w(x)dx): h ~ a}.
Verifying (2.2) and (2.3) we deduce that
P+h = h+ = max{h,O}
is the projection of h onto C+ for every h E 1-£. Also, one can check di-
rectly that h+ is actually the nonnegative function closest to h for arbitrary
weight w . •
(d (wed)
(h,l) = 1a h(x)w(x) dx = 10 hW-1(x) dx, (2.12)
and its greatest convex minorant Hw, with a nondecreasing derivative hw,
say.
We prove that h"- defined as h"- = hw W E C~ is the projection of h
onto C~ by checking (2.2) and (2.3). For the former one, we need the
following lemma (cf. Marshall and Olkin [55, Proposition A.2.(iii) , p. 444]).
Lemma 1 If G ~ H are functions of bounded variation on an interval
[A, D), which are equal at the endpoints, then
holds for every non decreasing function 9 for which both the integrals exist.
2.1 Elements of Hilbert Space Theory 15
PROOF (cf. Marshall and Proschan [56]). The statement is easily verified
for all indicator functions l[y,D)(X), A < y < D. Therefore this is true for
all positive combinations of the indicator functions, and, by the Lebesgue
monotone convergence theorem, for arbitrary nonnegative nondecreasing
functions as well. Therefore the reversed inequality holds for the nonnega-
tive nonincreasing functions. Finally, we represent an arbitrary nondecreas-
ing function as the difference of two nonnegative terms, the nonincreasing
and nondecreasing ones, and apply the above statements to both parts . •
(g, h) =
l o
W (d)
gW-1(y)hW-1(y) dy
=
l o
W (d)
gW-1(y) Hw(dy)
<
l o
W (d)
gW-1(y) Hw(dy)
= lad g(x)h'(x)w(x) dx
= (g,h,). (2.15)
and we have
16 2. Basic Notions
Thus the equality holds for the integrals over the whole of the open set
W-l({Hw < Hw}). For the remaining part W- 1 ({Hw = Hw}) the con-
clusion holds, since Hw = Hw implies that h" = h there. Summing up,
we have
(h,h,,) = (h",h,,),
which is the desired conclusion . •
The construction of the L2-projection onto the family of monotone func-
tions under uniform weighting was presented in Moriguti [58]. For the gen-
eral case we refer to Rychlik [85]. For simplicity, we treated elements of
L2-spaces as functions rather than equivalence classes up to almost sure
equality, and we follow this convention later. For example, h E C+ gener-
ally means almost sure nonnegativity, and monotonicity can be precisely
defined by comparisons of integrals over the intervals of the same weight.
In Examples 2 and 3 we were able to determine projections for arbitrary h,
but there are no general rules for constructing projections onto other con-
vex cones. Then we apply arguments suitable for specific functions h and
cones. Usually, we first try to describe the shape of the projection in a para-
metric way, and then precisely determine the parameters. We finally note
that more often the characterization (2.2) and (2.3) is used for determining
projections of principal interest. In the problems under study we first try to
find the projection by minimizing the distance of the Hilbert space point,
representing a functional, to a cone. The solution of the auxiliary problem
satisfies (2.2) and (2.3) in particular which are needed for estimating values
of the functional over the cone.
(2.16)
(2.17)
where
F- 1 (x) = sup{y: F(y):::; x}, 0:::; x < 1, (2.18)
2.2 Statistical Linear Functionals 17
Observe that the functions F-l- J.'F form the convex cone of nondecreasing
functions integrating to O. Our purpose is to determine sharp bounds for
normalized statistical functionals represented as T(F-l ) /mF and T(F- l -
J.'F)/UF for general and restricted classes of quantile functions. Now we
present exemplary linear functionals of statistical importance acting on
quantile functions in L2([0, 1), dx). Considering bounds on narrower classes
of distributions, it is convenient to transform the quantile functions so that
other L2-spaces are studied.
Quantiles F-l(P) of order 0 < p < 1. They characterize distributions by
describing levels that divide respective populations into subsets contain-
ing desired proportions of elements. Quantiles are often used for defining
critical levels of tests and interval estimates. In order to evaluate F-l(P)
in terms of moments, we represent it as a limit of continuous linear L2_
functionals
F-l(P) = q'"p
lim -I-lq F-l(x) dx = lim (F-l, _I_I(p q»),
q- p p q'"p q- p ,
(2.20)
= Fj:nF(x), (2.21)
say. Since
f 1m ()
X = Fj:n
' (X ) = n (n -1)
j _ 1 X j-l( 1 - X )n- ,
j (2.22)
18 2. Basic Notions
we have
EFXj:n 1 +00
-00 x Fj:nF{dx)
= 11 F- 1 {x)/i:n{X) dx
Note that Fj:n and /i:n are the distribution and density functions, re-
spectively, of the jth order statistic of the standard uniform Li.d. sam-
ple of size n. Order statistics are directly used for estimating quantiles
of order j /n and describing the lifetime of the j -out-of-n reliability system
which contains n independent identical elements and operates until at least
n + 1 - j of its elements do. We can also study the expectations of linear
combinations of order statistics (so called L-statistics)
(2.24)
Formula (2.25) is valid for arbitrary coefficients C1, • •• , Cn, and distribution
function F with a finite expectation. The supremum is attained for some
distributions in 'Pn(F). A detailed characterization of the distributions as
well as arguments leading to (2.25) are presented in Section 5.1. Maximizing
(2.25) over a family of marginals F, we first determine the F providing the
extreme value of (F-1, ge), and then take the joint distributions in 'Pn(F) ,
for which the expectation of the L-statistic is actually equal to (F-1, ge) for
the specified F. By definition, ge is a nondecreasing step function with n-l
jumps at most, located at some points of the form j In, j = 1, ... ,n - 1. In
particular, for arbitrary 1 ~ j ~ n,
sup
PeP.. (F)
Eplj:n =
n
n
+1-
.1 1
J . (j-1)/n
F-1(x) dx
= (F- 1, n +~ _ j l[(j-1)/n,l) )
say (cf. also Caraux and Gascuel [19], Rychlik [76]). One can more easily de-
termine projections of simple step functions presented in (2.25) and (2.27)
than of polynomials appearing in respective formulae (2.24) and (2.23) for
the LLd. case. This explains the surprising fact that we have more results
and of simpler forms for arbitrarily dependent samples than for standard
independent observations. Note finally that the projection method allows
us to measure sensitivity of L-statistics upon dependence by evaluating
(2.28)
Record values. Record values in numerical sequences are ones that exceed
all the preceding ones. For a random sequence Xi, j ~ 1, record values,
and respective record occurrence times L n , n ~ 1, are random increasing
sequences. By convention, we assume
(2.29)
Due to another convention, the first value of record occurs at time 1 and
equals Xl. Like extreme order statistics, records are applied in estimating
strength of materials, predicting natural disasters, sport achievements, and
the like. They were first studied by Chandler [20]. Comprehensive studies of
records can be found in Ahsanullah [1] and Arnold et al. [8]. IT observations
are independent identically distributed and the distribution does not have
an atom at its right support endpoint, then the sequence of records is
infinite almost surely. Formulae (2.29), (2.30), and (2.31) are well defined
for arbitrary original sequences, but we assume further that Xi, j ~ 1,
are independent and have an identical continuous distribution function F,
say. It is obvious that if a current value of a record is given, then the
conditional distribution of the next one is identical with the distribution of
the parent variable under the condition that it exceeds the actual record
value. In particular, for records R~ of an i.i.d. sequence xy, j ~ 1, with
the standard exponential distribution function
we have
(2.33)
2.2 Statistical Linear Functionals 21
and
EFRn = EvF-1V{R~)
= (>0 F- 1 {1 _ e- xn e- z dx
Z)
10 n!
= 11
F-1{x)fn{x) dx
= (F-1,fn) (2.34)
with
fn{x) = [-In{l- x)]n In!,
which is a desired inner product representation.
kth record values. An increasing sequence of record values arises from
a nondecreasing sequence of sample maxima (Xn :n ), n ~ I, by crossing
out all repetitions. For arbitrary fixed k, the sequence of kth greatest order
statistics (Xn+1-k:n), n ~ k, is nondecreasing as well. By analogy, we can
define occurrence times and values of kth records in the following way
L(k)
0 = k, Itok) = X 1:k, (2.35)
L(k)
n = min{j > L~k21 : Xj > X L(Io)
n-l
H-k:L(Io) },
n-l
(2.36)
R(k)
n = X L~Io)+1_k:L~Io), n~1. (2.37)
The kth records were introduced by Dziubdziela and Kopocinski [26]. There
is also another convention of defining kth record occurrence times that
consists in subtracting k - 1 from L~k) defined in (2.35) and (2.36). This
implies that we start counting records only when the first k observations
are carried out. In particular, the first value of the kth record occurs at
moment 1 then. However, we are concerned with record values here, which
are not affected by the particular definitions of occurrence times. In contrast
with standard records, a random variable X L(Io) observed at the kth record
time L~k) for k ~ 2 does not necessarily become the kth record value
immediately. Generally, we have
p{R(k)
nH
> yIR(k)
n
= x) = [1- F{Y)] k
1 _ F{x) , Y
> x, (2.38)
{cf. (2.32)), which means that the distribution of the next value of the kth
record, when the current one is known, is the same as the distribution of
22 2. Basic Notions
the minimum of k original variables Xj under the condition that they are
greater than the current record. Relation (2.38) shows that distributions of
nth values of kth records from the i.i.d. sample with distribution function F
and standard first records of the i.i.d. sample with distribution function
=
1
EFR~k) EV(F1:kF)-lV(R:;)
= 00
F- 1(1 _ e-z/k) xn e-Z dx
o n!
= 11 F-1(X)f~k)(x) dx
where
k,~ 1, n ~ 0, (2.40)
is the density function of the nth value of the kth record of the LLd. stan-
dard uniform sequence, and fA 1) = fn defined in (2.34). Formula (2.38) is
true for arbitrary F. However, for the discontinuous F function F-1 F1-;-iV
is not strictly increasing and may transform a record in the exponential se-
quence into a number equalizing a previous score. Therefore (2.39) is true
only for continuous F.
EFS~k) = 11 F- 1(x)fAk)(1 - x) dx
= 1
o
1 kn+1
F-1(x) _ _ (_lnx)nxk-1 dx.
n!
2.2 Statistical Linear Functionals 23
D ( ) _ F(x) - F(y)
.c'11I x - I _ F(y) , x ~ y,
= F-l(F(y) + [1 - F(Y)]X)/;-i:n-i(X) dx
= r
10
l
F-l( )f... . (x -
F(Y)) l[F(II),l)(X) dx
X }-I.n-I 1 - F(y) 1 - F(y)
(2.41)
for 1 :$ i < j :$ n. Similar results are obtained for the reversed conditioning.
The conditional distribution of the X j :n when Xk:n = z is given for k > j is
identical with the distribution of the jth order statistic from the sample of
size k -1 with the common distribution function of Xl under the condition
that Xl :$ z which equals
F lz( x ) = F(x)
F(z) , x:$ z.
H both Xi:n = Y and Xk:n = z for i < j < k are known, then the condi-
tional distribution of X j :n coincides with that of Xj-i:k-l-i of the doubly
truncated population with distribution function
(see also Arnold et al. [7, Section 2.4]). Respective analogues of (2.41) are
EF(Xj:n!Xk:n = z)
11°
= EFlzXj:k-l
EF(Xj:n!Xi:n = y, Xk:n = z)
= EFI,Xj-i:k-l-i
1°
III
1 -1 ( X - F(y) ) l[F(y),F(z»(X)
= F (X)!i-i:k-l-i F(z) _ F(y) F(z) _ F(y) dx (2.43)
= 1- [l-F(X)f
1- F(y)
= 1 - [1 - Fly(x))k
= Fu(Fly)(x).
2.3 Restricted Families of Distributions 25
= EFl:kFlyRn-m-1
= EF R(k)
III n-m-l
r
10
1 F-1( )ik) (x - F(Y)) l[F(y),I)(X) d
x n-m-l 1 _ F(y) 1 _ F(y) x
r 1 k n- m (1 - x)k-l
10 F- 1(x) (n - m - 1)![1 - F(y)]k
1 - F(Y)] n-m-l
x [ In 1- x l[F(y),I) (x)dx. (2.44)
p1Z()= V- 1F(min{x,z})
x V-IF(z)'
Accordingly, elementary calculations lead us to the linear functional repre-
sentation
EF(RmIRn = z)
[Link] :n -
11
1
=
rF
10
1 -1
(x)fm:n-l
( -In(I-X)) l[o,F(z)) (x)
-In(1 _ F(z)) -(1 _ x) In(1 _ F(z)) dx.
In this case we have another pair of location and scale parameters: the left
endpoint of support a F = 0 and the square root of the second moment m F,
respectively.
It is convenient to study the upper halves of p-1 - [Link], which form the
cone
and extend the functions to the whole unit interval using (2.49). This im-
plies the modification of functionals which consists in symmetrizing them
about 1/2. Indeed, by (2.49), this yields
= r1
J1/2
2[p-1(X) _ J.L]hS(x) dx (2.51)
for
h 8 (x) = [h(x) - h(l- x-)]/2. (2.52)
Observe that the norm of p-1 - [Link] in L 2 ([1/2, 1), 2dx) is aF.
2.3 Restricted Families of Distributions 27
°
with C~ u (C~ u) and C; u (C~ u) denoting intersections of (2.53) «2.54),
respecti~ely) "Vrith linear-~ubsp~es of functions integrating to (vanishing
at 0, respectively). The apparently awkward notation is justified below. We
say that F(x) succeeds the standard uniform distribution function
U(x) = x, ° ~ x ~ 1,
in the convex order (written as F tc U) if F- 1U = F-1 is convex on [0,1).
The reversed relation F :jc U is defined by the convexity of U- 1 F = F on
the support of F, being equivalent with the concavity of F-1. By conven-
tion, we call distributions satisfying F tc U and F :jc U the decreasing
and increasing density distributions, respectively. The convex order, de-
fined in van Zwet [100], is a partial order of (absolutely) continuous distri-
bution functions, invariant under location and scale transformations (see
Dharmadhikari and Joag-dev [25, Theorem 9.1, p. 217]). Therefore we can
generalize (2.53) and (2.54) as
C{.,w = {g E L 2([aw,dw),w(x)dx) : 9 - nondecreasing
and convex}, (2.55)
C~w = 2
{g E L ([aw,dw),w(x)dx) : 9 - nondecreasing
and concave}, (2.56)
by taking compositions F-1 W for an arbitrarily fixed distribution func-
tion W with support [a,d) = [aw,d w ), and density w. We also introduce
convex cones C~cw (C~cw) and C~cw (C~cw) by adding conditions
ld g(x)w(x) dx = 0,
g(a) = 0,
respectively, to definition (2.55) «2.56), respectively). These definitions are
justified by properties of compositions of F- 1 with W. In particular, we
have
Function AF(X), called the failure (hazard) rate, describes the infinitesi-
mal probability of failure at a short period of time following x under the
condition of surviving until x. Condition F ~c V defines the family of dis-
tributions with increasing failure rate (IFR, for short). Likewise, F tc V
coincides with nonincrease of AF, and defines the family of decreasing failure
rate (DFR) distributions. We specify here bounds on statistical functionals
restricted to IFR and DFR distributions. Since U -<c V, every DFR distri-
bution has a decreasing density, and the increase of density implies that of
the failure rate.
van Zwet [100] (see also Lawrence [49]) defined a counterpart ~. of
convex order for symmetric distributions, and called it s-order. We have
F ~8 W for F and W symmetric about [Link] and [Link], respectively, if
F-lW(x) is concave for x ~ [Link]. This is equivalent with convexity of
F-l W on the lower half of the support of W, and concavity-convexity
of W- l F on the support of F. Combining constructions leading to (2.50),
(2.55), and (2.56), we define the convex cones
AF(X) = -In[l -
x
F(x)] =.!.
x
r AF(Y) dy
10 (2.61)
2.3 Restricted Families of Distributions 29
P(x)
x
=.! r f(y) dy
x Jo
is nondecreasing (nonincreasing, respectively). Accordingly, the relation de-
fines the family of life distributions with increasing (decreasing) density on
the average: although f may be multimodal, the larger values in [aF, dF)
are more (less) probable than the smaller ones. The star order is scale in-
variant. In order to make it invariant with respect to translations as well,
we generalize the definition as follows: P !:::* (~*)W iff aF,aW are finite
and [F-IW(X) - p-lW(aw)]/(x - aw) is nondecreasing (nonincreasing)
on [aw, dw). The definition enables to establish mean-variance bounds on
statistical functionals by projecting them onto convex cones of p- 1 W - [Link]
described by the generalized star relations:
_.
c~ w = {g E L2([aw,dw),w(x)dx) : g(x), g(x) - g(aw) are
x-aw
nondecreasing and
ld g(x)w(x) dx = (g,l) = 0
for all 9 E C~ is assumed. Let C( denote the extension of C~ by dropping the
integral condition. Note that each C( is translation invariant (i.e., fulfills
the assumption of Lemma 2 with go = 1). Therefore (h, 1) = 0 implies that
PI' h E C~ and coincides with p2 h. Otherwise we replace h by
(h,l)
ho =h- (1,1/ =h - (h,I), (2.64)
is equivalent to
F k- 1 ~ F- 1 in L2([0, l),dx).
More intuitively, we can express the relation in terms of random variables.
If Fk, k ~ 1, are distribution functions of a given family, and F generated
by the projection method does not belong to the family, then Fk, k ~ 1,
attains in the limit the bound attained by F if
that is, F;l(X) converges to F-1(X) in the mean square for a standard
uniform random variable X.
For instance, (2.53) contains the quantile functions of distributions with
decreasing density, and, possibly, an atom at the left endpoint ap of sup-
port. Suppose that the solution of the projection problem, dual to a problem
of evaluating a linear statistical functional over decreasing density distribu-
tions, has an atom. This can be replaced by absolutely continuous (uniform,
say) components with increasing concentration about ap so that we obtain
decreasing density distributions with quantile functions tending to that of F
in L2([0, 1), dx). Moreover, these absolutely continuous approximations can
be modified so as to preserve desired moments of limiting F.
3
Quantiles
The main results of this chapter come from Rychlik [90]. The bounds for
quantiles of general distributions were obtained by Moriguti [58]. Those for
symmetric and symmetric unimodal distributions may also be concluded
from the Chebyshev and Gauss inequalities, respectively. Vysochanskii and
Petunin [102] presented a refinement of the Gauss inequality for unimodal
distributions. Further generalizations can be found in Dharmadhikari and
Joag-dev [25, Section 1.5]. We also notice that the Markov inequality yields
F- 1 (P) ~ fJ.F
I-p
-1(p) mF
F ~ (I _ p)1/2 .
1 d
g(x) H(dx) ~ 1 d
g(x)li(x) dx ~
[
1d
g2(x) dx l d
li2 (x) dx
]1/2
(3.1)
{ -1 if o ~ x <p,
li(x) = ~, if p ~ x < 1,
P
Illil1 2 ,
I-p
3.1 General and Symmetric Distributions 35
and
F-l(X) - J1. h(x)
a = Ilhll·
Moriguti [58] also showed that
F- 1 (q) _ F-l(P) (1 1)1/2
--='------"'--'-
aF
< -- + -
- 1- q p , O<p<q<l. (3.3)
Special cases of (3.2) and (3.3) for the median and interquartile distance
yield
for a < b < c ::; d and M > O. Then for every 9 E C( w there exists
ga.{3 E C(
_c
w defined as -
(3.6)
If
JbC(X - a)w(x) dx
fC ( )
Jb W x dx
~
a
d( ) ( )l
x - a w x dx, (3.8)
then p( wh = 1.
Othe1wise there exists a unique fJ* < b that solves
and
a* = a*(fJ*) = h.
[ d(x - fJ*)w(x) dx
]-1
>0 (3.10)
such that
/' () _ ( _ (x - fJ*)+ (3.11)
Pt,c wh x -a. x-fJ*)+ - d' •
Jmax{{3.,a}(y - fJ*)w(y) dy
Precisely, if relation "~" holds in (3.9) with fJ replaced by a, then fJ* ~ a,
and (3.11) is linear on [a, d). In the opposite case, fJ* > a and the projection
is actually a two-piece broken line.
PROOF. Our purpose is to minimize
D(a,fJ) = Ilga{3-hW
= a2 1d (x - fJ)2w(x) dx
r
max{{3,a}
2a 1
- Jbc w(x) dx Jb (x - fJ)w(x) dx + Jbc w(x) dx (3.12)
with respect to a and fJ. For fixed fJ < b, we easily find optimal
JbC(x - fJ)w(x) dx
a* = a* (fJ) = d > 0, (3.13)
J: w(x) dx Jmax{{3,a}(x - fJ)2w(x) dx
which plugged into (3.12) gives
38 3. Quantiles
By differentiation,
x [l b
C
w(x) dx
d
r
lmax{/3,a}
(x - j3)2w(X) dx
- 1r d
max{/3,a}
(x - j3)w(x) dx r (x -
1b
j3)w(x) dxl.
Since the factor in the first line is positive, the sign of the derivative is
identical with that of the last two which is denoted by K(j3). Observe that
lmax{/3,a}
r (x - j3)w(x) dx
- lbr w(x) dx lmax{/3,a} d
evaluated at j3 = b satisfies
K'(b)
d = Ew(Xlb < X < c) - Ew(XIX > b) < O. (3.16)
fb w(x) dx fbc w(x) dx
l
Finally,
K"(j3) = -w(j3)
c
(x - j3)w(x) dx ~0 (3.17)
under the convention that w(j3) = 0 for j3 < aw. By (3.17) and (3.16),
K'(j3) is constant for j3 < a and nonincreasing to K'(b) < 0 for a ~ j3 < b.
If K'(a) ~ 0, which means that (3.8) holds, then K'(j3) ~ 0 and K(j3)
is nonincreasing and, by (3.15), positive for all j3 ~ b. Therefore (3.14) is
minimized at j3 = -00, which implies that the projection is a constant.
This amounts to (h,l) = 1, which proves the first assertion.
If K'(a) > 0, then K(j3) changes its sign once from minus to plus at
some j3 < b, at which (3.14) is minimized. Equations (3.9) and K(j3) = 0
are equivalent, and allow us to rewrite (3.13) as (3.10). Note that K(a) 2:: 0
3.2 Distributions with Monotone Density and Failure Rate 39
implies that the solution to (3.9) satisfies (3 $ a, and that the projection is
linear. In the opposite case a broken line with break at (3 > a is obtained.
This proves the final claim of Lemma 5. •
Using (3.11) and writing (3
'-*
= max{(3.,a}, we obtain
IIPgcwh -IW = IIPSw h 11 2 - 1
= -.
I; (x - (3.)2W(X) dx - [J;-. (x - (3.)w(x) dX]2
[J:. (x - (3.)w(x) dx]2 (3.18)
P~cwh(x) -1
IlPgcw h -III
(x - (3.)+ - I; (y - (3.)w(y) dy
= {d -. d }1/2' (3.19)
It!.. (y - (3.)2 W(Y) dy - [It!.. (y - (3.)w(y) dy]2
Letting c ~ b, we reduce the right-hand sides of (3.8) and (3.9) to b - a
and b - (3, respectively. Analyzing the resulting formulae, we derive bounds
for quantiles. We introduce some notation before presenting them. For
a random variable X with distribution function W, and arbitrary (3 E (a, d),
denote the expectation and variance of X under the condition of exceeding
level (3 by
I;xw(x)dx
J.'w ((3) = Ew(XIX > (3) = d ' (3.20)
1/3 w(x) dx
= Varw(XIX > (3)
I; x 2w(x)dx - I-''tv((3)
=
I; w(x) dx (3.21)
We also define a distance between the pth quantile and the conditional
mean in respective standard deviation units
8 ((3) = W-1(p) - I-'w((3) (3.22)
w ow ((3) .
Theorem 4 (P tc W) IfW- 1 (p) $ I-'w, then p-l(p) $ I-'F.
If J.'w < W-l(p) $ J.'w + u~/(I-'w - aw), then
p-l(p) - I-'F < r ( ) _ W- 1 (p) - I-'w
_ uw aw - , (3.23)
UF UW
and the equality holds for the location-scale transformations of W, that is,
for
P(x) = W ( I-'w + uw- X-I-')
u- . (3.24)
40 3. Quantiles
Finally, if W-l(P) > [Link] + a'tv/([Link] - aw), then there exists a unique
solution (3* = (3*(P) E (aw, W-l(P)) to equation
[W- 1 (P) - [Link] «(3)][[Link] «(3) - (3] = a'tv«(3), (3.25)
and
which becomes an equality for the uniform distribution on [I' - .;aa, I' +
.;aa].
If 2/3 < p < 1, then
we get
F- 1(P) - f,LF < (_2__ 1)1/2 (3.31)
O"F - 'Yv(P)
This is an equality if F is the combination of an atom at f,L-0"['Y/(2-'Y)]1/2,
and the exponential distribution with location h(2-'Y)j1/2f,L-'Y0" and scale 0",
with respective probabilities 1 - 'Y and 'Y.
We see that bounds (3.29) and (3.31) tend to infinity if p /' 1, and the
same holds generally for (3.26). It follows from the fact that for b = W- 1(P)
large enough, f3 = f3(b) satisfies
I;
x(x - f3)w(x) dx
b = -'--d.,------- (3.32)
1/3 (x - f3)w(x) dx
(cf. (3.9) for f3 > a). Therefore b(f3) /' d, as f3 /' d, and the same holds for
the inverse. Furthermore, 1 - W(f3) ~ 0, whereas the nominator of A~ w
remains bounded below from zero, as b /' d. _0
(3.33)
for some a > 0, and a ::; [3 ::; b ::; 'Y ::; c with a(-y - [3) ::; M such that
°
that g+ provides a better approximation. The latter two imply that this
can be further improved by adding a constant C > such that
(g+ + C, 1) = (h,I).
We can repeat the truncation and translation operations several times until
we eventually obtain a function positive at c.
42 3. Quantiles
Then we assume that g( e) > 0, and separately study three cases: g(b) :5 0,
g(a) :5 0 < g(b), and g(a) > O. In the first one, we construct the line
passing through (b,O) and tangent to the graph of gin (b, e], and truncate
it at level gee). The resulting function has the desired form, runs above 9
and below h in [a, e), and the reversed relations hold in [e, d). In the second
case, we take the linear function secant to 9 at (3 E [a, b) satisfying g((3) = 0
and b. The truncation defined above provides the conclusions of the former
case. Finally, if g(a) > 0, then it suffices to replace point ((3,0) by (a, 0) in
the construction of the second case .•
D(a, (3,,) = l d
[a(min{x,,} - (3) + 1- h(x)]2w(x)dx (3.34)
with respect to a and (3. However, this is a simple matter, because (3.34)
is a convex quadratic function in both arguments. Therefore we have
Lemma 7 Under the hypotheses of Lemma 5, we have
p{
-<
wh(x) = a("{Hmin{x,,} - (3("{)] + 1 (3.35)
L<min{:Il,'Y}w(:Il)d:ll _ (3( )
1< w(:Il)d:ll '
a(,) = (3.37)
= Varwmin{X,,}
= l d
[min{x,,}]2 W(X) dx - ti~("{), (3.39)
W-1(P) - tiw(,)
= (3.40)
OW("{)
3.2 Distributions with Monotone Density and Failure Rate 43
g. = limP{ wh E C~ w.
c'-"b -" -"
_.
then Pt wh = 1. Otherwise
(3. = (3.53)
46 3. Quantiles
= la
d[ l[b C) (X)
f3 + a(x - a)l[b,d)(X) - Ibc ~(x) dx
]2 w(x) dx
= f32 - 2f3 - It ;)
b W x dx
r
+ a 2 d(x - a)2w(x) dx
Jb
Ew«X - a)21X > b) - Ew{X - alb < X < c)Ew{X - alX > b)
> Varw{XIX > b) > O.
Moreover, by (3.55), f3 < 1 < 1/[f: w{x) dxl, and so (3.53) satisfies both
constraints deduced from the geometric arguments presented in Lemma 8.
3.3 Distributions with Monotone Density and Failure Rate on the Average 47
JbC(x-a)w(x) dx J,d( ) ( )
I.c w(x) dx - b X - a w x dx
[Link]-III b 1/2' (3.57)
{Jbd(X-a)2 w(x)dX - [Jbd(x-a)w(x)dxf}
p~.wIt(x)-l _ (x-a)I[b,d)(x)-Jbd(y-a)w(y)dy· ( )
IIP~. wit-III - { d d 2}1/2' 3.58
- Jb (y-a)2 w(y)dy - [Jb (y-a)w(y)dy]
Letting c '\t b, and using the notation of Theorem 4, we state its analogue
for F ~* W.
Theorem 6 (F ~* W) Using (3.20) and (3.21) with f3 = b = W- 1(P),
write
17=17w(a,b) = Ew[(X - a)I[b,d)(X))
with
8~(P) = 12t?~(0,p) = 1 + 6p2 - 4p 3 - 3p4. (3.64)
Relation (3.63) becomes an equality for the mixture of an atom at J1. -
V30'(1- p2)/8u(P) with probability p, and the uniform distribution on [J1. +
V30'(p2 + 2p - 1)/8u(P), J1. + V30'(1 + p2)/8u(P)] with probability 1 - p.
lim
p,-1(p)_
k
(
J1.= ~
)1/2
k-too 0' 1- p
The claim can be justified by some informal arguments. To fix the ideas,
we consider the approximation of two-point distribution F defined in The-
orem 2 by Fk ~* U, that is, ones with nondecreasing Fk{x)/{x - ak) on
respective supports [ak, dk). Note that F is constant beside of its two jumps,
and the starshaped Fk have to increase at rate 1/{x - ak) at least on their
supports. However, letting ak -t -00, we relax conditions on the increase
rate of Fk itself. Accordingly, such Fk can be starshaped and approximate
well any two-point distribution at and beside jump points. The formal proof
of Theorem 7 is constructive. It is also possible to provide special construc-
tions of sequences satisfying the statement for specific W with unbounded
densities about respective left support ends.
1
F- W{x) -
(jF
JLF -gp
_ () _
x - {i!.l--
~
I-p'
p' if a ~ x < b = W-1(P),
if b ~ x < d.
(3.68)
Our proof consists in constructing a sequence gk E C~ w, k -t 00, that
converges to (3.68) in L2([aw, dw), w(x)dx). For sufficiently large k, we
define piecewise linear continuous nondecreasing starshaped functions
k4 (x - a) - k - J¥-, if a ~ x ~ a + k- 3 ,
-J!=1!. p , if a + k- 3 ~ X ~ b,
k~
b-a - k- J!=1!.
p ,
if b< x
-
<
-
b + b-a
kv'p(l-p) ,
Observe that
= 1 a+k- 3
[k 4 (x - a) - k]2W(X) dx
because
and
Relations
[Fk"lW(x) - /l-F]/UF = 91e(X)
define a sequence of Fie -<* W with common mean /l- and variance u 2 for
which Fk"l(p) /" Jp/(l- p) .•
for q '\t ponto e:,.c w as defined in (2.60). If p < 1/2, then hB $ O. Its
projection onto the family of nonnegative functions e+ ::> e:,.cw is P+h B =
o which actually belongs to e:,. w, and so P:" whB = O. This sequence of
elaborated arguments leads us-to the trivial c~nclusion
Otherwise we have
Lemma 10 Let
hB = Ml[b,c) E L 2 ([[Link],dw),2w(x)dx)
for /l-w < b < c $ dw and M > O. Then for every 9 E e:,.c w there exists
9a{3 E etc w defined as -
(3.74)
= (3.75)
Otherwise
I;. (x - f3.)2 w(x) dx
(3.76)
d ] 2'
2 [ I{3.(x - f3.)w(x) dx
(x - 13.)+
= 1/2 . (3.77)
[2 I;. (y - f3.)2 w(y) dy]
Letting e ":>I b, recalling (3.20) through (3.22), and introducing
~w = EwlX - [Link] =2 i d
(x - f.L)w(x) dx (3.78)
which becomes the equality for the uniform distribution on interval [I' -
vau, I' + vau].
If 5/6 $ p < 1, then
C11
tIo)
54 3. Quantiles
ones. Certainly, the general bounds are merely valid for wider classes of
life distributions (e.g., NBU, HNBUE and C-class defined in Klefsjo [46])
which are generated by more general orders than the star order. On the
other hand, there are still interesting open problems of accurate bounds on
quantiles of distributions succeeding a given one in these general orders.
Bounds for expectations of order statistics from general i.i.d. samples, due
to Moriguti [58], are presented in Section 4.1. Bounds of Sections 4.2 and 4.4
for populations with decreasing density and failure rate, and symmetric uni-
modal ones, were obtained by Gajek and Rychlik [33]. Results of Section 4.3
for restricted families determined by the star order come from Rychlik [89].
Section 4.5, partially based on Okolewski and Rychlik [65], is devoted to
the study of quantile estimation bias in various nonparametric families of
distributions.
We do not discuss here bounds for the sample maximum and range from
discrete populations obtained by LOpez-Blazquez [50, 51]. We also merely
mention sharp inequalities due to Papadatos [69] for the expectations of
order statistics and their differences of nonnegative samples expressed in
terms of the population mean. Blom [17] and van Zwet [100] developed
another method of evaluating the expectations of order statistics from re-
stricted families of distributions defined by means of stochastic orders using
quantiles of the parent distribution. The method is based on the Jensen in-
equality. Papadatos [67, 68] derived some attainable bounds on the variance
of order statistics from general and symmetric distributions, respectively,
measured in the standard deviation units of the parent distribution. Anal-
ogous evaluations for covariances of order statistics are known, but only
small samples of size n = 2 and 3 were treated (see, e.g., Ma [52] and
Papadatos [70]). Bounds and approximations for moments of order statis-
tics were reviewed in David [22, Chapter 4], Arnold and Balakrishnan [5,
Chapters 4,5], and Rychlik [84].
56 4. Order Statistics of Independent Samples
with
1 1/(n-1}
n- 1 x-JL )]
F(x) = [
~ 1+
(
(2n _ 1)1/2 U
_ (2n - 1)1/2 < X - JL < (2n _ 1)1/2 (4.3)
n-1 - U - ,
n n
EF LCj(Xj:n - JLF)/UF ~ LCj(fj:n -1) (4.4)
j=l j=l
°
continuously increase for a E [0, (j - 1)/(n - 1)] and so do , 0 (1), ranging
from to 'U-1)/(n-1)(I) > Fj:n (l) = 1. The line 10 • for some 0.* E (0, (j-
1)/(n-l» such that 10 • (1) = 1 becomes a part of the lower convex envelope
of Fj:n on [0.*,1], and the remaining part coincides with Fj:n . Therefore,
the greatest convex minorant Fj:n of Fj:n has form
p. ( ) _ { Fj:n(x),
3: nX -
if °
~ x ~ 0.*,
/i:n(a*)(x - 1) + 1, if 0.* ~ x < 1,
(4.6)
we calculate
F(x) = {° ,
f:-1 (1 +
):n
BZ-I:')
".'
'J71
il
~
~". -<_.!.
".
B'
_.!.B -< Z-I:' < f;,n(a.)-1
B' (4.12)
1, 71 Z-I:' > Ji:n(a.)-1
'J ". - B .
_ 2 n 2{n - 1) [ n 2n - 3 {n - 2)2n-1]
Ilh:nll = (2n _ 3){2n _ 1) 1 - (n _ 1)4(n-1) . (4.13)
0, if ~:::; 0,
F{x) = { fj-;~ (lIhnll~), if 0:::; ~ < f(I'ii~~il)' (4.16)
1 if .L > hn(a.)
, m - IIhnll'
for 2 :::; j :::; n - 1, and
0, if ~:::; 0,
F{x) = { [(2n _ 1)-1/2~]1/(n-1) , if 0:::; ~ < {2n - 1)1/2, (4.17)
1, if ~ ~ (2n - 1)1/2
for the sample maximum. Both (4.16) and (4.17) are life distributions.
Density function hn is decreasing and so projection p/' hn = 1 f/. C+,
but the constant can be approximated in L2{[0, 1), dx) by functions
This implies that the trivial bound for the sample minimum
(4.18)
EFXn:n - J.I. = n r
11/2
[F-1(X) - J.I.][x n- 1 - (1- X)'~-1] dx
[I.;, dxt'
r
< n [ZR-1 - (1- z)R- 11' u/v'2
(4.21)
(cf. (4.2» was obtained by Moriguti [57] through use ofthe Schwarz inequal-
ity. The extreme distributions for which the equality holds have quantile
functions equal to x n- 1 - (1- x)n-1 up to affine transformations.
Bounds for general order statistics of symmetric populations are derived
by projecting differences Sj:n onto the family of nondecreasing functions
in L2([1/2, 1), 2dx). For this purpose we first analyze variability of the
differences, using the auxiliary results of Gajek and Rychlik [33, Lemma 3,
p.167].
Lemma 12 Consider the function defined in (4.20) for x E [1/2,1) and
(n + 1)/2 < j < n. This is nonnegative, and equal to 0 at 1/2 and 1, and
increasing-decreasing. For j ::; n - 2 function (4.20) is concave-convex if
j ::; [n + (3n - 5)1/2]/2, and convex-concave-convex othennise. Also, Sn-1:n
is concave for n ::; 7, and convex-concave for n ~ 8.
Now we merely apply the first statement of the lemma. Further properties
are needed in Section 4.4. Observe first that, by symmetry, Sj:n ::; 0 for
j ::; (n + 1)/2, and so
'-(2nnif'
(2!-2W n - 2i )
liPs /i:nW = n [F2j-1:2n-1 (a*) + F2n-2jH:2n-1 (a*)
(4.26)
1, 17 - 2B '
c~
_c
w = {g E L 2 ([aw,dw),w(x)dx) : 9 - nondecreasing, convex,
g(aw) = O}. (4.28)
In the sequel we frequently abstract from the specific form of h and assume
the following
P:cufn:n
_ = fn:n E C~ u,_c
h(x), if a ~ x ~ p,
haP(x) ={ h(P) + a(x - P), if P ~ x < d. (4.34)
Then we have
g(x) :$ 1000(x) :$ h(x), if 0:$ x :$ 6.,
h(x) :$ 1011(x) :$ g(x), if 6.:$ x < d,
and the latter holds in a left neighborhood of 0, as well. By convexity of g,
we have
1011(x) :$ g(x), a:$ x :$ o.
This relation, combined with g( 'Y) = h('Y), implies the existence of (3 E b, oj
such that 1011 ((3) = h((3) and 1611 ~ h' ((3). Therefore
l~ = h'(-y) =
> h(6.1 ~(-y) ~ g(6.1 =~('Y) ~ g'('Y+). (4.37)
If 'Y > a, then, by (4.37), we have h'('Y) < g'(-y-), and h < 9 on a left
neighborhood of 'Y. Summing up, under condition g(-y+) < h(-y+), three
cases are possible:
64 4. Order Statistics of Independent Samples
under condition ~ < d (cf. Case (iv)). For x < ,,(, we proceed as in
Case (iii). We have
P:
_c
wh(x) = a.(x - a)
4.2 Life Distributions with Decreasing Density and Failure Rate 65
whose slope 0:. can be easily determined. Note that for fixed {3 E [a, b],
= 2K({3)L(,8) , (4.42)
say. We minimize D (0:. ({3),,8) by determining the set
/C = {,8 E (a, b) : K({3) ~ O},
and analyzing sign changes of L in /C. If /C = (a,l\;] for some I\; E (a, b]
(which actually holds in special cases W = U, V), then L({3) is positive,
negative, or negative-positive, with zero at some .oX E /C, and so optimal
{3. = a, 1\;, and .oX in the respective cases. Specific forms of
does not exceed the number of sign changes of the sequence ao, ... , am. The
first and last signs of (4.43) are identical with the signs of the first and last
nonzero elements of ao, ... , am, respectively.
This is equivalent to saying that Bk,m, °
$ k $ m, for fixed m form
a Chebyshev system on [0,1] (see Karlin and Studden [45, Theorem 1.4.1]).
In addition, differentiation and integration operators defined in (4.41) and
(4.42), acting on compositions
and .8. = .8.(j,n) being the smaller of the smallest positive zeros of poly-
nomials
4.2 Life Distributions with Decreasing Density and Failure Rate 67
j (n - j + 3)'
Ku(x) = 3 L(j + 1 - k)!k:n+2(X) - (n _ j)! . !;-l:n+2(X)
k=1
. 3 (n - j + 2)!
+ (n - J -"2) (n _ j)! hn+2(X), (4.48)
j-l 3j - k
Lu(x) = L(n+ 1- -2-)!k:n+2(X)
k=1
(n - j)(n - j -I)! () (4.49)
- 4 j:n+2 X .
(4.53)
for
(21- 2)(2"-2;)
= n '-(2"" if' F2j-l:2n-l(')'.. )
+ (1-1'.. )[2a~+2a .. hnh'.. )+ !J,n(')''')] , (4.54)
(4.55)
68 4. Order Statistics of Independent Samples
0, i'Il £
m - < 0,
F(x) = { fJ-;~ (Brii) , il 0< £ <
!;,n(-r.)
'I - m - B '
1 - (1 - 'Y.) exp ( il
'I
£
m -
>
!;,n("Y.)
B .
(4.58)
Formulae (4.45) through (4.47), and (4.53) through (4.55) are concluded
from (4.34) and (4.41). Moreover, (4.48) and (4.49) as well as (4.56) and
(4.57) differ from K and L defined in (4.42) by positive functional factors
common for all j, and original variables are replaced by V (x) in the lat-
ter case. All bounds are achieved by absolutely continuous distributions.
The left-hand parts of (4.50) and (4.58) are the inverses of polynomials
equal to /j,n up to scale factors, and cannot be written explicitly. They
have uniform and exponential right tails, respectively. Bound (4.15) was
derived directly from the Schwarz inequality. Using the following integral
approximations of harmonic series
n+ 1 n + 1/2
In n+ 1 -J. < I'V;on < In n +1/2- J. ,
we deduce that condition 1'V;,n ::; 2 providing bound (4.52) is true for
j ::; (1 - e- 2 )(n + 1/2) and false for j ~ (1 - e- 2 )(n + 1). The difference
between both estimates is (1 - e- 2 )/2 ~ 0.43233 which implies that for
given n the condition has to be directly checked for one j at most. Bounds
of Theorem 12 are tighter than those of Theorem 11, because a narrower
class of distributions was treated there. This is confirmed by numerical
comparisons presented in Gajek and Rychlik [33, Table I, p. 162].
4.3 Distributions with Monotone Density and Failure Rate on the Average 69
h((3), if a ~ x ~ (3,
hf3a(x) = { h(x), if (3 ~ x ~ a, (4.59)
h(al:::~(f3) (x - a) + h((3) , if a ~ x < d,
are also possible here. However, starshapedness is not apparent here. For
fixed (3 E [a, c), let
sf3( "')
.... -_ h(a) - h((3) , a =f. a E [(3, c), (4.61)
a-a
denote the slopes of lines
passing through B = (a, h(fJ)) and points A = (a, h(a)) , fJ $ a < c, of the
graph of h (cf. the last line of (4.59)). Write
8a (a) = h'(a) ~ 0
for completeness. We first prove the following.
Lemma 16 For every fJ E [a,c) there exists a(fJ) E B(fJ) = [max{fJ,b},c)
such that (4.61) increases on (fJ,a(fJ)) and decreases on (a(fJ),c).
In consequence, (4.59) is actually starshaped if a E [fJ,a(fJ)]. Observe that
a(fJ) is the point where lBA.{,8) is tangent to h and therefore can be deter-
mined by
h(a) - h(fJ) = h'(a)(a - a). (4.63)
The former inequality holds, because the middle term defines the line that
is tangent to h at a«(3) and majorizes h. The latter is a consequence of
the fact that the line is the smallest starshaped function in [a«(3) , d) that
passes through Band ..4«(3) = (a(.B) , h(a«(3))). The above arguments show
that hf3ii:(f3) actually satisfies (4.60).
However, h -:f:. hf3ii:(f3) ~ h and, due to Lemma 2, the approximation can
be further improved by a downward translation of hf3ii:(f3). We consider the
case g(a«(3)) < h(a«(3)) now. Define
Also, we have
g(x) ::; h«(3) ::; g(x), x E [a, (3].
Consequently, hf3a defined as h«(3), h, and IBD in [a, (3j, [(3, aj, and [a, d),
respectively, lies closer to h than g. Since (3 ::; a < a«(3) , we see that
hf3a E cC w·
If (3 ~·a, it may happen that lBD has no sign changes in [a,o). Then
IBD(a) = h(a) = 0 and lBD < h in (a,o), and lBD > h in (0, d). The
last relation is a consequence of concavity and ultimate decrease of h. If
lBD = 0::; h, referring to Lemma 2 we decrease the distance to h by adding
72 4. Order Statistics of Independent Samples
a positive constant I~D. H IBD has a positive slope, we can take a line IBD
running through (8,g(8)) with a slightly smaller slope. Observe that IBD
lies closer to h than the original one in [a, 8] and [8, d). Both modifications
lead to linear functions that cross h once at some a < ii«(3). Therefore we
are in position to apply the construction of the previous paragraph, which
ends the proof. •
Below we determine the optimal parameters. For a ~ (3 ~ c and (3 ~ a ~
c, a '" a, set
k(p) = i d
[1 - h(x)] (x - a)w(x) dx ~ 0, (4.68)
_.
such that Pt' wh = ha•.a., defined as in (4.59).
Since h is strictly increasing from h(a) = 0 to h(c) = suph > 1 (cf. (4.29)
and (4.30)), p is actually well defined.
PROOF OF LEMMA 17. By Lemma 15, we should minimize
(cf. (4.61) and (4.62)) with respect to two parameters (3 E [a, c) and a :f;
0: E [max{(3,b},o((3)) c (a,c). Fixing (3 and differentiating (4.71) with
respect to 0:, we obtain
8D((3,0:)
80:
= 0:: a [h'(o:) -l~A] id[lBA(X) - h(x)]w(x) dx
8K((3,0:)
J::l
vO:
_ h'(o:) -l~A
-
0: - a
I d( _ )2 ()d
0<
X a wx x>,0
which implies that (4.72) is the product of a positive function and increasing
K((3, .). Since the integrand is positive for 0: = 0((3), we have K((3, 0((3)) >
O. If 0: = (3 > a, then lBA = h((3). It follows that K((3) < 0 as (3 '\.- a and
K((3) > 0 for (3 /" c.
We can summarize the behavior of (4.71) as follows. If (3 is small enough
then (4.72) is negative for 0: close to (3, and changes its sign at an 0:.((3) E
((3,0((3)) where (4.64) vanishes and the unique minimum of D((3,·) is at-
tained. If (3 ;::: S satisfying K(S) = 0 then (4.72) is positive for all 0: > (3.
Then D((3,·) is minimized at 0:.((3) = (3 which gives a constant approxi-
mation h{3f3 = h((3). It remains to choose (3 E [a,e) such that ((3,0:*((3))
minimizes (4.71), where 0:*((3) > (3 satisfies K((3,o:*((3)) = 0 for (3 < Sand
0:*((3) = (3 for (3 ;::: S.
By Lemma 2, a necessary condition for that is
L((3,o:*((3)) = l d
[h{3o<.(f3)(x) - h(x)]w(x) dx = 0 (4.73)
(cf (4.65)). It is clear that (4.67) strictly increases from negative L(a) to
positive L(c) = L(c,o:*(c)). We have L($) = 0 if
We show that L((3, 0:* ((3)) is also increasing, when 0:*((3) > (3 is determined
by (4.69). Consider
74 4. Order Statistics of Independent Samples
dL((3, a*((3»
d(3
(4.74)
Plugging
8K(p,a.(p»
00. ((3) = 813
8K(p,a. (13»
d(3
8a.
= h'((3) f:.(p) [X - a*((3)](x - a)w(x) dx >0
[h'(a*((3» -lkA.(p)l f:.(p) (X - a)2w(x) dx
[ [d (x_a)W(X)dx]2}. (4.75)
ia.(p)
The last two lines are positive by the Schwarz inequality and so is the whole
expression in the curly brackets. Because h' ((3) and the denominator are
positive as well, the same holds for (4.75).
/J
We are thus led to the following conclusions. H ~ /J then ((3*, a. ((3.» =
/J)
(/J, is the unique pair satisfying necessary condition (4.73) for minimiz-
ing(4.71). This gives the first statement of Lemma 17. To see that ~ /J /J
coincides with (4.68), we note that (4.66) satisfies
and
[('(f3) = h'(f3) i d
(x - a)w(x) dx > 0.
Therefore [( is incr~asing and has a sin&.le ze!o at p. The same holds for L
which vanishes at f3. Hence conditions f3 ~ f3 and (4.68) are equivalent. If
p < p, then
L(P) = lim_L(f3,a.(f3» > 0,
f3/'f3
and we can make L(f3, a. (f3» smaller by decreasing f3. There is some f3. E
[a, P) that satisfies (4.73), because the opposite contradicts the existence
of the solution. Taking a. = a.(f3.), we see that (4.73) and (4.70) are
identical, and (4.69) holds by the definition of a.(f3.) . •
(4.76)
we have
fd . [1- hnW(x)](x - a)w(x)dx ~ 0, (4.77)
lW-l(f3)
then EFXj:n ~ J1.F for all F b W.
Otherwise there exists a pair (a., f3.), aw ~ f3. < W- 1 «j -1)/(n - 1»,
aw < a. E [f3., W-l«j -1)/(n -1))), determined by equations
hnW(a) - h nW(f3)
a-a
Ia
d(x _ a)2w(x)dx
a-a a
I
+hnW(a) - hn W (f3) d(x _ a)w(x)dx
such that
EFXj:n - J1.F
0'
<
-
B = B O~.w (.}, n) (4.80)
F -
76 4. Order Statistics of Independent Samples
for
B2 = (fj:n W(!3.W [1 13
• w(x) dx + i~ w(x) dx 1
3-(2nnif
(2!-2)(2n-~i)
+ n [F2j-1:2n-1 W(a.) - F2j-1:2n-1 W(!3.)]
!
The equality in (4.80) is attained by
0, it X-IJ<_l-hnW(.B.)
~ q B'
r-1(Bx-1J + 1) it _l-li: nW (i3.) < X-IJ
F(x)= 3:n q ,
~ B - q
< _l-lj:n W(o:.)
W (a + (0:. hn
-a)[B9+l- hn W(i3.)])
W(o:.)-hn W(i3.) ,
it
~
X-IJ
q -
B
> _ 1-
B
'
Ij:n W(o:.)
.
(4.82)
Thivial bounds identical with general ones for the sample minimum are
consequences of
P~.w(h -1) = P: wh -1 = 0
- -'
under (4.68), rewritten as (4.77). These apparently hold for small order
statistics. Equations (4.78) and (4.79) follow from (4.69) and (4.70). Plug-
ging in W = U, V we specify (4.77) to (4.82).
Proposition 8 (decreasing density on the average) If for given 2 ~
j ~ n - 1, and ~ defined in (4.76)
(4.83)
= 0, (4.84)
= 0, (4.85)
and then
(4.86)
4.3 Distributions with Monotone Density and Failure Rate on the Average 77
for
0,
:-l(BX-/f + 1)
f 1:n u '
F(x) =
a.[B T H-/j,n(.8.)]
Ij,n(a.)-lj,n({3.) ,
1,
L
A
a
(1- e-{3 - e: ) fi,nV((3) +(1 +~) e- a fi,nV(a)
-1 + Fj:n V((3) - Fj:n V(a) = 0, (4.91)
such that
EFXj:n - /-LF <
_
B _ BO
-
(. )
)-- V ),n , (4.92)
UF -'
where
78 4. Order Statistics of Independent Samples
'-(2"" iT
(2!-2W"-~;)
+ n [F2j-l:2n-l V(.8*) - F2j-l:2n-l V(a*)]. (4.93)
I
Bound (4.92) becomes the equality for
B '
1 a.[B9+1-I;,,, V(f3.)] zf x-J.! > _1-/;:n V(a.)
- exp (- 1;,,, V(a.)-I;,,, V(f3.) ), q - B •
(4.94)
Distribution functions (4.50) and (4.58) are similar to (4.88) and (4.94),
respectively. The essential difference is that the latter ones have jumps of
height .8., V (.8.), respectively. If
[1 - /i:n V(.8*)]aF
[Link] = BO (.)
tov),n
, (4.95)
where .8k \.t a, and ak are sufficiently large. Alternatively, we can also an-
alyze the norm convergence of hk W- 1 in £2([0,1), dx). For more reasoning
details we refer the reader to Rychlik [89]. Summarizing, we have
4.4 Symmetric Unimodal Distributions 79
EFXj:n - [Link]
O'F
< V3 (~ _
- n +1
1) '
(4.98)
where the equality holds for F uniformly distributed on [J.t - V30', J.t + ..;30'].
80 4. Order Statistics of Independent Samples
(4.99)
where
and (4.100), (4.101), and (4.103) are derived from general formulae by
elementary calculations.
Note that (4.103) has a density symmetric about [Link], a finite support
with uniform ends, and a (principally infinite) peak at the center. The
statement of Theorem 15 is weaker than those of Theorems 11 and 12 in
that we were not able to determine explicitly the pairs (j, n) for which
L(a) > 0, and the resulting optimal bounds are determined by the minimal
distributions W in the class with respect to the order.
In Table 4.1 numerical evaluations of mean-variance bounds are pre-
sented for the jth order statistics, 10 ~ j ~ 20, of i.i.d. samples of size
n = 20, coming from general (G), symmetric (S), symmetric unimodal
(SUN) populations, and those with decreasing density and failure rate on
the average (DDA and DFRA, respectively). For j ~ 9, all bounds are triv-
ial except for the first case. The values of the third column were presented
in Gajek and Rychlik [33, Table III]. Numerical bounds for the DDA and
DFRA samples of size 15 can be found in Rychlik [89].
82 4. Order Statistics of Independent Samples
(4.104)
Both (4.104) and (4.105) are attained by the (limiting) two-point distribu-
tion
P(x) =
(4.109)
4.5 Bias of Quantile Estimates 83
EFX2 :n - F-l(P)
-----..;;;....;.. ~
-B
_-- -B/'(2
_ )
,n,p, (4.110)
O'F
where
1
0, it X-J.! hn(a.)
'J u B'
f -l(-B~) if - hn(a.) < z-J.! < _ hn(2/n)
F( ) 2:n - u ' B - u - B'
(4.113)
X = ~, it
'J
_ hnJ2/n)
l!.
<
-
< nF2:n]2/n)
X-J.!
u (n-2)l!. '
1, it X-J.! > nF2:n (2/n)
'J u - (n 2)l!.
Note that
B(j,n,p) - B(j,n,p) ~ [P(1- p)tl/2. (4.114)
Precisely, we have equality in (4.114) under conditions of (4.105), and strict
inequality holds only in the exceptional cases treated in the last two state-
ments. This implies that the bias oscillation presumably does not depend
on the sample size, is smaller for the central quantiles, and increases to
infinity on the tails. The only positive effect of the sample increase is that
the absolute deviation of the bias decreases to 1/[4p(1 - p)j1/2, because
Fj:n(P) -+ 1/2 by the de Moivre-Laplace theorem. However, for no quan-
tile the bias oscillation tends to zero. This is obvious, because we take into
account the distributions with nonunique pth quantiles (estimation prob-
lems in such cases are discussed in Feldman and Tucker [28]). Removing
such distributions would not help here either, because the rate of conver-
gence of Xj:n to F- 1 (P) depends on the slope increase of Fat F-l(P), and
it is impossible to determine uniform rates without imposing conditions on
the slope (cf., e.g., Zielmski [104]). In fact, all the bounds of Theorem 16
are optimal for the class of strictly increasing F. However, the conclusions
on attainability should be formulated more carefully then.
The proof of Theorem 16 is based on the representation
84 4. Order Statistics of Independent Samples
and Lemma 3. The upper bound (4.104) is the L2- norm of the derivative
of the greatest convex minorant of Fj :n - 1[p,1), and the centered quantile
pieces on [O,p] and [p, 1] with slopes -[l-Fj:n (P)]/p and [l-Fj:n(p)]/(l-p),
respectively. An easy computation leads us to the final claim.
In order to get the lower bounds it suffices to find the greatest convex mi-
norant for 1[p,1) - Fj:n. Distribution function Fl:n is concave, and therefore
the greatest convex minorant of 1[p,1) (x) - Fj:n(x) coincides with -Fl:n(X)
and the straight line F 1:n (P)(x -1)/(1- p) on [O,p) and [p, 1), respectively.
For 2 :5 j < n, the problem is that function 1[p,1) (x) - Fj:n(x) = -Fj:n(x)
°
for :5 x :5 p = j / n is decreasing concave-convex with the inflection point
(j - l)/(n - 1) < p. Two cases are possible: either the straight line join-
ing (0, -Fj:n(O)) = (0,0) with (p, -Fj:n(P)) lies entirely below the graph
of -Fj:n, or the line crosses the graph there. In the former case, the line
becomes a part of the greatest convex minorant. In the latter one, the mino-
rant consists of the line passing through (0,0), and tangent to the graph at
a point 0:* E ((j -1)/(n-1),p), and -Fj:n(x) itself on [o:*,p]. In both cases,
the line passing through (p, -Fj:n(P)) and (1, 1[p,1)(1) - Fj:n(l)) = (1,0)
is the remaining part of the minorant. The problem of settling which case
actually holds is equivalent to checking the sign of the expression
°
We have b(j,n) :5 in the first case, which leads to (4.105), and (4.108)
holds otherwise. A thorough study carried out in Okolewski and Rych-
lik [65] shows that (4.115) is positive for j = 2 with n ;::: 7 only. This
is based on analysis of sign changes of sequences defined as integrals of
a fixed function with sequences of totally positive functions (see Karlin [43]
and [44, Chapter 1] for details). We preserved the dual notation j/n and p
in Theorem 16 on purpose. The reason is that the formulae remain true if
°
p ~ jfn as well. For instance, both (4.104) and (4.105) hold for all j, n with
< p = (j - l)/(n - 1) < 1. In fact, the former is true for all 1 < j < n
with p ;::: (j - l)/(n - 1) and so is the latter for all 1 < j < n with
p:5 (j - l)/(n - 1).
Now we restrict ourselves to the distributions following a given W in the
convex and star orders. We present the upper bounds only. As above, all
are true for arbitrary j,n with p ;::: (j - l)/(n - 1). The auxiliary dual
problem to solve is to maximize the functional
EFXj:n - F-1(P)
over C~cw and c~. w' For convenience, we replace the constant subtrahend
in the orackets by a variable one to make the difference constant. This does
not affect the limit
EFXj:n - F-1(p)
(4.119)
g+ = max{g,O} E Cf':w
lies closer to h, and, by Lemma 2, so does
86 4. Order Statistics of Independent Samples
and define
6 = inf{x~: g(x) ~ h(x)}, (4.120)
setting 6 = d if d < 00 and g(d-) :5 h(d-).
If 6 > c, the straight line leo secant to 9 at e and 6 has a positive slope
g(6) - g(e)
°= 6 -e) ,
satisfies
Jd (x-'"'I.)h(x)w(x)dx
0.=0.('"'1.) = 'Y. 2,(4.123)
J~ (x-'"'I.)2w(x)dx- [J~ (x-'"'I.)w(x)dx]
d
{3.={3.h.) = o.h.)l (x-'"'I.)w(x)dx. (4.124)
'Y.
4.5 Bias of Quantile Estimates 87
(3.(a,,) = l d
[a(x-,)l["!,d)(x)-h(X)]W(X)dX
= a i d
(x - ,)w(x) dx
=
D(a, (3. (a, ,), ,)
= a.(b) ld(X-b)W(X)dX,
88 4. Order Statistics of Independent Samples
EFXj :n - F- 1(P)
Theorem 17 (F tc W) Set
J2 = J~(b) =
=
Varw(X - b)+
-1 AX-JL
F(x) = W(W (P) + 1] + t9-U -)I[tL- i)/J,+oo) (x). (4.131)
A
Note that
is tight and attained by the mixture of the atom at J.' - 0'[(1- p)/(l +p)Jl/2,
and the exponential distribution with scale parameter 0'/ (1- p2)1/2 starting
from the atom, with respective weights p and 1 - p.
a = g(6) + (3 > O.
6-a
Then the function (4.134) is the smallest possible one for arguments x E
(b,6), and the greatest one for x E (6, d) among all starshaped functions
starting from (a, -(3) and passing through (6,g(6)). We have thus proved
(4.121) and (4.122), with (4.119) replaced by (4.134), which is precisely the
assertion of the lemma. •
90 4. Order Statistics of Independent Samples
for
t(x - a)h(x)w(x) dx
(4.135)
Jcd(x - a)2w(x) dx - [t(x - a)w(x) dX] 2'
The proofis similar to that of Lemma 19. We first find the solution (4.136)
to the minimization problem of 119at3 - hW with respect to 13. Then we
look for Ct* minimizing 119at3.(a)-y - hll 2 and find (4.135), which is positive
by arguments analogous to those presented in the last part of the proof of
Lemma 19.
Taking h(x) = hn W(x) for x ~ c, and passing to the limit, we get
and the normalized function 9a.t3. /119a •.8.11 identical to (3.58). This implies
that the distribution functions which attain the bounds in Theorem 18,
and Propositions 12 and 13 coincide with those attaining the bounds of
Theorem 6, and Propositions 5 and 6, respectively.
Theorem 18 (P ~* W) Under the notation (3.60), b = W-1(P), and
EFX j :n - p-1 (P) ::; B~. w(j, n,p) = Tl;i:( (a~~) . (4.139)
(J'F - wa,
Proposition 12 (decreasing density on the average) If (j -l)/(n-
1) ::; p E (0,1), then the following inequality is sharp
EFX j :n - p-1(P) -0
< Bt.u(j, n,p)
2v'3j 1 - Pj+1:n+1 (P)
= (4.140)
n +1 9u(P)
4.6 Open Problems 91
Notations 8u(p) and 8v (p) , that appear in (4.140) and (4.141), were defined
in (3.64) and (3.67), respectively. It is worth pointing out that the upper
bias deviations in large samples are significantly different for the classes
of distributions determined by the convex and star orders. To see this,
suppose that p is fixed and take a sequence of order statistics X j : n , j =
j(n),n -+ 00, such that p ~ (j -l)/(n -1) -+ p. If F?::c W, then (4.130)
holds and, since X j :n -+ b = W-l(p) almost surely, (4.128) tends to O.
On the other hand, (4.129) does not depend on the sample size and is
positive. In consequence, the bound in (4.130) tends to 0, and so do those
of (4.132) and (4.133). If F ?::. W, then (4.138) tends to (b- a)/2 > 0, and
(3.60) remains fixed positive. Therefore the right-hand sides of all (4.139)
to (4.141) have positive limits.
Table 4.2 contains values of extreme upper deviations of order statistics
estimates Xj:n. 1 $ j < n = 20, of quantiles of order p = j In for the fam-
ilies of general (G), decreasing density (~O) and failure rate (OFR), and
decreasing density, failure rate on the average (OOA and OFRA, respec-
tively). In general populations, the estimates are more stable for the central
quantiles than for the extreme ones which confirms the theoretical analysis
of bias oscillation (cf. (4.114». Otherwise the deviations are increasing in j.
It follows from the fact that all the families are more concentrated on the
left, and the upper quantiles are more dispersed and thus more difficult to
be estimated.
TABLE 4.2. Sharp uniform variance bounds on upper bias deviations of esti-
mators Xj:20 of pth quantiles, p = j/20, for various families of distributions
(independent case).
lik [87]. Especially, it was shown that general bounds (5.13) are attained
by the IFRA distributions. In fact, the claim can be extended to families
of distributions defined by F j. W for general W. We also write explicitly
the bounds for F ~. U that have decreasing densities on the average. The
mean-variance bounds for order statistics with symmetric unimodal and
U-shaped distributions of parent variables, described in Section 5.5, were
obtained in Gajek and Rychlik [32] and Rychlik [84], respectively. We call
a symmetric distribution U-shaped if it has nonincreasing and nondecreas-
ing density on the lower and upper halves of its support. Bias of quantile
estimation in dependent samples is studied in Section 5.6. Section 5.7 deals
with the extreme deviation of expected order statistics under violating the
independence assumption.
It is worth pointing out that Papadatos [71] established sharp uniform
bounds on distribution functions of order statistics under the assumption
that the maxima of the subsamples of a given size have a specified common
distribution. The special cases were the independent and arbitrarily depen-
dent identically distributed samples discussed in this book. Balakrishnan
et al. [14] presented the best possible mean-variance bounds on the ex-
pectations of order statistics from the samples taken without replacement
from finite populations. As the size of a population increases to infinity, the
bounds reduce to those of the Li.d. samples from arbitrary populations (cf.
Section 4.1). On the other hand, in the case of exhaustive drawing with-
out replacement, these results coincide with ones for arbitrarily dependent
samples with arbitrary common marginal distributions, presented in Sec-
tion 5.2, which are actually attained by exhaustive sampling models. For
some quantile bounds on order statistics of dependent observations and
respective L-statistics, we refer the reader to Rychlik [84, Section 5].
for Ge being the greatest convex function on the unit interval satisfying
(2.26). Recalling Lemma 1, we obtain
Ep LCjYj:n
j=l
= (00 y (t. C;G',.(dy»)
< 1+
-00
00
y GeF(dy)
= 11 F- 1(x) Ge(dx)
= 11 F-1(x)ge(X) dx
= (F-1,ge), (5.5)
holds. We have
for all 1 ~ i ~ K" and 1 ~ j ~ >. such that ki - 1 ~ lj-l < lj ~ ki. Relation
C = K uniquely determines the distributions of all order statistics that
solve (5.4) and attain equality in (5.5). Then
G-.3: n (X ) -_ nF(x) - k i - 1
,
ki - ki- 1
for ki - 1 < j ~ ki . If C "I K, there are also other solutions. For instance,
for the sample mean we have
Gc(x) = x,
K = {O,n},
C = {O, ... ,n},
and (5.6) simply means that each }j:n should belong to the domain of Y1 .
Indeed,
1 n
Ep- 2:}j:n =
11
F-1(x)dx = [Link]
n j=l 0
relation
P{Yj-l:n ~ F-1{{j - l)jn) ~ Yj:n = Yn :n ) = 1 (5.7)
combined with (5.1) and (5.2) are respective conditions for equality. The
stochastically largest (Le., uniformly smallest) distribution function
G .. F{ ) = [nF{x) + 1 - j]+
,.n X n + 1 -3.
of the jth order statistic is uniquely determined. There are various ways of
constructing dependent samples with the same marginal distribution and
the extreme distribution of Yj:n. The simplest one is a random rearrange-
e: X)
ment of
Y F-
e: 1)
1:n = ... = Yj-l:n = 1 1
< F-
e:
1
< F- 1 1 + [1 - j : 1] X)
= Yj:n = ... = Yn :n
for some X uniformly distributed on [0,1]. For the sample minimum, (2.27)
and (5.7) imply the trivial claims EFY1:n ~ J.l.F, becoming the equality
for identical observations Y1 :n = Yn :n = Y1 • Condition (5.7) excludes the
possibility of constructing an absolutely continuous joint distribution of the
sample with the stochastically maximal jth order statistic except for the
sample maximum.
In the first paper in this field of research Mallows [53] constructed a den-
sity function of the sample with identical uniform marginals that provide
the maximal expectation of the sample maximum. Lai and Robbins [47]
extended the construction to arbitrary possibly nonidentical marginal dis-
tributions. Lai and Robbins [48] and Tchen [98] constructed infinite se-
quences of variables with identical and arbitrary distributions, respectively,
such that all sample maxima are stochastically maximal. Bounds (2.27) for
general order statistics of identically distributed samples were proved in-
dependently in Caraux and Gascuel [19] and Rychlik [76]. In the former,
some inequalities for nonidentically distributed observations were presented
(with conditions for equality established in Rychlik [82]). In the latter, the
problem of constructing sequences with stochastically extreme order statis-
tics was also discussed. Asymptotic properties of sequences of stochastically
extreme maxima and other order statistics were studied in Lai and Rob-
bins [48] and Rychlik [77], respectively.
100 5. Order Statistics of Dependent Observations
In the rest of this chapter we present tight bounds for expected order and
L-statistics of dependent samples for various families of marginal distribu-
tions. We also indicate the marginals for which the bounds are attained.
Formally, in each case one should say that these are attained by the joint
distributions specified for arbitrary marginals and fixed L-statistics by (5.1)
and (5.2) with (5.6) (replaced by (5.7) for single order statistics in partic-
ular), and the given marginal. However, we drop repeating the reference to
the construction of the joint probability and confine ourselves to describing
the optimal marginal.
(5.9)
_ 1 n n r1
d=-Ldj =LCj=Gc(l)= 10 gc(x)dx. (5.10)
n j=1 j=1 0
P dk. - ii)
(y, =J.L+U-- ki - ki-1
1 = , (5.12)
C n
= 101[F-1(X) - J.L][ge(X) - d] dx
ii
is based on (2.25), (5.10), (2.10), and (5.8). The Schwarz inequality provides
the sharp bound here, because Be - is nondecreasing, and so
( U) j-l
P Y1=J.L- C =-;;:-=I-P(Yi=J.L+uC). (5.14)
Bounds (5.13) for the sample maximum and general j were presented in
Arnold [2] and Gascuel and Caraux [34], respectively. We proceed now to
present analogous results for symmetrically distributed random variables.
Below in Theorem 20 we omit the formal presentation of the distribution
that attains the bound, because this needs introducing a rather complicated
notation.
102 5. Order Statistics of Dependent Observations
n
EF "lj·[Link]
L.J Cj· ~ C 8( c ) =
[
1"
-2
nL.J (dj-dnH- j )2
]1/2
, (5.15)
.
3=1 (TF n j= L( n+3)/2 J
where L·J denotes the floor of a number, is tight. Equality is attained by
a unique (up to location-scale tmnsformations) marginal symmetric distri-
bution supported on n points at most whose probabilities are multiplicities
of l/n.
PROOF. Inequality (5.15) follows from
EF
j=l
n
2: Cj(lj;n - J.L) = I 1/2
1
[F-1(X) - [Link]~(x) dx (5.16)
with
gc(x) - gc(1- x-)
g~(x) = 2
n
= "L.J dj - d2nH - j 1[(j-1)/n,j/n) ( x ) , (5.17)
j=L(n+3)/2J
1/2 ~ x < 1, and the Schwarz inequality. Note that (5.17) is a nonnegative
non decreasing piecewise constant function with L(n + 1)/2J values at most.
Thus its antisymmetric extension onto [0,1/2) is also a nondecreasing step
function with n values at most and jumps at some points j In, 1 ~ j ~
n - 1. The extension coincides with the quantile function of the extreme
distribution providing equality in (5.15) up to an affine transformation.
This justifies the latter claim of Theorem 20.•
is sharp and becomes the equality for the three-point marginal distribution
P(Y1=J.L) = 2Ii-1_!1
n 2'
(5.20)
(5.21)
5.2 General and Symmetric Distributions 103
Inequality (5.19) and its special case for j = n can be found in Gascuel
and Caraux [34] and Arnold [2], respectively. Here we merely mention the
second moment bounds for general L-statistics
< II(gc)+llmF
EFlj:n ~ ( ~
n+ -)
.)1/2 mF (5.23)
t
j=1
Cj (Yj:n
S
- y) < [
~tcP;-cP
3=1
]
1/2
, (5.27)
t Cj{Yj:n - y) < [
2~ t (dj - dn+1_j)2
]
1/2
, (5.28)
j=1 S j=L(n+3)/2J
1/2
Ln
CjYj:n < [
~ t,{dj)~ ]
, (5.29)
j=1 m
for
J~({3) =
=
Varw(X - (3)+
X - JL)
F(x) =W (
[3* + T}W([3*) + t1w([3*)-u-
A A
(5.38)
which becomes the equality for F being the uniform distribution function
on fJ.t - V3u, JL + V3u].
n)
Otherwise
EFYj:n - JLF < ~ (9 j - 9 - 1/2 (5.39)
UF - 3 n+ 1- j ,
and the equality holds for the mixture of the Dirac distribution at JL -
3u [en + 1 - j)/(9j - 9 - n)]1/2 and the uniform distribution on
[JL -
3 (n + 1 - j )
u 9j _ 9 _ n
1/2 3j - 3 + n
,JL + U [en + 1 - j)(9j - 9 - n)]1/2
1
with respective probabilities 3(j - 1)/(2n) - 1/2 and 3(n + 1 - j)/(2n).
Proposition 15 (decreasing failure rate) If (j - 1)/n ~ 1 - e- 1 ~
0.63212, then
EFYj:n - JLF < I n
(5.40)
UF -
n
n+ 1 -J.,
and the equality holds in (5.40) for the exponential distribution with location
JL - U and scale u.
Otherwise for
'Y ='Yv(j, n) = ( 1 - j -1)
-n- e,
we have
EFYj:n - JLF
UF
< C =d
- toV
(. ) = [
J,n
2 -1]
'Yv(j,n)
1/2
, (5.41)
(5.42)
5.3 Distributions with Monotone Density and Failure Rate 107
Note that (5.42) is the mixture of the exponential distribution with location
J.t - ulO and scale ul(-YO), with probability;, and the Dirac measure at
J.t - u I 0 with probability 1 - ;. This is a D FR life distribution if J.t = u I 0 .
For small order statistics (5.40) is attained by a life exponential distribution
with J.t = u.
JOw we consider analogous bounds for F ~c W. First we try to describe
P-jcwh for
n
h(x) = n+1_ j l[W-1((j-l)/n),dw) (x)
1
= 1 - W(b) l[b,d)(X).
_ _
a. - a.(,8.) - d d}
f: {/3
.
w(x) dx
• (5.45)
fb w(x) dx famln ., (,8. - x)w(x) dx
108 5. Order Statistics of Dependent Observations
Precisely, if
d~
,Q
= a
t X2W(X) dx - taxw(x)bdx I: xw(x) dx '
t w(x) dxbad (5.46)
Ia w(x) dx Ia xw(x) dx - Ia xw(x) dx
fJ
(Le., the left-hand side of (5.44) is not greater than the right-hand side for
(3 = d), then (3. = /J, and (5.43) is linear. This may happen if d is finite
and b is close to d. A probabilistic interpretation of Lemma 23 is given in
Theorem 22. For a < (3 < d, we define
· [· r
192 = 19~ «(3) = Varw«(3 - X)+
such that
Distribution function (5.54) jumps from W(.8*) > (j -1)/n to 1 at its right
endpoint. H W = U, both cases are possible. Since dv = +00, the first one
is only applicable for the IFR distributions.
Proposition 16 (increasing density) If (j - 1)/n < 2/3, then
EFlj:n - JLF < [(j - 1)(8n - 9j + 9)]1/2
(5.55)
UF - 3(n + 1 - j) ,
which becomes the equality for the combination of the uniform distribution
on
4n - 3j + 3 [ 9j - 9 ] 1/2]
[ JL - U [(j _ 1)(8n _ 9j + 9)]1/2' JL + U 8n - 9j + 9 '
and an atom at the right end of the interval, with respective coefficients
3(j -1)/(2n) and 1- 3(j -1)/(2n).
If (j - 1)/n ~ 2/3, then (5.38) holds with the equality for the uniform
distribution on lJL - J3u, JL + J3u].
Proposition 17 (increasing failure rate) For {3* > In[n/(n + 1 - j)]
defined by the equation
1 - (1 + {3)e-.8 n +1- j n
ln (5.56)
e- - +
.81{3= ·
3- l 1·
n+-3
and
(5.57)
we have
EFlj:n - JLF < j - 1 (Jv({3*) (5.58)
UF - n + 1- j e-.8· -1 + {3*.
The equality in (5.58) holds if
0,
F(x)= { l-exp(-I+e-.8·-(Jz~,,),
1,
(5.59)
This is the exponential distribution with location JL - u(1 - e- fJ .) / (J and
scale U /(J, right truncated at JL + u(e-.8· -1 + {3*)/8. The jump probability
is e- fJ ••
Likewise, for determining respective second moment bounds for life dis-
tributions F tc W and F ~c W, we need to find best approximations
of
n
h(x) = n + 1 _ j 1[(i-1)/n,1) W(x)
110 5. Order Statistics of Dependent Observations
In the first case, under the conditions of the latter statement of Theorem 21,
the optimal parameter 13 defined in (5.33) is positive which implies that
PSwh E C~"w provides the solution to our problem. Then we have
fbdxw(x) dx JLw(b)
(5.62)
= mw '
d
fb w(x) dx
[
fod
x 2w(x) dx
] 1/2
respectively.
Theorem 23 (F !:c W) If for b = W-1«j -1)jn)
(5.63)
EFYj:n JLw(b)
--<
- - -, (5.64)
mF mw
F(x) =W (m :"F)' w
we have
EF}j:n < VW({3.)
(5.66)
mF - J.,tW({3.) - (3. '
which becomes the equality for
Again, we point out analogies between Theorems 21 and 23. Under equiv-
alent conditions (5.30) and (5.63), we obtain analogous bounds (5.31) and
(5.64), respectively. In the opposite case, the bounds are determined by
the same parameter {3. and related by (5.60). Moreover, both the extreme
distributions (5.37) and (5.67) are location-scale modifications of W with
identical jump W ({3.) at the left support end.
Proposition 18 (decreasing density) If (j - l)/n ~ 1/3, then
EF}j:n
mF
~ J3
2
(1 + 1) ,
j -
n
(5.68)
ld ld
If
d xw(x) dx > x 2w(x) dx, (5.73)
Otherwise
(5.75)
(5.77)
say, and
ptwh(x) 1 . {x }
(5.78)
IlPtwhll = :y mm a*' 1 .
If (5.73) is not true (which is possible for dw < 00 only), then the projection
is actually linear and satisfies
(5.79)
ptwh(x) x
(5.80)
IIP~whll
_c
Theorem 24 (F ~c W) If
2
n mw d (5.81)
n + 1 - j JLw(b) < w
for b = W- 1 «j - l)jn), then there exists a unique b < a* < dw solving
W(a)tLw(a)[a - tLw(a)J = W(a)a~(a) + [(j - l)jnJatLw(b) (5.82)
5.3 Distributions with Monotone Density and Failure Rate 113
= ;;lv(a.) =
ld
;;2 Ew(min{X,a.})2
= lo
a•
x 2w(x) dx + a~
a.
w(x) dx (5.83)
we have
EFYj:n < n ;;w(a.)
(5.84)
mF - n+ 1-j a.
The equality in (5.84) holds for
....L<~
mF D, (5.85)
....L>~
mF - D'
!(1 - e- a )
a
- e- a = (1 _ j - 1) (1 -In
n
n
n+1-J
.) (5.89)
114 5. Order Statistics of Dependent Observations
and
(5.90)
we ha1Je
EFY;:n < n iiv(a.}
(5.91)
mF - n+1-j a.
The equality holds iff
0,
F(z} = { 1 - exp (-ii ':F) (5.92)
1,
which is the equality iff Y 1 is uniformly distributed on [J-t - V3u, J-t + v'3u].
116 5. Order Statistics of Dependent Observations
(5.100)
is sharp, and this becomes the equality for the combination of uniform dis-
tribution on
B() B ()
gj:n X = n +j-1
1- jgn+1-j:n x ,
and the same holds for the respective projections and their norms. This
explains the relation between the bounds in both theorems and the fact
that the distributions with extreme expectations of jth smallest and largest
order statistics are identical. We can immediately conclude similar bounds
for some L-statistics.
In contrast with the independent case, all the bounds are nontrivial for
any order statistic except that of the sample minimum. Table 5.1 con-
tains numerical values of mean-variance bounds for order statistics from
dependent samples of size n = 20 coming from the following families of
distributions: general (G), symmetric (S), symmetric unimodal (SUN) and
U-shaped (SUS), distributions with decreasing density and failure rate (DD
and DFR, respectively), and with decreasing density and failure rate on the
average (DDA and DFRA, respectively). In fact, all the bounds depend on
TABLE 5.1. Sharp uniform mean-variance bounds on expectations of order statistics from dependent samples
of size 20 for various families of distributions.
j G S SUN SUS DD DFR DDA DFRA ID IFR
2 0.22942 0.16644 0.15692 0.08660 0.08660 0.05129 0.09028 0.05250 0.21558 0.21583 Q1
3 0.33333 0.24845 0.23424 0.17321 0.17321 0.10536 0.18543 0.10999 0.31208 0.31285 g,
4 0.42008 0.32219 0.30376 0.25981 0.25981 0.16252 0.28232 0.17247 0.39167 0.39325
5 0.5 0.39529 0.37269 0.34641 0.34641 0.22314 0.37897 0.23995 0.46398 0.46673
6 0.57735 0.47141 0.44444 0.43341 0.43301 0.28768 0.47458 0.31254 0.53287 0.53725
7 0.65465 0.55328 0.52164 0.52398 (")
.....
.....
~
118 5. Order Statistics of Dependent Observations
value (j-1)/n only, because they were derived by means of projecting func-
tions of (j-1)/n. Therefore all estimates for Y;:n hold also true for any Yk:m
provided that (k - l)/m =
(j - l)/n. In particular, the moment bounds
for the sample medians Yn :2n+1 do not change under increase of the sample
size. Moreover, one can check that all the bounds for restricted families of
distributions presented in this chapter increase continuously in (j - l)/n.
Consequently, Table 5.1 may provide fair approximations of bounds for Yk:m
from various families of parent distributions with (k-1)/m ~ (j -1)/20 for
some 2 ::; j ::; 20. Some tables for second moment bounds were presented in
Gajek and Rychlik [33]. Also, the third column (SUN) of Table 5.1 comes
from that paper. The numerical results for distributions with monotone fail-
ure rate and monotone failure rate on the average were earlier presented in
Rychlik [87].
p-l(P) + sup
PE'P" (F)
Ep(-Xj:n) = r P-l(x)(l(p,I) - G_e)(dx),
10
1
with
(j/'(j,n,p) = n (1_p)I/2
n+1-j p
The upper bound in (5.101) is attained by the two-point marginal distribu-
tion (4.106). The lower one is attained by the three-point distribution
P (X=~-U [1_P~0P_1)/nr/2) =
j-1
n
j-1
P(X =~) = p-~,
(j _ l)/n ] 1/2)
P (X
[
= ~+U1_p+ (j -l)/n = 1-p.
I!/'(j,n,p) j - 1 ] 1/2 1
C/'(j,n,p) = Fj:n(P) [ p(n + j _ 1 _ np) ~ 2' (5.102)
B-/' (j,n,p) n+1-j 1
(j/'(j,n,p) = [1 - Fj:n(P)] n(l _ p) ~ 2' (5.103)
o < - r (x -
J-y. -y.)hj:n,p,q(x)w(x) dx
=
n(l - q + p) + 1 - j
(n + 1 _ j)(q _ p)
r
J-y. (x - -y.)w(x) dx
as q '\t p and c = W- 1 (q),'Y* '\t b = W- 1 (P). It follows that the bound for
F tc W has the form of (4.126) with hnW(x) replaced by gj:nW(x) =
n1[b,d)(x)j(n + 1- j). For F t* w, we have (4.137) with hnW replaced
by gj:n W. In both cases the normalized projection does not depend on the
specific form of h. Therefore the upper bounds in the dependent case are
attained by the same marginal distributions that provide the analogous
bounds for independent observations. Note that the same holds when F is
general.
Theorem 30 (F tc W) If (j -l)jn ~ p < 1 and b = W- 1 (P), then, with
notation (4.127) and (4.129), we have
EpYj:n - F- 1 (P)
< n 1/w(b)
--=........::..;.;.;.....--=-'- (5.105)
- n + 1 - j 'i?w (b)
A '
EpYj:n - F- 1 (P) -0 ( )
< C t • v j,n,p
Up
-0
Ct.v(j,n,p)
1 -p
}
1~
n {
= n+1-j 1 + p[l-1n(l-p)J2
(5.110)
B~.u(j, n,p) p
-+ 2(1 + p)'
(5.112)
C~.u(j, n,p)
B~. v(j, n,p) -In(l-p)
-0 ) -+ 2[1 -In(l - p)]'
(5.113)
Ct.v(j,n,p
TABLE 5.2. Sharp uniform variance bounds on upper bias deviations of estima-
tors lj:20 of pth quantiles, p = i/20, for various families of distributions (depen-
dent case).
j G DD DFR DDA DFRA
1 4.35890 1.57425 0.95119 1.71534 0.99748
2 3.15789 1.51700 0.95214 1.75672 1.04201
3 2.64497 1.47348 0.95925 1.77758 1.08589
4 2.35294 1.44088 0.96058 1.78338 1.12918
5 2.16506 1.41737 0.96825 1.77967 1.17204
6 2.03670 1.40175 0.97840 1.77127 1.21477
7 1.94681 1.39329 0.99127 1.76200 1.25780
8 1.88422 1.39159 1.00716 1.75471 1.30171
9 1.84257 1.39655 1.02647 1.75159 1.34727
10 1.81818 1.40836 1.04973 1.75437 1.39544
11 1.80907 1.42749 1.07763 1.76456 1.44744
12 1.81444 1.45479 1.11111 1.78373 1.50480
13 1.83450 1.49150 1.15142 1.81369 1.56953
14 1.87044 1.53947 1.20024 1.85672 1.64422
15 1.92450 1.60128 1.25988 1.91583 1.73235
16 2 1.68034 1.33333 1.99488 1.83837
17 2.10042 1.78017 1.42374 2.09772 1.96708
18 2.22222 1.89832 1.52944 2.22108 2.11701
19 2.29416 1.97386 1.60128 2.29389 2.22208
5.7 Extreme Effect of Dependence 123
by many authors (we refer merely to monographs of Huber [41] and Hampel
et al. [37] and references given there). For standard i.i.d. parametric mod-
els, many violations of marginals were studied. Dependence-robustness for
location models was analyzed in Rychlik [80]. The projection method pro-
vides tools for determining expectation sensitivity of arbitrary L-statistics
against dependence of observations
n n
t
PE'P.. (F) j=1 j=1
= 11o
p-l(X) [Gc(X) -
3=1
cj/;:n(x)] dx, (5.114)
(5.115)
is actually increasing and thus p? h1: n = h1:n . Another simple solution
emerges for the sample maximum. Then
_nxn-l, if 0<x<1_1
hn:n(x) = { n (1- x n - 1) , i f 1_1<x<1
- n' (5.116)
n - - ,
has antiderivative
if 0 < x < 1 - 1
H (x) - { _xn
' - - n' (5.117)
n:n - n(x -1) + 1- x n , if 1- ~ :::; x:::; 1,
which is decreasing on [0, I-lin], increasing on [1- lin, 1], and concave
on both intervals. Therefore
and
P?hn:n(x) = { -( - r-)n
(I 1 t- 1
'
if 0<x<1_1
- n' (5.119)
n1- n , if 1- 1n <x<1.
- -
124 5. Order Statistics of Dependent Observations
Lemma 25 Set
If hj:n«j -1)/(n -1)) < 0, then there exist unique r E «j -1)/(n -1), 1)
such that hj:n(r) = 0 and s E [(j - 1)/(n - 1),r) such that either s =
(j - 1)/(n -1) if Sj:n((j -1)/(n -1)) ::; hj:n((j -1)/(n -1)) or s is the
solution to Sj:n(x) = hj:n(x) otherwise.
If
hj:n ( nj-l)
_ 1 < 0 and Sj:n(s) ::; Sj:n (j-l)
-n- , (5.122)
then
p/' hj:n{x) = {Sj:n(s), II:ff 0::; x < s, (5.123)
hj:n(x), s ::; x ::; 1.
If
For the prevailing number of cases, we have (5.122) with s > (j -1)/(n-l).
Then (5.123) is continuous, and can be written as
EP Y 1:n - EpX1:n
sup
PEPn(P) ap
= A = A(I)
L..l. L..l.
n - 1
,n = ----:-
(2n - 1)1/2
(5.126)
F(x) = 1- -;
[1 ( n- x-p,
1 - (2n _ 1)1/2 -a-
1 )] 1/(n-1)
,
(5.128)
for
(5.129)
t1 2 =
n
j _ 1 F j :n
2 (j~
- 1) + ( t-
j -
~ 1) [n
+ n 1 _ j - /i:n{t)
]
2
n 2{1- t) 2n
+ (n + 1 _ j)2 n +1_ j [1 - Fj:n{t)]
2i-2 2.. -2i)
+ n ·-1 .. _ .. - · [1 - F2j-l:2n-l (t)]. (5.131)
EXj :n - EYj:n
A -(.
u. J,n
)
= sup
PE1'.. (F) UF
-In
. f EYj:n - EXj :n
FE:Fp UF
EF-Y';+l-j:n - EF-X;+l_j:n
= sup
F-E:F UF
= t1(n + 1- j,n). (5.133)
5.8 Open Problems 127
Accordingly, we conclude that the lower bound on the deviation of the ex-
pectation of the jth smallest order statistic under dependence is identical
with the upper bound on the deviation of the jth largest one. Similar ar-
guments explain the coincidence of (5.126) with the mean-variance bound
for the maximum in the i.i.d. sample {see (4.2)), and the fact that these
bounds are attained by distribution functions of variables mutually sym-
metric about J.t (see (4.3) and (5.127)). This is a consequence of relations
sup EFY1:n - EFX1:n = [Link] - EFX1:n = EF-X:L;'n - [Link]-' (5.134)
PE'P.. (F)
and both the suprema in (5.135), taken over all P E 'Pn{F) are achieved
by the same distribution function.
Numerical calculations based on Theorem 32 show that for general margi-
nal distributions central order statistics are more robust under dependence
than the extreme ones. We can also see that each jth smallest order statistic
is mOre sensitive than the respective jth greatest one. In fact, Rychlik [88]
measured dependence-robustness of order statistics in terms of scale pa-
rameters generated by central absolute moments of order 1 ~ p ~ 00. Nu-
merical analysis shows similarity of conclusion for cases p = 1 and p = 2,
which were presented above. For p = +00 (i.e., for bounded observations),
the conclusions are diametrically opposite. An interesting fact to note here
is that (5.135) holds true for the classes of marginal distributions with finite
pth moment for arbitrary 1 ~ p ~ 00.
2. From the results for single order statistics we can immediately con-
clude bounds for some L-statistics. For instance, we have
1
L ¥i:n
k
= sup
PE'P,.(F)
EF k
+ -
1
J
.
l=j
(5.136)
j 1
L
k-1
= sup EF [ -Y}:n +- ¥i:n
PE'P,.(F) n n 1=j+1
n+l-k ]
+ n Yk:n, (5.137)
=
Therefore we obtain bounds for spacings Y}:n - y}-l:n in standard de-
viation units by multiplying all respective mean-variance bounds for
single jth order statistics of various families of sample distributions
by factor n/U - 1). In all these cases, the bounds are attained by
the same joint distributions. This is not so for the differences of other
pairs of order statistics. Then, due to (2.25) and (2.26), we have
PROOF. Fix n, and take n random variables Ui, which are uniformly dis-
tributed on intervals [(i-l)/n, i/n], 1 ~ i ~ n. Here Ui may be independent
or arbitrarily dependent, for example generated for a single uniform random
variable Uo by means of transformations
.
U . -_ i - I
n
+ Uo
, 1~ i ~ n.
1 ~ i ~ n. (6.1)
Let (Y1 , •.. , Yn ) coincide with a randomly chosen vector Ui. We easily
check that lj, 1 ~ j ~ n, are uniformly distributed on [0,1]. Now extend
the sequence by adding i.i.d. random variables lj, j > n, with the same
uniform distribution, and independent of (Y1 , .•. , Yn). If (Yi, ... , Yn ) = Ui
for 1 ~ i < n, then L1 = i + 1, and
R1 = Un E[1 - ~,1] .
In order to complete the proof we notice that otherwise L1 > nand
R1 E (Un, 1] C [1 - ~, 1] .•
6.2 General and Symmetric Distributions 133
EFRn-ILF
< -1
{I I
[-In(1-x)] 2n dx-1
f/2
: 10 hn(x)dx-1f/2
rTF (n!)2 0
= [C) I
= [(~) - l f / 2 (6.2)
(1 - x)f~k)(x) = 1- F~k)(x)
k
Ln -:-d-In(1
j
= (1 - x)k - x)]j, (6.4)
j=O J.
where f~k)(x) and F~k)(x) are the density function defined in (2.40) and
the respective distribution function. Then
(6.5)
for
(6.7)
EFRn
mF
< c:) 1/2, n ~ 1, (6.8)
EFR~k)
< 11!~k)lI, n ~ 1, k ~ 2 (6.9)
mF
6.2 General and Symmetric Distributions 135
0, if ~ $ 0,
F(x) = { (f~k»)-1(IIfi,k)II~), if 0 $ ~ < f,1~~)i,)' (6.11)
1 if a: > f~k)().)
, m - IIl~k II '
respectively. Both (6.10) and (6.11) define life distributions. Formally, in-
equalities (6.8) and (6.9) hold true for n = 0, and they amount to 1 then.
They are attained by degenerate random variables concentrated at mF > O.
Obviously, this can be approximated by life random variables in the mean
square with an arbitrarily desired accuracy.
Since
Sn(x) = fn(x) - fn(l- x), x E [~, 1) ,
is the difference of increasing and decreasing functions with identical values
at 1/2, this is strictly increasing from 0 to 00. Evaluating expectations of
the first records in symmetrically distributed sequences we simply use the
Schwarz inequality
2D2 = /1 s;(x) dx
C:) -(n~)2 10
1/2
1
= Inn x Inn(1 - x) dx. (6.13)
which has a smooth symmetric density with the infinite support and sym-
metry center J1..
Raqab [74] also considered kth records in symmetrically distributed pop-
ulations using the greatest convex minorant approach. He proposed numeri-
cal evaluations that substantially improved bounds for general populations
(6.5) and (6.6), and nonsharp ones derived in Grudzien and Szynal [35]
based on direct application of the Schwarz inequality. However, Raqab's
paper lacks a theoretical justification for use of specific constructions of
the greatest convex minorant and precise description of the scope of appli-
cability. Therefore we do not present details here.
136 6. Records and kth Records
V(n - 1), if k = 2,
b = { V (n(2k-s)+[n 2+4 n(k-1)(k-2)]1/2) 'f k> 3 (6.16)
2(k-1)(k-2) ,1 _,
6.3 Life Distributions with Decreasing Density and Failure Rate 137
1 f~)V(y)
00
dy = T L f!::)V(x).
-z
m=O
n
L amfM.)(x),
n
fa(x) = x E (0,1),
m=O
does not exceed the number of sign changes of the sequence a = (ao, ... , an).
The first and the last signs of the sum are identical to the signs of the first
and last nonzero elements of a, respectively.
Note that it suffices to prove the claim for combinations of
·m+1
g(xjj,m + 1) = fM.)V(x)e- = J m!Z x me- jZ l[o,oo)(x) (6.17)
°
every positive integer i, and arbitrary nondecreasing sequences < Xl <
... < Xi, < VI < ... < Vi, we have
°
densities g(Xjj, v), v > 0, is strictly totally positive of order ooj that is, for
138 6. Records and kth Records
(see, e.g., Karlin [44, Chapter 1]). The variation diminishing property for
linear combinations of density functions of totally positive families is stated
in Karlin and Studden [45, Corollary 4.1].
A thorough analysis of (4.41) and{4.42) enables us to establish bounds
described in the following two theorems. We point out apparent analogies
with Theorems 11 and 12. Here we take into account k ~ 2 and n ~ 1, as
the remaining cases k = 1 and n = 0 were discussed above.
Theorem 35 (decreasing density) If (1 l/kt+l : ; 3, then
+
E F R(k)
mF
n <v'3 [
-
1-( _k_ )n+l]
k+1'
(6.18)
E R(k)
F
mF
n ::; D = D~cu(k,n)
-
= 11(J~k»o< •.[Link], (6.19)
lr
for
11(J~k»o<•.8.W = k(2k ~ n
+l C:)Fi~k-l)(,8*)
+ (1- ,8.)[f~k)(,8.W + (1- ,8.)2a.f~k){,8.)
+ {I - ,8.)3a ;/3, (6.20)
(6.21)
where ,8. < b (the first inflection point of f~k) defined in (6.16») is the
smaller of the smallest positive zeros of combinations
- kfn(k)
- 2 (x) +(k- 2 5) f n(k)
- 1 (x ) , (6.22)
6.3 Life Distributions with Decreasing Density and Failure Rate 139
Lu(x) = ~~
2k m=O
[3 (_k_)n-m -1]
k+ 1
,(Ie)(x)
m
_
4k(k
(k - l)(k - 2) lie) ( )
+ 1) n-l X • (6.23)
(6.24)
Theorem 36 (decreasing failure rate) If n $ 2k - 1, then
EFmale ) n+ 1
--"--- < - - (6.25)
mF - V2k'
which becomes the equality for the exponential distribution with the scale
parameter mF/V2.
Otherwise, under notations (4.34) and (4.41), we have
for
+
k(_k_) (2n)
n
2k - 1
2n+l F.(21e-l) ( )
2n 'Y. , (6.27)
n +1 (Ie)
= 2k(1 - 'Y.) [1 - Fn+l h'.)]
where 'Y. < V«n - .fii)/(k - 1)) is the minimum of the smallest positive
zeros of
Kv(x) = ~ n - m,(Ie)( ) _
~ 2k2 m x
1- k 3 ,(Ie) ( )
k2 n-2 X
m=O
+ (k - 1)2(2k - 1) ,(Ie) ( ) (6.29)
2k2 n-l X ,
140 6. Records and kth Records
L ()
VX
= ~
L...J
2k-n+m,(k)( )_ (k-1)2,(k) ( )
2k2 m X 2k2 n-l X . (6.30)
m=O
The equality holds in (6.26) if
0, if .at. < 0 (10)
m -,
F(x) = (f~k»-l
{ (D;i) , if O::;;i::; in A-r·) ,
1 - (1 - , .. ) exp ( - D:;_!:Ie)(-y.»), if ;i ~ i~Ie~-r.).
(6.31)
We see that the extreme expectations of kth record values are attained by
the uniform and exponential distributions in Theorems 35 and 36, respec-
tively, if n is small. This is true for n = 1 and all k in particular, which
is an immediate consequence of the linearity of respective projections im-
plied by lacking a convex piece in the left parts of i1k) (x) and i1 k)V(x).
Distribution functions F$;!) (x) of the jth record values appearing in (6.20),
(6.21), (6.27), and (6.28) have alternative representations
F$;!){x) = G(-ln(l- x)jj,m + 1),
where the right-hand side denotes the gamma distribution functions with
densities defined in (6.17), composed with the standard exponential quan-
tile function. Absolutely continuous distribution functions (6.24) and (6.31)
have the same form up to a multiplicative factor on the left, and uniform
are crucial for determining its projection onto various convex cones. We
easily check that each CPn integrates to 0, starts from the origin, decreases
to
cp(1 _ e-n+l) = _ (n - l)n-1 ,
n!
and increases to +00 at 1 passing through the horizontal axis at 1 - e- n .
Its antiderivative «Pn needed for the Moriguti [58] projection, is therefore
concave decreasing, convex decreasing, and convex increasing in [0,1 -
e-n+l], [1 - e-n+l, 1 - e-n+l], and [1 - e- n , 1], respectively. This vanishes
°
at and 1, and is negative in the middle. Thus we deduce that its greatest
convex minorant <I>n is linear in [0, a*] for some a* E [1_e- n+l, 1- e-n+l],
that is determined by equation
XCPn(X) = «pn(x) = Fn(x) - Fn-1(X) = fn(x) (6.34)
(cf. (6.4)), and coincides with «Pn elsewhere. Finally,
o -,
P CPn(x) = «pn(x) = cpn(max{a*,x}).
Theorem 37 (general distributions) For n 2:: 2 we have
where
F(x) = CPn-1 ( x
~-O"- -11-) l[I'-O"rpn(a.)/~,oo)(x). (6.38)
for 1/2 < x < 1. This is necessary for application of the greatest
convex minorant construction. Furthermore, we could ask for respec-
tive stricter bounds for symmetric unimodal distributions. Note that
these are unknown in the case of first records as well.
2. What are the mean-variance bounds for kth records of Li.d. sam-
ples with decreasing density and failure rate? For distributions with
increasing density and failure rate neither second moment nor mean-
variance bounds are known. A more general problem consists in estab-
lishing analogous results for distributions preceding and succeeding
a given one in the convex order.
3. What are the respective evaluations for the record values coming from
populations with distributions determined by the star order relations
with a fixed one?
4. Results of Section 6.4 should be completed by considering increments
of kth records R~k) - R~k21 in various classes of parent distributions.
It is of interest to compare kth record increments for different k.
5. We have practically no evaluations for the lower records. Generally,
they cannot be concluded from analogous results determined for the
upper records.
6. Establishing bounds on EF(Rn - Xn:n) in various families of distri-
butions we evaluate the rate of increase of the gap between the ele-
ments of nondecreasing sequences of sample maxima and their strictly
increasing subsequence. Analyzing EF(~k) - X n+1-k:n), we derive
similar evaluations for the kth largest order statistics.
7
Predictions of Order and Record
Statistics
)_ 1
1 [F-l(X) - y]2 dx
2 _ E ((X _ )21X _F_--,l(~y)~~.,.....,...._ _ (7.1)
mFI~ - F Y >Y - 1 _ F(y) ,
where y is the value of a previous observation. We have
Jo1[F-l(x) - y]2 dx
2
mpo
III
< 1- F(y)
EF(X _y)2
= 1- F(y)
u} + (IJ.F - y)2
= 1- F(y)
= M;o =M,
III
(7.2)
that is, F has a jump at the left endpoint y of its support. It occurs that the
optimal constants for bounds expressed in terms of m ~o are identical with
. III
those measured in larger units M ~o • These have a more intuitive meaning
III
depending on moments and value of the parent distribution at a fixed ob-
servation point, and allow us to determine distributions attaining bounds
uniquely up to the three mentioned parameters. The former are indepen-
dent of the distribution at points preceding y, and so admit an ambiguity
in the description of extremal distributions. Therefore we choose represen-
tations of bounds in terms of M ~o -units. Results for general distributions,
III
and ones following the uniform and exponential distributions in the convex
order, are stated in Sections 7.1 and 7.2, respectively. Formal proofs will
be published in Rychlik [91].
o F(x + y) - F(y)
1'111(x) = F11I(x + y) = 1 _ F(y) , x ~ 0, (7.4)
(7.5)
form the convex cone C+ of nondecreasing elements of L2 ([0, 1), dx) starting
from 0. Combining (7.3) and (7.4) with (2.23), and further applying the
projection method, the Schwarz inequality, and (7.2), we obtain
7.1 General Distributions 147
For i +2 ~ j ~ n - 1, we have
if ~ < 0,
°
- M
·· .(~)
if < ~ < f J-"n-. 1=1>
Ilfi-i:n-ill'
.. . ( .2i.::.2. )
if ~ > f J-"n-. 1-p
M - IIfi-i:n-ili
(7.10)
Finally,
holds and becomes the equality for the combination of degenerate and power
distributions
0, if z'M" < 0,
F(x)= { p+(I-p)[ 1. Z-"]l/(n-i-l) i~ 0< Z-" <v'2n-2i-l
v'2n-21-1 M ' :/ - M =-_:::-;----::-'
1, if z'M" ~ v'2n - 2i - 1.
(7.12)
All the distribution functions (7.8), (7.10), and (7.12) have jumps at y.
Also, (7.10) has an atom at the right end of its support, and (7.12) does not.
We emphasize the fact that basic relation (7.3) and, accordingly, issuing
inequalities (7.7), (7.9), and (7.11) hold true for continuous distributions F.
They are the best possible, but attainable merely in the limit by continuous
distributions tending to (7.8), (7.10), and (7.12), respectively, in the sense
described earlier (see, e.g., the last two paragraphs of Section 2.3). The
same reservations concern the other results of this chapter.
EF(R(k)
n
- R(k)IR(k)
m m
= y) = EFoIv R(k)
n-m-l> n > m, (7.13)
which becomes the equality for the mixture of Dirac and Weibull distribu-
tions
(F1l/S) 1 + 1) 3 - 11,
o -1 (x) = g[(x O~x~1.
The composition
(7.18)
whose convexity is implied by that of p-l W(x). Note that the reversed im-
plication is false, because we can extend FI~ and WI~ on the left to y and z,
respectively, in an arbitrary way. Compositions (7.18) form the convex cone
+ 2 ( w(z +x) )
CtcWI~ C L [0, dw - z), 1 _ W(z) dx .
and
(cf. (7.1)). Equalities in (7.19) and (7.20) hold if P-l(X) is constant for
x E [p, 1). Convexity and nondecrease of P-1W(x) impose the single value
of the quantile function on the whole domain. The conditional expectations
have a practical sense if y is the only support point of P, and they are equal
to 0, but all the evaluations are meaningless then.
Some results of this section can be immediately concluded from the pre-
vious one. For instance, (7.12) is a decreasing density distribution function
7.2 Distributions with Decreasing Density and Failure Rate 151
which implies that (7.11) is the best possible bound for the restricted class.
In the remaining cases 3 ~ i + 2 ~ j < n we have
EF(Xj:n - Xi:nIXi:n = y)
= 1 1 P
- [F-l(x + p) - y]!i-i:n-i (1 ~ p) /'~ p
~ 10f1-P[F-l(X + p) - y]P: u!i-i:n-i
_c
(-1
x )
p -
dx
I-p
< lIP: u!i-i:n-illm~oIII
_c .
< lIP:
_c
u!i-i:n-iIlM~o.
III
(7.21)
P: U O !i-i:n-i(X) = P:
_c Ip _c
u!i-i:n-i (-1
x ),
_ P
r~
and the identity
1 1 P
- ~j-i:n-i (1 ~ p) -g (1 ~ p) 1 p= 1 1
[!i-i:n-i(x)-g(X Wdx.
The second inequality in (7.21) follows from the Schwarz inequality com-
bined with changes of variables, and (7.2) implies the last one. Equality in
(7.21) is obtained if
p,+ f..
tcU 3-I:n-t I-p
.( z )
F-1(x+p)-y = + M~o, o ~ X < 1, (7.22)
II P>-c u !i-i:n-ill III
F-1(x) = y, o ~ x <po (7.23)
Recalling Theorem 11, we can now state the following assertions.
Theorem 40 (order statistics, decreasing density) If i +2 ~ j <
min{(2n + i + 2)/3, n - I}, then
which becomes the equality for F being the convex combination of the Dirac
measure at y and the uniform distribution on [y, y + yaM] with respective
coefficients p and 1 - p.
152 7. Predictions of Order and Record Statistics
where the right-hand side is defined by (4.46) with 0. = 0.(13.) and 13. =
13.(j - i, n - i) uniquely determined by (4.47) through (4.49). Inequality
(7.25) becomes the equality for
< 0,
1
??::::1l.
M
o < ??::::1l. < !i-im-.({3.)
71 - M IIUj-':n-.)".[Link]'
-
~J 1;-':n-.({3.) < x-y < 1;-':n-.({3.)+OI.(I-{3.) (7.27)
IIUj-':n-.)".p.11 - M - IIUj-':n-.)".p. II '
x-y > !j-i:n-i ({3. )+01. (1-{3.)
M - IIUj-i:n-;}".II.11 '
respectively.
In the case F ~c V, the projection method is essentially exploited for pre-
dicting sample maxima as well. By arguments similar to ones used in (7.21),
for i + 2 ~ j ~ n we deduce that
EF(Xj :n - Xi:nIXi:n = y)
= 1 00
[F-l(p + (1 - p)V(x)) - yjli-i:n-i V(x)v(x) dx
F- 1 (p+(1-p)V(x)) -y =
P: +v Ii-i:n-i V(x) Mpo,
c
x~O,
IIPtc vli-i:n-NII Iy
F- 1 (x) = y, O~x<p
(7.29)
7.2 Distributions with Decreasing Density and Failure Rate 153
(7.30)
(7.31)
(7.32)
respectively.
Note that (7.16) has decreasing density and failure rate. This means that
general bounds (7.15) are the optimal ones for the predictions of the values
of the first records from populations with decreasing density and failure
rate. Therefore it remains to study differences of nonsuccessive values of
kth records for k ~ 2. In order to obtain desired conclusions, we are reduced
to replacing !;-i:n-i(X) by f~".lm-l(x) in (7.21) and (7.28), and references
to Theorems 35 and 36.
Theorem 42 (kth records, decreasing density) If (l+l/kt- m ::; 3,
we have
(7.33)
This is the equality for the combination of the exponential distribution with
location parameter y and scale M/V2, and an atom at y.
Otherwise
(k) (k)1 (k) )
(
EF Rn - Rm Rm = y <
-
II(f(k) V) -1 II (7.36)
M FO n-m-l ct. V (-y.) ,
III
where the bound is defined by (6.27), and the parameters 0. = 0.(1'.) and
1'. = 'Y.(k,n-m-l) are determined from (6.28) to (6.30). Inequality (7.36)
becomes the equality for F defined in (7.31) and (7.32) with /i-i:n-i replaced
by f~k}.m-l'
Distribution functions (7.26), (7.31), and their modifications described in
Theorems 42 and 43 have unbounded derivatives in the right neighborhoods
of the left endpoints of their supports. The shape conditions entail unique
extensions of them to the left: only a jump at y and no mass of the left are
admitted. This means that the analogous bounds expressed in mFO -units
III
are attained by the same distributions only. This is not so in the remaining
cases. It is worth noticing that bounds (7.24) and (7.33) as well as (7.29)
and (7.35) are attained by the same elements of families with decreasing
density and failure rate, respectively.
2. Solve the analogous problem for the distributions defined by the star
order relations.
4. One can try to retrieve missing values of order statistics using (2.43).
However, the problem is that the quantile functions of doubly trun-
cated distributions Ell:do not form a convex cone, and one cannot
simply use our projection method here.
of sequences of L-statistics
~
n
Cj,nXj:n = ~
n
h ~ - h (i- -1)]
[(i) ;- Xj:n
158 8. Further Research Directions
1}j = k + (n - j)(m + 1) ~ 1, 1 $ j $ n,
have expectations
EFX(j n m k)=
, "10
r F- 1(x) (j-1)!(m+1) (1_x)l1i- 1[1_(1_x)m+1]j- 1dx
1 j.
TIi=l1},
for m i: 1, and
(8.1)
for positive x and y. The reversed inequality defines F ~+ W. The su-
peradditive order implies the star one and so describes larger classes of
8. Further Research Directions 159
distributions. In fact, relation (8.1) is applied for the life distributions and
should be satisfied for all 0 = aw ~ x + y < dw, but it can be modi-
fied by subtracting F-1W(aw) = F-I(O) from all the terms in (8.1) when
aw "# O. Especially for W = V, composition
defines the hazard function of F, and (8.1) describes the new worse than
used (NWU) distributions. The interpretation of the subadditivity of the
hazard function is that surviving time x + y by a single device is more
likely than under replacement by a new one in meantime x. The reversed
relation defines the new better than used (NBU) distributions. Since mean-
variance bounds presented here for F :5. W cannot be improved in general
populations, the same holds for the larger classes F:5+ W, and the NBU
distributions in particular. However, a difficult problem is to determine
bounds for F !:+ W, because the superadditivity does not allow a natural
graphical interpretation. Verifying the property of a function at fixed x + y,
we should study its values at all pairs x and y. Since
n
9j:n(X) =
n + 1 -). l[(j-I)/(n,l) (x)
is superadditive for (j - lIn ~ 1/2, we have the conclusion of Corollary 1
for the subadditive distribution functions F !:+ U then. We conjecture
that the bound on quantiles and order statistics of dependent samples of
the NWU populations are attained by Poisson distributions.
The method based on the greatest convex minor ants provides the best
L2-approximations of general functions by monotone ones. This was fruit-
fully exploited in determining sharp mean-variance and second moment
bounds on various statistical functionals in general and symmetric popu-
lations. For other families of distributions, the projection heavily depends
on properties of functions being projected. The reason is that there are no
known procedures which allow us to determine projections of general func-
tions onto corresponding convex cones of quantile functions and their mod-
ifications. What we strongly need here is, for example a general method of
projecting onto the family of convex functions. It is not clear here whether
one should try to generalize the Moriguti [58] method of convex minor ants
for this problem, or develop quite a different approach.
We obtain usually nons harp bounds on the statistical functionals defined
in Section 2.2 in terms of
and IIF-I - JtFllp, p > 1, by using the Holder inequality instead of the
Schwarz one. These are scale parameters expressed in terms of pth roots
160 8. Further Research Directions
Th(F- 1) = 11 F-l(X)h(x) dx
< 11 F-l(x)h(x) dx
(8.3)
(8.4)
Finally, it is of interest if projections onto convex sets that are not neces-
sarily convex cones provide meaningful evaluations of functionals. IT so, we
could significantly extend the class of functionals for which respective eval-
uations hold by adding ones that act directly on the distribution functions,
densities, and other characteristics of distributions.
References
[3] Arnold, B.C. (1985), p-Norm bounds on the expectation of the max-
imum of possibly dependent sample, J. Multivar. Anal. 17, 316-332.
[6] Arnold, B.C. and R.A. Groeneveld (1974), Bounds for deviations be-
tween sample population statistics, Biometrika 61, 387-389.
[24] David, H.A., H.O. Hartley, and E.S. Pearson (1954), The distribution
of the ratio, in a single normal sample, of range to standard deviation,
Biometrika 41, 482-493.
[25] Dharmadhikari, S. and K. Joag-dev (1988), Unimodality, Convexity,
and Applications, Academic, New York.
[26] Dziubdziela, W. and B. Kopocinski (1976), Limiting properties of the
kth record values, Zastos. Mat. 15, 187-190.
[27] Fahmy, S. and F. Proschan (1981), Bounds on differences of order
statistics, Amer. Statist. 35,46-47.
[28] Feldman, D. and H.G. Thcker (1966), Estimation of non-unique quan-
tiles, Ann. Math. Statist. 37, 451-457.
[29] Franco, M. and J.M. Ruiz (1996), On characterization of continuous
distributions by conditional expectation of record values, Sankhya
A 58, 135-14l.
[30] Franco, M. and J.M. Ruiz (1999), Characterization based on condi-
tional expectations of adjacent order statistics: A unified approach,
Proc. Amer. Math. Soc. 127, 861-874.
[31] Gajek, L. and A. Okolewski (2001), Projection method for moment
bounds on record statistics from restricted families, submitted for
publication.
[32] Gajek, L. and T. Rychlik (1996), Projection method for moment
bounds on order statistics from restricted families. I. Dependent case,
J. Multivar. Anal. 57, 156-174.
[33] Gajek, L. and T. Rychlik (1998), Projection method for moment
bounds on order statistics from restricted families. II. Independent
case, J. Multivar. Anal. 64, 156-182.
[34] Gascuel, O. and G. Caraux (1992), Bounds on expectations of order
statistics via extremal dependences, Statist. Probab. Lett. 15, 143-148.
[35] Grudzien, Z. and D. Szynal (1985), On the expected values of kth
records and associated characterizations of distributions, in: Proba-
bility and Statistical Decision Theory, Vol. A (F. Konecny, J. Mo-
gyor6dy, and W. Wertz,eds.) Reidel, Dordrecht, 119-127.
[36] Gumbel, E.J. (1954), The maxima of the mean largest value and of
the range, Ann. Math. Statist. 25, 76-84.
[37] Hampel, F.R., E.M. Ronchetti, P.J. Rousseeuw, and W.A. Stahel
(1986), Robust Statistics. The Approach Based on Influence Func-
tions, Wiley, New York.
166 References
[38] Hartley, H.O. and H.A. David (1954), Universal bounds for mean
range and extreme observation, Ann. Math. Statist. 25,85-99.
[39] Hawkins, D.M. (1971), On the bounds of the range of order statistics,
J. Amer. Statist. Assoc. 66, 644-645.
(40) Huang, M.L. and P. Brill (1999), A level crossing quantile estimation
method, Statist. Probab. Lett. 45, 111-119.
[41] Huber, P.J. (1981), Robust Statistics, Wiley, New York.
[44) Karlin, S. (1968), Total Positivity, Vol. I, Stanford Univ. Press, Stan-
ford, CA.
[45) Karlin, S. and Studden, W.J. (1966), Tchebyshev Systems: With Ap-
plications in Analysis and Statistics, Interscience, New York.
[46] Klefsjo, B. (1983), A useful ageing property based on the Laplace
transform, J. Appl. Prob. 20, 615-626.
[47] Lai, T.L. and H. Robbins (1976), Maximally dependent random vari-
ables, Proc. Nat. Acad. Sci. U.S.A. 73,286-288.
[48] Lai, T.L. and H. Robbins (1978), A class of dependent random vari-
ables and their maxima, Z. Wahrsch. Verw. Gebiete 42,89-111.
[49] Lawrence, M.J. (1975), Inequalities for s-ordered distributions, Ann.
Statist. 3, 413-428.
[50] Lopez-Blazquez, F. (1998), Discrete distributions with maximum
value of the maximum, J. Statist. Plann. Inference 70, 201-207.
[51) Lopez-Blazquez, F. (2000), Bounds for the expected value of spacings
from discrete distributions, J. Statist. Plann. Inference 84, 1-9.
[52) Ma, C. (1992), Variance bound of function of order statistic, Statist.
Probab. Lett. 13, 25-27.
[53] Mallows, C.L. (1969), Extrema of expectations of uniform order
statistics, SIAM Rev. 11,410-411.
[54] Mallows, C.L. and D. Richter (1969), Inequalities of Chebyshev type
involving conditional expectations, Ann. Math. Statist. 40, 1922-
1932.
References 167
[60] Nagaraja, H.N. (1981), Some finite sample results for the selection
differential, Ann. Inst. Statist. Math. 33, 437-448.
[63] Nair, K.R. (1948), The distribution of the extreme deviate from the
sample mean and its studentized form, Biometrika 35, 118-144.
[72] Plackett, R.L. (1947), Limits of the ratio of mean range to standard
deviation, Biometrika 34, 120-122.
[74] Raqab, M.Z. (1997), Bounds based on greatest convex minorants for
moments of record values, Statist. Probab. Lett. 36, 35-41.
[77] Rychlik, T. (1992), Weak limit theorems for stochastically largest or-
der statistics, in: Order Statistics and Nonparametrics. Theory and
Applications (I.A. Salama and P.K. Sen, eds.), North-Holland, Ams-
terdam, 141-154.
exponential distribution, 20, 40, kth lower record value, 22, 143
41, 44, 48, 67, 89, 106, kth (upper) record occurrence ti-
109, 111, 114, 139, 142, me, 21
153 kth (upper) record value, 21, 22,
131-143, 148, 149, 153,
failure rate, 28 154
Vol. 87: J. Milller, Lectures on Random Voronoi V~l.106: ~d Niederreiter and Peter Jau-Shyong
Tessellations. vii, 134 pages, 1994. ShIUe (Editors), Monte Carlo and Quasi-Monte Carlo
Methods in Scientific Computing. xiv, 372 pages, 1995.
Vol. 88: J. E. Kolassa, Series Approximation Methods
Vol. 107: Masafumi Akahira and Kei Takeuchi Non-
in Statistics. Second Edition. ix, 183 pages, 1997.
i
Regular Statistical Estimation. vii, 183 pages, 995.
Vol. 8?: P. Cheeseman, and R W. Oldford (Editors),
Vo~. 108: 'Yesley L. Schaible (Editor), Indirect
Selectmg Models from Data: AI and Statistics lV. xii,
487 pages, 1994. EstImators m US Federal Programs. viii 195 pages
1995. "
Vol.. ~O: A. Csenki, Dependability for Systems with a
PartitIOned State Space: Markov and Semi-Markov Vol. 109: Helmut Rieder (Editor), Robust Statistics
Theory and Computationallmplementation. x, 241 Data Analysis, and Computer Intensive Methods. xiv,
pages, 1994. 427 pages, 1996.
Vol. 91: J.D. Malley, Statistical Applications of Vol. 110: D. Bosq, Nonparametric Statistics for
Jordan Algebras. viii, 10 I pages, 1994. Stochastic Processes. xii, 169 pages, 1996.
Vol. 92: M. Eerola, Probabilistic Causality in Vol. III: Leon Willenborg and Ton de Waal
Longitudinal Studies. vii, 133 pages, 1994. Statistical Disclosure Control in Practice. xiv', 152
pages, 1996.
Vol. 93: Bernard Van Cutsem (Editor), Classification
and Dissimilarity Analysis. xiv, 238 pages, 1994. Vol. 112: Doug Fischer and Hans-J. Lenz (Editors),
Learning from Data. xii, 450 pages, 1996.
Vol. 94: Jane F. Gentleman and G.A. Whitmore
(Editors), Case Studies in Data Analysis. viii, 262 Vol. 113: Rainer Schwabe, Optimum Designs for
pages, 1994. Multi-Factor Models. viii, 124 pages, 1996.
Vol. 95: Shelemyahu Zacks, Stochastic Visibility in Vol. 114: C.C. Heyde, Yu. V. Prohorov, R. Pyke, and
Random Fields. x, 175 pages, 1994. S.T. Ra.c~ev (Edit~rs), A~ens Conference on Applied
Probablhty and TIme Series Analysis Volume I:
Vol. 96: Ibrahim Rahimov, Random Sums and Applied Probability in Honor of J.M. Gani. viii, 424
Branching Stochastic Processes. viii, 195 pages, 1995. pages, 1996.
Vol. 97: R. Szekli, Stochastic Ordering and Vol. 115: P.M. Robinson and M. Rosenblatt (Editors),
Dependence in Applied Probability. viii, 194 pages Athens Conference on Applied Probability and Time
lW~ , Series Analysis Volume n: Time Series Analysis in
Memory ofE.J. Hannan. viii, 448 pages, 1996.
Vol. 98: Philippe Barbe and Patrice Bertail The
Weighted Bootstrap. viii, 230 pages, 1995: Vol. 116: Genshiro Kitagawa and Will Gersch,
Smoothness Priors Analysis of Time Series. x, 261
Vol. 99: C.C. Heyde (Editor), Branching Processes: pages, 1996.
Proceedings of the First World Congress. viii, 185
pages, 1995. Vol. 117: P~ul Glasserman, Karl Sigman, and David
D. Yao (EdItOrs), Stochastic Networks. xii, 298,1996.
Vol. 100: Wlodzimierz Bryc, The Normal
Distribution: Characterizations with Applications. viii, Vol. 118: Radford M. Neal, Bayesian Learning for
139 pages, 1995. Neural Networks. xv, 183, 1996.
Vol. 101: H.H. Andersen, [Link], D. Sl!Irensen, Vol. 119: Masanao Aoki and Arthur M. Havenner
and P.S. Eriksen, Linear and Graphical Models for the Applications of Computer Aided Time Series '
Multivariate Complex Normal Distribution. x, 184 Modeling. ix, 329 pages, 1997.
pages, 1995.
Vol. 120: Maia Berkane, Latent Variable Modeling
Vol. 102: A.M. Mathai, Serge B. Provost, and Takesi and Applications to Causality. vi, 288 pages, 1997.
Hayakawa, Bilinear Forms and Zonal Polynomials. x,
378 pages, 1995. Vol. 121: Constantine Gatsonis, James S. Hodges,
Robert E. Kass, Robert McCulloch, Peter Rossi and
Vol; 103: Anestis Antoniadis and Georges Oppenheim Nozer D. Singpurwalla (Editors), Case Studies in
(EdItors), Wavelets and Statistics. vi, 411 pages, 1995. Bayesian Statistics, Volume III. xvi, 487 pages, 1997.
Vol. 104: Gilg U.H. Seeber, Brian J. Francis, Reinhold Vol. 122: .Timothy G. Gregoire, David R Brillinger,
Ha~[Link], and G~riele Steckel-Berger (Editors),
Peter J. Dlggle, Estelle Russek-Cohen, William G.
Stattsttcal Modellmg: 10th International Workshop, Warren, and Russell D. Wolfinger (Editors), Modeling
Innsbruck, July 10-14th, 1995. x, 327 pages, 1995. Longitudinal and Spatially Correlated Data. x, 402
pages, 1997.
Vol. 123: D.Y. Lin and T.R. Fleming (Editors), Vol. 142: Gyilrgy Terdik, Bilinear Stochastic Models
Proceedings of the First Seattle Symposium in and Related Problems of Nonlinear Time Series
Biostatistics: Survival Analysis. xiii, 308 pages, 1997. Analysis: A Frequency Domain Approach. xi, 258
pages, 1999.
Vol. 124: Christine H. Muller, Robust Planning and
Analysis of Experiments. x, 234 pages, 1997. Vol. 143: Russell Barton, Graphical Methods for the
Design of Experiments. x, 208 pages, 1999.
Vol. 125: Valerii V. Fedorovand Peter Hackl, Model-
Oriented Design of Experiments. viii, 117 pages, 1997. Vol. 144: L. Mark Berliner, Douglas Nychka, and
Timothy Hoar (Editors), Case Studies in Statistics and
Vol. 126: Geert Verbeke and Geert Molenberghs, the Atmospheric Sciences. x, 208 pages, 2000.
Linear Mixed Models in Practice: A SAS-Oriented
Approach. xiii, 306 pages, 1997. Vol. 145: James H. Matis and Thomas R. Kiffe,
Stochastic Population Models. viii, 220 pages, 2000.
Vol. 127: Harald Niederreiter, Peter Hellekalek,
Gerhard Larcher, and Peter Zinterhof (Editors), Monte Vol. 146: Wim Schoutens, Stochastic Processes and
Carlo and Quasi-Monte Carlo Methods. xii, 448 Orthogonal Polynomials. xiv, 163 pages, 2000.
pages, 1997.
Vol. 147: Jiirgen Franke, Wolfgang HlIrdle, and
Vol. 128: L. Accardi and C.C. Heyde (Editors), Gerhard Stahl, Measuring Risk in Complex Stochastic
Probability Towards 2000, x, 356 pages, 1998. Systems. xvi, 272 pages, 2000.
Vol. 129: Wolfgang Hlirdle, Gerard Kerkyacharian, Vol. 148: S.E. Ahmed and Nancy Reid, Empirical
Dominique Picard, and Alexander Tsybakov, Bayes and Likelihood Inference. x, 200 pages, 2000.
Wavelets; Approximation, and Statistical
Applications. xvi, 265 pages, 1998. Vol. 149: D. Bosq, Linear Processes in Function
Spaces: Theory and Applications. xv, 296 pages,
Vol. 130: Bo-Cheng Wei, Exponential Family 2000.
Nonlinear Models. ix, 240 pages, 1998.
Vol. 150: Tadeusz Calinski and Sanpei Kageyarna,
Vol. 131: Joel L. Horowitz, Semiparametric Methods Block Designs: A Randomization Approach, Volume
in Econometrics. ix, 204 pages, 1998. I: Analysis. ix, 313 pages, 2000.
Vol. 132: Douglas Nychka, Walter W. Piegorsch, and Vol. 151: HAkan Andersson and Tom Britton,
Lawrence H. Cox (Editors), Case Studies in Stochastic Epidemic Models and Their Statistical
Environmental Statistics. viii, 200 pages, 1998. Analysis. ix, 152 pages, 2000.
Vol. 133: Dipak Dey, Peter Miiller, and Debajyoti Sinha Vol. 152: David Rlos Insua and Fabrizio Ruggeri,
(Editors), Practical Nonparametric and Semiparametric Robust Bayesian Analysis. xiii, 435 pages, 2000.
Bayesian Statistics. xv, 408 pages, 1998.
Vol. 153: Parimal Mukhopadhyay, Topics in Survey
Vol. 134: Yu. A. Kutoyants, Statistical Inference for Sampling. x, 303 pages, 2000.
Spatial Poisson Processes. vii, 284 pages, 1998.
Vol. 154: Regina Kaiser and Agustin Maravall,
Vol. 135: Christian P. Robert, Discretization and Measuring Business Cycles in Economic Time Series.
MCMC Convergence Assessment. x, 192 pages, 1998. vi, 190 pages, 2000.
Vol. 136: Gregory C. Reinsel and Raja P. Velu, Vol. 155: Leon Willenborg and Ton de Waal,
Multivariate Reduced-Rank Regression. xiii, 272 Elements of Statistical Disclosure Control. xvii, 289
pages, 1998. pages, 2000.
Vol. 137: V. Seshadri, The Inverse Gaussian Vol. 156: Gordon Willmot and X. Sheldon Lin,
Distribution: Statistical Theory and Applications. xi, Lundberg Approximations for Compound
360 pages, 1998. Distributions with Insurance Applications. xi, 272
pages, 2000.
Vol. 138: Peter Hellekalek and Gerhard Larcher
(Editors), Random and Quasi-Random Point Sets. xi, Vol. 157: Anne Boomsma, Marijtje A.J. van Duijn,
352 pages, 1998. and Tom A.B. Snijders (Editors), Essays on Item
Response Theory. xv, 448 pages, 2000.
Vol. 139: Roger B. Nelsen, An Introduction to
Copulas. xi, 232 pages, 1999. Vol. 158: Dominique Ladiray and Benoit Quenneville,
Seasonal Adjustment with the X-II Method. xxii, 225
Vol. 140: Constantine Gatsonis, Robert E. Kass, pages, 200 1.
Bradley Carlin, Alicia Carriquiry, Andrew Gelman,
Isabella Verdinelli, and Mike West (Editors), Case Vol. 159: Marc Moore (Editor), Spatial Statistics:
Studies in Bayesian Statistics, Volume IV. xvi, 456 Methodological Aspects and Some Applications. xvi,
pages, 1999. 282 pages, 2001.
Vol. 141: Peter Muller and Brani Vidakovic (Editors),
Bayesian Inference in Wavelet Based Models. xi, 394 Vol. 160: Tomasz Rychlik, Projecting Statistical
pages, 1999. Functionals. ix, 169 pages, 2001.