Modelling and Assessment of Survival Probability of Shock Model with Two Kinds of Shocks

Abstract

The study aims at modelling and assessment of survival probability of a component experiencing two kinds of shocks namely, damage shock and fatal shock. Shocks are occurring randomly in time as events of a Poisson Process. The two cases of fixed/random threshold of components are studied. Survival probabilities of proposed models are derived. Maximum likelihood estimators (MLEs) of survival probabilities are obtained using the data from life testing experiments. Fisher information and asymptotic distribution of MLEs of parameters are obtained when a constant threshold is considered. Computation and comparison of estimators of two cases (constant threshold and random threshold) are made through simulation studies. The study recommends the consideration of threshold as a random variable.

Share and Cite:

Munoli, S. and , Suhas (2019) Modelling and Assessment of Survival Probability of Shock Model with Two Kinds of Shocks. Open Journal of Statistics, 9, 484-493. doi: 10.4236/ojs.2019.94033.

1. Introduction

Study of reliability can be broadly classified into two major aspects i.e. reliability modelling and reliability assessments. Modelling reliability aims at development of reliability model under certain assumptions. Here one discusses about configuration of associated components, methods of reliability enhancement, the conditions under which the system functions at its best, etc. Here the point of interest is “system survival” and compensatory measures to increase system survival.

On the counter part, reliability assessment involves processes such as identification of various modes of system failure, accordingly defining life testing experiments, assigning parameter values to life distributions, repair distributions, threshold distributions and inferential procedures to evaluate the modelled reliability. To assess reliability, one needs data pertaining to system functioning or failure and is obtained through life testing experiments.

Shock models have been studied by several authors. Here are a few references of contributions towards these models. Esary et al. (1973) studied the life distribution properties of a device subjected to a sequence of shocks occurring randomly in time according to a homogeneous Poisson process [1]. Abdel-Hameed and Proschan (1973) extended the results of [1] by assuming shocks that occur according to a nonhomogeneous Poisson process [2]. Barlow and Proschan (1975) discussed different varieties of general and cumulative shock models [3]. Abdel-Hameed and Proschan (1975) and Klefsjo (1981) considered general case in which shocks occur according to a nonstationary pure birth process [4] [5]. Ross (1981) presented a generalized Poisson shock model using increasing failure rate average (IFRA) stochastic processes [6]. Posner and Zuckerman (1986) proposed a replacement model for a semi-Markov shock model with additive damage [7]. Anderson (1987) considered a general shock model in which the time intervals between shocks have infinite expectations [8]. Savits (1988) studied nonhomogeneous Poisson shock processes and the effect of different mean functions for the shock processes on the distributional properties of the joint component lifetimes [9]. Gut (1990) considered theory for stopped two-dimensional random walks to describe cumulative shock models [10]. Skoulakis (2000) studied system reliability subjected to shocks generated by renewal point process [11]. Mallor and Santos (2003) introduced a new shock model by considering the magnitude of the shock and interarrival time between shocks, which generalizes some of the classical shock models [12]. Munoli and Bhat (2011) derived the reliability function of non-accumulating damage shock model with successive shocks causing greater damage. They obtained Maximum Likelihood and Bayes estimators of reliability functions using the data from a type II censored sample without replacement life testing experiment [13]. Lindskog and McNeil (2003) studied Poisson shock models in the context of insurance loss modelling and credit risk modelling [14].

In most of the works cited above, it is observed that shocks are of one type, causing damage to the component (damage shock). But there are some catastrophic shocks whose impact is very high and the component fails at the occurrence of such shocks (fatal shock). The present study endeavours to model and assess survival probability of a component experiencing two kinds of shocks namely, damage shock and fatal shock. The two cases of fixed/random threshold of components are considered. Section 2 deals with the model description and discussions on real life examples. Assessment of Survival function is explained in Section 3. Section 4 deals with validation of derived results through simulation studies. The case of threshold being a random variable is considered in Section 5. Comparison of estimators obtained in two cases of fixed and random threshold is outlined in Section 6. Results and findings are also discussed in the same Section.

2. Modelling Survival Probability

A component or system is experiencing shocks occurring randomly in time as events of a Poisson process with intensity λ , λ > 0 . Shocks are of two types. One is fatal shock, which causes failure of the system or component. Another is damage shock, which causes some amount of damage to the component. Damages are non-accumulating. Damages (X) are assumed to be independent and identical exponential random variables with parameter θ , θ > 0 . The system or component fails either due to experiencing a fatal shock or whenever the amount of damage due to a damage shock exceeds its threshold u. If the damage does not exceed its threshold, the component functions as good as new one. The probability of occurrence of fatal shock is ( 1 p ) and that of damage shock is p.

