Abstract

The maximum likelihood method is the most widely used estimation method. On the other hand, it can produce substantial bias, and an approximate confidence interval based on the maximum likelihood estimator cannot be valid when the sample size is small. Because the sizes of the record values are considerably smaller than the original sequence observed in the majority of cases, a method appropriate for this situation is required for precise inference. This paper provides the exact confidence intervals for unknown parameters and exact predictive intervals for the future upper record values by providing some pivotal quantities in the two-parameter Rayleigh distribution based on the upper record values. Finally, the validity of the proposed inference methods was examined from Monte Carlo simulations and real data.

1. Introduction

The cumulative distribution function (cdf) and probability density function (pdf) of the random variable (RV), , with the Rayleigh distribution are given, respectively, by where is the location parameter and is the scale parameter. The Rayleigh distribution is used because the life of the model theory reliability plays an important role in modeling the life of the random phenomenon. Moreover, it is used in many applications, such as reliability, life tests, and survival analysis because its failure rate is a linear function of time. Therefore, this distribution has been studied by many authors in the case, where samples are censored due to a range of reasons. Dyer and Whisenand [1] examined the properties of the -optimum best linear unbiased estimators (BLUEs) of the scale parameter in the Rayleigh distribution and provided an approximate -optimum BLUE based on order statistics. Sinha and Howlader [2] derived the highest posterior density (HPD) credible interval for the scale parameter and the reliability function in a Rayleigh distribution. Ali Mousa and Al-Sagheer [3] obtained the maximum likelihood estimators (MLEs) and Bayes estimators for , , and the reliability function of the Rayleigh distribution obtained based on the progressively Type-II censored data. Raqab and Madi [4] discussed the Bayesian predictive methods for the total time on test using doubly censored data with a Rayleigh distribution and the scale parameter and applied the methods to a real data set that represented the deep-groove ball bearing failure times. For the same real data, Kim and Han [5] applied a Bayesian inference method based on the conjugate prior of the scale parameter of the Rayleigh distribution under general progressive censoring and S. Dey and T. Dey [6] applied this by providing point and interval estimation methods for the scale parameter of the Rayleigh distribution under progressive Type-II censoring with binomial removal. This paper considered a two-parameter Rayleigh distribution based on the upper record values that are used extensively to build statistical modeling arising in many real-life situations involving weather, sports, economics, and life tests. The record values are described as follows.

Let , be a sequence of independent and identically distributed (iid) RVs from a continuous probability distribution. If for all , then is an upper record value. The indices at which the upper record values occur are given by the record times , , where , , with . Chandler [7] first studied the record values and their basic properties. Ahsanullah [8] provided detailed descriptions of the general theory and applications for the well-known probability distributions based on the records. Seo and Kim [9] provided inference methods to estimate unknown parameters and predicted future upper record values from the extreme value distribution using both frequentist and Bayesian approaches. Note that the sizes of the record values are actually considerably smaller than the observed original sequence in the majority of case; a method appropriate for this situation is required for precise inference. The maximum likelihood method is the most extensively used estimation method. On the other hand, the approximate confidence intervals (CIs) based on the asymptotic normality of the MLE can yield inappropriate results because and are supported by and , respectively. Moreover, the asymptotic normality of the MLE requires the suitable regularity conditions and it is difficult to prove that the regularity conditions are satisfied when the record values are observed from the two-parameter Rayleigh distribution. This paper constructs exact CIs for unknown parameters of the Rayleigh distribution based on the upper record values by providing some pivotal quantities, which are much more efficient than the maximum likelihood method in terms of computation cost. Another aim of this paper is to construct exact predictive intervals (PIs) for the future upper record values based on the past upper record values from the Rayleigh distribution because it is very important to correctly predict in many fields such as earthquakes, flood, and rainfall.

The remainder of the paper is structured as follows. Section 2 provides some pivotal quantities and derives the exact CIs for unknown parameters and PIs for the future upper record values in the Rayleigh distribution based on the upper record values. Section 3 assesses the validity of the proposed method through Monte Carlo simulations and real data. Section 4 concludes the paper.

2. Inference Based on Pivotal Quantity

The likelihood function for is given by Arnold et al., [10] as

Let be the first upper record values from the two-parameter Rayleigh distribution. The likelihood function based on record values is given by

The MLEs and can be found by solving the following likelihood equations for and simultaneously:

On the other hand, the MLEs cannot be expressed in closed form and their exact distributions are difficult to derive. Alternatively, by the asymptotic normality of the MLE, the approximate CIs for and can be obtained as where denotes the upper point of the standard normal distribution and the variances and are the diagonal elements of the asymptotic variance-covariance matrix obtained by inverting the Fisher information matrix for unknown parameters : under certain regularity conditions. Nevertheless, it can provide inappropriate results because the supports of and do not coincide with that of the normal distribution and the record values are rarely observed, as mentioned before.

2.1. Confidence Interval

This subsection develops inference methods based on the pivotal quantities to construct exact CIs for unknown parameters and PIs for the future upper record values. Note that the proposed method is much easier to calculate than the maximum likelihood method. The following provides some pivotal quantities.

Let

are the upper record values with a standard exponential distribution. From this result, the following spacing can be obtained: which are the iid RVs from the standard exponential distribution (see Arnold et al., [10]). Based on the spacing, a pivotal quantity having a distribution with degrees of freedom and a pivotal quantity can be derived as having the distribution with degrees of freedom. Because they have independent RVs, the following pivotal quantity is obtained: which has a distribution with and degrees of freedom. An exact CI for based on the pivotal quantity can be constructed asfor any , where is the upper percentile of the distribution with and degrees of freedom.

Moreover, because has the distribution with degrees of freedom, an exact CI for based on the pivotal quantity can be constructed aswhere is the upper percentile of the distribution with degrees of freedom. Note that because the precise CI (13) depends on the nuisance parameter , this paper shows how to address the nuisance parameter based on a generalized pivotal quantity, and an exact CI for is proposed based on the generalized pivotal quantity.

Let be the unique solution of , where has a distribution with and degrees of freedom. The unique solution can then be given by

Moreover, let be the RV from the distribution with degrees of freedom. The generalized pivotal quantity from the pivotal quantity is given by

Here, the samples can be obtained by generating the RVs and . s are ordered as . Therefore, an exact CI for based on the generalized pivotal quantity can be constructed:where denotes the largest integer less than or equal to . In Section 3, the proposed CIs are examined in terms of the coverage probability (CPs) to determine if they are valid CIs.

2.2. Predictive Interval

This subsection develops a method for predicting the future upper record values based on the observed upper record values by providing a pivotal quantity. Let be a future upper record value. The conditional density function of , given , defined by Ahsanullah [11], is given byfrom the Markov property of the record values. Assuming that the observed upper record values, , arise from the Rayleigh distribution with the pdf (2), the conditional density function (17) is written as

Let

Because the Jacobian of transformation is the density function of is given by which is the pdf of the distribution with degrees of freedom. Suppose that and are known. An exact PI based on the pivotal quantity for the future upper record value is obtained as

When and are unknown, they can be substituted with by and in PI (22) based on the fact that is the generalized pivotal quantity for constructing the exact CI for . In the same way, the generalized pivotal quantity is given by and an exact PI for based on the generalized pivotal quantity can be constructed as follows:

3. Application

This section assesses the proposed methods through a Monte Carlo simulation and presents a real data set.

3.1. Simulation Study

The proposed exact CIs (12) and (16) are assessed in terms of their CPs and average lengths (ALs). The upper record values were first generated from the standard Rayleigh distribution with and for different , and the CIs (12) and (16) were calculated based on the generated samples by using the provided methods in Section 2.1. The CPs and ALs of the exact CIs were obtained over 10,000 simulations. These values are reported in Table 1.

Table 1 shows that the CPs matched their corresponding nominal levels even in a small sample size and that all ALs decrease with increasing sample size.

3.2. Real Data

To illustrate the proposed inference procedure, the survival times in (days) of a group of lung cancer patients (from Lawless [12, p. 319]) were considered as follows:

From the data, the observed upper record values were 6.96, 9.30, 10.18, 11.94, and 12.94. Soliman and Al-Aboud [13] showed that the Rayleigh distribution fits the observed record values well. These record values are employed to obtain the proposed CIs (12) and (16). Moreover, the exact PIs for the future upper record values were computed, as listed in Table 2.

4. Concluding Remarks

This paper proposes methods for inferencing the exact CIs for unknown parameters (, ) in the Rayleigh distribution based on the upper record values and exact PIs for the future upper record values by providing some pivotal quantities. Because the proposed exact CI (13) and PI (22) depend on the nuisance parameters, this study proposed generalized pivotal quantities and to solve the drawback. The proposed methods were more computationally convenient than the maximum likelihood method. Moreover, the proposed exact CIs provide very good performance even in small sample sizes. If the location parameter of the Rayleigh distribution is of interest, the exact CI (12) should be used because it does not have any nuisance parameter.

Competing Interests

The authors declare that they have no competing interests.