In this paper a two–parameter weighted Garima distribution which includes one–parameter Garima distribution has been proposed for modeling real lifetime data. Statistical properties of the distribution including shapes of probability density function, moments and moment related measures, hazard rate function, mean residual life function and stochastic orderings have been discussed. The estimation of its parameters has been discussed using the method of maximum likelihood. Application of the proposed distribution has been discussed.

Keywords: garima distribution, moments, hazard rate function, mean residual life function, stochastic ordering, maximum likelihood estimation, goodness of fit

Fisher¹ firstly introduced the idea of weighted distributions to model ascertainment biases which were later reformulated by Rao² in a unifying theory for problems where the observations fall in non–experimental, non–replicated and non–random. When a researcher collects observations in the nature according to certain stochastic model, the distribution of the collected observations will not have the original distribution unless every observation has been given an equal chance of being included. Suppose the original observation $x_0$ comes from a population having a probability density function (pdf.), $f_o(x,{\theta }_1)$, where ${\theta }_1$ may be a parameter vector; and observation $x$ is collected according to a probability re–weighted with a weight function $w(x,{\theta }_2)>0$,${\theta }_2$ being a new parameter vector, then $x$ comes from a population having pdf

$f\left(x;{\theta }_1,{\theta }_2\right)=k.w\left(x;{\theta }_2\right)f_o\left(x;{\theta }_2\right)$          (1.1)

where $k$ is a normalizing constant. Recall that such types of distributions are known as weighted distributions. A weighted distribution having weight function $w\left(x,{\theta }_2\right)=x$ is called length–biased distribution. Extensive discussions on some general probability models leading to weighted probability distributions and their applications and the occurrence of $\stackrel{⌢}{\theta }$ in the problems in sampling have been discussed by Patil & Rao.^3,4

Shanker⁵ proposed a lifetime distribution named Garima distribution for modeling behavioral Science data and defined by its pdf and cumulative distribution function (cdf)

$f_1\left(x;\theta \right)=\frac{\theta }{\theta +2}\left(1+\theta +\theta x\right)e^{-\theta x};x>0,\theta >0$             (1.2)

$F_1\left(x,\theta \right)=1-\left[1+\frac{\theta x}{\theta +2}\right]e^{-\theta x};x>0,\theta >0$            (1.3)

The first four raw moments (moments about origin) of Garima distribution obtained by Shanker⁵ are

$${\mu }_1{}^{\prime }=\frac{\theta +3}{\theta \left(\theta +2\right)},{\mu }_2{}^{\prime }=\frac{2\left(\theta +4\right)}{{\theta }^2\left(\theta +2\right)},{\mu }_3{}^{\prime }=\frac{6\left(\theta +5\right)}{{\theta }^3\left(\theta +2\right)},{\mu }_4{}^{\prime }=\frac{24\left(\theta +6\right)}{{\theta }^4\left(\theta +2\right)}$$

The central moments (moments about the mean) of Garima distribution obtained by Shanker⁵ are given by

$$\begin{array}{l}{\mu }_2=\frac{{\theta }^2+6\theta +7}{{\theta }^2{\left(\theta +2\right)}^2}\\ {\mu }_3=\frac{2\left({\theta }^3+9{\theta }^2+21\theta +15\right)}{{\theta }^3{\left(\theta +2\right)}^3}\\ {\mu }_4=\frac{3\left(3{\theta }^4+36{\theta }^3+134{\theta }^2+204\theta +111\right)}{{\theta }^4{\left(\theta +2\right)}^4}\end{array}$$

Statistical properties including shapes for different values of parameter, hazard rate function, mean residual life function, stochastic ordering, mean deviations, order statistic, Bonferroni and Lorenz curves, Renyi entropy measures and stress–strength reliability of Garima distribution have been discussed in Shanker.⁵ Estimation of parameter using both the method of moment and the method of maximum likelihood along with application of Garima distribution has been explained in Shanker.⁵ The suitability, superiority and application of Garima distribution for modeling behavioral science data over exponential and Lindley distribution, introduced by Lindley,⁶ have also been discussed by Shanker.⁵ Ghitany et al.,⁷ have detailed study on statistical properties, estimation of parameter and application of Lindley distribution to model waiting time in a bank and established that Lindley distribution gives much closer fit than exponential distribution. Further, Ghitany et al.,⁸ introduced a two–parameter weighted Lindley distribution (WLD) defined by its pdf and cdf

