Mustafa Kamal

Department of Basic Sciences,College of Science and Theoretical Studies,Saudi Electronic University,Dammam,32256,Saudi Arabia

ABSTRACT In general,simple subsystems like series or parallel are integrated to produce a complex hybrid system.The reliability of a system is determined by the reliability of its constituent components.It is often extremely diきcult or impossible to get specific information about the component that caused the system to fail.Unknown failure causes are instances in which the actual cause of system failure is unknown.On the other side,thanks to current advanced technology based on computers,automation,and simulation,products have become incredibly dependable and trustworthy,and as a result,obtaining failure data for testing such exceptionally reliable items have become a very costly and time-consuming procedure.Therefore,because of its capacity to produce rapid and adequate failure data in a short period of time,accelerated life testing(ALT)is the most utilized approach in the field of product reliability and life testing.Based on progressively hybrid censored(PrHC)data from a three-component parallel series hybrid system that failed to owe to unknown causes,this paper investigates a challenging problem of parameter estimation and reliability assessment under a step stress partially accelerated life-test(SSPALT).Failures of components are considered to follow a power linear hazard rate(PLHR),which can be used when the failure rate displays linear,decreasing,increasing or bathtub failure patterns.The Tempered random variable(TRV)model is considered to reflect the eあect of the high stress level used to induce early failure data.The maximum likelihood estimation(MLE)approach is used to estimate the parameters of the PLHR distribution and the acceleration factor.A variance covariance matrix(VCM)is then obtained to construct the approximate confidence intervals(ACIs).In addition,studentized bootstrap confidence intervals(ST-B CIs)are also constructed and compared with ACIs in terms of their respective interval lengths(ILs).Moreover,a simulation study is conducted to demonstrate the performance of the estimation procedures and the methodology discussed in this paper.Finally,real failure data from the air conditioning systems of an airplane is used to illustrate further the performance of the suggested estimation technique.

KEYWORDS Step stress partially accelerated life test;progressive hybrid censoring;data masking;power linear hazard rate distribution;hybrid system

1 Introduction

Computers,cellphones,and many other systems of competitive innovation and modernization are getting increasingly complex,with numerous subsystems and sub-assemblies in each.Furthermore,these subsystems and subassemblies are made up of several components,making the life testing procedure for such systems more complicated.The study of reliability and life testing of systems,subsystems,or components is entirely depend on lifetime data,which is a combination of two important pieces of information.The first is to determine the product’s failure time,and the second is to determine the reasons for its failure.The failure times of the system can easily be recorded,but the cause of the failure is not always recognized due to a variety of factors,including a lack of adequate diagnosis,time and expense constraints for comprehensive failure analysis,and numerous component failures with destructive consequences.As a result,in the literature of reliability and life testing analysis,such data in which the true cause of system failure is unknown and only a minimum random subset of the reasons that are responsible for system failure can be recognized is referred to as masked data.See[1,2]for more details.

The overall quality of today’s products due to the existing advanced technology based on computers,automation and simulation has been improved drastically,which makes them extremely reliable and trustworthy.Consequently,gathering failure data for testing of such extremely reliable products using ordinary reliability tests(ORTs)has become a very costly and time-consuming process,making the use of ORTs impractical.Hence,ALTs are a more advanced approach for obtaining fast failure data by testing items under higher stress than normal,and then,a life stress relationship is used to get the product’s life characteristics under normal usage settings.According to[3],stress loading in ALT may be performed in a variety of ways,although the most commonly used stress loadings are constant,step,and progressive stress loadings.Many scholars so far have looked at the ALT models,including[4–10].Assuming a lognormal lifetime distribution,Li et al.[5]proposed two types of Bayesian accelerated acceptance sampling plans for illustrating product reliability based on the product’s operating characteristic curve under Type-I censoring.The first plan addresses both producer and customer risks at the same time,whereas the second exclusively considers consumer risk.Rahman et al.[9]used MLE methods for estimating the parameters of Burr-X life distribution parameters assuming that failure under arithmetically increasing stress levels of CSALT forms a geometric process.Ma et al.[10]proposed an optimum hybrid accelerated test plan under many experimental design restrictions by combining ALT with accelerated degradation testing and modeling the degradation process with an inverse Gaussian process.

