Mathematical Problems in Engineering
Volume 2021 (2021), Article ID 9957944, 8 pages
https://doi.org/10.1155/2021/9957944

Research Article

Estimation of Constant Stress Partially Accelerated Life Test for Fréchet Distribution with Type-I Censoring

Abdullah Ali H. Ahmadini,¹ Wali Khan Mashwani,² Rehman Ahmad Khan Sherwani,³ Shokrya S. Alshqaq,¹ Farrukh Jamal,⁴ Miftahuddin Miftahuddin,⁵ Kamran Abbas,⁶ Faiza Razaq,⁶ Mohammed Elgarhy,⁷ and Sanaa Al-Marzouki⁸

¹Department of Mathematics, Faculty of Science, Jazan University, Jazan, Saudi Arabia
²Institute of Numerical Sciences, Kohat University of Science & Technology, Kohat 26000, Pakistan
³Collerge of Statistical and Actuarial Sciences, University of the Punjab, Lahore, Pakistan
⁴Department of Statistics, The Islamia University of Bahawalpur, Bahawalpur 61300, Pakistan
⁵Department of Statistics, Faculty of Mathematics and Sciences, Syiah Kuala University, Banda Aceh, Indonesia
⁶Department of Statistics, University of Azad Jammu and Kashmir, Muzaffarabad, Pakistan
⁷The Higher Institute of Commercial Sciences, Al Mahalla Al Kubra, Algarbia 31951, Egypt
⁸Statistics Department, Faculty of Science, King Abdul Aziz University, Jeddah 21551, Saudi Arabia

Correspondence should be addressed to Miftahuddin Miftahuddin; miftah@unsyiah.ac.id

Received 3 April 2021; Revised 22 April 2021; Accepted 30 April 2021; Published 12 May 2021

Academic Editor: Dr. Dilbag Singh

Copyright © 2021 Abdullah Ali H. Ahmadini et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Modern reliability engineering accelerated life tests (ALT) and partially accelerated life tests (PALT) are widely used to obtain the timely information on the reliability of objects, products, elements, and materials as well as to save time and cost. The ALTs or PALTs are useful in determining the failed manners of the items at routine conditions by using the information of the data generated from the experiment. PALT is the most sensible method to be used for estimating both ordinary and ALTs. In this research, constant stress PALT design for the Fréchet distribution with type-I censoring has been investigated due to a wide applicability of the Fréchet distribution in engineering problems especially in hydrology. The distribution parameters and acceleration factor are obtained by using the maximum likelihood method. Fisher's information matrix is used to develop the asymptotic confidence interval estimates of the model parameters. A simulation study is conducted to illustrate the statistical properties of the parameters and the confidence intervals by using the R software. The results indicated that the constant stress PALT plan works well. Moreover, a numerical example is given to exemplify the performance of the proposed methods.

1. Introduction

The continuous progress in manufacturing design makes the products and materials highly reliable [1]. It is becoming more and more difficult to get information under normal use conditions about the lifespan of objects and materials with higher reliability at the testing time, and it becomes very expensive and time-consuming under such conditions [2]. The sample is examined for these objects at more serious circumstances than usual [3]. These situations belong to stresses, which may be in the form of voltage, force, temperature, vibrations, moisture, etc [4–7]. To develop a new life testing method, these problems have motivated the researcher to get information timely on the reliability of components, objects, and products.

Particularly, when test items are run at both normal and higher than normal stress conditions, the more suitable test to be applied is the partially accelerated life test (PALT) [8–10]. Furthermore, it is the most important method to be utilized for evaluating the acceleration factor and thus extrapolating the accelerated test data to use condition. The PALT combines both accelerated and ordinary life tests. Thus, PALT is reasonable for estimating the acceleration factor (β > 1), “which is the ratio of the hazard rate at the accelerated condition to that at normal use conditions” [11].

The commonly used methods in PALT are constant stress PALT (CSPALT) [4, 12, 13], step stress PALT (SSPALT) [5, 14–16], and progressive stress PALT (PSPALT) [17–19]. In CSPALT, products are examined at either accelerated condition or at normally used condition before the test is being terminated.