$f_2(x;\theta ,\beta )=\frac{{\theta }^{\beta +1}}{\left(\theta +\beta \right)\Gamma \left(\beta \right)}{x}^{\beta -1}\left(1+x\right)e^{-\theta x};x>0,\theta >0,\beta >0$       (1.4)

$F_2(x;\theta ,\beta )=1-\frac{\left(\theta +\beta \right)\Gamma \left(\beta ,\theta x\right)+{\left(\theta x\right)}^{\beta }{e}^{-\theta x}}{\left(\theta +\beta \right)\Gamma \left(\beta \right)};x>0,\theta >0,\beta >0$        (1.5)

where $\Gamma \left(\beta ,\theta x\right)$ is the upper incomplete gamma function defined as

$\Gamma \left(\beta ,z\right)={\displaystyle \underset{z}{\overset{\infty }{\int }}{e}^{-y}}{y}^{\beta -1}dy;y\ge 0,\beta >0$       (1.6)

It can be easily shown that Lindley distribution is a particular case of WLD at $\beta =1$. Shanker et al.,⁹ proposed an extension of WLD named a three–parameter weighted Lindley distribution which includes Lindley distribution and WLD as particular cases.

Shanker¹⁰ has obtained discrete Poisson–Garima distribution (PGD), discussed its statistical properties, estimation of parameter and applications for count data from biological sciences. Shanker & Shukla^11,12 have also proposed size–biased Poisson–Garima distribution (SBPGD) and Zero–truncated Poisson–Garima distribution (ZTPGD), along with estimation of their parameter and applications for data which structurally excludes zero counts.

The organizations of the rest of the paper are as follows: Section 2 deals with two–parameter weighted Garima distribution (WGD) which includes one–parameter Garima distribution proposed by Shanker⁵ along with the shapes of the pdf and the cdf of WGD. Section 3 deals with statistical constants and associated measures of WGD including coefficient of variation, skewness, kurtosis, index of dispersion. Section 4 deals with the reliability properties of WGD including hazard rate function, mean residual life function and stochastic ordering. Section 5 deals with the maximum likelihood method for the estimation of parameters of the distribution. Section 6 deals with the goodness of fit of the proposed distribution with a real lifetime data over other one parameter and two–parameter lifetime distributions. Finally, the conclusions of the paper have been presented in section 7.

The pdf of the weighted Garima distribution (WGD) can be expressed as

$f_3\left(x;\theta ,\beta \right)=K{x}^{\beta -1}{f}_o\left(x;\theta \right);x>0,\theta >0,\beta >0$ (2.1)

where, K is the normalizing constant and $f_o(x;\theta )$is the pdf of Garima distribution given in (1.2). Thus the pdf of WGD can be obtained as

$f_3\left(x;\theta ,\beta \right)=\frac{{\theta }^{\beta }}{\left(\theta +\beta +1\right)}\frac{x^{\beta -1}}{\Gamma \left(\beta \right)}\left(1+\theta +\theta x\right)e^{-\theta x};x>0,\theta >0,\beta >0$  (2.2)

where

$\Gamma \left(\beta \right)={\displaystyle \underset{0}{\overset{\infty }{\int }}{e}^{-y}}{y}^{\beta -1}dy;y>0,\beta >0$ is the complete gamma function.

We say that $X$ follows WGD with parameters $\theta$ and $\beta$ if its pdf is given by (2.2) and we denote it by $X~\text{WGD}\left(\theta ,\beta \right)$. It can be easily verified that Garima distribution and size–biased Garima distribution (SBGD) are particular cases of WGD at $\beta =1$ and $\beta =2$, respectively. The behavior of the pdf of WGD for different combinations of parameters $\theta$ and $\beta$ are shown in Figure 1.

Figure 1 Behavior of the pdf of WGD for different combinations of the parameters $\theta$ and $\beta$.

The cdf of the WGD can be expressed as

$F_3\left(x;\theta ,\beta \right)=1-\frac{{\left(\theta x\right)}^{\beta }{e}^{-\theta x}+\left(\theta +\beta +1\right)\Gamma \left(\beta ,\theta x\right)}{\left(\theta +\beta +1\right)\Gamma \left(\beta \right)};x>0,\theta >0,\beta >0$    (2.3)

