THE LENGTH-BIASED WEIGHTED WILSON HILFERTY DISTRIBUTION AND ITS APPLICATIONS Текст научной статьи по специальности «Математика»

Shivendra Pratap Singh, Surinder Kumar and Naresh Chandra Kabdwal , No 2

Vnliimf 19 Turip 2024

LENGTH-BIASED WEIGHTED WILSON HILFERTY DISTRIBUTION... Volu^ne 19' T ' 2024

THE LENGTH-BIASED WEIGHTED WILSON HILFERTY DISTRIBUTION AND ITS APPLICATIONS

Shivendra Pratap Singh

Department of Statistics, Babasaheb Bhimrao Ambedkar University, Lucknow, India

[email protected]

Surinder Kumar

Department of Statistics, Babasaheb Bhimrao Ambedkar University, Lucknow, India

[email protected]

Naresh Chandra Kabdwal*

Department of Mathematics and Statistics, Banasthali Vidyapith, Rajasthan, India

[email protected]

Abstract

In this article, we propose a new length-biased weighted form of Wilson Hilferty distribution named as Length-Biased Weighted Wilson Hilferty Distribution. The various Statistical properties of the proposed distribution like, reliability function, hazard rate function, reverse hazard rate function, moment generating function, quantile function, the coefficient of variation etc. are considered to understand its nature. Furthermore, we have used the method of maximum likelihood for estimation of the parameters of proposed distribution. Also, we obtain the Shannon's entropy, stochastic ordering, Lorenz and Bonferroni curves. The performance of the proposed distribution is compared with competitive distributions using two real data sets.

Keywords: Wilson Hilferty distribution, length-biased weighted Wilson Hilferty distribution, hazard function, reversed hazard function, maximum likelihood estimation.

1. Introduction

The weighted distribution arises when the observations are recorded from random process, the probability of recorded observations are not equal, and instead they are recorded according to some weighted function. The concept of weighted distributions was first given by [2]. Subsequently, [3] introduced a general form for model specification and data interpretation problems and identified that many situations can be modelled by weighted distributions.

Let T is as non-negative random variable with the probability density function (pdf) f(t), then the weighted distribution is given by

M (!)

Shivendra Pratap Singh, Surinder Kumar and Naresh Chandra Kabdwal , No 2

Vnliimf 19 Turip 2024

LENGTH-BIASED WEIGHTED WILSON HILFERTY DISTRIBUTION..._Volu^"P ' Tцnp, 2024

on the support of T, where w(t) > 0 and H = J w(t)f(t)dt is considered as a normalizing factor that forces gw(t) to integrate to unity. When we replace w(t) = t (i. e. the length of units) in equation (1), we get a special case of the weighted distribution called length biased distribution see, [5]. A r. v. T, is said to have a length biased weighted distribution if its pdf is defined as

gL{t) = w^m (2)

Where, n = J wL(t) Wj(t) f(t)dt and Wj(t) > 0, provided thet wL(t) = t.

Weighted distribution in general and length biased distributions are specifically very useful and convenient for the analysis of life time data. Weighted distributions are commonly used in study related to reliability, biomedicine, ecology, analysis of family data, branching process and various other fields of research, see [6], [8], [10] and [12].

Many length-biased weighted distributions with applications in different fields have been presented in the literature, see [9] discussed the characterization of inverse Gaussian and gamma distribution through their length biased distributions, [13] introduced a new class of weighted exponential distribution, [14], [15] discussed the length-biased weighted generalized Rayleigh distribution and also studied the length biased weighted Weibull distribution for rainfall data in India, a weighted Lindley distribution for survival data is given by [16], the length-biased lognormal distribution with application in the analysis of oil field exploration data is discussed by [17], [18] proposed the length-biased weighted exponential and Rayleigh distribution and its properties, different methods of estimation of parameters are applied for weighted exponential distribution by [19] which introduced by [13], [20] presented the length-biased weighted Lomax distribution and application with cancer data, [21] proposed inverted weighted exponential distribution and its properties, [24] discussed the length-biased exponential distribution for Bayesian reliability estimation, [23] proposed the length-biased weighted Lindley distribution, [25] proposed weighted exponentiated inverted exponential distribution and its properties, time and failure censoring schemes for Marshall Olkin alpha power extended Weibull distribution is presented by [26]. Recently, [28] proposed power weighted Sujatha distribution and application to survival times of patients head and neck cancer data, [27] proposed a weighted intervened exponential distribution as a lifetime model.

