7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
7.2
CiteScore
3.7
Impact Factor
Generic selectors
Exact matches only
Search in title
Search in content
Post Type Selectors
Search in posts
Search in pages
Filter by Categories
ABUNDANCE ESTIMATION IN AN ARID ENVIRONMENT
Case Study
Editorial
Invited review
Letter to the Editor
Original Article
REVIEW
Review Article
SHORT COMMUNICATION
View/Download PDF

Translate this page into:

Original article
26 (
3
); 200-204
doi:
10.1016/j.jksus.2013.11.003

Recurrence relations for the moments of order statistics from doubly truncated modified Makeham distribution and its characterization

King Saud University, College of Science, Department of Statistics and Operations Research, P.O. Box 2455, Riyadh 11451, Saudi Arabia
Cairo University, Faculty of Economics & Political Science, Department of Statistics, Giza, Egypt
Al-Azhar University, Faculty of Science, Department of Mathematics, Nasr City, 11884 Cairo, Egypt

*Corresponding author at: King Saud University, College of Science, Department of Statistics and Operations Research, P.O. Box 2455, Riyadh 11451, Saudi Arabia. Tel.: +966 114676335; fax: +966 114673036 ismail@ksu.edu.sa (Ali A. Ismail)

Disclaimer:
This article was originally published by Elsevier and was migrated to Scientific Scholar after the change of Publisher.

Available online 22 November 2013

Peer review under responsibility of King Saud University.

Abstract

In this study a general form of recurrence relations of continuous function for doubly truncated modified Makeham distribution is obtained. Recurrence relations between single and product moments of order statistics from doubly truncated modified Makeham distribution are given. Also, a characterization of modified Makeham distribution from the right and the left is discussed through the properties of order statistics.

Keywords

Modified Makeham distribution
Bathtub hazard rate function
Single and product moments
Order statistics
Recurrence relations
Characterization

Acronyms and Notations

MMD

modified Makeham distribution

h(x)

hazard rate function

f(x)

density function

F(x)

cumulative distribution function

b

MMD shape parameter (b > 0)

a

MMD scale parameter (a > 0)

1

1 Introduction

Order statistics arise obviously in many real life applications. The moments of order statistics have assumed considerable interest in recent years and have been tabulated quite widely for several distributions. For an extensive survey, see for example (Arnold et al., 1992; Hendi et al., 2006). Many authors have presented some recurrence relations satisfied by single and product moments of order statistics see for example (Khan et al., 1983a,b; Athar and Islam, 2004; among others). A number of recurrence relations satisfied by these moments of order statistics are available in the literature. Balakrishnan and Malik (1985) derived some identities involving the density functions of order statistics. These identities are useful in checking the computation of the moments of order statistics. Balakrishnan and Malik (1986) established some recurrence relations of order statistics from the linear-exponential distribution. Balakrishnan et al. (1988) reviewed several recurrence relations and identities for the single and product moments of order statistics from some specific distributions. Recently, Bekçi (2009) considered some recurrence relations for the moments of order statistics from uniform distribution.

Makeham distribution is an important life distribution and has been commonly used to fit actuarial data (see, Marshall and Olkin, 2007). According to Scollnik (1995), Makeham distribution of mortality represents a failure law, where the hazard rate is a mixture of non-aging failure distribution and the aging failure distribution with exponential increase in failure rates. For a description on the genesis and applications of Makeham distribution one may refer to Makeham (1860).

Feng et al. (2008) used the least squares type estimation to estimate the parameters of Makeham distribution. Aboutahoun and Al-Otaibi (2009) discussed relations between single and product moments of Makeham distribution.

The distribution presented in this paper is called modified Makeham distribution (MMD). This distribution has frequently been used to describe human mortality and to establish actuarial tables. Its hazard rate function is a product of a monotonically decreasing and a monotonically increasing function. The hazard rate function is given by

(1.1)
h ( x ) = b a x a b 1 exp x a b , where a is the scale parameter and b is the shape parameter. It is noted that for b < 1 the hazard rate function has a bathtub curve. This is a very desirable characteristic of the reliability models. It is important to note that the scale parameter a can be taken equal to 1 with no loss in generality. This simplifies many of the expressions.

Hazard rate function is the conditional probability that failure of the device will occur during a small interval [t;t + Δt] given that device has survived to time t. The bathtub hazard rate function is typical for the lifetime description of most technical devices. At the beginning the failure rate decreases, after a while it becomes stable, and then it increases due to natural wear out processes.

The cumulative distribution function of the MMD is of the following form:

(1.2)
F ( x ) = 1 exp 1 exp x a b . Derivative of the cumulative distribution function (1.2) gives density function of the following form:
(1.3)
f X ( x ) = b a x a b 1 exp x a b exp 1 exp x a b .
For basic characteristics of this distribution see (Kosznik-Biernacka, 2006, 2007, 2011).

In this paper, some recurrence relations satisfied by single and product moments of order statistics from doubly truncated MMD are established. Recurrence relations for the non-truncated case are given as a special case. A characterization of this distribution based on the properties of the order statistics has also been considered.

The doubly truncated probability density function of the MMD is given by

(1.4)
f ( x ) = 1 P Q b a x a b 1 exp x a b exp 1 exp x a b , Q 1 x P 1 . where 1 P = exp 1 exp P 1 a b and 1 Q = exp 1 exp Q 1 a b

Let P 2 = 1 P P Q and Q 2 = 1 Q P Q .

From (1.4), we have

(1.5)
1 F ( x ) = P 2 + a b a x b 1 e x a b f ( x ) . Let X1:n ⩽ X2:n ⩽ ⋯ ⩽ Xn:n be the order statistics from continuous distribution functions (c.d.f) F1(x) and probability density function (p.d.f) f1(x). The p.d.f of order statistics Xr:n and the joint p.d.f of two order statistics Xr:n and Xs:n are f r : n ( x ) = n ! ( r 1 ) ! ( n r ) ! [ F 1 ( x ) ] r 1 [ 1 F 1 ( x ) ] n r f 1 ( x ) , f r , s : n ( x , y ) = n ! ( r 1 ) ! ( s r 1 ) ! ( n s ) ! [ F 1 ( x ) ] r 1 [ F 1 ( y ) F 1 ( x ) ] s r 1 [ 1 F 1 ( y ) ] n s f 1 ( x ) f 1 ( y ) , x < y ; 1 r < s n . Then for any monotone continuous function ϕ(x), we have
(1.6)
α r : n = E [ ϕ ( x ) ] = n ! ( r 1 ) ! ( n r ) ! Q 1 P 1 ϕ ( x ) [ F ( x ) ] r 1 [ 1 F ( x ) ] n r f ( x ) dx
and for any measurable joint function ϕ(x,y), we have
(1.7)
α r , s : n = E [ ϕ ( x , y ) ] = n ! ( r 1 ) ! ( s r 1 ) ! ( n s ) ! Q 1 P 1 x P 1 ϕ ( x , y ) [ F ( x ) ] r 1 [ F ( y ) F ( x ) ] s r 1 [ 1 F ( y ) ] n s f ( x ) f ( y ) dydx .
Also, assume that

μ r : n ( k ) = E ( X r : n k ) , μ r , s : n ( j , k ) = E ( X r : n j X s : n k ) , μ 0 : n ( k ) = Q 1 k , k = 1 , 2 , ; n = 0 , 1 , . , μ n : n 1 ( k ) = P 1 k , E ( e ( X 0 : n / a ) b ) = e ( Q 1 / a ) b and E ( e ( X n : n 1 / a ) b ) = e ( P 1 / a ) b .

Now, to derive recurrence relations for the moments of order statistics from modified Makeham distribution, let us consider the following theorems as shown in the next sections.

2

2 Recurrence relations for single moments

Theorem 1

Let Xr:n ⩽ Xr+1:n, (1 ⩽ r ⩽ n) be an order statistics, Q1 ⩽ Xr:n ⩽ P1, n ⩾ 1 and for any measurable function ϕ(x), then the recurrence relation for the single moments of order statistics is

(2.1)
α r : n = P 2 α r : n 1 + Q 2 α r 1 : n 1 + a b nb E [ X 1 b e ( X / a ) b ϕ ( x ) ] .

Proof

From (1.6), we have

(2.2)
α r : n α r 1 : n 1 = n 1 r 1 Q 1 P 1 ϕ ( x ) [ F ( x ) ] r 1 [ 1 F ( x ) ] n r + 1 dx . In view of (1.5) and (2.2), we have
(2.3)
α r : n α r : n 1 = n 1 r 1 Q 1 P 1 ϕ ( x ) [ F ( x ) ] r 1 [ 1 F ( x ) ] n r P 2 + a b a x b 1 e ( x / a ) b f ( x ) dx = ( n 1 ) ( n r ) P 2 ( α r : n 1 α r 1 : n 2 ) + a b nb E [ ϕ ( x ) X 1 b e ( X / a ) b ] ,
(2.4)
But ( n r ) α r : n + r α r + 1 : n = n α r : n 1 .
(2.5)
Or ( n r ) α r 1 : n 1 + ( r 1 ) α r : n 1 = ( n 1 ) α r 1 : n 2 .
Now, from (2.3) and (2.5) we obtain α r : n = P 2 α r : n 1 + ( 1 + P 2 ) α r 1 : n 1 + a b nb E [ ϕ ( x ) X 1 b e ( X / a ) b ] . Since 1 + P2 = Q2, then the theorem is proved.  □

