On the Estimation of Parameter of Weighted Sums of Exponential Distribution

The random variableZ n,α = Y 1 + 2 α Y 2 + ⋅ ⋅ ⋅ + n α Y n , with α ∈ R and Y 1 , Y 2 , . . . being independent exponentially distributed random variables with mean one, is considered. Van Leeuwaarden and Temme (2011) attempted to determine good approximation of the distribution of Z n,α . The main problem is estimating the parameter α that has the main state in applicable research. In this paper we show that estimating the parameter α by using the relation between α and mode is available. The mean square error values are obtained for estimating α by mode, moment method, and maximum likelihood method.


Introduction
The exponential distribution is one of the most applicable distributions in survival models and phenomena with memoryless property.Random size sample from the exponential distribution and considering the discrete distribution are the basis for creating new distributions.Proschan [1] showed that combinations distributions with constant failure rate have decreasing hazard function.So, in recent years, new distributions are introduced based on generalization and correction of the exponential distribution.For example, Adamidis et al. [2] find a new bivariate distribution, exponential-geometric distribution, with decreasing failure rate that is used in survival models.This distribution combines with exponential and geometric distributions.In the case that parameter of exponential distribution has decreasing-increasing failure rate the new distribution is called exponential-weibull that in fact is a generalized of exponential-geometric distribution introduced in 2010 by Silva et al. [3].In all these studies trough estimating the parameters algorithm EM has been used.
Consider  , =  1 + 2   2 + ⋅ ⋅ ⋅ +     as the random variable with  ∈ R and  1 ,  2 , . . .being independent exponentially distributed random variables with mean one.In special case distribution function of the random variable  , is a simple alternating series.For the case  = 1 accurs in various contexts, such as linear combinations of order statistics, noise in radio receivers, and run models, Kingman and Volkov [4].For all  ≥ −1/2 the random variable  , follows the central limit theorem and leads to random variable normal standard.The result of normal approach of  ( , < ) would be useful for  → ∞ in terms of X's, which are near to mean  , .For  = −1 the random variable  , is independent from  and identity distribution and perfectly shows that the random variable has exponential distribution.One example of extreme value theory is that  , has Gumbel distribution as  → ∞.Because  = −2 the random variable  , is Kolmogorov distribution.In general case, the distribution function is a series with alternating sign that Van Leeuwarden and Temme [5] derived an approximation uniform expansion for P ( , < ) by applying an extended version of the saddle point method.
In this paper in order to estimate the parameter , first we use the exact density of linear combination of independent exponential distribution.Determining the state of distribution for different values of  would lead us to exact value of it.Some results show the simulations between distributions in this family.Also to reach the estimation of parameter  we suggest sample mode.In Section 1 we review 2 Chinese Journal of Mathematics the density function of  , and its characteristics.We will show mode and sum result of it in Section 2. For estimating mode there are numerous formula [6].But here, we used the following formula: where  Mo is a lower bound in mode class,  1 and  2 are absolute difference frequency of mode class before and after it, respectively, and c is a length of mode class.In statistical table, mode class is a class that has high frequency.At the end, using trough mean square error we compare three estimators for , according the sample mode, the moment sample, and maximum likelihood estimator.

Density of 𝑍 𝑛,𝛼
The distribution function of the random variable  , is where Then we can write In special case, for  = 2 the distribution function is as follows: For  ≤ −3 and different values of , the figures of distribution will be shown in Figures 1(a) and 1(b).
Based on the shape, it is clear that distributions are fair proximity to  2, ; that is,  , ≈  2, ; ∀ ≥ 2,  ≤ −3. (6) This point can limit the estimating parameter , so we can estimate it accurately.With some non complex calculation of  2, ()/ = 0, we can obtain the mode: Next theorem is a good result to show the effect of  in distribution function based on mode value.
In Table 1 the distribution function values of  2, in mode point by considering different values of  evaluated and showed in, as it's clearer, results repeated for negative points.Due to Theorem 1, correspondence negative values can be obtained.
Then we can see lim The formula in (7) showed that the relations of mode and mean with parameter  lead us to better estimator.Since we can estimate mode and mean trough sample data, we can find two estimators for parameter .Solution of the equation that known and unknown sides are sample and distribution mode, respectively, lead us to following answer: where lambert (⋅) is a defined function in MAPLE software.The other estimator is obtained from the mean sample by MoM method and expectation of distribution.Let  be the sample mean of  2, 's; then αMoM = ln ( − 1) If  < 1, the moment estimator is not achieved.

Simulating and Comparing the Estimators
The simulation that was performed with a sample size of 100 is considered to have been repeated 100 times.In Table 2,

Table 1 :
The distribution function values by considering different values .

Table 2 :
Comparing MSE values by three methods., the phrase inside the log function that may be negative instead of some parameters , MSE values for all the three estimators method calculated.By numeric method, we compare three estimators, MLE, MoM, and estimator based on distribution mode.As we expected, based on Table2, MSE is more preferred than the two other estimators.While for some values  we cannot submit MoM estimator, MSE for sample mod estimator prefers MoM.Perhaps the use of other estimator modes improved alpha estimation.