However,there are situations when these sorts of relationships are not possible.Because of the prevalence of this issue in ALTs,PALTs are considered to be preferable option and are usually implemented in two ways.The first is constant stress PALT,while the second is SSPALT.In SSPALT,a sample of components or systems is first tested under normal usage circumstances for a certain amount of time,and then systems or components that have not failed are allocated to testing at accelerated conditions until all items fail or a predetermined censorship scheme is met.SSPALT analysis has been considered by many authors since it was first suggested by[11,12]as a TRV model.Bhattacharyya et al.[13]developed a tampered failure rate model using the TRV model.Bai et al.[14–16]are others who considered SSPALT using the same concept of the TRV model for different distributions and censoring schemes.Assuming the TRV model,Zhang et al.[17]discussed the MLEs of the unknown parameters for the extended Weibull distribution.Ismail et al.[18,19]investigated the MLEs of the parameters of the Weibull and Burr Type-XII distributions,respectively,and compared the results based on two different PrHC schemes.Mahmoud et al.[20–22]deals with SSPALT based on an adaptive type PrHC scheme to obtain estimates of parameters using MLEs and Bayesian estimates(BEs)of generalized Pareto,two-parameter exponentiated Weibull and Lindley distributions,respectively.

A considerable number of studies have been carried out on parameter estimation using mask data based on ORTs since it was first introduced by[1,2].Considering different failure distributions,Guess et al.[23–33]utilizes the MLE technique for estimating model parameters for a single component or a series or parallel systems of two or three components,whereas Reiser et al.[34–43]considered BE technique based on different priors.As it was discussed earlier,many real-life systems or machines these days are made of hybrid structures which are a combination of series and parallel subsystems.More complex systems can have many hybrid subsystems which are generally connected in parallel-series,series-parallel,series-parallel-series,and parallel-seriesparallel configurations and so on.Peng et al.[44]developed a Bayesian approach for system reliability evaluation and prediction in which pass-fail data,lifespan data,and deterioration data are integrated coherently at multiple system levels.Yang et al.[45]developed an Adaptive Bayesian Melding reliability evaluation technique for analyzing and assessing the reliability of a hierarchical system with imperfect prior knowledge.They also expanded the concept to a broader multi-level hierarchical structure by employing a more effective method of pooling inconsistent priors.Yang et al.[46]proposed a Bayesian reliability approach for complex systems with dependent life metrics and developed a likelihood decomposition method to convert the overall likelihood into a product of explicit and implicit evidence-based likelihood functions.As of now,only a limited number of research papers considering hybrid systems based on ORTs have been published.Considering three component parallel-series and series-parallel systems,Wang et al.[47]obtained the MLEs of the parameters based on masked data assuming constant and linear hazard rate of independent components.Sha et al.[48]considered a hybrid system of three dependent components and obtained the MLEs based on mask data assuming a bivariate exponential distribution.Cai et al.[49]considered the same system as in[48]and derived the MLEs based copula function under mask cause of failure.Recently,Rodrigues et al.[50,51]used more complex structures based on four or five components and obtained MLEs and BEs under incomplete data.

Only a few studies on ALTs so far that focused on hybrid systems and masked data have been published in available research.Reference[52]under SSALT,describes the procedure of obtaining MLEs of the parameters of the Weibull distribution for a series system based on an expectation minimization algorithm assuming symmetric masking.Considering masked interval data in SSALT,Fan et al.[53]obtained estimates of the parameters of the exponential distribution.Xu et al.[54,55]described the general Bayesian analysis of the series system masked failure data for the log-location-scale and Weibull distributions respectively under SSALT.Assuming the same exponential hazard rate for a four components hybrid system,Shi et al.[56,57]obtained the MLEs of unknown parameters of the model based on the masked data under SSPALT and CSPALT respectively.Shi et al.[58]considered two different hybrid systems of three components and then obtained MLEs of the modified Weibull distribution.

To the best of our knowledge so far,no article has been published that considers SSALT for PLHR distribution for a hybrid system under masked data.Much of the sources listed above considered hazard rates that are monotonic in nature,but there are cases where the hazard rate is not monotonic.The primary goal of this work is to describe the SSPALT using a more flexible PLHR that can be used when failure rates indicate non-monotonic characteristics for hybrid systems.To illustrate the considered estimation procedure under SSPALT using PrHC masked data,a three-component hybrid system is considered.The rest of the paper is organized as follows,Section 2 addresses the formulation of the SSPALT model for PrHC masked data from a hybrid system and some useful assumptions.The MLEs and corresponding ACIs and ST-B CIs of the parameters and acceleration factor are discussed in Section 3.In Section 4,we conducted a simulation study to demonstrate the performance of the estimation procedures and the methodology discussed in this paper.Section 5 addresses the suggested approach’s real-life data applicability.Finally,we conclude our study in Section 6.

2 Design and Assumptions of the Model

2.1 PLHR Distribution