where $\Gamma \left(\beta ,\theta x\right)$ is the upper incomplete gamma function defined in (1.6). Behavior of the cdf of the WGD for different combinations of the parameters $\theta$ and $\beta$ are shown in Figure 2.

Figure 2 Behavior of the cdf of WGD for different combinations of the parameters $\theta$ and $\beta .$

The rth raw moments (moment about origin), ${\mu }_r{}^{\prime }$,of WGD (2.2) can be derived as

$${\mu }_r{}^{\prime }=E\left(X^r\right)={\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^r{f}_2\left(x;\theta ,\beta \right)dx={\displaystyle \underset{0}{\overset{\infty }{\int }}{x}^r\frac{x^{\beta -1}}{\left(\theta +\beta +1\right)}\frac{{\theta }^{\beta }}{\Gamma \left(\beta \right)}\left(1+\theta +\theta x\right)e^{-\theta x}dx}}$$

$=\frac{\left(\theta +\beta +r+1\right)\Gamma \left(\beta +r\right)}{{\theta }^r\left(\theta +\beta +1\right)\Gamma \left(\beta \right)};r=1,2,3,...$              (3.1)

Thus the first four raw moments of WGD are obtained as

$${\mu }_1{}^{\prime }=\frac{\beta \left(\theta +\beta +2\right)}{\theta \left(\theta +\beta +1\right)}$$

$${\mu }_2{}^{\prime }=\frac{\beta \left(\beta +1\right)\left(\theta +\beta +3\right)}{{\theta }^2\left(\theta +\beta +1\right)}$$

$${\mu }_3{}^{\prime }=\frac{\beta \left(\beta +1\right)\left(\beta +2\right)\left(\theta +\beta +4\right)}{{\theta }^3\left(\theta +\beta +1\right)}$$

$${\mu }_4{}^{\prime }=\frac{\beta \left(\beta +1\right)\left(\beta +2\right)\left(\beta +3\right)\left(\theta +\beta +5\right)}{{\theta }^4\left(\theta +\beta +1\right)}$$

Using relationship ${\mu }_r=E{\left(X-{\mu }_1{}^{\prime }\right)}^r={\displaystyle \sum _{k=0}^r\left(\begin{array}{l}r\\ k\end{array}\right)}{\mu }_k{}^{\prime }{\left(-{\mu }_1{}^{\prime }\right)}^{r-k}$ between central moments and raw moments, the central moments of WGD are

$${\mu }_2=\frac{\beta \left\{{\theta }^2+\left(2\beta +4\right)\theta +\left({\beta }^2+3\beta +3\right)\right\}}{{\theta }^2{\left(\theta +\beta +1\right)}^2}$$

$${\mu }_3=\frac{2\beta \left\{{\theta }^3+\left(3\beta +6\right){\theta }^2+\left(3{\beta }^2+9\beta +9\right)\theta +\left({\beta }^3+4{\beta }^2+6\beta +4\right)\right\}}{{\theta }^3{\left(\theta +\beta +1\right)}^3}$$

$${\mu }_4=\frac{3\beta \left\{\begin{array}{l}\left(\beta +2\right){\theta }^4+\left(4{\beta }^2+16\beta +16\right){\theta }^3+\left(6{\beta }^3+34{\beta }^2+58\beta +36\right){\theta }^2\\ +\left(4{\beta }^4+28{\beta }^3+68{\beta }^2+72\beta +32\right)\theta \\ +\left({\beta }^5+8{\beta }^4+25{\beta }^3+38{\beta }^2+29\beta +10\right)\end{array}\right\}}{{\theta }^4{\left(\theta +\beta +1\right)}^4}$$

It can be verified that at $\beta =1$, the moments about origin and the moments about mean of WGD reduces to the corresponding moments of Garima distribution.

The coefficient variation (C.V.) coefficient of skewness$\left(\sqrt{{\beta }_1}\right)$, coefficient of kurtosis $\left({\beta }_2\right)$ and index of dispersion $\left(\gamma \right)$ of WGD are thus given as

$$C.V=\frac{\sigma }{{{\mu }^{\prime }}_1}=\frac{\sqrt{{\theta }^2+\left(2\beta +4\right)\theta +\left({\beta }^2+3\beta +3\right)}}{\sqrt{\beta }\left(\theta +\beta +2\right)}$$