Notation 1

For the non-truncated case, by putting Q1 = 0, P1 = ∞, Q2 = 0, P2 = 1, we get α r : n = α r : n 1 + a b nb E [ X 1 b e ( X / a ) ϕ ( x ) ] .

Special cases:

  • Putting ϕ(x) = xb, then we obtain the simple recurrence relations between single moments of order statistics

    (2.6)
    μ r : n ( b ) = P 2 μ r : n 1 ( b ) + Q 2 μ r 1 : n 1 ( b ) + a b n E ( e ( X / a ) b ) .

  • Putting ϕ ( x ) = e ( x / a ) b , we obtain

    (2.7)
    E ( e ( X r : n / a ) b ) = P 2 E ( e ( X r : n 1 / a ) b ) + Q 2 E ( e ( X r 1 : n 1 / a ) b ) + 1 n ,

  • If we put r = 1 and n = 1 in (2.6) and (2.7), then we get

    (2.8)
    μ 1 : 1 ( b ) = P 2 P 1 b Q 2 Q 1 b + a b E ( e ( X / a ) b ) , where E ( e ( X / a ) b ) = e Q 1 a b e P 1 a b ln u e 1 u du , which can be easily computed numerically.

    and

    (2.9)
    E ( e ( X r : n / a ) b ) = P 2 e ( P 1 / a ) b + Q 2 e ( Q 1 / a ) b ) + 1 .

  • Putting r = 1, n = n − r in (2.7), we obtain

    (2.10)
    E ( e ( X 1 : n r / a ) b ) = P 2 e ( P 1 / a ) b ) + Q 2 e ( Q 1 / a ) b ) + 1 n r .

  • By putting r = 1, n = 2 in (2.10), we get

    (2.11)
    E ( e ( X 1 : 2 / a ) b ) = P 2 e ( P 1 / a ) b ) + Q 2 e ( Q 1 / a ) b ) + 1 .

3

3 Recurrence relations for product moments

Theorem 2

Let Xr:n ⩽ Xr+1:n,r=1, 2, … , n − 1 be an order statistics from a random sample of size n, and for r  <  s > n, then the recurrence relation for the product moments of order statistics is

(3.1)
α r , s : n α r , s 1 : n = nP 2 n s + 1 ( α r , s : n 1 α r , s 1 : n 1 ) + a b ( n s + 1 ) b E [ ϕ ( x , y ) . Y 1 b . e ( Y / a ) b ] , where ϕ ( x , y ) = ϕ ( x , y ) y .

Proof

From (1.7), we have

(3.2)
α r , s : n α r , s 1 : n = n ! ( r 1 ) ! ( s r 1 ) ! ( n s + 1 ) ! Q 1 P 1 x P 1 ϕ ( x , y ) [ F ( x ) ] r 1 [ F ( y ) F ( x ) ] s r 1 [ 1 F ( y ) ] n s + 1 f ( x ) dydx . But
(3.3)
[ 1 F ( y ) ] = P 2 + a b a y b 1 e ( y / a ) b f ( y ) .
Now, by substituting from (3.3) in (3.2), then
(3.4)
α r , s : n α r , s 1 : n = n ! ( r 1 ) ! ( s r 1 ) ! ( n s + 1 ) ! × Q 1 P 1 x P 1 ϕ ( x , y ) [ F ( x ) ] r 1 [ F ( y ) F ( x ) ] s r 1 [ 1 F ( y ) ] n s P 2 + a b a y b 1 e ( y / a ) b f ( y ) f ( x ) dydx .
Hence, the theorem is proved.  □

Notation 2