Additionally, censoring is very ordinary in life testing experiments. It frequently happens that the experiment is censored when the experimenter may not be able to examine the lifespan of all units put on the test because of time restrictions and other limitations on the data collection. It is usually used when a distribution of exact lifetimes is known only for a part of the test items, and the remainder of the lifetimes is known only to exceed certain values under a life test. The two most common types of censoring are type-I censoring (time censoring) and type-II censoring (failure censoring).

For a brief overview of PALT, literature is abundant on designing PALT. Bai and Chung [8, 9] discussed both the problem of estimation and optimally designing constant stress PALT for exponential and lognormal distributions. Ismail [20] has applied maximum likelihood (ML) and Bayesian methods for evaluating the parameters of Pareto distribution of the second kind and the acceleration factor. Ismail et al. [21] considered the case of CSPALT when under type-I censoring two stress levels was included. Similarly, Cheng and Wang [22], Ismail [23], Kamal et al. [13], Wang [24, 25], and Hassan et al. [12] applied PALT and acceleration factors using different lifetime distributions; PALT and CSPALT have been studied via Bayesian approach [19, 26–28].

In this paper, the CSPALT plan using type-I censoring with the assumption that the lifetimes of the test items at use condition follow Fréchet distribution which is considered. Maurice Fréchet was a French mathematician who had identified one possible limit distribution for the largest order statistic in 1927 [29]. The probability density function (PDF) of Fréchet distribution is

where α is the shape parameter and θ is the scale parameter. This distribution has a heavy upper tail and is bounded on the lower side (t > 0). The cumulative distribution function (CDF) is

The reliability function of Fréchet distribution is

The Fréchet distribution is applied to extreme events and linked to the modeling of several real-world phenomena, including human lifetimes, flood and seismic analyses, radioactive emissions, and maximum one-day rainfalls.

The rest of the study is arranged as follows. In Section 2, the assumption and test procedures for Fréchet distribution with constant stress using type-I censoring are discussed. The maximum likelihood estimates (MLEs) of the distribution parameters and acceleration factor are presented in Section 3. Fisher information matrix and asymptotic confidence interval are discussed in Section 4. A simulation study is presented in Section 5. Finally, Section 6 presents the results, conclusion, and recommendations for future work.

2. Model Description

The following notations are used for model description:

(i)n = total number of objects in a PALT(ii)τ = censoring time of CSPALT(iii)t = lifespan of an object at normal use condition(iv)X = lifespan of an object at accelerated condition(v)β = acceleration factor (β > 1) and defined as β = T/X(vi) = observed lifetime of object at usual use condition(vii) = observed lifetime of unit at accelerated condition(viii) ≤ … ≤  ≤ τ = ordered failure time at use condition(ix) ≤ … ≤  ≤ τ = ordered failure time at accelerated condition

2.1. Test Procedure of CSPALT

In CSPALT, total test units “n” are randomly divided into two samples of sizes “nr” and n(1 − r), respectively, where “r” is the sample proportion. In the first sample, “nr” items are randomly chosen among “n” items, which are assigned to accelerated conditions and the other n(1 − r) are assigned to regular use conditions. Each test object of every sample is run until the censoring time “τ” and the test condition is not changed.

2.2. Assumptions

The following assumptions are also made for CSPALT.

(i)Under normal-use condition, lifespan of an object supports Fréchet distribution.(ii)The lifespan of an object at accelerated condition is attained by X = β⁻¹T, where β > 1, which is an acceleration factor. Therefore, the PDF at accelerated condition is

(iii)The lifespan of objects allocated to the normal condition are independently and identically distributed (i.i.d) random variables.(iv)The lifespan of objects assigned to accelerate are i.i.d. random variables.(v)The lifetimes and are mutually independent.

3. Parameter Estimation

The ML method is one of the most important and widely used methods in statistics. The main reasons are that the ML method is very vigorous and provides the estimates of parameters with good statistical properties such as consistency, asymptotic unbiased, asymptotic efficiency, and asymptotic normality.

Let and be the indicator functions under the use and accelerated conditions, respectively, such that