$$\sqrt{{\beta }_1}=\frac{{\mu }_3}{{\mu }_2^{3/2}}=\frac{2\left\{{\theta }^3+\left(3\beta +6\right){\theta }^2+\left(3{\beta }^2+9\beta +9\right)\theta +\left({\beta }^3+4{\beta }^2+6\beta +4\right)\right\}}{\sqrt{\beta }{\left\{{\theta }^2+\left(2\beta +4\right)\theta +\left({\beta }^2+3\beta +3\right)\right\}}^{3/2}}$$

$${\beta }_2=\frac{{\mu }_4}{{\mu }_2^2}=\frac{3\left\{\begin{array}{l}\left(\beta +2\right){\theta }^4+\left(4{\beta }^2+16\beta +16\right){\theta }^3+\left(6{\beta }^3+34{\beta }^2+58\beta +36\right){\theta }^2\\ +\left(4{\beta }^4+28{\beta }^3+68{\beta }^2+72\beta +32\right)\theta \\ +\left({\beta }^5+8{\beta }^4+25{\beta }^3+38{\beta }^2+29\beta +10\right)\end{array}\right\}}{\beta {\left\{{\theta }^2+\left(2\beta +4\right)\theta +\left({\beta }^2+3\beta +3\right)\right\}}^2}$$

$$\gamma =\frac{{\sigma }^2}{{{\mu }^{\prime }}_1}=\frac{{\theta }^2+\left(2\beta +4\right)\theta +\left({\beta }^2+3\beta +3\right)}{\theta \left(\theta +\beta +1\right)\left(\theta +\beta +2\right)}$$

Behaviors of coefficient of variation (C.V), coefficient of skewness $\left(\sqrt{{\beta }_1}\right)$, coefficient of kurtosis $\left({\beta }_2\right)$ and index of dispersion $\left(\gamma \right)$ of WGD have been prepared for different values of $\theta$ and $\beta$ are presented in Tables 1, 2, 3, 4.

It is obvious from Table 1 that for a given $\beta$, the C.V. increases as the $\theta$ increases, whereas for a given$\theta$, the C.V. decreases as the value of $\beta$ increases. It is obvious from Table 2 that for a given $\theta$($\beta$), $\left(\sqrt{{\beta }_1}\right)$ decreases (increases) as the $\beta$($\theta$) increases. It is obvious from Table 3 that for a given $\theta$($\beta$), the coefficient of Kurtosis (${\beta }_2$) decreases (increases) as the $\beta$($\theta$) increases.

It is obvious from Table 4 that for a given $\beta$, the index of dispersion decreases as $\theta$ increases. Similarly, for a given$\beta$, the index of dispersion decreases as $\beta$ increases.

In this section three important reliability properties namely hazard rate function, mean residual life function and stochastic ordering of WGD has been discussed

a. Hazard rate function

The survival (reliability) function of WGD can be expressed as

$S\left(x;\theta ,\beta \right)=1-F_3\left(x;\theta ,\beta \right)=\frac{{\left(\theta x\right)}^{\beta }{e}^{-\theta x}+\left(\theta +\beta +1\right)\Gamma \left(\beta ,\theta x\right)}{\left(\theta +\beta +1\right)\Gamma \left(\beta \right)}$ (4.1.1)

The hazard (or failure) rate function, $h(x)$ of WGD is thus expressed as

$h\left(x\right)=\frac{f_3\left(x;\theta ,\beta \right)}{S\left(x;\theta ,\beta \right)}=\frac{x^{\beta -1}{\theta }^{\beta }\left(1+\theta +\theta x\right)e^{-\theta x}}{{\left(\theta x\right)}^{\beta }{e}^{-\theta x}+\left(\theta +\beta +1\right)\Gamma \left(\beta ,\theta x\right)};x>0,\theta >0,\beta >0$    (4.2)

The behavior of $h(x)$ of WGD for different combinations of $\theta$ and $\beta$ are shown in Figure 3.

Figure 3 Behavior of the hazard function of WGD for different combinations $\theta$ and $\beta .$

Mean residual life function

The mean residual life function $\mu \left(x\right)=E\left(X-x|X>x\right)$ of the WGD can be derived as