Let us consider threshold of the component as a known constant. The survival probability of the component at time t is given by:

S 1 ( t ) = k = 0 e ( p λ t ) ( p λ t ) k k ! ( 1 e u θ ) k e ( 1 p ) λ t (1)

The first part of above expression under summation represents the probability that the component experiences k number of damage shocks during ( 0 , t ) , second part is the probability that the damages due to all k shocks are less than its threshold u. The third component is the probability that during ( 0 , t ) the component do not experience a fatal shock. During ( 0 , t ) , the component may experience 0 shock or 1 shock or 2 shocks and so on … and hence the summation with k = 0 , 1 , 2 , .

Taking the terms independent of k outside the summation and further simplification leads to

S 1 ( t ) = e λ t [ 1 p ( 1 e u θ ) ] (2)

S 1 ( 0 ) = 1 , S 1 ( ) = 0 , S 1 ( t ) t , implying that the survival probability of the component is non-increasing function of time t.

The proposed model can be used to describe following real life examples:

Example 1: Heart disease is the leading cause of death for both men and women. Heart is a muscular organ. Most common heart diseases are heart attack and cardiac arrest. But heart attacks don’t always kill instantly. Heart attacks occur when blood flow to the heart muscle is temporarily blocked, starving the muscle tissue of oxygen which causes scarring and damage to heart muscle (damage shock with amount of damage tolerable and the person survives with this heart attack). For a heart attack to quickly lead to death, the damage due to the heart needs to be great enough to cause the heart to beat irregularly and eventually stop entirely (damage due to a damage shock exceeding threshold). On the contrary, cardiac arrest occurs suddenly and often without warning. Cardiac arrest is abrupt loss of heart function, breathing and consciousness. It results from an electric disturbance in the heart that disrupts its pumping action, stopping blood flow to different organs and can lead to death (fatal shock). In this scenario, one will be interested to study the probability that the patient survives or fails with heart attack (damage shock) or with what probability the patient succumbs to cardiac arrest (fatal shock). By studying this, physicians, nutritionist can suggest remedial measures for: avoiding heart attack; if not avoidable, reducing the amount of damage and to make the cardiac arrest (fatal shock) a null event.

Example 2: Credit-scoring systems aid the decision of whether to grant credit to an applicant or not. Traditionally, this is done by estimating the probability that an applicant will default. In recent years the aim has been changing towards choosing the customers of higher profit. In this case it is important to know that when a customer will default (fatal shock). It is possible that if the time to default is long, the acquired interest will compensate or even exceed losses resulting from default. Another factor that affects profitability is the case in which customers close their account early or pay off the loan early by switching to another lender. Depending on when the actual repayment occurred, the lender will lose a proportion of the interest on the loan [15]. One can consider the number of times partial repayments made as damage shocks and the amount repaid is as damage caused due to damage shock. On the other hand, if a customer pays off the loan early then it will be considered as amount of damage due to a damage shock exceeding the threshold (outstanding loan amount). In this case, one will be interested to investigate the probability that a customer will default or the probability that a customer pays off the loan early.

3. Assessment of the Survival Probability

Suppose r components each with threshold u having life distribution H ( t ) = 1 S 1 ( t ) are subjected to life testing experiment and the experiment is conducted until all of them fail. Out of r components, let r 1 components fail due to damage shock i.e. damage exceeding the threshold u and r 2 ( = r r 1 ) components fail due to occurrence of fatal shock. It is assumed that damage due to fatal shock is not observable. Let the ith component fail at n i t h shock, i = 1 , 2 , , r . Let X i j , i = 1 , 2 , , r ; j = 1 , 2 , , n i be a random variable (r.v.) representing the amount of damage due to jth damage shock to ith component. X i j ’s are assumed to be independent exponential random variables with parameter θ , θ > 0 .

Whenever a component fails due to a damage shock, it is assumed that the damage due to a shock at which the component fails is not observable but is known to exceed the component’s threshold u. Let t i j be the time epoch at which jth shock has occurred to ith component, j = 1 , 2 , , n i ; i = 1 , 2 , , r . The inter arrival times ( t i j t i j 1 ) are exponential random variables with parameter p λ . For r 1 components that fail due to damage exceeding threshold, the joint distribution of r.v.s n i , t i 1 , t i 2 , , t i n i , X i 1 , , X i n i 1 is