In reliability and life testing experiments/studies,hazard rate is the most important function since it plays a very important role in characterizing the aging process of the systems and hence in classifying failure time distributions,more details can be seen in[59].Commonly used hazard rates are constant,linear and power hazard rates.However,distributions derived from these hazard rates are very useful and popular in reliability and life testing theory when the failure rate depicts monotonic properties,but they cannot be used to fit non-monotonic failure rates.Therefore,in this paper,a more flexible PLHR distribution introduced by[60]that is derived from the combination of linear and power hazard rates is considered.Different shapes of hazard rate and density function of PLHR distribution with different values of parameters are given in Figs.1 and 2,respectively.The Probability density function(PDF),cumulative distribution function(CDF),reliability function(RF)and the hazard rate(HR)of PLHR distribution are given respectively by the following equations:

Figure 1:(Continued)

Figure 1:Shapes of the HR function of PLHR distribution with different values of the parameters

Figure 2:Shapes of the PDF of PLHR distribution using different values of the parameters

2.2 Basic Assumptions

In this paper,the following assumptions are made in order to describe the SSPALT model for a hybrid structure under PrHC masked data:

1.Xξ,ξ=1,2 represents the two stress levels,i.e.,X1and X2are normal and accelerated stress level,respectively,used in SSPALT.

2.The system under consideration is a series parallel system containing three independent componentsj=1,2,3 which is described by Fig.3 as follows:

Figure 3:A series parallel hybrid system of three components

3.The lifetimeTξi,ξ=1,2;i=1,2,...,nof systemitested under SSPALT are i.i.d.at both normal and accelerated test conditions.

4.The data masking mechanism is statistically independent of the various stress conditions used in the experiment and the actual cause of failure of the system.

5.The lifetime ofjth component at normal stress X1,in the considered hybrid system follows PLHR distribution and the corresponding PDF,CDF,SF and HR are given respectively by Eqs.(1)–(4).

6.Let the total lifetime T of a component in hybrid system under SSPALT is explained by TRV model and can be written as follows:

whereXrepresents lifetime of the component at X1,τis the time when items are switched from X1to X2and is known as stress change time and?>1 is used to reflect the effect of stress change and is known as acceleration factor.Now lifetime of the component at X2by using Eqs.(1)–(5)can be obtained as

Definition 2.1:Masking Probability

Suppose there arensystem that are tested in the experiment andtξirepresents observed value ofTξi,whereξ=1,2 represents the stress levels.Suppose also thattξij,i=1,2,3,...,n;j=1,2,3 are realizations of the life timeTξijof thejthcomponent of the systemi.Letωξi?Ωξibe a masked event corresponding to thejthcomponent of the systemithat is one of the possible causes for system failure.Ifωiinclude only single component,e.g.,ωi={1},then it can be said that the failure cause of the system is exact,otherwise cause is called masked[1,47,56].LetCξijbe thejthcomponent that is the exact cause of failure of systemiwith masked eventωξi?Ωξi,then the masking probability can be defined as follows:

Now,using the Assumption 4,above expression for masking probability can be written as follows:

2.3 Formulation of the SSPALT with PrHC Mask Data

PrHC scheme was first proposed and applied by[61]in traditional life tests.In SSPALT,suppose experiment is first started by randomly choosingnsimilar systems of three components described in the Fig.3 to be assigned to test at the normal stress level X1with some predefined progressive values of random removals patternse1,e2...,em,stress change timeτand a test termination timet0.Now,test will progress by removinge1systems randomly at the timet1of first failure.Similarly,e2systems will be removed randomly at the timet2of second failure and so on.If,n1is the number of total observed failures occurred at normal use condition X1anden1is total systems removed at X1before stress change timeτ,then total remaining surviving systems will be(n?nu?en1).At this pointτof the experiment,all the remaining(n?n1?en1)survival systems are removed from the test and are assigned to test at accelerated test condition X2(>X1)with hybrid test termination timet0?=min(tm,m,n,t0)following the same procedure as it was on stress X1.Ifmthfailuretm,m,nis observed before the pre-set time constraintt0,then we will stop the test by removing all the remainingem=n?m?(e1+e2+···+em?1)test systems from the experiment.Ifmthfailuretm,m,nis not observed beforet0,then test will be terminated by removing all the remaininged=n?d?(e1+e2+···+ed?1)test systems,wheredis the number system failures beforet0.Letn2be the total number of observed failures at accelerated test condition X1during time interval[τ,t0],for more details see[62,63].In many complicated systems these days,the cause that exactly responsible for system failure is not known.Let Ωξi={{1},{2},{3},{1,2},{1,3},{2,3},{1,2,3}}be the set of all the possible events that can be a reason for the system failure andωξibe a subset(observed value)of the Ωξi.If observed subsetωξicontains exactly one element,then cause of the system failure is known.In contrast,if observed subsetωξicontains more than one element,then the exact reason of system failure is not known(or Masked).In this scenario,following two cases of failure data are observed in general:

Now,using Eqs.(10),(12)and the PrHC masked data given in Eq.(16),we can write the likelihood function for the hybrid system under SSPALT as follows:

where,ti=tξ,i,m,n,r=m,r?=0 for first data set andr=d=n1+n2,r?=n?d?(e1+e2+···+ed?1)for second data set given in Eq.(12).

3 Estimation of Model Parameters

3.1 Point Estimates

3.2 Interval Estimates

3.2.1 ACIs

3.2.2 Bootstrap CIs

This subsection deals with the bootstrap re-sampling approach for constructing parameters CIs in which the original data is assumed to be a population and then several samples are generated using the original data to create the CIs,for more details see[64–66].We use the following algorithm to construct ST-B CIs[56]:

4 Simulation Study

Here in this section of the paper,we are going to perform a simulation study to investigate and compare the performance of estimates for the hybrid system under SSPALT for PrHC masked data by utilizing the Monte Carlo simulation technique.The performance of MLEs assessed with their respective mean square errors(MSEs)and relative absolute biases(RABs).95% ACIs are also constructed and their performance is investigated in terms of their respective ILs.First,we generated the hybrid progressive censored data from the considered distribution under SSPALT following some of the steps given in[58,64,67]and the complete steps of the algorithm are given as:

9.Replicate Steps 2–9,10,000 times,and obtain average MLEs,MSEs and RABs of the parameters.

10.Compute the 95% ACIs for the parameters of PLHR distribution and acceleration factor.

11.For different values of[τ,n,m,t0,(e1,e2,...,em)],?,γ,κand test schemes,the above procedure from Steps 1 to 10 should be repeated.

In our study,we used four different censoring schemes generated from six different combinations ofnandmwith different fixed values ofτandt0.The censoring schemes,initially fixed values of parameters and the values ofτ,n,m,t0that we used in this paper are given in Table 1.Averages of MLEs based on 10,000 replications described in simulation study and their respective MSEs and RABs are reported in Tables 2 and 3.95% ACIs and ST-B CIs along with their expected interval length are reported in Tables 4 and 5.

Table 1:PrHC schemes,combinations of(n,m) and(τ,t0) with fixed values of parameters ?,γ,κ

Table 2:The MLEs of parameters of with their respective MSEs and RABs with Schemes 1 and 2

Table 2(continued).(n,m)?γκ MLEMSERABMLEMSERABMLEMSERAB CS 1 with(?=1.6,γ =1.1,κ=1.2) and(τ =1.0,t0=1.25)(60,35)1.6260.5290.1721.0760.2490.0581.1260.1780.028(60,45)1.9950.0080.0991.0450.770.4111.4950.6990.327(70,35)1.7310.3810.0841.1810.3970.1331.2310.3260.086(70,45)1.6270.5280.1711.0770.250.0581.1270.1790.029(80,35)1.6260.5290.1721.0760.2490.0581.1260.1780.028(80,45)1.650.4950.1491.10.2820.0731.150.2120.039 CS 1 with(?=1.6,γ =1.1,κ=1.2) and(τ =1.25,t0=1.5)(60,35)1.6270.5280.1711.0770.250.0581.1270.1790.029(60,45)1.730.3820.0841.180.3960.1331.230.3260.086(70,35)1.6260.5290.1721.0760.2490.0581.1260.1790.025(70,45)1.6270.5270.1711.0770.2510.0581.1270.180.029(80,35)1.6240.5320.1741.0740.2460.0561.1240.1760.027(80,45)1.6470.4990.1511.0970.2790.0711.1470.2080.038 CS 2 with(?=1.6,γ =1.1,κ=1.2) and(τ =0.75,t0=1.0)(60,35)1.6270.5270.1711.0770.2510.0581.1270.180.029(60,45)1.6270.5280.1721.0770.250.0581.1270.1790.028(70,35)1.5430.6460.270.9930.1320.0181.0430.0610.004(70,45)1.6240.5320.1751.0740.2450.0561.1240.1750.027(80,35)1.6260.530.1731.0760.2480.0571.1260.1770.028(80,45)1.6250.530.1731.0750.2480.0571.1250.1770.028 CS 2 with(?=1.6,γ =1.1,κ=1.2) and(τ =1.0,t0=1.25)(60,35)1.6280.5270.1711.0780.2510.0581.1280.180.029(60,45)1.6270.5270.1711.0770.2510.0581.1270.180.029(70,35)1.6250.530.1731.0750.2480.0571.1250.1770.028(70,45)1.6480.4980.1511.0980.2790.0711.1480.2090.038(80,35)1.6130.5470.1861.0630.2310.051.1130.160.023(80,45)1.6250.5310.1731.0750.2470.0571.1250.1760.028 CS 2 with(?=1.6,γ =1.1,κ=1.2) and(τ =1.25,t0=1.5)(60,35)1.6280.5260.171.0780.2510.0591.1280.1810.029(60,45)1.6480.4970.151.0980.2810.0721.1480.210.038(70,35)1.5440.6450.270.9940.1320.0181.0440.0620.004(70,45)1.6260.5290.1721.0760.2480.0571.1260.1780.028(80,35)1.6520.4920.1471.1020.2850.0741.1520.2150.04(80,45)1.5440.6450.270.9940.1320.0181.0440.0620.004