For the non-truncated case, by putting Q1 = 0, P1 = ∞, Q2 = 0,P2 = 1, we obtain α r , s : n α r , s 1 : n = n n s + 1 ( α r , s : n 1 α r , s 1 : n 1 ) + a b ( n s + 1 ) b E [ ϕ ( x , y ) . Y 1 b . e ( Y / a ) b ] . Special cases:

  1. If we choose ϕ(x,y) = xjyb, then we get

    (3.5)
    μ r , s : n ( j , b ) μ r , s 1 : n ( j , b ) = nP 2 n s + 1 ( μ r , s : n 1 ( j , b ) μ r , s 1 : n 1 ( j , b ) ) + a b n s + 1 . E ( X j e ( Y / a ) b ) , which is the recurrence relation between product moments of order statistics for the doubly truncated MMD

  2. If we choose ϕ(x,y) = yb and s=r+1 in above theorem and noting that μ r , r + 1 : n ( 0 , b ) = μ r + 1 : n ( b ) and μ r , r : n ( j , b ) = μ r : n ( j + b ) ,see ( Khan et al., 1983b ) , we get

    (3.6)
    μ r + 1 : n ( b ) μ r : n ( b ) = nP 2 n r μ r + 1 : n 1 ( b ) μ r + 1 : n 1 ( b ) + a b n r . E ( e ( Y / a ) b ) , which is the recurrence relation between single moments of order statistics for the doubly truncated MMD.

Notation 3

for the non-truncated case, if Q1 = 0, P1 = ∞, Q2 = 0,P2 = 1, then we have

(3.7)
μ r , s : n ( j , b ) = μ r , s 1 : n ( j , b ) n n s + 1 μ r , s : n 1 ( j , b ) μ r , s 1 : n 1 ( j , b ) + a b n s + 1 . E ( X j e ( Y / a ) b ) .

4

4 Characterization of MMD

The unconditional p.d.f. of (s − r)th order statistics in a sample of size (n − r)th order statistics is given by

(4.1)
f ( x s : n | x r : n = x ) = ( n r ) ! ( s r 1 ) ! ( n s ) ! [ F ( y ) F ( x ) ] s r 1 [ 1 F ( y ) ] n s f ( y ) [ 1 F ( x ) ] n r , x y , where f(xs:nxr:n = x) is the conditional p.d.f. of Xs:n given Xr:n = x and the sample drawn from f ( y ) [ 1 F ( x ) ] , x y on the left truncated at x.

Therefore, for the left truncated at x, it can be seen that Q 1 = x , Q = F ( x ) , P 1 = , P = 1 , Q 2 = 1 and P 2 = 0 . Also, if the parent distribution truncated from the right at y, then f ( x r : n | x s : n = y ) = ( s 1 ) ! ( r 1 ) ! ( s r 1 ) ! [ F ( x ) ] r 1 [ F ( y ) F ( x ) ] s r 1 f ( x ) [ F ( y ) ] s 1 , x y . Thus in the case of right truncation at x, we have Q 1 = 0 , Q = 0 , P 1 = x , P = F ( x ) , Q 2 = 1 F ( x ) and P 2 = 1 F ( x ) F ( x ) . In this section we obtain a characterization of MMD via the properties of order statistics as shown from the following theorems.

Theorem 3

If F(x) is a distribution function, 0 < F(x) < 1, then X is modified Makeham if and only if E ( e ( X r + 1 : n / a ) b | X r : n = x ) = ( n r ) e ( x / a ) b + 1 n r .

Proof

The necessary condition is proved from (2.10) by using the properties of left truncation. To prove the sufficient condition, we have n r [ 1 F ( x ) ] n r x e ( y / a ) b [ 1 F ( y ) ] n r 1 f ( y ) dy = 1 n r + e ( x / a ) b . That is,

(4.2)
( n r ) x e ( y / a ) b [ 1 F ( y ) ] n r 1 f ( y ) dy = [ 1 F ( x ) ] n r 1 n r + e ( x / a ) b . Differentiating (4.2) with respect to X, we obtain ( n r ) e ( x / a ) b [ 1 F ( x ) ] n r 1 f ( x ) = ( n r ) [ 1 F ( x ) ] n r 1 1 n r + e ( x / a ) b f ( x ) + [ 1 F ( x ) ] n r e ( x / a ) b b a x a b 1 , which gives h ( x ) = b a x a b 1 exp x a b . Then, the theorem is proved.  □

Without loss of generality, when a = b = 1, we have n r [ 1 F ( x ) ] n r x e y [ 1 F ( y ) ] n r 1 f ( y ) dy = 1 n r + e x . That is, ( n r ) x e y [ 1 F ( y ) ] n r 1 f ( y ) dy = [ 1 F ( x ) ] n r 1 n r + e x . Differentiating the above equation with respect to X, we obtain ( n r ) e x [ 1 F ( x ) ] n r 1 f ( x ) = ( n r ) [ 1 F ( x ) ] n r 1 1 n r + e x f ( x ) + [ 1 F ( x ) ] n r e x , which gives f ( x ) 1 F ( x ) = e x . The solution of above differential equation is

