Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring

Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with... Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring div.banner_title_bkg div.trangle { border-color: #A76330 transparent transparent transparent; opacity:0.75; /*new styles start*/ -ms-filter:"progid:DXImageTransform.Microsoft.Alpha(Opacity=75)" ;filter: alpha(opacity=75); /*new styles end*/ } div.banner_title_bkg_if div.trangle { border-color: transparent transparent #A76330 transparent ; opacity:0.75; /*new styles start*/ -ms-filter:"progid:DXImageTransform.Microsoft.Alpha(Opacity=75)" ;filter: alpha(opacity=75); /*new styles end*/ } div.banner_title_bkg div.trangle { width: 215px; } #banner { background-image: url('http://images.hindawi.com/journals/as/as.banner.jpg'); background-position: 50% 0;} Hindawi Publishing Corporation Home Journals About Us Advances in Statistics About this Journal Submit a Manuscript Table of Contents Journal Menu About this Journal · Abstracting and Indexing · Advance Access · Aims and Scope · Article Processing Charges · Articles in Press · Author Guidelines · Bibliographic Information · Contact Information · Editorial Board · Editorial Workflow · Free eTOC Alerts · Publication Ethics · Reviewers Acknowledgment · Submit a Manuscript · Subscription Information · Table of Contents Open Special Issues · Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Full-Text ePUB Linked References How to Cite this Article Advances in Statistics Volume 2015 (2015), Article ID 319051, 13 pages http://dx.doi.org/10.1155/2015/319051 Research Article Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring M. M. Mohie EL-Din , 1 S. E. Abu-Youssef , 1 Nahed S. A. Ali , 2 and A. M. Abd El-Raheem 2 1 Department of Mathematics, Faculty of Science, Al-Azhar University, Cairo, Egypt 2 Department of Mathematics, Faculty of Education, Ain Shams University, Cairo, Egypt Received 22 October 2014; Accepted 20 January 2015 Academic Editor: Shuo-Jye Wu Copyright © 2015 M. M. Mohie EL-Din 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 Based on progressive censoring, step-stress partially accelerated life tests are considered when the lifetime of a product follows power generalized Weibull distribution. The maximum likelihood estimates (MLEs) and Bayes estimates (BEs) are obtained for the distribution parameters and the acceleration factor. In addition, the approximate and bootstrap confidence intervals (CIs) of the estimators are presented. Furthermore, the optimal stress change time for the step-stress partially accelerated life test is determined by minimizing the asymptotic variance of MLEs of the model parameters and the acceleration factor. Simulation results are carried out to study the precision of the MLEs and BEs for the parameters involved. 1. Introduction In reliability analysis, it is not easy to collect lifetimes on highly reliable products with very long lifetimes, because very few or even no failures may occur within a limited testing time under normal conditions. For this reason, accelerated life tests (ALTs) or partially accelerated life tests (PALTs) are one of the most common approaches that are used to obtain enough failure data, in a short period of time. In ALTs all test units are subjected to higher than usual levels of stress, to induce early failures. In PALTs units are tested at both accelerated and use conditions. The information obtained from the test performed in the accelerated or partially accelerated test is used to estimate the failure behavior of the units under normal conditions. The stress loading in ALTs can be applied in different ways. Commonly used methods are constant-stress and step-stress. Nelson [ 1 ] discussed the advantages and disadvantages of each of such methods. In constant-stress ALT, each unit is run at constant high stress until either failure occurs or the test is terminated. In step-stress ALT, the stress on each unit is not constant but is increased step by step at prespecified times or upon the occurrence of a fixed number of failures. When a test involves two levels of stress with the first level as the normal condition and has a fixed time point for changing stress referred to as a step-stress partially ALT (SSPALT). PALTs were studied under step-stress scheme by several authors; for example, see Goel [ 2 ], DeGroot and Goel [ 3 ], Bhattacharyya and Soejoeti [ 4 ], Bai and Chung [ 5 ], Ismail and Aly [ 6 ], and Abdel-Ghani [ 7 ]. In ALTs or PALTs, tests are often stopped before all units fail. The estimate from the censored data is less accurate than the estimate from complete data. However, censored data is more than offset by the reduced test time and expense. The most common censoring scheme is type-II censoring. Consider units are placed on life test, and the experimenter terminates the experiment after a prespecified number of units fail. In this scenario, only the smallest lifetimes are observed. In conventional type-II censoring schemes do not allow to remove units at points other than the terminal point of the experiment. A generalization of type-II censoring is the progressive type-II censoring. It is a method which enables an efficient exploitation of the available resources by continual removal of a prespecified number of surviving test units at each failure time. On other hand, the removal of units before failure may be intentional to save time and cost or when some items have to be removed for use in another experiment. This paper will concentrate on SSPALTs under progressive type-II censoring. It can be described as follows. Consider an experiment in which units are placed on a life testing experiment. At the time of the first failure, units are randomly removed from the remaining surviving units. Similarly, at the time of the second failure, units from the remaining units are randomly removed. The test continues until the th failure occurs at which time, all the remaining, units are removed. If , then , which is the complete sample situation. If , then which corresponds to conventional type-II censoring. A recent account on progressive censoring schemes can be found in the book by Balakrishnan and Aggarwala [ 8 ]. The paper is organized as follows. In Section 2 , a description of the model, test procedure, and its assumptions are presented. In Section 3 , the MLEs of the SSPALT model parameters are derived. The BEs of model parameters using Lindley’s approximation and MCMC method are obtained in Section 4 . In Section 5 , the approximate and bootstrap confidence bounds for the model parameters are constructed. In Section 6 , estimation of optimal stress change time is obtained. Section 7 contains the simulation results. Conclusion is made in Section 8 . 2. Model Description 2.1. Power Generalized Weibull Distribution The power generalized Weibull (PGW) distribution is an extension of Weibull distribution. It was introduced by Bagdonavičius and Nikulin [ 9 ] as a baseline distribution for the accelerated failure time model. It not only contains distributions with unimodal and bathtub hazard shape but also allows for a broader class of monotone hazard rate. The PGW distribution is specified by the probability density function (pdf): the corresponding survival function is and the corresponding hazard rate function is given by Particular cases of the power generalized Weibull distribution are as follows. (1) If , PGW distribution tends to Weibull distribution. (2) If and , PGW distribution tends to exponential distribution. 2.2. Assumptions and Test Procedure The following assumptions are used throughout the paper in the framework of SSPALT. (1) identical and independent items are put on a life test. (2) The lifetime of each unit has PGW distribution. (3) The test is terminated at the time of the th failure, where is prefixed . (4) Each of the items is run under normal use condition. If it does not fail or remove from the test by a prespecified time , it is put under accelerated condition. (5) At the time of the th failure, a random number of the surviving items , are randomly selected and removed from the test. Finally, at the time of the th failure, the remaining surviving items are removed from the test and the test is terminated. (6) Let be the number of failures before time at normal condition, and let be the number of failures after time at stress condition, then, the observed progressive censored data are where and . (7) The tampered random variable (TRV) model holds. It was proposed by DeGroot and Goel [ 3 ]. According to tampered random variable model the lifetime of a unit under SSPALT can be written as where is the lifetime of the units under normal condition, is the stress change time, and is the acceleration factor . (8) From the TRV model in ( 5 ), the pdf of two parameters PGW distribution under SSPALT is given by 3. Maximum Likelihood Estimation In this section, the MLEs of the model parameters are obtained. Let , be the observed values of the lifetime obtained from a progressive censoring scheme under SSPALT, with censored scheme . The maximum likelihood function of the observations takes the following form: where From ( 6 ) in ( 7 ), we get where . The log-likelihood function may then be written as and thus we have the likelihood equations for , , and , respectively, as Now, we have a system of three nonlinear equations in three unknowns , , and . It is clear that a closed form solution is very difficult to obtain. Therefore, an iterative procedure such as Newton Raphson can be used to find a numerical solution of the above nonlinear system. 4. Bayes Estimation In this section, the square error loss (SEL) function is considered to obtain BEs of the model parameters , , and . Unfortunately, in many cases the BEs cannot be expressed in explicit forms. So, approximate BEs are obtained under noninformative prior (NIP) and informative prior (IP) using Lindley’s approximation and Markov chain Monte Carlo (MCMC) method. 4.1. Noninformative Prior Assume that the parameters , , and are independent and the NIP for each parameter is as follows: then, the joint NIP of the parameters is given by The joint posterior density function of the parameters , , and can be written from ( 9 ) and ( 15 ) as Based on SEL function, the Bayes estimator of the function of the parameters is Unfortunately, we cannot compute this integral explicitly. Therefore, we adopt two different procedures to approximate this integral; such procedures are Lindley’s approximation and MCMC method. 4.1.1. Bayes Estimation Using Lindley’s Approximation for NIP In this sub-subsection, the approximate BEs , , and under SEL function using Lindley’s approximation are obtained. According to Lindley in [ 10 ], any ratio of the integral is of the form where is a function of , and only; is log of likelihood function; is log of joint prior of , and . Then, can be evaluated as where are the MLEs of the parameters , , , , , , , , , , , . is the variance-covariance matrix of unknown vector . Form the prior distribution in ( 15 ) and ( 19 ), the values of the BEs of various parameters are 4.1.2. Bayesian Estimation Using MCMC Method for NIP In this sub-subsection, MCMC method is considered to generate samples from the posterior distribution and then compute the BEs of , , and under SSPALT using progressive type-II censoring. From the joint posterior density function in ( 16 ), the conditional posterior distributions of , , and can be written, respectively, as The conditional posterior distributions of , , and in ( 23 ), ( 24 ) and ( 25 ) cannot be reduced analytically to well known distribution, but the plot of them shows that they are similar to normal distribution. So, to generate random samples from this distribution, we use the Metropolis method with normal proposal distribution; see Metropolis et al. [ 11 ]. The following algorithm is proposed to generate , , and from the posterior distribution and then obtain the BEs. Algorithm 1. Step 1 . Start with , , and . Step 2 . Set . Step 3 . Generate from proposal distribution . Step 4 . Calculate the acceptance probability Step 5 . Generate . Step 6 . If , accept the proposal distribution and set . Otherwise, reject the proposal distribution and set . Step 7 . To generate do the Steps 2–6 for not . Step 8 . To generate do the Steps 2–6 for not . Step 9 . Set . Step 10 . Repeat Steps 3–9 times. Step 11 . Obtain the BEs of , , and using MCMC under SEL function as 4.2. Informative Prior Assume that the parameters and are dependent have prior density , and is independent of them with NIP as ; then, the joint IP of the parameters is given by The joint posterior density function of the parameters , , and can be written from ( 9 ) and ( 28 ) as Based on SEL function, the Bayes estimator of the function of the parameters is Unfortunately, this integral cannot be reduced to a closed form. Therefore, we adopt two different procedures to approximate this integral; such procedures are Lindley’s approximation and MCMC method. 4.2.1. Lindley’s Approximation of Bayes Estimation for IP In this sub-subsection, Lindley’s approximation is used to obtain BEs of , , and under SEL function. From ( 19 ) and ( 28 ), the BEs of various parameters are 4.2.2. Bayesian Estimation Using MCMC Method for IP In this sub-subsection, the BEs of , , and under SSPALT using progressive type-II censoring for IP case are obtained. From the joint posterior density function in ( 29 ), the conditional posterior distributions of , , and can be written, respectively, as The conditional posterior distributions of , , and in ( 32 ), ( 33 ) and ( 34 ) cannot be reduced analytically to well known distribution, but the plot of them shows that they are similar to normal distribution. So, the Metropolis method is used to generate random samples from this distribution, with normal proposal distribution. The following algorithm is used to generate , , and from the posterior distribution in the case of IP and then obtain the BEs. Algorithm 2. Step 1 . Start with , , and . Step 2 . Set . Step 3 . Generate from proposal distribution . Step 4 . Calculate the acceptance probability Step 5 . Generate . Step 6 . If , accept the proposal distribution and set . Otherwise, reject the proposal distribution and set . Step 7 . To generate do the Steps 2–6 for not . Step 8 . To generate do the Steps 2–6 for not . Step 9 . Set . Step 10 . Repeat Steps 3–9 times. Step 11 . Obtain the BEs of , , and using MCMC under SEL function as 5. Interval Estimation In this section, the approximate and bootstrap confidence intervals (CIs) of the parameters , , and are derived. 5.1. Approximate Confidence Intervals In this subsection, the approximate confidence intervals of the parameters are derived based on the asymptotic distributions of the MLEs of the elements of the vector of unknown parameters . It is known that the asymptotic distribution of the MLEs of is given by Miller [ 12 ]: where , is the variance-covariance matrix of the unknown parameters . The approximate 100 two sided confidence intervals for , , and are, respectively, given by where is the 100 th percentile of a standard normal distribution. 5.2. Bootstrap Confidence Intervals In this subsection, confidence intervals based on the parametric bootstrap method for the unknown parameters , , and using percentile interval are derived; for more details see Efron and Tibshirani [ 13 ]. The following algorithm is implemented to obtain a bootstrap sample. Algorithm 3. Step 1 . From an original data, , compute the MLEs of the parameters , , and . Step 2 . Use , , and to generate a bootstrap sample with same . Step 3 . As in Step 1 based on compute the bootstrap sample estimates , , and of , , and , respectively. Step 4 . Repeat Steps 1–3 times and arrange each estimate in ascending order to obtain the bootstrap sample , , and . Then, the 100 percentile bootstrap confidence intervals for , , and are, respectively, given by 6. Estimation of Optimal Stress Change Time In this section, the optimal change stress time is found by minimizing the asymptotic variance of MLEs of the model parameters and the acceleration factor. The asymptotic variance of , , and is given by the trace of the inverse of the Fisher information matrix. NMinimize option of Mathematica is used to find the time which minimize the asymptotic variance of MLEs. Table 7 represents the values of for different values of , , and . 7. Simulation Studies In this section, simulation studies are conducted to investigate the performances of the MLEs and BEs in terms of their mean square errors (MSEs) for different choices of , , and ; the results were concluded in Table 5 . Also, the approximate and percentile bootstrap CIs are computed; the results were concluded in Table 6 . The progressive censoring schemes used in the Monte Carlo simulation study are given in Table 1 . The estimation procedure is performed according to the following algorithm. Table 1: The progressive censoring schemes used in the Monte Carlo simulation study. Algorithm 4. Step 1 . Specify the values of , , , and . Step 2 . Specify the values of the parameters , , and . Step 3 . For given values of the prior parameters and generate from and from . Step 4 . Use the model given by ( 6 ) to generate progressively censored data for given , ; the set of data can be considered as where and . Step 5 . Use the progressive censored data to compute the MLEs of the model parameters. The Newton Raphson method is applied for solving the nonlinear system ( 11 ) to obtain the MLEs of the parameters. Step 6 . Compute the BEs of the model parameters relative to SEL function for NIP case and IP case based on Lindley’s approximation. Step 7 . Compute the BEs of the model parameters relative to SEL function for NIP case and IP case based on MCMC algorithm, with , . Step 8 . Compute the approximate confidence bounds with confidence levels for the three parameters of the model. Step 9 . Compute the bootstrap confidence intervals for the model parameters, using Algorithm 3 . Step 10 . Replicate the Steps 3–9, times. Step 11 . Compute the average values of the mean square errors (MSEs) associated with the MLEs and BEs of the parameters. Step 12 . Do Steps 1–11 with different values of , , and . 7.1. Illustrative Example In this subsection, we present an example to illustrate the estimation procedure and the two considered approximate and bootstrap CIs for the parameters , , and . In this example, we simulate a sample of size , based on population parameter values , , and . The stress change time is chosen to be equal to . Under progressive type-II censoring, , , and ; then the failure times of SSPALT were listed in Table 2 . Table 3 presents MSEs of the MLEs and BEs of the model parameters. Table 4 includes approximate and percentile bootstrap CIs of the model parameters and their lengths. Table 2: SSPALT simulation data with true values , , , and . Table 3: MSEs of the estimates of , , and with population parameter values: ( , , and ), , , , and , , . Table 4: 95% approximate and bootstrap CIs and their lengths of the parameters , , and with population parameter values , , , and . Table 5: MSEs of the estimates of , , and with true values , , , and . Table 6: Approximate CIs, bootstrap CIs, and their lengths of the parameters , , and with true values , , , and . Table 7: The optimal change stress time with true values , , and . 8. Conclusion In this paper, we considered the statistical inference procedure for the unknown parameters of the PGW distribution and the acceleration factor , when the data are progressive type-II censored from SSPALTs. We obtained the MLEs and BEs of the considered parameters and studied their performance through their MSEs. We also constructed approximate and bootstrap CIs for the parameters. From the results in Tables 5 – 7 , we observe the following. (1) The MSEs of MLEs and BEs of the considered parameters decrease as the sample size increases. (2) The BEs of the considered parameters obtained from both Lindley’s approximation and MCMC method give more accurate results through the MSEs than MLEs. (3) The BEs of , obtained from MCMC method give more accurate results through the MSEs than the BEs obtained from Lindley’s approximation. (4) For the large samples , the BEs of acceleration factor obtained from MCMC method give more accurate results through the MSEs than the BEs of acceleration factor obtained from Lindley’s approximation. (5) For the small samples , the BEs of acceleration factor obtained from Lindley’s approximation give more accurate results through the MSEs than the BEs of acceleration factor obtained from MCMC method. (6) The bootstrap CIs give more accurate results than the approximate CIs since the lengths of the former are less than the lengths of latter, for different sample sizes, observed failures, and censoring schemes, in most cases. (7) The optimal stress change time decreases as the sample size increases. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgments The authors would like to thank the two anonymous referees for their constructive comments on an early version of this paper. References W. Nelson, Accelerated Testing: Statistical Models, Test Plans and Data Analysis , Wiley, New York, NY, USA, 1990. P. K. Goel, “Some estimation problems in the study of tampered random variables,” Tech. Rep. 50, Department of Statistics, Carnegie-Mellon University, Pittspurgh, Pa, USA, 1971. M. H. DeGroot and P. K. Goel, “Bayesian estimation and optimal designs in partially accelerated life-testing,” Naval Research Logistics Quarterly , vol. 26, no. 2, pp. 223–235, 1979. View at Publisher · View at Google Scholar · View at MathSciNet G. K. Bhattacharyya and Z. Soejoeti, “A tampered failure rate model for step-stress accelerated life test,” Communications in Statistics—Theory and Methods , vol. 18, no. 5, pp. 1627–1643, 1989. View at Publisher · View at Google Scholar · View at MathSciNet 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. View at Publisher · View at Google Scholar · View at Scopus A. A. Ismail and H. M. Aly, “Optimal planning of failure-step stress partially accelerated life tests under type-II censoring,” Journal of Statistical Computation and Simulation , vol. 80, no. 12, pp. 1335–1348, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus M. M. Abdel-Ghani, “The estimation problem of the log-logistic parameters instep partially accelerated life tests using type-I censored data,” The National Review of Social Sciences , vol. 41, no. 2, pp. 1–19, 2004. N. Balakrishnan and R. Aggarwala, Progressive Censoring: Theory, Methods, and Applications , Birkhäuser, Boston, Mass, USA, 2000. View at Publisher · View at Google Scholar · View at MathSciNet V. Bagdonavičius and M. Nikulin, Accelerated Life Models: Modeling and Statistical Analysis , Chapman & Hall/CRC, 2002. D. V. Lindley, “Approximate Bayes methods,” Trabajos de Estatadistica , vol. 31, pp. 223–237, 1980. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, “Equation of state calculations by fast computing machines,” The Journal of Chemical Physics , vol. 21, no. 6, pp. 1087–1092, 1953. View at Publisher · View at Google Scholar · View at Scopus R. G. Miller, Survival Analysis , Wiley, New York, NY, USA, 1981. B. Efron and R. J. Tibshirani, An Introduction to the Bootstrap , vol. 57 of Monographs on Statistics and Applied Probability , Chapman and Hall, London, UK, 1993. View at Publisher · View at Google Scholar · View at MathSciNet (function (i, s, o, g, r, a, m) { i['GoogleAnalyticsObject'] = r; i[r] = i[r] || function () { (i[r].q = i[r].q || []).push(arguments) }, i[r].l = 1 * new Date(); a = s.createElement(o), m = s.getElementsByTagName(o)[0]; a.async = 1; a.src = g; m.parentNode.insertBefore(a, m) })(window, document, 'script', '//www.google-analytics.com/analytics.js', 'ga'); ga('create', 'UA-8578054-2', 'auto'); ga('send', 'pageview'); http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Advances in Statistics Hindawi Publishing Corporation

Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring

Download PDF
14 pages

Loading next page...
 
/lp/hindawi-publishing-corporation/estimation-in-step-stress-accelerated-life-tests-for-power-generalized-hRcCH39sIl

References

References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.

Publisher
Hindawi Publishing Corporation
Copyright
Copyright © 2015 M. M. Mohie EL-Din et al.
ISSN
2356-6892
Publisher site
See Article on Publisher Site

Abstract

Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring div.banner_title_bkg div.trangle { border-color: #A76330 transparent transparent transparent; opacity:0.75; /*new styles start*/ -ms-filter:"progid:DXImageTransform.Microsoft.Alpha(Opacity=75)" ;filter: alpha(opacity=75); /*new styles end*/ } div.banner_title_bkg_if div.trangle { border-color: transparent transparent #A76330 transparent ; opacity:0.75; /*new styles start*/ -ms-filter:"progid:DXImageTransform.Microsoft.Alpha(Opacity=75)" ;filter: alpha(opacity=75); /*new styles end*/ } div.banner_title_bkg div.trangle { width: 215px; } #banner { background-image: url('http://images.hindawi.com/journals/as/as.banner.jpg'); background-position: 50% 0;} Hindawi Publishing Corporation Home Journals About Us Advances in Statistics About this Journal Submit a Manuscript Table of Contents Journal Menu About this Journal · Abstracting and Indexing · Advance Access · Aims and Scope · Article Processing Charges · Articles in Press · Author Guidelines · Bibliographic Information · Contact Information · Editorial Board · Editorial Workflow · Free eTOC Alerts · Publication Ethics · Reviewers Acknowledgment · Submit a Manuscript · Subscription Information · Table of Contents Open Special Issues · Special Issue Guidelines Abstract Full-Text PDF Full-Text HTML Full-Text ePUB Linked References How to Cite this Article Advances in Statistics Volume 2015 (2015), Article ID 319051, 13 pages http://dx.doi.org/10.1155/2015/319051 Research Article Estimation in Step-Stress Accelerated Life Tests for Power Generalized Weibull Distribution with Progressive Censoring M. M. Mohie EL-Din , 1 S. E. Abu-Youssef , 1 Nahed S. A. Ali , 2 and A. M. Abd El-Raheem 2 1 Department of Mathematics, Faculty of Science, Al-Azhar University, Cairo, Egypt 2 Department of Mathematics, Faculty of Education, Ain Shams University, Cairo, Egypt Received 22 October 2014; Accepted 20 January 2015 Academic Editor: Shuo-Jye Wu Copyright © 2015 M. M. Mohie EL-Din 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 Based on progressive censoring, step-stress partially accelerated life tests are considered when the lifetime of a product follows power generalized Weibull distribution. The maximum likelihood estimates (MLEs) and Bayes estimates (BEs) are obtained for the distribution parameters and the acceleration factor. In addition, the approximate and bootstrap confidence intervals (CIs) of the estimators are presented. Furthermore, the optimal stress change time for the step-stress partially accelerated life test is determined by minimizing the asymptotic variance of MLEs of the model parameters and the acceleration factor. Simulation results are carried out to study the precision of the MLEs and BEs for the parameters involved. 1. Introduction In reliability analysis, it is not easy to collect lifetimes on highly reliable products with very long lifetimes, because very few or even no failures may occur within a limited testing time under normal conditions. For this reason, accelerated life tests (ALTs) or partially accelerated life tests (PALTs) are one of the most common approaches that are used to obtain enough failure data, in a short period of time. In ALTs all test units are subjected to higher than usual levels of stress, to induce early failures. In PALTs units are tested at both accelerated and use conditions. The information obtained from the test performed in the accelerated or partially accelerated test is used to estimate the failure behavior of the units under normal conditions. The stress loading in ALTs can be applied in different ways. Commonly used methods are constant-stress and step-stress. Nelson [ 1 ] discussed the advantages and disadvantages of each of such methods. In constant-stress ALT, each unit is run at constant high stress until either failure occurs or the test is terminated. In step-stress ALT, the stress on each unit is not constant but is increased step by step at prespecified times or upon the occurrence of a fixed number of failures. When a test involves two levels of stress with the first level as the normal condition and has a fixed time point for changing stress referred to as a step-stress partially ALT (SSPALT). PALTs were studied under step-stress scheme by several authors; for example, see Goel [ 2 ], DeGroot and Goel [ 3 ], Bhattacharyya and Soejoeti [ 4 ], Bai and Chung [ 5 ], Ismail and Aly [ 6 ], and Abdel-Ghani [ 7 ]. In ALTs or PALTs, tests are often stopped before all units fail. The estimate from the censored data is less accurate than the estimate from complete data. However, censored data is more than offset by the reduced test time and expense. The most common censoring scheme is type-II censoring. Consider units are placed on life test, and the experimenter terminates the experiment after a prespecified number of units fail. In this scenario, only the smallest lifetimes are observed. In conventional type-II censoring schemes do not allow to remove units at points other than the terminal point of the experiment. A generalization of type-II censoring is the progressive type-II censoring. It is a method which enables an efficient exploitation of the available resources by continual removal of a prespecified number of surviving test units at each failure time. On other hand, the removal of units before failure may be intentional to save time and cost or when some items have to be removed for use in another experiment. This paper will concentrate on SSPALTs under progressive type-II censoring. It can be described as follows. Consider an experiment in which units are placed on a life testing experiment. At the time of the first failure, units are randomly removed from the remaining surviving units. Similarly, at the time of the second failure, units from the remaining units are randomly removed. The test continues until the th failure occurs at which time, all the remaining, units are removed. If , then , which is the complete sample situation. If , then which corresponds to conventional type-II censoring. A recent account on progressive censoring schemes can be found in the book by Balakrishnan and Aggarwala [ 8 ]. The paper is organized as follows. In Section 2 , a description of the model, test procedure, and its assumptions are presented. In Section 3 , the MLEs of the SSPALT model parameters are derived. The BEs of model parameters using Lindley’s approximation and MCMC method are obtained in Section 4 . In Section 5 , the approximate and bootstrap confidence bounds for the model parameters are constructed. In Section 6 , estimation of optimal stress change time is obtained. Section 7 contains the simulation results. Conclusion is made in Section 8 . 2. Model Description 2.1. Power Generalized Weibull Distribution The power generalized Weibull (PGW) distribution is an extension of Weibull distribution. It was introduced by Bagdonavičius and Nikulin [ 9 ] as a baseline distribution for the accelerated failure time model. It not only contains distributions with unimodal and bathtub hazard shape but also allows for a broader class of monotone hazard rate. The PGW distribution is specified by the probability density function (pdf): the corresponding survival function is and the corresponding hazard rate function is given by Particular cases of the power generalized Weibull distribution are as follows. (1) If , PGW distribution tends to Weibull distribution. (2) If and , PGW distribution tends to exponential distribution. 2.2. Assumptions and Test Procedure The following assumptions are used throughout the paper in the framework of SSPALT. (1) identical and independent items are put on a life test. (2) The lifetime of each unit has PGW distribution. (3) The test is terminated at the time of the th failure, where is prefixed . (4) Each of the items is run under normal use condition. If it does not fail or remove from the test by a prespecified time , it is put under accelerated condition. (5) At the time of the th failure, a random number of the surviving items , are randomly selected and removed from the test. Finally, at the time of the th failure, the remaining surviving items are removed from the test and the test is terminated. (6) Let be the number of failures before time at normal condition, and let be the number of failures after time at stress condition, then, the observed progressive censored data are where and . (7) The tampered random variable (TRV) model holds. It was proposed by DeGroot and Goel [ 3 ]. According to tampered random variable model the lifetime of a unit under SSPALT can be written as where is the lifetime of the units under normal condition, is the stress change time, and is the acceleration factor . (8) From the TRV model in ( 5 ), the pdf of two parameters PGW distribution under SSPALT is given by 3. Maximum Likelihood Estimation In this section, the MLEs of the model parameters are obtained. Let , be the observed values of the lifetime obtained from a progressive censoring scheme under SSPALT, with censored scheme . The maximum likelihood function of the observations takes the following form: where From ( 6 ) in ( 7 ), we get where . The log-likelihood function may then be written as and thus we have the likelihood equations for , , and , respectively, as Now, we have a system of three nonlinear equations in three unknowns , , and . It is clear that a closed form solution is very difficult to obtain. Therefore, an iterative procedure such as Newton Raphson can be used to find a numerical solution of the above nonlinear system. 4. Bayes Estimation In this section, the square error loss (SEL) function is considered to obtain BEs of the model parameters , , and . Unfortunately, in many cases the BEs cannot be expressed in explicit forms. So, approximate BEs are obtained under noninformative prior (NIP) and informative prior (IP) using Lindley’s approximation and Markov chain Monte Carlo (MCMC) method. 4.1. Noninformative Prior Assume that the parameters , , and are independent and the NIP for each parameter is as follows: then, the joint NIP of the parameters is given by The joint posterior density function of the parameters , , and can be written from ( 9 ) and ( 15 ) as Based on SEL function, the Bayes estimator of the function of the parameters is Unfortunately, we cannot compute this integral explicitly. Therefore, we adopt two different procedures to approximate this integral; such procedures are Lindley’s approximation and MCMC method. 4.1.1. Bayes Estimation Using Lindley’s Approximation for NIP In this sub-subsection, the approximate BEs , , and under SEL function using Lindley’s approximation are obtained. According to Lindley in [ 10 ], any ratio of the integral is of the form where is a function of , and only; is log of likelihood function; is log of joint prior of , and . Then, can be evaluated as where are the MLEs of the parameters , , , , , , , , , , , . is the variance-covariance matrix of unknown vector . Form the prior distribution in ( 15 ) and ( 19 ), the values of the BEs of various parameters are 4.1.2. Bayesian Estimation Using MCMC Method for NIP In this sub-subsection, MCMC method is considered to generate samples from the posterior distribution and then compute the BEs of , , and under SSPALT using progressive type-II censoring. From the joint posterior density function in ( 16 ), the conditional posterior distributions of , , and can be written, respectively, as The conditional posterior distributions of , , and in ( 23 ), ( 24 ) and ( 25 ) cannot be reduced analytically to well known distribution, but the plot of them shows that they are similar to normal distribution. So, to generate random samples from this distribution, we use the Metropolis method with normal proposal distribution; see Metropolis et al. [ 11 ]. The following algorithm is proposed to generate , , and from the posterior distribution and then obtain the BEs. Algorithm 1. Step 1 . Start with , , and . Step 2 . Set . Step 3 . Generate from proposal distribution . Step 4 . Calculate the acceptance probability Step 5 . Generate . Step 6 . If , accept the proposal distribution and set . Otherwise, reject the proposal distribution and set . Step 7 . To generate do the Steps 2–6 for not . Step 8 . To generate do the Steps 2–6 for not . Step 9 . Set . Step 10 . Repeat Steps 3–9 times. Step 11 . Obtain the BEs of , , and using MCMC under SEL function as 4.2. Informative Prior Assume that the parameters and are dependent have prior density , and is independent of them with NIP as ; then, the joint IP of the parameters is given by The joint posterior density function of the parameters , , and can be written from ( 9 ) and ( 28 ) as Based on SEL function, the Bayes estimator of the function of the parameters is Unfortunately, this integral cannot be reduced to a closed form. Therefore, we adopt two different procedures to approximate this integral; such procedures are Lindley’s approximation and MCMC method. 4.2.1. Lindley’s Approximation of Bayes Estimation for IP In this sub-subsection, Lindley’s approximation is used to obtain BEs of , , and under SEL function. From ( 19 ) and ( 28 ), the BEs of various parameters are 4.2.2. Bayesian Estimation Using MCMC Method for IP In this sub-subsection, the BEs of , , and under SSPALT using progressive type-II censoring for IP case are obtained. From the joint posterior density function in ( 29 ), the conditional posterior distributions of , , and can be written, respectively, as The conditional posterior distributions of , , and in ( 32 ), ( 33 ) and ( 34 ) cannot be reduced analytically to well known distribution, but the plot of them shows that they are similar to normal distribution. So, the Metropolis method is used to generate random samples from this distribution, with normal proposal distribution. The following algorithm is used to generate , , and from the posterior distribution in the case of IP and then obtain the BEs. Algorithm 2. Step 1 . Start with , , and . Step 2 . Set . Step 3 . Generate from proposal distribution . Step 4 . Calculate the acceptance probability Step 5 . Generate . Step 6 . If , accept the proposal distribution and set . Otherwise, reject the proposal distribution and set . Step 7 . To generate do the Steps 2–6 for not . Step 8 . To generate do the Steps 2–6 for not . Step 9 . Set . Step 10 . Repeat Steps 3–9 times. Step 11 . Obtain the BEs of , , and using MCMC under SEL function as 5. Interval Estimation In this section, the approximate and bootstrap confidence intervals (CIs) of the parameters , , and are derived. 5.1. Approximate Confidence Intervals In this subsection, the approximate confidence intervals of the parameters are derived based on the asymptotic distributions of the MLEs of the elements of the vector of unknown parameters . It is known that the asymptotic distribution of the MLEs of is given by Miller [ 12 ]: where , is the variance-covariance matrix of the unknown parameters . The approximate 100 two sided confidence intervals for , , and are, respectively, given by where is the 100 th percentile of a standard normal distribution. 5.2. Bootstrap Confidence Intervals In this subsection, confidence intervals based on the parametric bootstrap method for the unknown parameters , , and using percentile interval are derived; for more details see Efron and Tibshirani [ 13 ]. The following algorithm is implemented to obtain a bootstrap sample. Algorithm 3. Step 1 . From an original data, , compute the MLEs of the parameters , , and . Step 2 . Use , , and to generate a bootstrap sample with same . Step 3 . As in Step 1 based on compute the bootstrap sample estimates , , and of , , and , respectively. Step 4 . Repeat Steps 1–3 times and arrange each estimate in ascending order to obtain the bootstrap sample , , and . Then, the 100 percentile bootstrap confidence intervals for , , and are, respectively, given by 6. Estimation of Optimal Stress Change Time In this section, the optimal change stress time is found by minimizing the asymptotic variance of MLEs of the model parameters and the acceleration factor. The asymptotic variance of , , and is given by the trace of the inverse of the Fisher information matrix. NMinimize option of Mathematica is used to find the time which minimize the asymptotic variance of MLEs. Table 7 represents the values of for different values of , , and . 7. Simulation Studies In this section, simulation studies are conducted to investigate the performances of the MLEs and BEs in terms of their mean square errors (MSEs) for different choices of , , and ; the results were concluded in Table 5 . Also, the approximate and percentile bootstrap CIs are computed; the results were concluded in Table 6 . The progressive censoring schemes used in the Monte Carlo simulation study are given in Table 1 . The estimation procedure is performed according to the following algorithm. Table 1: The progressive censoring schemes used in the Monte Carlo simulation study. Algorithm 4. Step 1 . Specify the values of , , , and . Step 2 . Specify the values of the parameters , , and . Step 3 . For given values of the prior parameters and generate from and from . Step 4 . Use the model given by ( 6 ) to generate progressively censored data for given , ; the set of data can be considered as where and . Step 5 . Use the progressive censored data to compute the MLEs of the model parameters. The Newton Raphson method is applied for solving the nonlinear system ( 11 ) to obtain the MLEs of the parameters. Step 6 . Compute the BEs of the model parameters relative to SEL function for NIP case and IP case based on Lindley’s approximation. Step 7 . Compute the BEs of the model parameters relative to SEL function for NIP case and IP case based on MCMC algorithm, with , . Step 8 . Compute the approximate confidence bounds with confidence levels for the three parameters of the model. Step 9 . Compute the bootstrap confidence intervals for the model parameters, using Algorithm 3 . Step 10 . Replicate the Steps 3–9, times. Step 11 . Compute the average values of the mean square errors (MSEs) associated with the MLEs and BEs of the parameters. Step 12 . Do Steps 1–11 with different values of , , and . 7.1. Illustrative Example In this subsection, we present an example to illustrate the estimation procedure and the two considered approximate and bootstrap CIs for the parameters , , and . In this example, we simulate a sample of size , based on population parameter values , , and . The stress change time is chosen to be equal to . Under progressive type-II censoring, , , and ; then the failure times of SSPALT were listed in Table 2 . Table 3 presents MSEs of the MLEs and BEs of the model parameters. Table 4 includes approximate and percentile bootstrap CIs of the model parameters and their lengths. Table 2: SSPALT simulation data with true values , , , and . Table 3: MSEs of the estimates of , , and with population parameter values: ( , , and ), , , , and , , . Table 4: 95% approximate and bootstrap CIs and their lengths of the parameters , , and with population parameter values , , , and . Table 5: MSEs of the estimates of , , and with true values , , , and . Table 6: Approximate CIs, bootstrap CIs, and their lengths of the parameters , , and with true values , , , and . Table 7: The optimal change stress time with true values , , and . 8. Conclusion In this paper, we considered the statistical inference procedure for the unknown parameters of the PGW distribution and the acceleration factor , when the data are progressive type-II censored from SSPALTs. We obtained the MLEs and BEs of the considered parameters and studied their performance through their MSEs. We also constructed approximate and bootstrap CIs for the parameters. From the results in Tables 5 – 7 , we observe the following. (1) The MSEs of MLEs and BEs of the considered parameters decrease as the sample size increases. (2) The BEs of the considered parameters obtained from both Lindley’s approximation and MCMC method give more accurate results through the MSEs than MLEs. (3) The BEs of , obtained from MCMC method give more accurate results through the MSEs than the BEs obtained from Lindley’s approximation. (4) For the large samples , the BEs of acceleration factor obtained from MCMC method give more accurate results through the MSEs than the BEs of acceleration factor obtained from Lindley’s approximation. (5) For the small samples , the BEs of acceleration factor obtained from Lindley’s approximation give more accurate results through the MSEs than the BEs of acceleration factor obtained from MCMC method. (6) The bootstrap CIs give more accurate results than the approximate CIs since the lengths of the former are less than the lengths of latter, for different sample sizes, observed failures, and censoring schemes, in most cases. (7) The optimal stress change time decreases as the sample size increases. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgments The authors would like to thank the two anonymous referees for their constructive comments on an early version of this paper. References W. Nelson, Accelerated Testing: Statistical Models, Test Plans and Data Analysis , Wiley, New York, NY, USA, 1990. P. K. Goel, “Some estimation problems in the study of tampered random variables,” Tech. Rep. 50, Department of Statistics, Carnegie-Mellon University, Pittspurgh, Pa, USA, 1971. M. H. DeGroot and P. K. Goel, “Bayesian estimation and optimal designs in partially accelerated life-testing,” Naval Research Logistics Quarterly , vol. 26, no. 2, pp. 223–235, 1979. View at Publisher · View at Google Scholar · View at MathSciNet G. K. Bhattacharyya and Z. Soejoeti, “A tampered failure rate model for step-stress accelerated life test,” Communications in Statistics—Theory and Methods , vol. 18, no. 5, pp. 1627–1643, 1989. View at Publisher · View at Google Scholar · View at MathSciNet 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. View at Publisher · View at Google Scholar · View at Scopus A. A. Ismail and H. M. Aly, “Optimal planning of failure-step stress partially accelerated life tests under type-II censoring,” Journal of Statistical Computation and Simulation , vol. 80, no. 12, pp. 1335–1348, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus M. M. Abdel-Ghani, “The estimation problem of the log-logistic parameters instep partially accelerated life tests using type-I censored data,” The National Review of Social Sciences , vol. 41, no. 2, pp. 1–19, 2004. N. Balakrishnan and R. Aggarwala, Progressive Censoring: Theory, Methods, and Applications , Birkhäuser, Boston, Mass, USA, 2000. View at Publisher · View at Google Scholar · View at MathSciNet V. Bagdonavičius and M. Nikulin, Accelerated Life Models: Modeling and Statistical Analysis , Chapman & Hall/CRC, 2002. D. V. Lindley, “Approximate Bayes methods,” Trabajos de Estatadistica , vol. 31, pp. 223–237, 1980. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, “Equation of state calculations by fast computing machines,” The Journal of Chemical Physics , vol. 21, no. 6, pp. 1087–1092, 1953. View at Publisher · View at Google Scholar · View at Scopus R. G. Miller, Survival Analysis , Wiley, New York, NY, USA, 1981. B. Efron and R. J. Tibshirani, An Introduction to the Bootstrap , vol. 57 of Monographs on Statistics and Applied Probability , Chapman and Hall, London, UK, 1993. View at Publisher · View at Google Scholar · View at MathSciNet (function (i, s, o, g, r, a, m) { i['GoogleAnalyticsObject'] = r; i[r] = i[r] || function () { (i[r].q = i[r].q || []).push(arguments) }, i[r].l = 1 * new Date(); a = s.createElement(o), m = s.getElementsByTagName(o)[0]; a.async = 1; a.src = g; m.parentNode.insertBefore(a, m) })(window, document, 'script', '//www.google-analytics.com/analytics.js', 'ga'); ga('create', 'UA-8578054-2', 'auto'); ga('send', 'pageview');

Journal

Advances in Statistics – Hindawi Publishing Corporation

Published: Feb 17, 2015

There are no references for this article.