Table 3:The MLEs of parameters of with their respective MSEs and RABs with Schemes 3 and 4

Table 3(continued).(n,m)?γκ MLEMSERABMLEMSERABMLEMSERAB CS 4 with(?=1.6,γ =1.1,κ=1.2) and(τ =1.25,t0=1.5)(60,35)1.6520.4920.1471.1020.2850.0741.1520.2150.04(60,45)1.6270.5270.1711.0770.250.0581.1270.180.029(70,35)1.6250.5310.1731.0750.2470.0571.1250.1760.028(70,45)1.6260.5290.1721.0760.2490.0531.1260.1780.028(80,35)1.6130.4470.1861.0630.2310.051.1130.160.023(80,45)1.5740.3450.1270.9940.1320.0181.0440.0620.004

Table 4:Lower(L)and Upper(U)limits and corresponding ILs of 95% ACIs and ST-B CIs with Schemes 1 and 2

Table 4(continued).(n,m)?γκ ACI(L,U,IL)ACI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)BCI(L,U,IL)BCI(L,U,IL)CS 1 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.25,t0=1.5)(60,35)1.086,2.17,1.084 0.826,1.202,0.376 1.064,1.192,0.184 1.032,2.06,1.04 0.784,1.141,0.361 1.02,1.143,0.175(60,45)1.089,2.169,1.082 0.826,1.204,0.378 1.064,1.194,0.137 1.035,2.059,1.038 0.784,1.143,0.363 1.02,1.145,0.13(70,35)1.083,2.171,1.088 0.827,1.2,0.373 1.063,1.191,0.128 1.029,2.061,1.044 0.785,1.139,0.358 1.019,1.142,0.122(70,45)1.078,2.174,1.096 0.827,1.198,0.371 1.064,1.188,0.124 1.025,2.064,1.052 0.785,1.137,0.356 1.02,1.139,0.118(80,35)1.078,2.174,1.096 0.827,1.198,0.371 1.064,1.188,0.136 1.025,2.064,1.052 0.785,1.137,0.356 1.02,1.139,0.129(80,45)1.074,2.176,1.134 0.827,1.196,0.394 1.064,1.186,0.127 1.021,2.066,1.088 0.785,1.136,0.378 1.02,1.137,0.121 CS 2 with(?=1.6,γ =1.1,κ=1.2)and(τ =0.75,t0=1.0)(60,35)1.162,2.134,1.297 0.819,1.251,0.432 1.062,1.234,0.172 1.105,2.026,1.245 0.778,1.188,0.415 1.018,1.183,0.163(60,45)1.078,2.174,1.062 0.827,1.198,0.371 1.064,1.188,0.124 1.025,2.064,1.019 0.785,1.137,0.356 1.02,1.139,0.118(70,35)1.083,2.171,1.088 0.826,1.2,0.374 1.063,1.191,0.128 1.029,2.061,1.044 0.784,1.139,0.359 1.019,1.142,0.122(70,45)1.078,2.174,1.096 0.827,1.198,0.371 1.064,1.188,0.124 1.025,2.064,1.052 0.785,1.137,0.356 1.02,1.139,0.118(80,35)1.074,2.176,1.102 0.827,1.196,0.369 1.064,1.186,0.122 1.021,2.066,1.058 0.785,1.136,0.354 1.02,1.137,0.116(80,45)1.081,2.173,1.083 0.827,1.2,0.373 1.064,1.19,0.126 1.028,2.063,1.039 0.785,1.139,0.358 1.02,1.141,0.12 CS 2 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.0,t0=1.25)(60,35)1.027,2.199,1.172 0.832,1.168,0.336 1.063,1.163,0.165 0.976,2.088,1.125 0.79,1.109,0.322 1.019,1.115,0.157(60,45)1.162,2.134,0.972 0.818,1.252,0.434 1.062,1.234,0.157 1.105,2.026,0.933 0.777,1.189,0.417 1.018,1.183,0.149(70,35)1.17,2.13,0.969 0.818,1.256,0.438 1.062,1.238,0.172 1.112,2.022,0.93 0.777,1.193,0.42 1.018,1.187,0.163(70,45)1.975,1.979,1.108 0.682,2.515,0.233 0.587,2.367,0.148 1.877,1.879,1.063 0.648,2.388,0.224 0.563,2.269,0.141(80,35)1.072,2.178,1.106 0.828,1.195,0.367 1.064,1.186,0.127 1.019,2.068,1.061 0.786,1.135,0.352 1.02,1.137,0.121(80,45)1.078,2.174,1.096 0.827,1.198,0.371 1.064,1.188,0.124 1.025,2.064,1.052 0.785,1.137,0.356 1.02,1.139,0.118