f ( x ) = e x e 1 e x , which is the standard modified Makeham distribution.

Theorem 4

If F(x) is any distribution function, 0 < F(x) < 1, then X is modified Makeham if and only if E ( e ( X 1 : n / a ) b | X 2 : n = x ) = 1 F ( x ) F ( x ) e ( x / a ) b + 1 F ( x ) + 1 .

Proof

The sufficient condition is proved from (2.11) by using the properties of right truncation. To prove the necessity, we have 1 F ( x ) 0 x e ( y / a ) b f ( y ) dy = 1 F ( x ) F ( x ) e ( x / a ) b + 1 F ( x ) + 1 . Or

(4.3)
0 x e ( y / a ) b f ( y ) dy = [ 1 F ( x ) ] e ( x / a ) b + 1 + F ( x ) . Differentiating (4.3) with respect to x, we get

e ( x / a ) b f ( x ) = f ( x ) e ( x / a ) b [ 1 F ( x ) ] e ( x / a ) b b a x a b 1 + f ( x ) , which gives h ( x ) = b a x a b 1 e x a b . Thus, the theorem is proved.  □

5

5 Conclusion

In this study, general relations of order statistics for any continuous measurable function for doubly truncated MMD had been established. Further, a characterization of MMD based on the properties of the order statistics had been obtained. The recurrence relations for the single and product moments of order statistics are important in the theory of order statistics. The moments of order statistics were obtained by some other moments of order statistics. Through the properties of order statistics, the characterization results from left and right of MMD distribution were given.

References

  1. , , . Recurrence relations between moments of order statistics from doubly truncated Makeham distribution. Compt. Appl. Math.. 2009;28(3):277-290.
    [Google Scholar]
  2. , , , . A First Course in Order Statistics. New York: Wiley; .
  3. , , . Recurrence Relations for single and product moments of generalized order statistics from a General Class of Distribution. Metron. 2004;LXII(3):327-337.
    [Google Scholar]
  4. , , . Order statistics from the linear-exponential distribution, part I: increasing hazard rate case. Commun. Statist. Theor. Meth.. 1986;15(1):179-203.
    [Google Scholar]
  5. , , . Some general identities involving order statistics. Commun. Statist. Theor. Meth.. 1985;14(2):333-339.
    [Google Scholar]
  6. , , , . Recurrence relations and identities for moments of order statistics, II: specific continuous distributions. Commun. Statist. Theor. Meth.. 1988;17(8):2657-2694.
    [Google Scholar]
  7. , . Recurrence relations for the moments of order statistics from uniform distribution. Sci. Res. Essays. 2009;4(11):1302-1305.
    [Google Scholar]
  8. , , , . Estimation of parameters of the Makeham distribution using the least squares method. Math. Comput. Simul.. 2008;77:34-44.
    [Google Scholar]
  9. , , , . Order Statistics from Doubly Truncated Gompertz Distribution and its Characterizations. Egypt. Statist. J.. 2006;50(1):21-31.
    [Google Scholar]
  10. , , , . Recurrence relations between moments of order statistics. Naval Res. Logist. Quart.. 1983;30:419-441. Corrections, 32, 1985, 693
    [Google Scholar]
  11. , , , . Recurrence relations between product moments of order statistics. J. Statist. Plan. Inf.. 1983;8:175-183.
    [Google Scholar]
  12. , . Rodzina rozkladow Makehama. W Slupsku, PAP; .
  13. , . Makeham’s generalised distribution. Computat. Meth. Sci. Technol.. 2007;13(2):113-120.
    [Google Scholar]
  14. , . Sums of independent random variables of the modified Makeham distributions. Computat. Meth. Sci. Technol.. 2011;17(1-2):35-40.
    [Google Scholar]
  15. , . On the law of mortality and the construction of annuity tables. Assur. Magaz. J. Inst. Actuaries (Lond.). 1860;8:301-310.
    [Google Scholar]
  16. , , . Life Distributions: Structure of Nonparametric, Semiparametric, and Parametric Families. New York: Springer-Science; .
  17. , . Simulating random variates from Makeham’s distribution and from others with exact or nearly Log-Concave densities. Trans. Soc. Actuaries. 1995;47:409-437.
    [Google Scholar]
Show Sections