i = 1 r 1 ( p λ ) n i e p λ t n i θ n i 1 e θ j = 1 n i 1 x i j e u θ (3)

Similarly, for r 2 components that fail due to fatal shock at n i t h shock, the joint distribution of r.v.’s n i , t i 1 , t i 2 , , t i n i , X i 1 , , X i n i 1 is

i = 1 r 2 ( p λ ) n i 1 e p λ t n i 1 θ n i 1 e θ j = 1 n i 1 x i j ( 1 p ) λ e ( 1 p ) λ e ( t n i t n i 1 ) (4)

Combining the two cases, the joint distribution L 1 of all the involved r.v.’s is given by

L 1 = p n . r 2 λ n . ( 1 p ) r 2 θ n . r e p λ t .. e λ t . e θ ( x .. + r 1 u ) (5)

where,

n . = n 1 + n 2 , t . = i = 1 r 2 ( t n i t n i 1 ) ,

x .. = i = 1 r j = 1 n i 1 x i j , t .. = i = 1 r 1 t n i + 2 i = 1 r 2 t n i 1 i = 1 r 2 t n i .

Taking logarithm of sample likelihood ( L 1 ), equating partial derivatives of log ( L 1 ) with respect to θ , p and λ to zero, we get maximum likelihood estimators of θ , p and λ as

θ ^ = n . r x .. + r 1 u p ^ = t . ( n . r 2 ) t . n . + t .. r 2 λ ^ = n . t . + r 2 t .. t . ( t .. + t . ) } (6)

MLE S ^ 1 ( t ) of survival function S 1 ( t ) is obtained using the MLEs of parameters [16] in expression (2).

Asymptotic Distribution: To obtain the asymptotic distribution of θ ^ , p ^ , λ ^ , let us denote the Fisher Information Matrix of θ , p , λ as I ( θ , p , λ )

I ( θ , p , λ ) = [ n . r θ 2 0 0 0 b E ( t .. ) 0 E ( t .. ) n . λ 2 ] (7)

where,

E ( t .. ) = i = 1 r 1 n i p λ + 2 i = 1 r 2 n i 1 p λ i = 1 r 2 ( n i 1 p λ + 1 ( 1 p ) λ ) , (8)

b = n . r p 2 + r 2 ( 1 p ) 2 (9)

From the asymptotic properties of MLE under regularity conditions and multivariate central limit theorem we have,

[ ( θ θ ^ ) , ( p p ^ ) , ( λ λ ^ ) ] N 3 ( 0 , I 1 )

I 1 is inverse of Fisher Information Matrix “ I ( θ , p , λ ) ” and it is given by

I 1 = [ θ 2 n . r 0 0 0 n . b n . ( E ( t .. ) λ ) 2 E ( t .. ) ( 1 / λ 2 ) [ b n . ( E ( t .. ) λ ) 2 ] 0 E ( t .. ) ( 1 / λ 2 ) [ b n . ( E ( t .. ) λ ) 2 ] b ( 1 / λ 2 ) [ b n . ( E ( t .. ) λ ) 2 ] ] (10)

4. Simulation Study

Validation of proposed model and computation of estimators for different combinations of parameter values is done using Monte-Carlo Simulation. For the ith component, the random variables t i 1 , t i 2 , , t i n i , X i 1 , , X i n i 1 are generated as follows:

Step 1: A random number V i is generated from U ( 0 , 1 ) . If 0 < V i < p = p 0 , then it is considered that component has failed due to a damage shock.

Step 2: Initialize n i = 0 ; for λ = λ 0 , θ = θ 0 the exponential random variable X i 1 with parameter θ 0 is generated and the value of n i is incremented by one. The process of generation of exponential random variables and incrementation of n i is carried out until X i n i is greater than u.

Step 3: n i number of inter arrival times having exponential distribution with parameter p 0 λ 0 are generated. Addition of these inter arrival times results in t i n i . Thus we get n i , X i 1 , , X i n i 1 , t n i .

Step 4: If V i p 0 , it is considered that the component failure is due to a fatal shock. The process of generation of n i , X i 1 , , X i n i 1 is same as Step 2. n i 1 number of inter-arrival times having exponential distribution with parameter p 0 λ 0 are generated. Using these inter arrival times, t i n i 1 is obtained. One inter-arrival time ( t i n i t i n i 1 ) having exponential distribution with parameter ( 1 p 0 ) λ 0 is generated.