Table 4(continued).(n,m)?γκ ACI(L,U,IL)ACI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)BCI(L,U,IL)BCI(L,U,IL)CS 2 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.25,t0=1.5)(60,35)1.164,2.132,0.978 0.817,1.253,0.457 1.062,1.234,0.176 1.106,2.024,0.939 0.776,1.19,0.439 1.018,1.183,0.167(60,45)1.151,2.139,0.909 0.819,1.244,0.425 1.063,1.227,0.169 1.094,2.031,0.872 0.778,1.181,0.408 1.019,1.176,0.161(70,35)1.116,2.156,1.046 0.824,1.221,0.397 1.064,1.208,0.145 1.061,2.047,1.004 0.782,1.159,0.381 1.02,1.158,0.138(70,45)1.083,2.171,1.083 0.826,1.2,0.334 1.063,1.191,0.135 1.029,2.061,1.039 0.784,1.139,0.321 1.019,1.142,0.128(80,35)1.072,2.178,1.117 0.828,1.195,0.369 1.064,1.186,0.128 1.019,2.068,1.072 0.786,1.135,0.354 1.02,1.137,0.122(80,45)1.078,2.174,1.087 0.827,1.198,0.372 1.063,1.189,0.123 1.025,2.064,1.043 0.785,1.137,0.357 1.019,1.14,0.117

Table 5:Lower(L)and Upper(U)limits and corresponding ILs of 95% ACIs and ST-B CIs with Schemes 3 and 4

Table 5(continued).(n,m)?γκ ACI(L,U,IL)BCI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)CS 3 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.0,t0=1.25)(80,35)1.162,2.134,0.972 0.818,1.252,0.434 1.062,1.234,0.157 1.105,2.026,0.933 0.777,1.189,0.417 1.018,1.183,0.149(80,45)1.17,2.13,0.969 0.818,1.256,0.438 1.062,1.238,0.172 1.112,2.022,0.93 0.777,1.193,0.42 1.018,1.187,0.163 CS 3 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.25,t0=1.5)(60,35)1.082,2.205,1.12 0.829,1.189,0.353 1.057,1.173,0.194 1.029,2.094,1.075 0.787,1.129,0.339 1.013,1.124,0.184(60,45)1.162,2.134,0.972 0.818,1.252,0.434 1.062,1.234,0.157 1.105,2.026,0.933 0.777,1.189,0.417 1.018,1.183,0.149(70,35)1.17,2.13,0.969 0.818,1.256,0.438 1.062,1.238,0.172 1.112,2.022,0.93 0.777,1.193,0.42 1.018,1.187,0.163(70,45)1.168,2.164,0.988 0.821,1.241,0.394 1.057,1.217,0.164 1.11,2.055,0.948 0.78,1.178,0.378 1.013,1.167,0.156(80,35)1.151,2.139,0.909 0.819,1.244,0.425 1.063,1.227,0.169 1.094,2.031,0.872 0.778,1.181,0.408 1.019,1.176,0.161(80,45)1.085,2.204,1.095 0.829,1.191,0.364 1.058,1.174,0.121 1.031,2.093,1.051 0.787,1.131,0.349 1.014,1.125,0.115 CS 4 with(?=1.6,γ =1.1,κ=1.2)and(τ =0.75,t0=1.0)(60,35)1.085,2.204,1.117 0.829,1.191,0.32 1.058,1.174,0.199 1.031,2.093,1.072 0.787,1.131,0.307 1.014,1.125,0.189(60,45)1.178,2.158,0.968 0.818,1.249,0.453 1.056,1.222,0.174 1.12,2.049,0.929 0.777,1.186,0.435 1.012,1.171,0.165(70,35)1.165,2.133,0.996 0.818,1.254,0.426 1.063,1.235,0.163 1.107,2.025,0.956 0.777,1.191,0.409 1.019,1.184,0.155(70,45)0.929,2.359,1.63 0.862,1.028,0.266 1.036,1.052,0.156 0.883,2.24,1.564 0.818,0.976,0.255 0.993,1.008,0.148(80,35)1.074,2.176,1.12 0.827,1.196,0.369 1.064,1.186,0.132 1.021,2.066,1.075 0.785,1.136,0.354 1.02,1.137,0.125(80,45)1.075,2.177,1.113 0.828,1.197,0.369 1.065,1.187,0.147 1.022,2.067,1.068 0.786,1.137,0.354 1.021,1.138,0.14 CS 4 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.0,t0=1.25)(60,35)1.725,2.165,0.97 0.82,1.24,0.421 1.057,1.215,0.167 1.107,2.056,0.864 0.779,1.177,0.404 1.013,1.165,0.179(60,45)1.192,2.152,1.185 0.818,1.257,0.446 1.055,1.231,0.19 1.133,2.043,1.137 0.777,1.194,0.428 1.011,1.18,0.181(70,35)1.173,2.161,0.969 0.819,1.246,0.415 1.056,1.22,0.166 1.115,2.052,0.93 0.778,1.183,0.398 1.012,1.169,0.158(70,45)1.074,2.176,1.12 0.827,1.196,0.369 1.064,1.186,0.132 1.021,2.066,1.075 0.785,1.136,0.354 1.02,1.137,0.125(80,35)0.725,2.361,1.102 0.861,1.027,0.369 1.036,1.05,0.129 0.689,2.242,1.058 0.818,0.975,0.354 0.993,1.007,0.123(80,45)1.069,2.179,1.11 0.829,1.192,0.363 1.064,1.184,0.126 1.016,2.069,1.065 0.787,1.132,0.348 1.02,1.135,0.12

