Modelling and Assessment of Survival Probability of Shock Model with Two Kinds of Shocks ()
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
. 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
. 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
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:
(1)
The first part of above expression under summation represents the probability that the component experiences k number of damage shocks during
, 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
the component do not experience a fatal shock. During
, the component may experience 0 shock or 1 shock or 2 shocks and so on … and hence the summation with
Taking the terms independent of k outside the summation and further simplification leads to
(2)
, 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
are subjected to life testing experiment and the experiment is conducted until all of them fail. Out of r components, let
components fail due to damage shock i.e. damage exceeding the threshold u and
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
shock,
. Let
be a random variable (r.v.) representing the amount of damage due to jth damage shock to ith component.
’s are assumed to be independent exponential random variables with parameter
.
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
be the time epoch at which jth shock has occurred to ith component,
. The inter arrival times
are exponential random variables with parameter
. For
components that fail due to damage exceeding threshold, the joint distribution of r.v.s
is
(3)
Similarly, for
components that fail due to fatal shock at
shock, the joint distribution of r.v.’s
is
(4)
Combining the two cases, the joint distribution
of all the involved r.v.’s is given by
(5)
where,
Taking logarithm of sample likelihood (
), equating partial derivatives of
with respect to
and
to zero, we get maximum likelihood estimators of
and
as
(6)
MLE
of survival function
is obtained using the MLEs of parameters [16] in expression (2).
Asymptotic Distribution: To obtain the asymptotic distribution of
, let us denote the Fisher Information Matrix of
as
(7)
where,
(8)
(9)
From the asymptotic properties of MLE under regularity conditions and multivariate central limit theorem we have,
is inverse of Fisher Information Matrix “
” and it is given by
(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
are generated as follows:
Step 1: A random number
is generated from
. If
, then it is considered that component has failed due to a damage shock.
Step 2: Initialize
; for
the exponential random variable
with parameter
is generated and the value of
is incremented by one. The process of generation of exponential random variables and incrementation of
is carried out until
is greater than u.
Step 3:
number of inter arrival times having exponential distribution with parameter
are generated. Addition of these inter arrival times results in
. Thus we get
.
Step 4: If
, it is considered that the component failure is due to a fatal shock. The process of generation of
is same as Step 2.
number of inter-arrival times having exponential distribution with parameter
are generated. Using these inter arrival times,
is obtained. One inter-arrival time
having exponential distribution with parameter
is generated.
Steps 1 to 4 are repeated for
. The real life examples that are quoted in Section 2 justifies the consideration of r as 20, 25 (moderate sample) and
(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
and
are computed using which estimators of
are obtained and for given set of parameters, the survival probability
is computed for
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
of the component is exponential random variable with parameter
.
The survival probability of the component at time t is given by:
(11)
, 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
having life distribution
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
of all the involved r.v.’s is given by
(12)
where,
and the maximum likelihood estimators of
and
are obtained as
(13)
where
and
are to be obtained numerically using Newton-Raphson method. MLE
of survival function
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
case. For the ith component, the random variables
are generated as follows:
For given value of
an exponential random variable
is generated and the remaining procedures are same as in Step 1 to Step 4 of Section 4. On the same lines estimators of
are computed for same set of r and t values of Section 4.
6. Comparison of Estimators
simulations are carried out to compute mean square errors (MSEs) of
using the formula
(14)
Table 1 and Table 2 provide the relative efficiencies (E) of
as compared to
which is computed as the ratio of
to
[17].
Table 1. Survival probabilities, MSEs and Efficiencies for
.
Figures in bold represents
.
Table 2. Survival probabilities, MSEs and Efficiencies for
.
Figures in bold represent
.
7. Results and Conclusion
From the tables it is clear that survival functions
and
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,
overestimates the true value of survival function
and
underestimates the true value of survival function
, 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.
performs better than
, 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.