Steps 1 to 4 are repeated for r = 20 , 25 , 35 , 50 . The real life examples that are quoted in Section 2 justifies the consideration of r as 20, 25 (moderate sample) and r = 30 , 50 (large samples). In the example of heart disease, the number of patients will be moderately large and among them the cases of cardiac arrests (death due to fatal shock) will be small. To get the non-zero number of cases of failure due to fatal shocks, r needs to be moderately large. The example of credit scoring systems, the number of loan (credit) accounts in any financial company is large and hence, the use of values of r as 35, 50.

The statistics n . , x .. , t .. and t . are computed using which estimators of S 1 ( t ) are obtained and for given set of parameters, the survival probability S 1 ( t ) is computed for t = 0.5 ( 0.25 ) 2.0 units of time.

5. Modelling Survival Probability When Threshold (Ui) Is Random Variable

The threshold of the component or system is considered as random variable by many researchers. It is trivial assumption as the damage sustaining capacity depends on threshold and threshold of the component may vary due to raw materials used in its manufacturing, technology with which it is manufactured, inbuilt capacity, physical properties, random factors that influence its shock sustaining capacity, etc. It is assumed that threshold U i of the component is exponential random variable with parameter σ , σ > 0 .

The survival probability of the component at time t is given by:

S 2 ( t ) = k = 0 e ( p λ t ) ( p λ t ) k k ! e ( 1 p ) λ t ( θ θ + σ ) k = e ( 1 p ) λ t e ( p λ t ) k = 0 ( θ θ + σ ) k ( p λ t ) k k ! = e λ t ( 1 p θ σ + θ ) (11)

S 2 ( 0 ) = 1 , S 2 ( ) = 0 , S 2 ( t ) t , implying that the survival probability of the component is non-increasing function of time t.

5.1. Assessment of the Survival Probability

Suppose r components each with threshold U i , i = 1 , 2 , , r having life distribution H ( t ) = 1 S 2 ( t ) are subjected to life testing experiment and the experiment is conducted until all of them fail. Continuing the life testing experiment as in Section 3.

The joint distribution L 2 of all the involved r.v.’s is given by

L 2 = p n . r 2 λ n . ( 1 p ) r 2 θ n . r e p λ t .. e λ t . e θ x .. ( σ σ + θ ) r 1 σ r e u . (12)

where,

n . = n 1 + n 2 , t . = i = 1 r 2 ( t n i t n i 1 ) , x .. = i = 1 r j = 1 n i 1 x i j ,

t .. = i = 1 r 1 t n i + i = 1 r 2 t n i 1 i = 1 r 2 t n i , u . = i = 1 r u i .

and the maximum likelihood estimators of p , λ , σ and θ are obtained as

p ^ = t . ( n . r 2 ) t . n . + t .. r 2 λ ^ = n . t . + r 2 t .. t . ( t .. + t . ) σ ^ = ( r 1 + r n . ) θ + x .. θ 2 n . r x .. θ θ ^ = σ ( σ i = 1 r u i r ) r 1 + r σ i = 1 r u i } (13)

where σ ^ and θ ^ are to be obtained numerically using Newton-Raphson method. MLE S ^ 2 ( t ) of survival function S 2 ( t ) is obtained using the MLEs of parameters in expression (11).

5.2. Simulation Study

The estimators are computed using Monte-Carlo simulation for random threshold U i case. For the ith component, the random variables U i , t i 1 , t i 2 , , t i n i , X i 1 , , X i n i 1 are generated as follows:

For given value of σ = σ 0 an exponential random variable U i is generated and the remaining procedures are same as in Step 1 to Step 4 of Section 4. On the same lines estimators of S 2 ( t ) are computed for same set of r and t values of Section 4.

6. Comparison of Estimators

M = 10000 simulations are carried out to compute mean square errors (MSEs) of S ^ i ( t ) , i = 1 , 2 using the formula

M S E i ( S ^ i ( t ) ) = j = 1 M ( S i ( t ) S ^ i j ( t ) ) 2 M , i = 1 , 2 (14)

Table 1 and Table 2 provide the relative efficiencies (E) of S ^ 2 ( t ) as compared to S ^ 1 ( t ) which is computed as the ratio of M S E 1 ( S ^ 1 ( t ) ) to M S E 2 ( S ^ 2 ( t ) ) [17].

Table 1. Survival probabilities, MSEs and Efficiencies for p = 0.4 , λ = 0.8 , u = 1.2 , σ = 0.83 ,   θ = 0.6 .

Figures in bold represents M S E 2 ( S ^ 2 ( t ) ) .