Table 5(continued).(n,m)?γκ ACI(L,U,IL)BCI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)ACI(L,U,IL)BCI(L,U,IL)CS 4 with(?=1.6,γ =1.1,κ=1.2)and(τ =1.25,t0=1.5)(60,35)1.178,2.158,0.968 0.818,1.249,0.453 1.056,1.222,0.174 1.12,2.049,0.929 0.777,1.186,0.435 1.012,1.171,0.165(60,45)1.17,2.13,0.969 0.818,1.256,0.438 1.062,1.238,0.172 1.112,2.022,0.93 0.777,1.193,0.42 1.018,1.187,0.163(70,35)1.162,2.134,0.972 0.818,1.252,0.434 1.062,1.234,0.157 1.105,2.026,0.933 0.777,1.189,0.417 1.018,1.183,0.149(70,45)1.083,2.171,1.068 0.826,1.2,0.374 1.063,1.191,0.128 1.029,2.061,1.025 0.784,1.139,0.359 1.019,1.142,0.122(80,35)1.081,2.173,1.092 0.827,1.2,0.373 1.064,1.19,0.126 1.028,2.063,1.048 0.785,1.139,0.358 1.02,1.141,0.12(80,45)1.074,2.176,1.136 0.827,1.196,0.466 1.064,1.186,0.114 1.021,2.066,1.57 0.785,1.136,0.447 1.02,1.137,0.108

Figure 4:The plot of simulated samples and histogram of the parameters(?,γ,κ) respectively with Scheme 1

Fig.4 presents the plots of simulated samples and the histograms of the parameters to verify the convergence of the parameters.For this demonstration,plots are made using the Scheme 1 simulation results.The graphs indicate that the parameters are converging,and the histograms reveal that they are asymptotically normal over large number of simulation runs.As a result,the simulation study is consistent with the statistical assumptions for parameter estimation.The step-by-step process of the proposed estimation method is demonstrated through a flow chart given in Fig.5.

Figure 5:The flow chart to demonstrate the proposed estimation method

From the results listed in Tables 2 and 3,it can be observed that,for fixedτandt0,the MSEs and RABs of the model parameters are decreasing in most of the test schemes as a result of an increase in the values ofnandmand MLEs ??,?γand ?κof the model parameters?,γandκare getting more closer to their respective real values and this is reasonable because when sample data is large,we can expect more precise estimates.A decreasing pattern in the values of MSEs and RABs can also be observed asnandmincreases simultaneously in all cases for fixed values ofn,τandt0.For fixednandm,the MSEs and RABs are also getting smaller as the values ofτandt0getting larger and MLEs are approaching to their respective real values and this is also very much possible because with a larger experimental time,we can expect more failures and hence sample size will be larger in that case.Results reported in Tables 4 and 5 shows that expected lengths of both 95% ACIs and ST-B CIs are getting shorter with increasing values ofnandm,for fixed values ofτandt0in almost all cases except for?which is quite acceptable since?is not a parameter of the distribution under study instead it is acceleration factor.It is also observed that both 95% ACIs and ST-B CIs are getting shorter when the values ofnandmincreases simultaneously.As a comparison between 95% ACIs and ST-B CIs,it is observed that the ST-B CIs provides more narrower expected ILs in almost every case.Therefore,now based on the above findings,we can conclude that the proposed model and methods of estimation in this paper have performed very well,and hence all the statistical assumptions for fitting the model and regarding the estimation are satisfactory.