$\mu \left(x\right)=\frac{1}{S\left(x;\theta ,\beta \right)}{\displaystyle \underset{x}{\overset{\infty }{\int }}t{f}_3\left(t;\theta ,\beta \right)dt-x}$

$=\frac{{\left(\theta x\right)}^{\beta }\left(\theta +\beta +2\right)e^{-\theta x}+\left\{{\beta }^2+\beta \theta +2\beta -\theta x\left(\theta +\beta +1\right)\right\}\Gamma \left(\beta ,\theta x\right)}{\theta \left\{{\left(\theta x\right)}^{\beta }{e}^{-\theta x}+\left(\theta +\beta +1\right)\Gamma \left(\beta ,\theta x\right)\right\}}$

Clearly $\mu \left(0\right)=\frac{\beta \left(\theta +\beta +2\right)}{\theta \left(\theta +\beta +1\right)}={{\mu }^{\prime }}_1$. The behavior of $\mu (x)$of the WGD for different combinations of $\theta$ and $\beta$ are shown in Figure 4.

Figure 4 Behavior of $\mu (x)$ of the WGD for different cpmbinations of $\theta$ and $\beta .$

Figure 5 Fitted pdf plots of considered distribution for the given dataset.

c. Stochastic ordering

A random variable $X$ is said to be smaller than a random variable $Y$ in the

1. Stochastic order $\left(X{\le }_{st}Y\right)$ if $F_X(x)\ge F_Y(x)$ for all x
2. Hazard rate order $\left(X{\le }_{hr}Y\right)$if $h_X(x)\ge h_Y(x)$ for all x
3. Mean residual life order $\left(X{\le }_{mrl}Y\right)$ if $m_X(x)\le m_Y(x)$ for all x
4. Likelihood ratio order $\left(X{\le }_{lr}Y\right)$if $\frac{f_X(x)}{f_Y(x)}$ decreases in x.

The following important interpretations due to Shaked & Shanthikumar¹³ are well known for establishing stochastic ordering of distributions.

$$\begin{array}{l}X{\le }_{lr}Y\Rightarrow X{\le }_{hr}Y\Rightarrow X{\le }_{mrl}Y\\ \Downarrow \\ X{\le }_{st}Y\end{array}$$

The WGD is ordered with respect to the strongest ‘likelihood ratio’ ordering as established in the following theorem.

Theorem

Let $X~\text{WGD}({\theta }_1,{\beta }_1)$ and $Y~\text{WGD}({\theta }_2,{\beta }_2)$. If ${\theta }_1>{\theta }_2$ and ${\beta }_1={\beta }_2$$\left(\text{or}{\beta }_1<{\beta }_2\text{and}{\theta }_1={\theta }_2\right)$, then $X{\le }_{lr}Y$ and thus $X{\le }_{hr}Y$,$X{\le }_{mrl}Y$ and $X{\le }_{st}Y$.

Proof

We have

$$\ln \frac{f_X\left(x;{\theta }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_{2,}{\beta }_2\right)}=\ln \left(\frac{{\theta }_1^{{\beta }_1}\left({\theta }_2+{\beta }_2+1\right)\Gamma \left({\beta }_2\right)}{{\theta }_2^{{\beta }_2}\left({\theta }_1+{\beta }_1+1\right)\Gamma \left({\beta }_1\right)}\right)x^{{\beta }_{{}_1}-{\beta }_2}\left(\frac{1+{\theta }_1+{\theta }_{1}x}{1+{\theta }_2+{\theta }_{2}x}\right)e-\left({\theta }_1-{\theta }_2\right)x$$

Now $\ln \frac{f_X\left(x;{\theta }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_{2,}{\beta }_2\right)}=\ln \left(\frac{{\theta }_1^{{\beta }_1}\left({\theta }_2+{\beta }_2+1\right)\Gamma \left({\beta }_2\right)}{{\theta }_2^{{\beta }_2}\left({\theta }_1+{\beta }_1+1\right)\Gamma \left({\beta }_1\right)}\right)+\left({\beta }_1-{\beta }_2\right)\ln x+\ln \left(\frac{1+{\theta }_1+{\theta }_{1}x}{1+{\theta }_2+{\theta }_{2}x}\right)-\left({\theta }_1-{\theta }_2\right)x$

This gives