The likelihood function for can be expressed as

where and .

The total likelihood for is

It is more convenient to work with the log-likelihood function. The log-likelihood function of (8) is

Let and be the numbers of items that failed at normal and accelerated conditions, respectively. Similarly, let and be the numbers of items censored at normal and accelerated conditions, respectively:

MLEs of α, θ, and β are obtained by solving the system of equations , respectively:

The above equations cannot be written in the closed form. So, an iterative procedure can be used to get MLEs. Here, we use the BFGS Quasi-Newton Optimization method which is available in the R software. Furthermore, for interval estimates of model parameters, we need the Fisher information matrix.

3.1. Fisher Information Matrix

Here, it is hard to obtain an explicit solution to nonlinear equations, so an iterative method such as the Newton Raphson method [30] is used to find the solution of equations to obtain MLEs. Also, we get the asymptotic confidence intervals. Hence, the asymptotic variance of the MLE can be attained by the inverse of the observed Fisher information matrix:

The 100 (1 − α) percent asymptotic confidence interval for α, θ, and β can be written as

respectively.

4. Simulation Study

A simulation study is carried out to evaluate the performance of proposed estimators. It is performed by using the R software (3.2.3) for demonstrating the theoretical outcomes of the estimation problems. The performances of MLEs are evaluated through mean square error (MSE) and absolute bias of estimates. Also, 90% and 95% confidence limits are constructed for parameters and acceleration factors. For this purpose, several data sets are generated from two-parameter Fréchet distribution with sizes . Samples are generated by inverse transformation of the Fréchet distribution and the process is replicated 1000 times.

The simulation steps are summarized below:

(1)For estimation, the sample size is varied to see the effect of large and small samples are considered. The random samples are generated from Fréchet distribution by using transformation , where u is uniformly distributed, i.e., U ∼ uniform (0, 1).(2)The values for true parameters are taken as (α = 0.75, β = 1.2, θ = 1), (α = 1, β = 1.1, θ = 1), and (α = 0.8, β = 1.2, θ = 0.8), respectively. The selected censoring level is 0.7 ∗ n(1 − r) and 0.3 ∗ nr.(3)The “n” test item is split with sample proportion r = 0.3, such that “nr” items are assigned at the accelerated condition and n(1 − r) are allotted to normal use condition.(4)The estimates of the model parameters and their corresponding summary statistics are obtained by the CSPALT model.(5)For the different combination of stresses, the MLEs for α, β, and θ, MSE, and absolute bias of the estimates and ninety percent and ninety-five percent CIs are derived.(6)Tables 1–3 summarize the results of the ML estimates for α, β, and θ. The numerical results presented in Table 4 are based on 1000 simulation replications.(7)Figures 1–3 present the MSEs against different sample sizes for the MLEs of α, β, and θ.

Table 1: Average estimates of CSPLAT with α = 0.75, β = 1.2, θ = 1, and r = 0.3.

η
Parameters
Estimates
MSE
Absolute bias
95% CI

50
α
1.5870
0.8109
0.8370
−0.8226, 3.9966

β
1.1830
0.1229
0.0169
0.49494, 1.8711

θ
0.4159
0.3485
0.5840
−1.2117, 2.0436

100
α
1.4755
0.5836
0.7352
−0.5724, 3.5234

β
1.2635
0.0715
0.0531
0.7797, 1.7473

θ
0.4147
0.3465
0.5858
−1.2112, 2.0408

150
α
1.4603
0.53274
0.7103
−0.5569, 3.5019

β
1.2914
0.04893
0.0776
0.8076, 1.7464

θ
0.4136
0.34611
0.5864
−1.2172, 2.0429

200
α
1.4386
0.49254
0.68863
−0.5043, 3.3997

β
1.2992
0.04160
0.09918
0.8576, 1.7275

θ
0.4126
0.34663
0.58743
−1.2125, 2.0402

300
α
1.4221
0.46364
0.67210
−0.4592, 3.3088

β
1.3089
0.03655
0.10886
0.8753, 1.7414

θ
0.4128
0.34592
0.58724
−1.2159, 2.0416