[1] introduced a Wison Hilferty distribution. For some recent developments of the Wilson Hilferty destitution the readers may, see [22]. Its probability density function (pdf), and cumulative distribution function (cdf), respectively as

r(a)

(3)

via, it3)

where, a and p are the shape and scale parameters, respectively, T(a) = J™ te-1e-tdt and y(x,y) = Jo wy-1e-wdw the complete and lower incomplete gamma functions, respectively.

In this paper, we propose a new length-biased distribution, called length-biased weighted Wilson-Hilferty distribution. Rest of the paper is structured as follows, in Section 2, the proposed distribution is introduced and its properties and reliability characteristics are discussed. In Section 3, the method of maximum likelihood is discussed for estimating the model parameters. Stochastic ordering and entropy are discussed in Section 4. Bonferroni and Lorenz curves and random number generation & Quantiles are discussed in Section 5 and 6, respectively. The applications of two real data sets are presented in Section 7. Finally, the conclusion is summarized in section 8.

Shivendra Pratap Singh, Surinder Kumar and Naresh Chandra Kabdwal , No 2 (78)

Vnliimf 19 Turp 2024

LENGTH-BIASED WEIGHTED WILSON HILFERTY DISTRIBUTION..._Volu^np ' J ' 2024

2. Length-Biased Weighted Wilson Hilferty Distribution, its Properties and Reliability Characteristics

In this section, we develop the Length-Biased Weighted Wilson Hilferty distribution. For this proposed new distribution, we present the pdf, cdf, reliability function, hazard function, moments, skewness and discuss some properties.

Consider the weight function as Wj(t) =~tk ,k> 0, wL(t) = t . Hence using considered weight function, unit length and equation (3) into the equation (2), we obtain the density of Length-biased weighted Wilson Hilferty distribution (LBWWHD) of the form

,k+i

m ; t,a,p,k>0 (5)

Where, a, p and k are the shape, scale and weighted parameters, respectively.

Figure 1: pdf plots of LBWWHD

Figure 1, clearly shows that LBWWH distribution is positively skewed. The cdf of LBWWHD is given by

/ k+1 a i3\ vla+-, -rrt3)

nt) = \f3J; t,a,p,k>0 (6)

where, f* ts-1 e-tdt = y(s, x) is a lower incomplete gamma function. The reliability function R (t) is given by

R(t) = 1