$$\frac{d}{dx}\ln \left(\frac{f_X\left(x;{\theta }_1,{\beta }_1\right)}{f_Y\left(x;{\theta }_2,{\beta }_2\right)}\right)=\frac{{\beta }_1-{\beta }_2}{x}+\frac{{\theta }_1-{\theta }_2}{\left(1+{\theta }_1+{\theta }_{1}x\right)\left(1+{\theta }_2+{\theta }_{2}x\right)}-\left({\theta }_1-{\theta }_2\right).$$

Thus, for ${\beta }_1={\beta }_2$ and ${\theta }_1\ge {\theta }_2$, or $\left({\beta }_1<{\beta }_2\text{and}{\theta }_1\ge {\theta }_2\right)$,$\frac{d}{dx}\ln \left(\frac{f_x\left(x;{\theta }_1,{\beta }_1\right)}{f_y\left(x;{\theta }_2,{\beta }_2\right)}\right)<0$. This shows that $X{\le }_{lr}Y$and thus $X{\le }_{hr}Y$,$X{\le }_{mrl}Y$ and $X{\le }_{st}Y$. This shows flexibility of WGD over Garima distribution.

Consider $(x_1,x_2,...,x_n)$ be a random sample of size n from WGD (2.2). The natural log likelihood function can be expressed as

$\ln L=n\left[\beta \ln \theta -\ln \left(\theta +\beta +1\right)-\ln \left(\Gamma \left(\beta \right)\right)\right]+\left(\beta -1\right){\displaystyle \sum _{i=1}^n\ln \left(x_i\right)+{\displaystyle \sum _{i=1}^n\ln \left(1+\theta +\theta x_i\right)}}-n\theta \overline{x}$

The maximum likelihood estimates ($\widehat{\theta }$,$\stackrel{⌢}{\beta }$) of $(\theta ,\beta )$ is the solution of the following log likelihood equations.

$\frac{\partial \ln L}{\partial \theta }=\frac{n\beta }{\theta }-\frac{n}{\theta +\beta +1}+{\displaystyle \sum _{i=1}^n\frac{1+x_i}{1+\theta +\theta x_i}}-n\overline{x}=0$

$\frac{\partial \ln L}{\partial \beta }=n\ln \theta -\frac{n}{\theta +\beta +1}-n\psi \left(\beta \right)+{\displaystyle \sum _{i=1}^n\ln \left(x_i\right)}=0$

where $\overline{x}$ being the sample mean and $\psi (\beta )$ is the digamma function defined as $\psi \left(\beta \right)=\frac{d}{d\beta }\ln \Gamma \left(\beta \right)$. Since these two log likelihood equations are not in closed forms, they cannot be solved analytically. However, the MLE’s $\left(\widehat{\theta },\widehat{\beta }\right)$ of $\left(\theta ,\beta \right)$ can be computed directly by solving the natural log likelihood equation using Newton–Raphson method available in R–software till sufficiently close estimates of $\stackrel{⌢}{\theta }$ and $\stackrel{⌢}{\beta }$ are obtained. In this paper starting values of parameters $\theta$ and $\beta$ are $\theta =0.5$ and $\beta =1.5$, respectively.

In this section a real lifetime data has been considered for the goodness of fit of WGD over one–parameter and two–parameter life time distributions. The data is regarding the tensile strength, measured in GPa, of 69 carbon fibers tested under tension at gauge lengths of 20mm, available in Bader & Priest.¹⁴

1.312     1.314      1.479      1.552      1.700      1.803      1.861      1.865      1.944      1.958      1.966      1.997      2.006      2.021      2.027      2.055      2.063      2.098      2.140      2.179      2.224      2.240      2.253      2.270      2.272      2.274      2.301      2.301      2.359      2.382      2.382      2.426      2.434      2.435      2.478      2.490      2.511      2.514      2.535      2.554      2.566      2.570      2.586      2.629      2.633      2.642      2.648      2.684      2.697      2.726      2.770      2.773      2.800      2.809      2.818      2.821      2.848      2.880      2.954      3.012      3.067      3.084      3.090      3.096      3.128      3.233      3.433      3.585      3.585