Table 2: Average estimates of CSPLAT with α = 1, β = 1.1, θ = 1, and r = 0.3

η
Parameters
Estimates
MSE
Absolute bias
95% CI

50
α
2.1466
1.4601
1.12465
−1.1840, 5.5040

β
0.7112
0.46457
0.38876
−0.9333, 2.2885

θ
0.5080
0.24811
0.49196
−0.8467, 1.8718

100
α
1.9824
1.04411
0.98242
−0.8141, 4.7996

β
0.7513
0.28392
0.34868
−0.4538, 1.9765

θ
0.5161
0.23707
0.48393
−0.8336, 1.8634

150
α
1.9443
0.94151
0.94430
−0.7090, 4.5968

β
0.8099
0.23370
0.29013
−0.2890, 1.9258

θ
0.5156
0.23648
0.48444
−0.8350, 1.8635

200
α
1.9318
0.90238
0.93182
−0.6499, 4.4828

β
0.8066
0.21864
0.29344
−0.2114, 1.9076

θ
0.5160
0.23574
0.48405
−0.8345, 1.8628

300
α
1.9043
0.84035
0.90430
−0.6121, 4.4129

β
0.8315
0.18293
0.26852
−0.1070, 1.8217

θ
0.5166
0.23464
0.48344
−0.8333, 1.8619

500
α
1.8991
0.82207
0.89914
−0.5948, 4.3840

β
0.8488
0.16980
0.25116
−0.0358, 1.7549

θ
0.5145
0.23628
0.48553
−0.8374, 1.8635

Table 3: Average estimates of CSPLAT with α = 0.8, β = 1.2, θ = 0.8, and r = 0.3

η
Parameters
Estimates
MSE
Absolute bias
95% CI

50
α
1.7092
0.9542
0.9092
−0.9064, 4.3248

β
1.2032
0.2057
0.0032
0.3143, 2.0921

θ
0.3449
0.2111
0.4551
−0.9225, 1.6123

100
α
1.5806
0.6579
0.7806
−0.6258, 3.7871

β
1.2256
0.1198
0.0257
0.5454, 1.9059

θ
0.3519
0.2028
0.4480
−0.8932, 1.5971

150
α
1.5592
0.6060
0.7592
−0.5721, 3.6905

β
1.2563
0.1070
0.0563
0.6059, 1.9067

θ
0.3465
0.2070
0.4535
−0.9124, 1.6055

200
α
1.5407
0.5723
0.7407
−0.5344, 3.6158

β
1.2771
0.0829
0.0772
0.6928, 1.8614

θ
0.3498
0.2058
0.4502
−0.8998, 1.5994

300
α
1.5204
0.5324
0.7203
−0.4893, 3.5301

β
1.2906
0.0682
0.0905
0.7487, 1.8324

θ
0.3481
0.2049
0.4519
−0.9058, 1.6019

500
α
1.5081
0.5092
0.7081
−0.4623, 3.4784

β
1.2925
0.0534
0.0926
0.8045, 1.7806

θ
0.3467
0.2040
0.4502
−0.9103, 1.6038

Table 4: Interval estimates of the parameters.

Parameters
90% CI
95% CI

α
−0.6142, 3.7657
−0.2623, 3.4137

β
0.26319, 2.5764
0.4491, 2.3905

θ
−0.5521, 0.9953
−0.4277, 0.8710

Figure 1: Graph of MSE versus sample size when α = 0.75, β = 1.2, and θ = 1.

Figure 2: Graph of MSE versus sample size when α = 1, β = 1.1, and θ = 1.

Figure 3: Graph of MSE versus sample size when α = 0.8, β = 1.2, and θ = 0.8.

5. Numerical Example

For illustration, we assume that the values of Fréchet constant stress model parameters are known and use the simulated data based on n = 100, α = 0.8, β = 1.3, θ = 0.5, and r = 0.3, and at accelerated condition, . The MLEs of the model parameters and acceleration factor obtained by using the R software are . The inverse of the estimated Fisher information matrix is

The 90% and 95% interval estimates for the model parameters and acceleration factor are provided in Table 4.