Table 2. Survival probabilities, MSEs and Efficiencies for p = 0.6 , λ = 0.4 , u = 0.9 , σ = 1.1 , θ = 0.7 .

Figures in bold represent M S E 2 ( S ^ 2 ( t ) ) .

7. Results and Conclusion

From the tables it is clear that survival functions S 1 ( t ) and S 2 ( t ) decrease as t increases. MLEs of parameters maximize the likelihood of the joint distribution function. Further, using invariance property of MLE, it is easy to obtain the MLE of survival function. In case of random threshold, the MLEs of parameters are obtained numerically and hence asymptotic distribution of MLEs cannot be obtained which is a limitation of this study. MLE, S ^ 2 ( t ) overestimates the true value of survival function S 2 ( t ) and S ^ 1 ( t ) underestimates the true value of survival function S 1 ( t ) , the reason for which may be simulated random variables in case of simulation study/a chance error if we are considering the real data. As the value of r increases, MSEs decrease. S ^ 2 ( t ) performs better than S ^ 1 ( t ) , so the study advocates consideration of threshold as a random variable and is a realistic assumption also.

Acknowledgements

The second author is thankful to Department of Science and Technology—Innovation in Science Pursuit for Inspired Research (DST-INSPIRE) for financial support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Esary, J.D., Marshall, A.W. and Proschan, F. (1973) Shock Models and Wear Processes. The Annals of Probability, 1, 627-649.
https://doi.org/10.1214/aop/1176996891
[2] Abdel-Hameed, M.S. and Proschan, F. (1973) Nonstationary Shock Models. Stochastic Processes and their Applications, 1, 383-404.
https://doi.org/10.1016/0304-4149(73)90019-7
[3] Barlow, R.E. and Proschan, F. (1975) Statistical Theory of Reliability and Life Testing. Holt. Rinhart and Winston, Inc., New York.
[4] Abdel-Hameed, M.S. and Proschan, F. (1975) Shock Models with Underlying Birth Processes. Journal of Applied Probability, 12, 18-28.
https://doi.org/10.2307/3212403
[5] Klefsjo, B. (1981) Survival under the Pure Birth Shock Model. Journal of Applied Probability, 18, 554-560.
https://doi.org/10.2307/3213305
[6] Ross, S.M. (1981) Generalised Poisson Shock Models. The Annals of Probability, 9, 896-898.
https://doi.org/10.1214/aop/1176994318
[7] Posner, M.J.M. and Zuckerman, D. (1986) Semi-Markov Shock Models with Additive Damage. Advances in Applied Probability, 18, 772-790.
https://doi.org/10.2307/1427187
[8] Anderson, K.K. (1987) Limit Theorems for General Shock Models with Infinite Mean Inter Shock Times. Journal of Applied Probability, 24, 449-456.
https://doi.org/10.2307/3214268
[9] Savits, T.H. (1988) Some Multivariate Distributions Derived from a Non-Fatal Shock Model. Journal of Applied Probability, 25, 383-390.
https://doi.org/10.2307/3214446
[10] Gut, A. (1990) Cumulative Shock Models. Advances in Applied Probability, 22, 504-507.
https://doi.org/10.2307/1427554
[11] Skoulakis, G. (2000) A General Shock Model for a Reliability System. Journal of Applied Probability, 37, 925-935.
https://doi.org/10.1239/jap/1014843073
[12] Mallor, F. and Santos, J. (2003) Classification of Shock Models in System Reliability. Monografias del Seminario Matematico García de Galdeano, 27, 405-412.
[13] Munoli, S.B. and Bhat, S.V. (2011) Reliability Estimation in Shock Model When Successive Shocks Cause Greater Damage (Non-Accumulating Damages). Research Journal of Mathematics and Statistics, 3, 61-66.
https://doi.org/10.1155/2011/721962
[14] Lindskog, F., and McNeil, A.J. (2003) Common Poisson Shock Models: Applications to Insurance and Credit Risk Modeling. ASTIN Bulletin: The Journal of the IAA, 33, 209-238.
https://doi.org/10.1017/S0515036100013441
[15] Stepanova, M. and Thomas, L. (2002) Survival Analysis Methods for Personal Loan Data. Operations Research, 50, 277-289.
https://doi.org/10.1287/opre.50.2.277.426
[16] Deshpande, J.V. and Purohit, S.G. (2005) Life Time Data: Statistical Models and Methods. World Scientific Publishing Company, Singapore.
https://doi.org/10.1142/5988
[17] Kale, B.K. (1999) A First Course on Parametric Inference. Narosa Publishing House, New Delhi.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.