The goodness of fit of WGD has been compared with one parameter exponential, Lindley and Garima distributions and two–parameter Gompertz distribution, generalized exponential distribution (GED) introduced by Gupta & Kundu¹⁵, lognormal distribution and WLD. The pdf and cdf of Gompertz distributions, lognormal and GED are presented in Table 5. The ML estimates of parameters, $-2\ln L$, Akaike Information criteria (AIC), K–S statistics and p–value of distributions for the considered dataset are presented in Table 6. The AIC and K–S Statistics are calculated using the formulae: $AIC=-2\ln L+2k$ and $\text{K-S}=\underset{x}{\text{Sup}}\left|F_n\left(x\right)-F_0\left(x\right)\right|$, where $k$ = the number of parameters, n = the sample size, $F_n\left(x\right)$is the empirical (sample) cumulative distribution function, and $F_0\left(x\right)$ is the theoretical cumulative distribution function. The best distribution is the distribution corresponding to lower values of $-2\ln L$, AIC, and K–S statistics.

It is quite obvious from table 6 that WGD is competing well with two–parameter lifetime distributions and gives quite satisfactory fit.

A two–parameter weighted Garima distribution(WGD) which includes one parameter Garima distribution proposed by Shanker⁵ has been suggested for modeling lifetime data from engineering. Its statistical properties including shapes of probability density function for different combinations of parameters, coefficients of variation, skweness, kurtosis, and index of dispersion have been explained. Reliability measures including hazard rate function, mean residual life function and the stochastic ordering have been studied. Estimation of parameters has been discussed using maximum likelihood. The goodness of fit of WGD has been explained with a real lifetime data and the fit has been found to be quite satisfactory over exponential, Lindley and Garima, Gompertz, generalized exponential, lognormal and weighted Lindley distributions.

Authors are grateful to the editor–in–chief of the journal and the anonymous reviewer for fruitful comments on the paper.

Author declares there is no conflict of interest towards this manuscript.

1. Fisher RA. The effects of methods of ascertainment upon the estimation of frequencies. Ann Eugenics. 1934;6(1):13–25.
2. Rao CR. On discrete distributions arising out of methods of ascertainment. In: Patil GP, editor. Classical and Contagious Discrete Distributions. Calcutta: Statistical Publishing Society; 1965. p. 320–332.
3. Patil GP, Rao CR. The Weighted distributions: A survey and their applications. In applications of Statistics. PR Krishnaiah, editor. Amsterdam: North Holland Publications Co; 1977. p. 383–405,
4. Patil GP, Rao CR. Weighted distributions and size-biased sampling with applications to wild-life populations and human families. Biometrics. 1978;34:179–189.
5. Shanker R. Garima distribution and Its Application to model behavioral science data. Biometrics & Biostatistics International Journal. 2016;4(7):1–9.
6. Lindley DV. Fiducial distributions and Bayes’theorem. Journal of the Royal Statistical Society. Series B. 1958;20(1):102–107.
7. Ghitany ME, Atieh B, Nadarajah S. Lindley distribution and its Application. Mathematics Computing and Simulation. 2008;78(4):493–506.
8. Ghitany ME, Alqallaf F, Al-Mutairi DK, et al. A two-parameter weighted Lindley distribution and its applications to survival data. Mathematics and Computers in simulation. 2011;81(6):1190–1201.
9. Shanker R, Shukla KK, Mishra A. A three-parameter weighted Lindley distribution. Statistics in Transition new series. 2017;18(2):291–300.
10. Shanker R. The discrete Poisson-Garima distribution. Biometrics and Biostatistics International Journal. 2017;5(2):1–7.
11. Shanker R, Shukla KK. Size-biased Poisson-Garima distribution with applications. Biometrics and Biostatistics International Journal. 2017;6(3):1-6.
12. Shanker R, Shukla KK. Zero-truncated Poisson-Garima distribution with applications. Biometrics and Biostatistics Open Access Journal. 2017 b;3(1):1-6.
13. Shaked M, Shanthikumar JG. Stochastic Orders and Their Applications. New York: Academic Press; 1994.
14. Bader MG, Priest AM. Statistical aspects of fiber and bundle strength in hybrid composites. In; hayashi T, Kawata K, Umekawa S. editors, ICCM-IV, Tokyo: Progress in Science in Engineering Composites; 1982. p. 1129–1136.
15. Gupta RD, Kundu D. Generalized exponential distribution. Australian and Newzealand Journal of statistics. 1999;41(2):173–188.