The asymptotic variance-covariance matrix of , θ, and β are

The 95% asymptotic confidence interval for α, θ, and β are presented, respectively:

From the results, we observe that the point estimates are stable. We also examine that the results support theoretical findings of CSPALT for Fréchet distribution.

6. Conclusion

From Tables 1–3, it is observed that the MLEs converge to the true values of the parameters with an increase in sample size “n.” Also, we find that the MLEs have good statistical properties for a fixed value of parameters, as we increase sample size, the MSE and bias of the estimators are decreased. Moreover, from Tables 1–3, the MSEs of α decrease as we increase the sample size. It can be seen from results in Tables 1–3 that when there is an increment in sample size, the interval of the estimates decreases. The outcomes of simulation suggest that PALT is an appropriate model using Fréchet distribution which enables us to save cost and time without using a higher level of stress to all test items. From simulation results, it is concluded that the MLEs have good statistical properties.

7. Recommendation

Some pertinent suggestions are also given for future research which is tantamount to provide a pathway for future researchers in the field of PALT. This work can be extended for type-II censoring of CSPALT using Fréchet distribution.

Data Availability

The data used in this article are freely available upon request from the authors and citing this paper in your manuscripts.

Conflicts of Interest

The authors declare that they have no conflicts of interest to report regarding the present study.

References

C. A. Lloyd, T. Roulstone, and R. E. Lyons, “Transport, constructability, and economic advantages of SMR modularization,” Progress in Nuclear Energy, vol. 134, p. 103672, 2021.

X. Shi, P. Lu, and Y. Shi, “Reliability estimation for hybrid system under constant-stress partially accelerated life test with progressively hybrid censoring,” Recent Patents on Engineering, vol. 14, no. 1, pp. 82–94, 2020.

V. S. Yakovyna, M. M. Seniv, I. I. Symets, and N. B. Sambir, “Algorithms and software suite for reliability assessment of complex technical systems,” Radio Electronics, Computer Science, Control, vol. 4, pp. 163–177, 2020.

A. Abd El-Raheem, “Optimal plans and estimation of constant-stress accelerated life tests for the extension of the exponential distribution under type-I censoring,” Journal of Testing and Evaluation, vol. 47, no. 5, pp. 3781–3821, 2018.

H. M. Aljohani and N. M. Alfar, “Estimations with step-stress partially accelerated life tests for competing risks Burr XII lifetime model under type-II censored data,” Alexandria Engineering Journal, vol. 59, no. 3, pp. 1171–1180, 2020.

S. Dey and M. Nassar, “Classical methods of estimation on constant stress accelerated life tests under exponentiated Lindley distribution,” Journal of Applied Statistics, vol. 47, no. 6, pp. 975–996, 2020.

M. M. M. El-Din, S. E. Abu-Youssef, N. S. A. Ali, and A. M. A. El-Raheem, “Estimation in constant-stress accelerated life tests for extension of the exponential distribution under progressive censoring,” Metron, vol. 74, no. 2, pp. 253–273, 2016.

D. S. Bai and S. W. Chung, “Optimal design of partially accelerated life tests for the exponential distribution under type-I censoring,” IEEE Transactions on Reliability, vol. 41, no. 3, pp. 400–406, 1992.

D. S. Bai, S. W. Chung, and Y. R. Chun, “Optimal design of partially accelerated life tests for the lognormal distribution under type I censoring,” Reliability Engineering & System Safety, vol. 40, no. 1, pp. 85–92, 1993.

R. Zaman, P. Nasiri, and A. Shadrokh, “Statistical inference for the lomax distribution under partially accelerated life tests with progressively type-ii censoring with binomial removal,” Jordan Journal of Mathematics and Statistics (JJMS), vol. 13, no. 3, pp. 439–458, 2020.

A. A. Ismail and A. Al Tamimi, “Optimum constant-stress partially accelerated life test plans using type-I censored data from the inverse Weibull distribution,” Strength of Materials, vol. 49, no. 6, pp. 847–855, 2017.