5 Real Data Application

In this section,an actual data set is utilized to further illustrate the performance of the suggested estimation technique and to demonstrate the application of the PLHR distribution in practice in the field of reliability engineering.The R statistical programming language is used for computation.The dataset reported in Table 6 is an uncensored dataset consisting of real failure times of an airplane’s air conditioning systems(in hours)was first discussed by[68].

Table 6:Failure times of the air conditioning system of an airplane(in hours)

The Kolmogorov-Smirnov(K-S)goodness of fit test is used to fit the PLHR distribution to real data.The K-S test is used to compare a given data sample with a reference continuous probability distribution.It is based on the K-S distance,which is the absolute maximum distance between the sample’s empirical distribution function and the reference distribution’s cumulative distribution function and its corresponding p-value.In this current example,the K-S distance was determined to be 0.14146 with a p-value of 0.5856,which is greater than 0.05.Fig.6 displays the plots of the empirical CDFvs.fitted CDF of PLHR distribution and histogram of datavs.fitted PDF of PLHR distribution.Consequently,it is evident from the K-S distance,p-value and Fig.6 that the PHLR distribution and considered sample data tabulated in Table 6 both have the same probability distribution.Therefore,given data can be used as an illustration for our model.

Now under SSPALT,for illustrative purposes,let’s consider the stress change timeτto be 24,the test termination timet0to be 130,and the sample sizemto be 20 with respect to the progressive censoring scheme(e1=e2=...=e6=0,e7=1,e8=0,e9=e10=1,e11=0,e12=e13=1,e14=0,e15=1,e16=e17=0,e18=e19=0,e20=1).So,utilizing the data provided in Table 6 with 0% masking,we have the following PrHC data reported in Table 7 on normal and accelerated stress.

For PrHC data in Table 7 with 0% masking under SSPALT,the MLEs of the parameters with initial valuesγ=0.042003 andκ=?0.3054115 which are the estimates of the parameters based on the complete data.And the initial value of the acceleration factor?is set to be 1.2.The MLEs with their corresponding standard errors are reported in Table 9.

Table 7:PrHC data of air conditioning systems with 0% masking at each stress level

Figure 6:Plot of empirical CDF vs.PLHR CDF and histogram of data vs.fitted PDF of the PLHR distribution

Now,to demonstrate the effect of masking,at each stress level,20% of failed systems are chosen at random to be masked.So,utilizing the data provided in Table 7 with 20% masking,the data obtained after masking with the PrHC scheme under SSPALT at each stress level is reported in the Table 8 as follows:

Table 8:PrHC data of air conditioning systems with 20% masking at each stress level

For PrHC data in Table 8 with 20% masking under SSPALT,we choose the same initial values ofγ,κand?as it was in the case of 0% masking.The MLEs with their corresponding standard errors are obtained and reported in Table 9.

From Table 9,it can be observed that the ML estimates are more accurate with smaller standard errors when the sample size is large or the masking proportion is 0%,than the estimates under a smaller sample size or the masking level is 20%.

Table 9:MLEs of parameters with their corresponding standard errors based on real data set

6 Conclusions

In this article,the SSPALT model has been developed for PrHC data to analyse the lifetime of a hybrid system of tree components under mask causes of failure.Assuming that the failure of components independently follows PLHR distribution,estimates of the parameters of PLHR distribution and the acceleration factor are then obtained using the MLE technique.The performance of MLEs is investigated through their respective MSEs and RABs.95% ACIs and ST-B CIs are also constructed and their performance is investigated in terms of their respective ILs.A simulation study has also been conducted to investigate and compare the performance of estimates for the hybrid system under SSPALT for PrHC masked data by utilizing the Monte Carlo simulation technique.Additionally,a real-world data application for an airplane’s air conditioning systems was utilized to demonstrate the proposed approach.As a comparison between 95% ACIs and STB CIs,it is observed that the ST-B CIs provide narrower expected ILs in almost every case.Based on the findings,it can be concluded that the proposed model and method of estimation performed well.Hence all the statistical assumptions for fitting the model and regarding the estimation are satisfactory.As a future research project,the present study may be extended to more complex systems using the Bayesian estimation technique with different censored data.

Data Availability:The data used in this paper is available in the paper.

Funding Statement:The authors received no specific funding for this study.

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

Appendix A