( ,k + 1 a

y(a+—- J t3)

On using some basic concept of an upper incomplete gamma integral's, it reduces to

r(a+OH, °t3) R(t) = ^ 3 ' ; t,a,p,k>0

r[a+

where, Jxrats 1e tdt = T(s,x) is an upper incomplete gamma function.

Table 1: Reliability function R(t) of LBWWHD for a = 2^ and 0 = 3

t k = 1 k = 2 k = 3 k = 4 k = 5

0.5 0.9999827 0.9999965 0.9999993 0.9999999 1.0000000

0.6 0.9999268 0.9999822 0.9999958 0.999999 0.9999998

0.7 0.9997537 0.9999303 0.9999809 0.9999949 0.9999987

0.8 0.9993024 0.9997745 0.9999294 0.9999785 0.9999936

0.9 0.9982712 0.9993722 0.999779 0.9999244 0.9999748

1.0 0.9961542 0.9984503 0.9993945 0.99977 0.9999149

(7)

3

Table 2: Reliability function R(t) of LBWWHD for k = 2 and P = 3

t a = 1 a = 2 a = 3 a = 4 a = 5

0.5 0.9999044 0.9999965 0.9999999 1.0000000 1.0000000

0.6 0.9997166 0.9999822 0.9999989 0.9999999 1.0000000

0.7 0.999292 0.9999303 0.9999935 0.9999994 0.9999999

0.8 0.9984419 0.9997745 0.9999691 0.9999958 0.9999994

0.9 0.9968914 0.9993722 0.9998803 0.9999773 0.9999957

1.0 0.9942659 0.9984503 0.9996054 0.9999000 0.9999746

Table 3: Reliability function R(t) of LBWWHD for k = 2 and a = 3

t 1 II P = 2 P = 3 4 II P = 5

0.5 0.9993884 0.999997 0.9999999 1.0000000 1.0000000

0.6 0.9955998 0.9999748 0.9999989 0.9999999 1.0000000

0.7 0.979182 0.9998513 0.9999935 0.9999993 0.9999999

0.8 0.9297582 0.9993323 0.9999691 0.9999967 0.9999994

0.9 0.8219026 0.9975858 0.9998803 0.999987 0.9999977

1.0 0.6472319 0.9927078 0.9996054 0.9999557 0.9999921

From Table 1,2 & 3, we conclude that, for the different values of a, p and k the reliability of the

distribution decreases with increase in the value of t.

Figure 2: Reliability plots of LBWWHD

Figure 2, shows the reliability behavior of the LBWWHD for varying values of shape parameter a, scale parameter p and weighted parameter k. Reliability function behaves like decreasing function.

The hazard function is defined as

№

h(t) =

R(t)

_3 (P

(J)0

, k+1

'eM-jt3} ^(c

¿¡m

r(«+k-+-) №. 7t3)

On simplifying, we get

, k+1

h(t)

_3(ff+ 3 t3*+*eXp{-ft3}

T[a+■

k+1 a_

3 . p

*3)

t, a, p. k> 0

(8)

- (p=1.5,a = 1.5,k=2)

-- (p-2.5,a- 1.5,k-2)

.... (p- 3.5.a- 1.5,k-2)

(P=45.a = 1B.k=2)

--(p = 5.5,a = 1.5,k=2) / t'

/ , '

- (p-3,a=5,k-4)

- - (p = 3.a=4.k = 4)

---- (p = 3.a = 3.k=4) / ,

(p = 3,a = 2rk=4] / *

--(p = 3.ct=1.k = 4) / /

/ /

' '

' —

______

t

3

Figure 3: Hazard rate plots of LBWWHD

Figure 3, shows the behaviour of hazard function for distinct values of a, p and k. Clearly, it shows that the hazard function of LBWWH behaves increasing hazard rate. The reverse hazard rate is defined as

№

3

Rh(t) =

(fi

Rh(t) =

k+1 3 f-3a+k

F(t)

t3a+K exp {-j t3} f(a +

r( ,k + 1\ ( , k + 1 a >

On simplifying, we get

.k+1

Rh(t) = ^ f U+, a ; t,a,p,k>0

Y(a+— Jt3)

Theorem 2.1. For r = 0, 1, 2, 3, ... rth moment of random variable T is given by

„/o -f ft+r+1^

..' - F(Tn- P\r'3 r[a+—

Vr = t(I ) = (J r(a+i±1

(9)

Proof. It T is a random variable with pdf f(t) from equation (5), then the rth moment is

( k+i\

o(a\(a+3)t3a+k^{ a 3}

fctn- w- r

E(T) = Vr=l ' ZT—TTu

I , k+1\

= -7^C<(3a+k+r)exp{- J t3}dt (10)

Theorem follows on taking y = ( j t3), and using the gamma function in equation (10).

Lemma 2.1. If a random variable T follows Length-biased weighted Wilson Hilferty distribution then on substituting r = 1,2 in equation (10), we obtain the mean and variance, respectively.

E(T) = (^) V 3 '

,B,2/3r(a + 1+t)

E(T2) = (H

r(a+—)

and,

Variance

r ( a +

1/3r(a +

k + 2> 3

r(a +

k + 1>

3

r(a + 1 + 3) K

a +

k + 2} 3

a +

k + U 3

Lemma 2.2. If a random variable T follows Length-biased Weighted Wilson Hilferty distribution then the coefficient of variation (C.V) is given by

,1/2

(11)

Proof. Coefficient of variation is given by,

C.V.=

^var(T) E(T)

or

k+2

ria+^yUa + ^l

1/2

k + 1\

k + 2>

i/3r(a + k~^)

k + 1

r(a + k~+-)

r(a + 1 + ^)r(a + k-^)-{r(

a+

k+2

1/2

k+2

k+2 r(a+—)

Lemma 2.2, follows on using Lemma 2.1.

Table 4: Coefficients of LBWWHD for P = 3

a k Mean Variance CV Skewness Kurtosis

2 1 1.521644 0.10416361 0.212102093 0.03819147 2.895791

2 1.590099 0.10032135 0.199192288 0.03175884 2.908649

3 1.65319 0.09697548 0.188368465 0.02692476 2.918847

4 1.71185 0.09402637 0.179126206 0.0231893 2.927098

5 1.766776 0.09139999 0.17111638 0.02023536 2.933893

4 1 1.482214 0.05474587 0.157857336 0.01590133 2.944386

2 1.519149 0.05351908 0.152283885 0.01427804 2.94851

3 1.554379 0.05239432 0.147260131 0.01291154 2.95208

4 1.588086 0.05135774 0.14270163 0.01174866 2.955199

5 1.620426 0.05039801 0.138540736 0.01074952 2.957945

2

2

2

3

Table 5: Coefficients of LBWWHD for a = 3

f k Mean Variance CV Skewness Kurtosis

1 1.306386 0.05475973 0.179126208 0.0231893 2.927098

2 1.348303 0.05323016 0.171116302 0.02023536 2.933893

2 3 1.387782 0.05185586 0.164088321 0.01785365 2.939574

4 1.425148 0.05061155 0.157857366 0.01590133 2.944386

5 1.460661 0.04947741 0.152283919 0.01427804 2.94851

1 1.645943 0.08692566 0.179126239 0.0231893 2.927098

2 1.698755 0.08449762 0.171116343 0.02023536 2.933893

4 3 1.748496 0.08231604 0.164088291 0.01785365 2.939574

4 1.795574 0.08034083 0.157857365 0.01590133 2.944386

5 1.840318 0.07854049 0.152283878 0.01427804 2.94851

According to the Table 4 & 5, skewness decreases and kurtosis increases whenever the values of k increases.

Lemma 2.3. If a random variable T follows Length-biased Weighted Wilson Hilferty distribution then harmonic mean (H) is given by

i = ("\1/3 r(a+!3)

h W r( Proof. The harmonic mean (H) is defined as

3) , k+i\

(12)

1=e(J) = C1 fWt

(13)

Using equation(5), we get

k+1

1 _ rm1 H~J0 t

(fi

l3(j) ' t3'+kexp{-f t3}

k+1

r^-i1)

dt

r(*+±±i)

1\l

t-3a+k — 1

exp{-— t3l}dt

Lemma 2.3, follows on using the transformation y = (^ ) t3, and the gamma function.

Lemma 2.4. If a random variable T follows Length-biased Weighted Wilson Hilferty distribution then moment generating function (MGF) and characteristic function (CF) of T are respectively, given by

M.

to = ZZo — (~) r(a+^±1)

3

k+i\

(14)

Proof: On using equation (5) and Taylor's series expansion the Lemma 2.4, follows.

3

3

3. Parameter Estimation

In this section, we estimate the parameters of the LBWWHD by using the maximum likelihood technique. Let T1, T2... Tn be the random sample of size n follows the LBWWHD(a,ft, k), then the likelihood function given as

I k+l\ 3" (f f^ '

The log-likelihood function can be written as

logL(t) = nlog3 + n(a + loga -n(a + k+;) logp

-nlogr( a+k+-)-jZti' + (3a + k)Z log ^ (16)

Differentiating equations (16) partially with respect to a, ft and k then equate to zero, we get

normal equations on the following form

dlogL(t) = 0 * ft = (17)

'f^=0* nloga+^a + ^-nlogp-^'-n^ (r(a + k^))

+3T1logti = 0 (18)

aiogm = 0* iioga- Ziogf} -^(^(a+^+Z log tt = 0 (19)

dk ' ° '

where, ip(z) = ^T(z) = ^^ is a logarithmic derivative of gamma function. As it seems, from

equations (17), (18) and (19), the analytical solution of a, p and k are not available.

Consequently, we have to use to non-linear estimation of the parameters using iterative method.

4. Stochastic Ordering and Entropy

Let X and Y be two independent random variables follows LBWWHD with shape parameter a, weighted parameter k and the scale parameters and p2, respectively.

When fx(t) and fY(t) be the density functions of X and Y, then X less than Y in likelihood order (X <lr Y) if is an increasing function of t. Here,

k+1

(wf3''"kM-k r^)

fx(t) ri ^k + 1\ , xa+!k±1 , ,

(a + —] 3(1;) '^expf-^}

k + 1 . a+-

fr(t) ' eXPrW2t}

(kf^'-expi-i;"}

k+1

Ml= (SlY+^expM^-^t'} (20)

fx(t) \bJ Ui bJ J v J

Differentiating equation (20), with respect to t, we get

k+1

dt\fx(ty (p2) (Pl (Pl p2) }

hence, tXÏM)>0 when, (fi — f2)>0 i.e. ft > ft.

Therefore, X <lr Y for, ft > Pt. If % < ir Y, then the following ordering shall also holds for the LBWWHD ([11]) ,

X <ir Y ^ X <hr Y ^ X <mri Y X<stY

for ft >ft.

The idea of entropy has significant importance in several academic disciplines, including probability and statistics, physics, communication theory, and economics. Entropy is a measure that quantifies the level of variety, uncertainty, or unpredictability shown by a given system. The entropy of a random variable T may be defined as a quantitative measure of the level of uncertainty or variation associated with it. The Shannon's entropy defined by

S(x) = -E[\ogf(t)]

Using equation (5), we get

k + 1 . a+-

S(x) = —E log

3 (L) 3 t3^exp{—f t3)

kTT{

[ ( ^^ h

k+i\

= —log( ^ ) — (3« + k)E(log(t)) +1E(t3) (21)

By solving the value of E(log(t)) and E(t3) and put in equation (21), we get

k + 1 ,a+—

. iM+L —Qq + Jlog^+vL+Î+lÙ + ^+Tl4)

■yr(a+L+l)J ( ( 3 ') r(a+±±i)

5. Bonferroni and Lorenz Curves

Let's assume that the random variable T is a non-negative with a continuous and twice differentiable cumulative distribution function. The Bonferroni curve of the random variable T is defined as

1

B(p)=—\ tf(t)dt

VP Jo

where, p = F(t), q = F 1(p) and ß = E(t)

B(p) =

rla+

k+i\

\1/3„( k+2■• ) T(a+—

I ^

Jo

3(j) 3 t3a+kexp{-jt3}

r(a+-

dt

(22)

(§)1/3<a++)

" i-3a+k+1

I "t

0

exp

H'3}d<

3

P

P

a

a

3

3

P

P

Substituting, y = (O) t3 in equation (22), we get Bonferroni curve

/ k+2 a 3\

B(P)=Ya+;m1

P r(a+—)

(23)

Lorenz curve is defined as

1 r"

L(P)=-I tf(t)dt = pB(p) ß Jo

Using equation (23), we get

/ k+2 a 3\ L(p)=Y(a+k+2k+ill

r(a+^

(24)

6. Random Number Generation and Quantiles

Random numbers of LBWWH can be easily generate by using the following function

=

(ß\ .t k + 1 m112

-, 1-U

where, U~U(0,1) and Q 1(a, z) is inverse of regularized incomplete gamma function, where regularized incomplete gamma function is defined as Q(a,z) = • Quantiles are given by

^[(iK'HiM-q)]

1/2

(25)

By putting q = 0.5 in equation (25), we get the median of LBWWHD

to.5 =

(ßM

a + -

k + 1

3

, 0.5)

3

l

2

7. Applications

In this section, we have considered two real data sets to check the suitability of the proposed distribution. Further, we have compared the distribution with the Length-Biased weighted Lindley distribution (LBWLD), Length-Biased Susheela distribution (LBSDi) and length-Biased Suja distribution (LBSD2), for suitability of proposed distribution. For this, we have used Akaike information criterion (AIC), Bayesian information criteria (BIC), Akaike Information Criterion Corrected (AICC) and Hannan-Quinn Information Criterion (HQIC), respectively. The AIC, BIC, AICC and HQIC are defined as:

AIC = 2K - 2logL, BIC = Klogn - 2 logL,

AICC = AIC + 2 K(K+1) , and HQIC = 2tflog(log(n)) - 2 logi.

n—K— 1

where n is the sample size, K is the number of parameters, and L denotes the likelihood function. Any probability model having smaller value of AIC, BIC and —LogL being the best model to fit the data set.

Dataset 1: The first dataset is taken from [7] which represent the tensile strength, measured in GPa, of 69 carbon fibers tested under tension at gauge lengths of 20mm.

Table 6: MLE, AIC and BIC for gauge lengths of 20 mm

Distribution a P k Log L AIC BIC AICC HQIC

LBWWHD 0.540 3.014 6.259 -48.895 103.790 110.492 104.148 106.449

LBWLD -- -- -- -87.8984 179.796 184.265 179.973 181.569

LBSD1 -- -- -- -92.3037 188.607 193.075 188.78 190.380

LBSD2 -- -- -- -134.128 186.607 188.841 186.66 187.493

According to the results in Table 6, the LBWWHD has the smallest values of these statistics, followed by LBWLD, LBSD1 and LBSD2. Therefore, the suggested distribution is the best choice for the tensile strength data. The fitted pdfs and empirical cdfs plots of the four models are sketched in Figure 4. Therefore, we assert that the LBWWHD fitting successfully the empirical plots of the data set.

Figure 4: Estimated densities and cdf plot of the models based on the real dataset 1.

Dataset 2- The second dataset used by [7], which represent the tensile strength, measured in GPa, of 63 carbon fibers tested under tension at gauge lengths of 10mm.

Table 7: MLE, AIC and BIC for gauge lengths of 10 mm

Distribution a P k Log L AIC BIC AICC HQIC

LBWWHD 0.1587 1.730 7.402 -58.7320 123.464 129.893 123.857 125.992

LBWLD -- -- -- -93.4265 190.85 195.133 191.046 192.538

LBSD1 -- -- -- -97.9972 199.994 204.280 200.188 201.680

LBSD2 -- -- -- -65.9256 133.851 135.995 133.914 134.694

For the dataset 2, we infer from the Table 7, the LBWWHD has the lowest values of AIC, BIC, AICC and HQIC followed by LBWLD, LBSD1 and LBSD2. Therefore, we conclude that LBWWHD is the most suitable choice for this dataset among the considered distributions. The fitted pdfs and empirical cdf plots of the four models are presented in Figure 5 and see that the LBWWHD fitting successfully the empirical plots of the data set.

Figure 5: Estimated densities and cdf plot of the models based on the real dataset 2.

8. Conclusion

In the present article a new distribution known as the Length-biased Weighted Wilson Hilferty distribution, has been proposed. The distribution in discussion is characterized by three parameters called shape, scale and weighted parameter. Through the use of certain formulae, the properties and characteristics of this distribution such as its moments, failure rate, reliability function etc., comprehensively examined and, the parameters estimation and stochastic comparison is also done. The examination and subsequent comparison of the criteria for AIC, BIC, AICC and HQIC have been conducted in relation to the Length-biased weighted Lindley distribution, Length-biased Sushila distribution and Length-biased Suja distribution. The actual lifetime of two sets of data has been successfully modelled and the resulting fit has been determined to be satisfactory.

Author Contributions: All authors have equal contribution. All authors reviewed the results and approved the final version of the manuscript.

Funding: This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors

Conflict of interest: The Authors declares that there is no conflict of interest.

References

[1] Wilson E. B., Hilferty M. M. (1931). The distribution of chi-square. Proc Natl Acad. Sci USA.;17:684-688.

[2] Fisher, R. A. (1934). The effects of methods of ascertainment upon the estimation of frequencies. Annals of Eugenics, 6, 13-25.

[3] Rao, C. R. (1965). On discrete distributions arising out of methods of ascertainment. Sankhya: The Indian Journal of Statistics, Series A, 311-324.

[4] Patil, G. P. and Rao, C. R. The weighted distributions: a survey of their applications. In: applications of statistics, Ed. by P. R. Krishiaiah, 383-405. North Holland Publishing Company, Netherlands, 1977.

[5] Patil G. P and Rao, C. R. (1978): Weighted distributions and size-biased sampling with applications to wildlife populations and human families. Biometrics, 179-189.

[6] Stene, J. (1981). Probability distributions arising from the ascertainment and the analysis of data on human families and other groups. Statistical distributions in scientific work, applications in physical, Social and life Sciences, 6, 233-244.

[7] Bader, M. and Priest, A. (1982). Statistical Aspects of Fiber and Bundle Strength in Hybrid Composites. In: Hayashi, T., Kawata, S. and Umekawa, S., Eds., Progress in Science and Engineering Composites, ICCM-IV, Tokyo, 1129-1136.

[8] Gupta, R.C. and Keating, J.P. (1985). Relations for reliability measures under length biased sampling. Scan. J. Statist., 13, 49-56.

[9] Khattree, R. (1989). Characterization of inverse Gaussian-and gamma distribution through their length-biased distributions. IEEE Trans. Reliab., 38: 610-611.

[10] Gupta, R. C., and S. N. V. A. Kirmani (1990). The role of weighted distribution in stochastic modeling. Commun. Stat. Theory Methods, 19(9), 3147-3162.

[11] Shaked, M., Shanthikumar, J. G. Stochastic Orders and Their Applications. Boston: Academic Press 1994.

[12] Oluyede, B.O. and George, E.O. (2002). On Stochastic Inequalities and Comparisons of Reliability Measures for Weighted Distributions. Mathematical problems in Engineering, 8, 1-13.

[13] Gupta, R. D. and Kundu, D. (2009). A new class of weighted exponential distributions. Statistics, 46(6), 621-634.

[14] Das, K. K., Roy, T. D. (2011a). Applicability of length-biased weighted generalized Rayleigh distribution. Advances in Applied Science Research, 2(4), 320-327.

[15] Das, K. K., Roy, T. D. (2011b). On some length-biased weighted Weibull distribution. Advances in Applied Science Research, 2(5), 465-475.

[16] Ghitany, M. E., Alqallaf, F., Al-Mutairi, H.A. (2011). A two parameter weighted Lindley distribution and its applications to survival data. Mathematics and Computers in Simulation, 81, 11901201.

[17] Ratnaparkhi, M. V. and Naik-Nimbalkar, U. V. (2012). The length-biased lognormal distribution and its application in the analysis of data from oil field exploration studies. Journal of Modern Applied Statistical Methods, 11, 225-260.

[18] Al-Kadim, K. A., Ali Hussein, N. (2014). New proposed length-biased weighted exponential and Rayleigh distribution with application. Mathematical Theory and Modelling, 4(7).

[19] Alqallaf, F., Ghitany, M. E. and Agostinelli, C. (2015). Weighted exponential distribution: Different methods of estimations. Applied Mathematics & Information Sciences, 9(3), 1167.

[20] Ahmad, A., Ahmad, S. P. and Ahmad, A. (2016). Length-biased weighted Lomax distribution: statistical properties and application. Pakistan Journal of Statistics and Operation Research, 12, 245-255.

[21] Oguntunde, P.E., Ilori, K. A. and Okagbue, H. I. (2018). The Inverted weighted exponential distribution with applications. International Journal of Advanced and Applied Sciences, 5(11), 46-50.

[22] Ramos PL, Almeida MP, Tomazella VL, Louzada F. (2019). Improved Bayes estimators and prediction for the Wilson-Hilferty distribution. Anais da Academia Brasileira de Ciencias, 91(3).

[23] Atikankul, Y., Thongteeraparp, A., Bodhisuwan, W., Volodin, A. (2020). The length-biased weighted Lindley distribution with applications, Lobachevskii Journal of Mathematics, 41, 308319.

[24] Mathew, J. and George, S. (2020). Length biased exponential distribution as a reliability model: a Bayesian approach, Reliability: Theory & Applications, 3(58):84-91.

[25] Ilori, A. K., and Jolayemi, E. T. (2021). The weighted exponentiated inverted exponential distribution. International Journal of Statistics and Applied Mathematics, 6(1), 45-50.

[26] Almetwally, E. M. (2022). Marshall Olkin alpha power extended Weibull distribution: Different methods of estimation on type-I and type-II censoring. Gazi University Journal of Science, 35(1), 293-312.

[27] Bhat, V. A. and Pundir, S. (2023): Weighted intervened exponential distribution as a lifetime distribution, Reliability: Theory & Applications, 2(73), 253-266.

[28] Shanker, R. and Shukla, K. K. (2023). Power weighted Sujatha distribution with properties and application to survival times of patients of head and neck cancer data, Reliability:

Theory & Applications, 3(74), 568-581.