A. S. Hassan, S. M. Assar, and A. N. Zaky, “Constant-stress partially accelerated life tests for inverted Weibull distribution with multiple censored data,” International Journal of Advanced Statistics and Probability, vol. 3, no. 1, pp. 72–82, 2015.

M. Kamal, S. Zarrin, and A. Islam, “Constant stress partially accelerated life test design for inverted Weibull distribution with type-I censoring,” Algorithms Research, vol. 2, no. 2, pp. 43–49, 2013.

A. M. Abd-Elfattah, A. S. Hassan, and S. G. Nassr, “Estimation in step-stress partially accelerated life tests for the Burr type XII distribution using type I censoring,” Statistical Methodology, vol. 5, no. 6, pp. 502–514, 2008.

I. Alam and A. Ahmed, “Parametric and interval estimation under step-stress partially accelerated life tests using adaptive type-II progressive hybrid censoring,” Annals of Data Science, pp. 1–13, 2020.

Ç. Çetinkaya, “Estimation in step-stress partially accelerated life tests for the power lindley distribution under progressive censoring,” Gazi University Journal of Science, 2021.

T. A. Abushal and A. A. Soliman, “Estimating the Pareto parameters under progressive censoring data for constant-partially accelerated life tests,” Journal of Statistical Computation and Simulation, vol. 85, no. 5, pp. 917–934, 2015.

D. Han and T. Bai, “On the maximum likelihood estimation for progressively censored lifetimes from constant-stress and step-stress accelerated tests,” Electronic Journal of Applied Statistical Analysis, vol. 12, no. 2, pp. 392–404, 2019.

G. Prakash, “Bayes estimation in step-stress PALT on type-I progressive hybrid Rayleigh data,” Afrika Statistika, vol. 14, no. 4, pp. 2165–2178, 2019.

A. Ismail, The Test Design and Parameter Estimation of Pareto Lifetime Distribution under Partially Accelerated Life Tests, Cairo University, Giza, Egypt, 2004.

A. A. Ismail, A. A. Abdel-Ghaly, and E. H. El-Khodary, “Optimum constant-stress life test plans for Pareto distribution under type-I censoring,” Journal of Statistical Computation and Simulation, vol. 81, no. 12, pp. 1835–1845, 2011.

Y.-F. Cheng and F.-K. Wang, “Estimating the Burr XII parameters in constant-stress partially accelerated life tests under multiple censored data,” Communications in Statistics—Simulation and Computation, vol. 41, no. 9, pp. 1711–1727, 2012.

A. A. Ismail, “Estimating the parameters of Weibull distribution and the acceleration factor from hybrid partially accelerated life test,” Applied Mathematical Modelling, vol. 36, no. 7, pp. 2920–2925, 2012.

L. Wang, “Estimation of constant-stress accelerated life test for Weibull distribution with nonconstant shape parameter,” Journal of Computational and Applied Mathematics, vol. 343, pp. 539–555, 2018.

L. Wang and Y. Shi, “Estimation for constant-stress accelerated life test from generalized half-normal distribution,” Journal of Systems Engineering and Electronics, vol. 28, no. 4, pp. 810–816, 2017.

M. Kumar, A. Pathak, and S. Soni, “Bayesian inference for Rayleigh distribution under step-stress partially accelerated test with progressive type-II censoring with binomial removal,” Annals of Data Science, vol. 6, no. 1, pp. 117–152, 2019.

G. Prakash, “Bayes estimation under different censoring patterns on constant-stress PALT,” Austrian Journal of Statistics, vol. 47, no. 4, pp. 60–74, 2018.

S. Roy, “Bayesian accelerated life test plans for series systems with Weibull component lifetimes,” Applied Mathematical Modelling, vol. 62, pp. 383–403, 2018.

Y. Wu, S. Younas, K. Abbas, A. Ali, and S. A. Khan, “Monitoring reliability for three-parameter Frechet distribution using control charts,” IEEE Access, vol. 8, pp. 71245–71253, 2020.

M. Bladt, L. J. R. Esparza, and B. F. Nielsen, “Fisher information and statistical inference for phase-type distributions,” Journal of Applied Probability, vol. 48, pp. 277–293, 2011.