Non-Parametric Recurrent Events Data Analysis: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(Created page with '==Recurrent Events Data Analysis== Recurrent Events Data Analysis, also called Recurrence Data Analysis (RDA), can be used in various applied fields such as reliability, medici…')
 
(Redirected page to Recurrent Event Data Analysis)
 
(40 intermediate revisions by 4 users not shown)
Line 1: Line 1:
==Recurrent Events Data Analysis==
#REDIRECT [[Recurrent_Event_Data_Analysis]]
 
Recurrent Events Data Analysis, also called Recurrence Data Analysis (RDA),  can be used in various applied fields such as reliability, medicine, social sciences, economics, business and criminology.
 
Whereas in life data analysis (LDA) it was assumed that events (failures) were independent and identically distributed (iid), there are many cases where events are dependent and not identically distributed (such as repairable system data) or where the analyst is interested in modeling the number of occurrences of events over time rather than the length of time prior to the first event, as in LDA.
 
Weibull++ provides both parametric and non-parametric approaches to analyze such data.
 
• The non-parametric approach is based on the well-known Mean Cumulative Function (MCF). The Weibull++ module for this type of analysis builds upon the work of Dr. Wayne Nelson, who has written extensively on the calculation and applications of MCF [31].
 
• The parametric approach is based on the General Renewal Process (GRP) model, which is particularly useful in understanding the effects of the repairs on the age of a system. Traditionally, the commonly used models for analyzing repairable systems data are perfect renewal processes (PRP), corresponding to perfect repairs, and nonhomogeneous Poisson processes (NHPP), corresponding to minimal repairs. However, most repair activities may realistically not result in such extreme situations but in a complicated intermediate one (general repair or imperfect repair/maintenance), which are well treated with the GRP model.
 
===Non-Parameteric Recurrence Data Analysis===
====Introduction====
Non-parametric recurrence data analysis provides a nonparametric graphical estimate of the mean cumulative number or cost of recurrence per unit versus age. In the reliability field, the Mean Cumulative Function (MCF) can be used to: [31]
 
• Evaluate whether the population repair (or cost) rate increases or decreases with age (this is useful for product retirement and burn-in decisions).
 
• Estimate the average number or cost of repairs per unit during warranty or some time period.
 
• Compare two or more sets of data from different designs, production periods, maintenance policies, environments, operating conditions, etc.
 
• Predict future numbers and costs of repairs, such as, the next month, quarter, or year.
 
• Reveal unexpected information and insight.
 
====The Mean Cumulative Function (MCF) for Recurrence Data====
 
In non-parametric analysis of recurrent events data, each population unit can be described by a cumulative history function for the cumulative number of recurrences. It is a staircase function that depicts the cumulative number of recurrences of a particular event, such as repairs, over time. Figure CHM depicts a unit's cumulative history function.
 
The nonparametric model for a population of units is described as the population of cumulative history functions (curves). It is the population of all staircase functions of every unit in the population. At age t, the units have a distribution of their cumulative number of events. That is, a fraction of the population has accumulated 0 recurrences, another fraction has accumulated 1 recurrence, another fraction has accumulated 2 recurrences, etc. This distribution differs at different ages,  <math>t</math> , and has a mean  <math>M(t)</math>  called the mean cumulative function (MCF) . The  <math>M(t)</math>  is the pointwise average of all population cumulative history functions(see Fig.MCFDemo).
 
 
For the case of uncensored data, the mean cumulative function  <math>M{{(t)}_{i}}\ </math> values at different recurrence ages  <math>{{t}_{i}}</math>  are estimated by calculating the average of the cumulative number of recurrences of events for each unit in the population at  <math>{{t}_{i}}</math> .
When the histories are censored, the following steps are applied.
 
1st Step - Order All Ages: Order all recurrence and censoring ages from smallest to largest. If a recurrence age for a unit is the same as its censoring (suspension) age, the recurrence age goes first. If multiple units have a common recurrence or censoring age, then these units could be put in a certain order or be sorted randomly.
 
2nd Step - Calculate the Number,  <math>{{r}_{i}}</math> , of Units that Passed Through Age  <math>{{t}_{i}}</math> :   
 
<math>\begin{align}
  & {{r}_{i}}= & {{r}_{i-1}}\quad \quad \text{if }{{t}_{i}}\text{ is a recurrence age} \\
& {{r}_{i}}= & {{r}_{i-1}}-1\text{  if }{{t}_{i}}\text{ is a censoring age} 
\end{align}</math>
 
<math>N</math> is the total number of units and  <math>{{r}_{1}}=N</math>  at the first observed age which could be a recurrence or suspension.
 
3rd Step - Calculate  <math>MCF</math>  Estimate,  <math>{{M}^{*}}(t)</math>: For each sample recurrence age  <math>{{t}_{i}},</math>  calculate the mean cumulative function estimate as follows:       
 
<math>{{M}^{*}}({{t}_{i}})=\frac{1}{{{r}_{i}}}+{{M}^{*}}({{t}_{i-1}})</math>
 
where  <math>{{M}^{*}}(t)=\tfrac{1}{{{r}_{1}}}</math>  at the earliest observed recurrence age,  <math>{{t}_{1}}</math> .
 
====Example 1====
A health care company maintains five identical pieces of equipment used by a hospital. When a piece of equipment fails, the company sends a crew to repair it. The following table gives the failure and censoring ages for each machine, where the + sign indicates a censoring age.
 
 
<math>\begin{matrix}
  Equipment ID & Months  \\
  \text{1} & \text{5, 10 , 15, 17+}  \\
  \text{2} & \text{6, 13, 17, 19+}  \\
  \text{3} & \text{12, 20, 25, 26+}  \\
  \text{4} & \text{13, 15, 24+}  \\
  \text{5} & \text{16, 22, 25, 28+}  \\
\end{matrix}</math>
<math></math>
Estimate the MCF values and the 95% confidence limits ignoring the repair duration.
====Solution to Example 1====
The MCF estimate is obtained as follows:
 
 
<math>\begin{matrix}
  ID & Months ({{t}_{i}}) & State & {{r}_{i}} & 1/{{r}_{i}} & {{M}^{*}}({{t}_{i}})  \\
  \text{1} & \text{5} & \text{F} & \text{5} & \text{0}\text{.20} & \text{0}\text{.20}  \\
  \text{2} & \text{6} & \text{F} & \text{5} & \text{0}\text{.20} & \text{0}\text{.20 + 0}\text{.20 = 0}\text{.40}  \\
  \text{1} & \text{10} & \text{F} & \text{5} & \text{0}\text{.20} & \text{0}\text{.40 + 0}\text{.20 = 0}\text{.60}  \\
  \text{3} & \text{12} & \text{F} & \text{5} & \text{0}\text{.20} & \text{0}\text{.60 + 0}\text{.20 = 0}\text{.80}  \\
  \text{2} & \text{13} & \text{F} & \text{5} & \text{0}\text{.20} & \text{0}\text{.80+0}\text{.20 =1}\text{.00}  \\
  \text{4} & \text{13} & \text{F} & \text{5} & \text{0}\text{.20} & \text{1}\text{.00 + 0}\text{.20 = 1}\text{.20}  \\
  \text{1} & \text{15} & \text{F} & \text{5} & \text{0}\text{.20} & \text{1}\text{.20 + 0}\text{.20 =1}\text{.40}  \\
  \text{4} & \text{15} & \text{F} & \text{5} & \text{0}\text{.20} & \text{1}\text{.40 + 0}\text{.20 = 1}\text{.60}  \\
  \text{5} & \text{16} & \text{F} & \text{5} & \text{0}\text{.20} & \text{1}\text{.60 + 0}\text{.20 = 1}\text{.80}  \\
  \text{2} & \text{17} & \text{F} & \text{5} & \text{0}\text{.20} & \text{1}\text{.80 + 0}\text{.20 = 2}\text{.0}  \\
  \text{1} & \text{17} & \text{S} & \text{4} & {} & {}  \\
  \text{2} & \text{19} & \text{S} & \text{3} & {} & {}  \\
  \text{3} & \text{20} & \text{F} & \text{3} & \text{0}\text{.33} & \text{2}\text{.00 + 0}\text{.33 = 2}\text{.33}  \\
  \text{5} & \text{22} & \text{F} & \text{3} & \text{0}\text{.33} & \text{2}\text{.33 + 0}\text{.33 = 2}\text{.66}  \\
  \text{4} & \text{24} & \text{S} & \text{2} & {} & {}  \\
  \text{3} & \text{25} & \text{F} & \text{2} & \text{0}\text{.50} & \text{2}\text{.66 + 0}\text{.50 = 3}\text{.16}  \\
  \text{5} & \text{25} & \text{F} & \text{2} & \text{0}\text{.50} & \text{3}\text{.16 + 0}\text{.50 = 3}\text{.66}  \\
  \text{3} & \text{26} & \text{S} & \text{1} & {} & {}  \\
  \text{5} & \text{28} & \text{S} & \text{0} & {} & {}  \\
\end{matrix}</math>
 
Confidence Limits for the MCF
Upper and lower conifidence limits for  <math>M({{t}_{i}})</math>  are:
 
 
<math>\begin{align}
  & {{M}_{U}}({{t}_{i}})= & {{M}^{*}}({{t}_{i}}).{{e}^{\tfrac{{{K}_{\alpha }}.\sqrt{Var[{{M}^{*}}({{t}_{i}})]}}{{{M}^{*}}({{t}_{i}})}}} \\
& {{M}_{L}}({{t}_{i}})= & \frac{{{M}^{*}}({{t}_{i}})}{{{e}^{\tfrac{{{K}_{\alpha }}.\sqrt{Var[{{M}^{*}}({{t}_{i}})]}}{{{M}^{*}}({{t}_{i}})}}}} 
\end{align}</math>
 
 
where  <math>\alpha </math>  ( <math>50%<\alpha <100%</math> ) is  confidence level,  <math>{{K}_{\alpha }}</math>  is the  <math>\alpha </math>  standard normal percentile and  <math>Var[{{M}^{*}}({{t}_{i}})]</math>  is the variance of the MCF estimate at recurrence age  <math>{{t}_{i}}</math> . The variance is calculated as follows:
 
 
<math>Var[{{M}^{*}}({{t}_{i}})]=Var[{{M}^{*}}({{t}_{i-1}})]+\frac{1}{r_{i}^{2}}\left[ \underset{j\in {{R}_{i}}}{\overset{}{\mathop \sum }}\,{{\left( {{d}_{ji}}-\frac{1}{{{r}_{i}}} \right)}^{2}} \right]</math>
 
 
where  <math>r</math>  is defined in Eqn.(R),  <math>{{R}_{i}}</math>  is the set of the units that have not been suspended by  <math>i</math>  and  <math>{{d}_{ji}}</math>  is defined as follows:
 
 
<math>\begin{align}
  & {{d}_{ji}}= & 1\text{  if the }{{j}^{\text{th }}}\text{unit had an event recurrence at age }{{t}_{i}} \\
& {{d}_{ji}}= & 0\text{  if the }{{j}^{\text{th }}}\text{unit did not have an event reoccur at age }{{t}_{i}} 
\end{align}</math>
 
====Example 2====
Using the data in Example 1, estimate the 95% confidence bounds.
====Solution to Example 2====
Using Eq. MCF Var the following table of variance values can be obtained:
 
 
Using Eqn.(MCF Bounds) and  <math>{{K}_{5}}=1.644</math>  for a 95% confidence level, the confidence bounds can be obtained as follows:
 
 
<math>\begin{matrix}
  ID & Months & State & MC{{F}_{i}} & Va{{r}_{i}} & MC{{F}_{{{L}_{i}}}} & MC{{F}_{{{U}_{i}}}}  \\
  \text{1} & \text{5} & \text{F} & \text{0}\text{.20} & \text{0}\text{.032} & 0.0459 & 0.8709  \\
  \text{2} & \text{6} & \text{F} & \text{0}\text{.40} & \text{0}\text{.064} & 0.1413 & 1.1320  \\
  \text{1} & \text{10} & \text{F} & \text{0}\text{.60} & \text{0}\text{.096} & 0.2566 & 1.4029  \\
  \text{3} & \text{12} & \text{F} & \text{0}\text{.80} & \text{0}\text{.128} & 0.3834 & 1.6694  \\
  \text{2} & \text{13} & \text{F} & \text{1}\text{.00} & \text{0}\text{.160} & 0.5179 & 1.9308  \\
  \text{4} & \text{13} & \text{F} & \text{1}\text{.20} & \text{0}\text{.192} & 0.6582 & 2.1879  \\
  \text{1} & \text{15} & \text{F} & \text{1}\text{.40} & \text{0}\text{.224} & 0.8028 & 2.4413  \\
  \text{4} & \text{15} & \text{F} & \text{1}\text{.60} & \text{0}\text{.256} & 0.9511 & 2.6916  \\
  \text{5} & \text{16} & \text{F} & \text{1}\text{.80} & \text{0}\text{.288} & 1.1023 & 2.9393  \\
  \text{2} & \text{17} & \text{F} & \text{2}\text{.0} & \text{0}\text{.320} & 1.2560 & 3.1848  \\
  \text{1} & \text{17} & \text{S} & {} & {} & {} & {}  \\
  \text{2} & \text{19} & \text{S} & {} & {} & {} & {}  \\
  \text{3} & \text{20} & \text{F} & \text{2}\text{.33} & \text{0}\text{.394} & 1.4990 & 3.6321  \\
  \text{5} & \text{22} & \text{F} & \text{2}\text{.66} & \text{0}\text{.468} & 1.7486 & 4.0668  \\
  \text{4} & \text{24} & \text{S} & {} & {} & {} & {}  \\
  \text{3} & \text{25} & \text{F} & \text{3}\text{.16} & \text{0}\text{.593} & 2.1226 & 4.7243  \\
  \text{5} & \text{25} & \text{F} & \text{3}\text{.66} & \text{0}\text{.718} & 2.5071 & 5.3626  \\
  \text{3} & \text{26} & \text{S} & {} & {} & {} & {}  \\
  \text{5} & \text{28} & \text{S} & {} & {} & {} & {}  \\
\end{matrix}</math>
 
The analysis presented in Examples 1 and 2 can be obtained automatically in Weibull ++ using the Non-Parametric RDA Specialized Folio, as shown next.
 
 
Note: In the above Folio, the  <math>F</math>  refers to failures and  <math>E</math>  refers to suspensions (or censoring ages).
The results with calculated MCF values and upper and lower 95% confidence limits are shown next along with the graphical plot.
 
 
====Example 3====
The following table displays transmission repairs on a sample of 14 cars with manual transmission in a preproduction road test [31]. Here + denotes the censoring ages (how long a car has been observed).
 
 
<math>\begin{matrix}
  Car ID & Mileage  \\
  \text{1} & \text{27099+}  \\
  \text{2} & \text{21999+}  \\
  \text{3} & \text{11891, 27583+}  \\
  \text{4} & \text{19966+}  \\
  \text{5} & \text{26146+}  \\
  \text{6} & \text{3648, 13957, 23193+}  \\
  \text{7} & \text{19823+}  \\
  \text{8} & \text{2890, 22707+}  \\
  \text{9} & \text{2714, 19275+}  \\
  \text{10} & \text{19803+}  \\
  \text{11} & \text{19630+}  \\
  \text{12} & \text{22056+}  \\
  \text{13} & \text{22940+}  \\
  \text{14} & \text{3240, 7690, 18965+}  \\
\end{matrix}</math>
 
The car manufacturer seeks to estimate the mean cumulative number of repairs per car by 24,000 test miles (equivalently 5.5 x 24,000 = 132,000 customer miles) and to observe whether the population repair rate increases or decreases as a population ages.
====Solution to Example 3====
The data is entered into a Non-Parametric RDA Specialized Folio in Weibull++ as follows.
 
The results are as follows,
 
The results indicate that after 13,957 miles of testing, the estimated mean cumulative number of repairs per car is 0.5. Therefore, by 24,000 test miles, the estimated mean cumulative number of repairs per car is 0.5.
The MCF plot is shown next.
 
A smooth curve through the MCF plot has a derivative that decreases as the population ages. That is, the repair rate decreases as each population ages. This is typical of products with manufacturing defects.
 
===Parametric Recurrence Data Analysis===
 
Weibull++'s Parametric RDA Specialized Folio is a tool to model recurrent event data.  It can capture the trend, estimate the rate and predict the total number of recurrences. The failure and repair data of a repairable system can be treated as one type of recurrence data. Past and current repairs may affect the future failure process. For most recurrence events, time (distance, cycles, etc.) is a key factor. With the time, the recurrence rate may keep constant, increase or decrease. For other recurrence events, not only the time, but also the number of events can affect the recurrence process, e.g. the debugging process in software development.
The parametric analysis approach utilizes the General Renewal Process (GRP) model [28]. In this model, the repair time is assumed to be negligible so that the processes can be viewed as point processes. This model provides a way to describe the rate of occurrence of events over time such as in the case of data obtained from a repairable system. This model is particularly useful in modeling the failure behavior of a specific system and understanding the effects of the repairs on the age of that system. For example, consider a system that is repaired after a failure, where the repair does not bring the system to an as-good-as-new or an as-bad-as-old condition. In other words, the system is partially rejuvenated after the repair. Traditionally in as-bad-as-old repairs, also known as minimal repairs, the failure data from such a system would have been modeled using a homogeneous or non-homogeneous Poisson process (NHPP). On rare occasions, a Weibull distribution has been used as well, in cases where the system is almost as-good-as-new after the repair, also known as perfect renewal process (PRP). However, for the intermediate states after the repair, there has not been a commercially available model, even though many models have been proposed in literature. In Weibull++, the GRP model provides the capability of modeling such systems with partial renewal (general repair or imperfect repair/maintenance) and allows for a variety of predictions such as reliability, expected failures, etc.
 
====The GRP Model====
 
In this model, the concept of virtual age is introduced. Denote by  <math>{{t}_{1}},{{t}_{2}},\cdots ,{{t}_{n}}</math>  the successive failure times and let  <math>{{x}_{1}},{{x}_{2}},\cdots ,{{x}_{n}}</math>  represent the time between failures ( <math>{{t}_{i}}=\sum_{j=1}^{i}{{x}_{j}})</math> . Assume that after each event, actions are taken to improve the system performance. Let  <math>q</math>  be the action effectiveness factor. There are two GRP models.
 
Type I:
 
<math>{{v}_{i}}={{v}_{i-1}}+q{{x}_{i}}=q{{t}_{i}}</math>
 
 
Type II:
 
 
<math>{{v}_{i}}=q({{v}_{i-1}}+{{x}_{i}})={{q}^{i}}{{x}_{1}}+{{q}^{i-1}}{{x}_{2}}+\cdots +{{x}_{i}}</math>
 
where  <math>{{v}_{i}}</math>  is the virtual age of the system right after  <math>i</math> th repair. The Type I model assumes that the  <math>i</math> th repair cannot remove the damage incurred before the ith failure. It can only reduce the additional age  <math>{{x}_{i}}</math>  to  <math>q{{x}_{i}}</math> .  The Type II model assumes that at the  <math>i</math> th repair, the virtual age has been accumulated to  <math>{{v}_{i-1}}+{{x}_{i}}</math> . The  <math>i</math> th repair will remove the cumulative damage from both current and previous failures by reducing the virtual age to  <math>q({{v}_{i-1}}+{{x}_{i}})</math> .
The power law function is used to model the rate of recurrence, which is:
 
 
<math>\lambda (t)=\lambda \beta {{t}^{\beta -1}}</math>
 
 
The conditional  <math>pdf</math>  is:
 
 
<math>f({{t}_{i}}|{{t}_{i-1}})=\lambda \beta {{({{x}_{i}}+{{v}_{i-1}})}^{\beta -1}}{{e}^{-\lambda \left[ {{\left( {{x}_{i}}+{{v}_{i-1}} \right)}^{\beta }}-v_{i-1}^{\beta } \right]}}</math>
 
 
MLE method is used to estimate model parameters. The log likelihood function is [28]:
 
 
<math>\begin{align}
  & \ln (L)= & n(\ln \lambda +\ln \beta )-\lambda \left[ {{\left( T-{{t}_{n}}+{{v}_{n}} \right)}^{\beta }}-v_{n}^{\beta } \right] \\
&  & -\lambda \underset{i=1}{\overset{n}{\mathop \sum }}\,\left[ {{\left( {{x}_{i}}+{{v}_{i-1}} \right)}^{\beta }}-v_{i}^{\beta } \right]+(\beta -1)\underset{i=1}{\overset{n}{\mathop \sum }}\,\ln ({{x}_{i}}+{{v}_{i-1}}) 
\end{align}</math>
 
where  <math>n</math>  is the total number of events during the entire observation period.  <math>T</math>  is the stop time of the observation.  <math>T={{t}_{n}}</math>  if the observation stops right after the last event.
 
====Confidence Bounds====
In general, in order to obtain the virtual age, the exact occurrence time of each event (failure) should be available, see Eqns.(Type I) and (Type 2). However, the times are unknown until the corresponding events occur. For this reason, there are no closed-form expressions for total failure number and failure intensity, which are functions of failure times and virtual age. Therefore, in Weibull++, Monte Carlo simulation is used to predict values of virtual time, failure number, MTBF and failure rate. The approximate confidence bounds obtained from simulation are provided. The uncertainty of model parameters is also considered in the bounds.
=====Bounds of Cumulative Failure (Event) Numbers=====
 
The variance of the cumulative failure number  <math>N(t)</math>  is:
 
 
<math>Var[N(t)]=Var\left[ E(N(t)|\lambda ,\beta ,q) \right]+E\left[ Var(N(t)|\lambda ,\beta ,q) \right]</math>
 
The first term accounts for the uncertainty of the parameter estimation. The second term considers the uncertainty caused by the renewal process even when model parameters are fixed. However, unless  <math>q=1</math> ,  <math>Var\left[ E(N(t)|\lambda ,\beta ,q) \right]</math>  cannot be calculated because  <math>E(N(t))</math>  cannot be expressed as a closed-form function of  <math>\lambda ,\beta ,</math>  and  <math>q</math> . In order to consider the uncertainty of the parameter estimation,  <math>Var\left[ E(N(t)|\lambda ,\beta ,q) \right]</math>  is approximated by:
 
<math>Var\left[ E(N(t)|\lambda ,\beta ,q) \right]=Var[E(N({{v}_{t}})|\lambda ,\beta )]=Var[\lambda v_{t}^{\beta }]</math>
 
where  <math>{{v}_{t}}</math>  is the expected virtual age at time  <math>t</math>  and  <math>Var[\lambda v_{t}^{\beta }]</math>  is::
 
 
<math>\begin{align}
  & Var[\lambda v_{t}^{\beta }]= & {{\left( \frac{\partial (\lambda v_{t}^{\beta })}{\partial \beta } \right)}^{2}}Var(\hat{\beta })+{{\left( \frac{\partial (\lambda v_{t}^{\beta })}{\partial \lambda } \right)}^{2}}Var(\hat{\lambda }) \\
&  & +2\frac{\partial (\lambda v_{t}^{\beta })}{\partial \beta }\frac{\partial (\lambda v_{t}^{\beta })}{\partial \lambda }Cov(\hat{\beta },\hat{\lambda }) 
\end{align}</math>
 
By conducting this approximation, the uncertainty of  <math>\lambda </math>  and  <math>\beta </math>  are considered. The value of  <math>{{v}_{t}}</math>  and the value of the second term in Eqn.(VarN) are obtained through the Monte Carlo simulation using parameters  <math>\hat{\lambda },\hat{\beta },\hat{q},</math>  which are the ML estimators. The same simulation is used to estimate the cumulative number of failures  <math>\hat{N}(t)=E(N(t)|\hat{\lambda },\hat{\beta },\hat{q})</math>  .
 
Once the variance and the expected value of  <math>N(t)</math>  have been obtained, the bounds can be calculated by assuming  <math>N(t)</math>  is lognormally distributed as:
 
 
<math>\frac{\ln N(t)-\ln \hat{N}(t)}{\sqrt{Var(\ln N(t))}}\tilde{\ }N(0,1)</math>
 
 
The upper and lower bounds for a given confidence level  <math>\alpha </math>  can be calculated by:
 
 
<math>N{{(t)}_{U,L}}=\hat{N}(t){{e}^{\pm {{z}_{a}}\sqrt{Var(N(t))}/\hat{N}(t)}}</math>
 
 
where  <math>{{z}_{a}}</math>  is the standard normal distribution.
If  <math>N(t)</math>  is assumed to be normally distributed, the bounds can be calculated by:
 
 
<math>N{{(t)}_{U}}=\hat{N}(t)+{{z}_{a}}\sqrt{Var(N(t))}</math>
 
 
<math>N{{(t)}_{L}}=\hat{N}(t)-{{z}_{a}}\sqrt{Var(N(t))}</math>
 
 
In Weibull++, the  <math>N{{(t)}_{U}}</math>  is the smaller of the upper bounds obtained from Eqns. (LogNorm) and (Norm). The  <math>N{{(t)}_{L}}</math>  is set to the largest of the lower bounds obtained from Eqns. (LogNorm) and (Norm). This combined method can prevent the out-of-range values of bounds for some small  <math>t</math>  values.
Bounds of Cumulative Failure Intensity and MTBF
For a given time  <math>t</math> , the expected value of cumulative MTBF  <math>{{m}_{c}}(t)</math>  and cumulative failure intensity  <math>{{\lambda }_{c}}(t)</math>  can be calculated using the following equations:
 
 
<math>{{\hat{\lambda }}_{c}}(t)=\frac{\hat{N}(t)}{t};{{\hat{m}}_{c}}(t)=\frac{t}{\hat{N}(t)}</math>
 
The bounds can be easily obtained from the corresponding bounds of  <math>N(t).</math>
 
 
<math>\begin{align}
  & {{{\hat{\lambda }}}_{c}}{{(t)}_{L}}= & \frac{\hat{N}{{(t)}_{L}}}{t};\text{  }{{{\hat{\lambda }}}_{c}}{{(t)}_{L}}=\frac{\hat{N}{{(t)}_{L}}}{t};\text{  } \\
& {{{\hat{m}}}_{c}}{{(t)}_{L}}= & \frac{t}{\hat{N}{{(t)}_{U}}};\text{  }{{{\hat{m}}}_{c}}{{(t)}_{U}}=\frac{t}{\hat{N}{{(t)}_{L}}} 
\end{align}</math>
 
 
Bounds of Instantaneous Failure Intensity and MTBF
 
The instantaneous failure intensity is given by:
 
 
<math>{{\lambda }_{i}}(t)=\lambda \beta v_{t}^{\beta -1}</math>
 
where  <math>{{v}_{t}}</math>  is the virtual age at time  <math>t.</math>  When  <math>q\ne 1,</math>  it is obtained from simulation. When  <math>q=1</math> ,  <math>{{v}_{t}}=t</math>  from Eqns. #(Type I) and (Type 2).
 
The variance of instantaneous failure intensity can be calculated by:
 
 
<math>\begin{align}
  & Var({{\lambda }_{i}}(t))= & {{\left( \frac{\partial {{\lambda }_{i}}(t)}{\partial \beta } \right)}^{2}}Var(\hat{\beta })+{{\left( \frac{\partial {{\lambda }_{i}}(t)}{\partial \lambda } \right)}^{2}}Var(\hat{\lambda }) \\
&  & +2\frac{\partial {{\lambda }_{i}}(t)}{\partial \beta }\frac{\partial {{\lambda }_{i}}(t)}{\partial \lambda }Cov(\hat{\beta },\hat{\lambda })+{{\left( \frac{\partial {{\lambda }_{i}}(t)}{\partial v(t)} \right)}^{2}}Var({{{\hat{v}}}_{t}}) 
\end{align}</math>
 
 
The expected value and variance of  <math>{{v}_{t}}</math>  are obtained from the Monte Carlo simulation with parameters  <math>\hat{\lambda },\hat{\beta },\hat{q}.</math>  Because of the simulation accuracy and the convergence problem in calculation of  <math>Var(\hat{\beta }),Var(\hat{\lambda })</math>  and  <math>Cov(\hat{\beta },\hat{\lambda }),</math>  <math>Var({{\lambda }_{i}}(t))</math>  can be a negative value at some time points. When this case happens, the bounds of instantaneous failure intensity are not provided.
Once the variance and the expected value of  <math>{{\lambda }_{i}}(t)</math>  are obtained, the bounds have been calculated by assuming  <math>{{\lambda }_{i}}(t)</math>  is lognormally distributed as:
 
 
<math>\frac{\ln {{\lambda }_{i}}(t)-\ln {{{\hat{\lambda }}}_{i}}(t)}{\sqrt{Var(\ln {{\lambda }_{i}}(t))}}\tilde{\ }N(0,1)</math>
 
 
The upper and lower bounds for a given confidence level  <math>\alpha </math>  can be calculated by:
 
 
<math>{{\lambda }_{i}}(t)={{\hat{\lambda }}_{i}}(t){{e}^{\pm {{z}_{a}}\sqrt{Var({{\lambda }_{i}}(t))}/{{{\hat{\lambda }}}_{i}}(t)}}</math>
 
 
where  <math>{{z}_{a}}</math>  is the standard normal distribution.
If  <math>{{\lambda }_{i}}(t)</math>  is assumed to be normally distributed, the bounds can be calculated by:
 
 
<math>{{\lambda }_{i}}{{(t)}_{U}}={{\hat{\lambda }}_{i}}(t)+{{z}_{a}}\sqrt{Var(N(t))}</math>
 
 
<math>{{\lambda }_{i}}{{(t)}_{L}}={{\hat{\lambda }}_{i}}(t)-{{z}_{a}}\sqrt{Var(N(t))}</math>
 
 
In Weibull++,  <math>{{\lambda }_{i}}{{(t)}_{U}}</math>  is set to the smaller of the two upper bounds obtained from Eqns.#(LognorIF)and (normIF).  <math>{{\lambda }_{i}}{{(t)}_{L}}</math>  is set to the largest of the two lower bounds obtained from Eqns.#(LognorIF) and (normIF). This combination method can prevent the out of range values of bounds when  <math>t</math>  values are small.
 
For a given time  <math>t</math> , the expected value of cumulative MTBF  <math>{{m}_{i}}(t)</math>  is:
 
 
<math>{{\hat{m}}_{i}}(t)=\frac{1}{{{{\hat{\lambda }}}_{i}}(t)}\text{  }</math>
 
 
The upper and lower bounds can be easily obtained from the corresponding bounds of  <math>{{\lambda }_{i}}(t):</math> 
 
<math>{{\hat{m}}_{i}}{{(t)}_{U}}=\frac{1}{{{{\hat{\lambda }}}_{i}}{{(t)}_{L}}}</math>
 
 
<math>{{\hat{m}}_{i}}{{(t)}_{L}}=\frac{1}{{{{\hat{\lambda }}}_{i}}{{(t)}_{U}}}</math>
 
====Bounds of Conditional Reliability====
Given mission start time  <math>{{t}_{0}}</math>  and mission time  <math>T</math> , the conditional reliability can be calculated by:
 
 
<math>R(T|{{t}_{0}})=\frac{R(T+{{v}_{0}})}{R({{v}_{0}})}={{e}^{-\lambda [{{({{v}_{0}}+T)}^{\beta }}-{{v}_{0}}]}}</math>
 
 
<math>{{v}_{0}}</math>  is the virtual age corresponding to time  <math>{{t}_{0}}</math> . The expected value and the variance of  <math>{{v}_{0}}</math>  are obtained from Monte Carlo simulation. The variance of the conditional reliability  <math>R(T|{{t}_{0}})</math>  is:
 
 
<math>\begin{align}
  & Var(R)= & {{\left( \frac{\partial R}{\partial \beta } \right)}^{2}}Var(\hat{\beta })+{{\left( \frac{\partial R}{\partial \lambda } \right)}^{2}}Var(\hat{\lambda }) \\
&  & +2\frac{\partial R}{\partial \beta }\frac{\partial R}{\partial \lambda }Cov(\hat{\beta },\hat{\lambda })+{{\left( \frac{\partial R}{\partial {{v}_{0}}} \right)}^{2}}Var({{{\hat{v}}}_{0}}) 
\end{align}</math>
 
 
Because of the simulation accuracy and the convergence problem in calculation of  <math>Var(\hat{\beta }),Var(\hat{\lambda })</math>  and  <math>Cov(\hat{\beta },\hat{\lambda }),</math>  <math>Var(R)</math>  can be a negative value at some time points. When this case happens, the bounds are not provided.
The bounds are based on:
 
 
<math>\log \text{it}(\hat{R}(T))\tilde{\ }N(0,1)</math>
 
 
<math>\log \text{it}(\hat{R}(T))=\ln \left\{ \frac{\hat{R}(T)}{1-\hat{R}(T)} \right\}</math>
 
 
The confidence bounds on reliability are given by:
 
 
<math>R=\frac{{\hat{R}}}{\hat{R}+(1-\hat{R}){{e}^{\pm \sqrt{Var(R)}/[\hat{R}(1-\hat{R})]}}}</math>
 
 
It will be compared with the bounds obtained from:
 
<math>R=\hat{R}{{e}^{\pm {{z}_{a}}\sqrt{Var(R)}/\hat{R}}}</math>
 
 
The smaller of the two upper bounds will be the final upper bound and the larger of the two lower bounds will be the final lower bound.
====Example 4====
The following table gives the failure times of the air-condition unit of an aircraft. The observation is ended by the time of the last failure [3].
 
 
<math>\begin{matrix}
  \text{50} & \text{329} & \text{811} & \text{991} & \text{1489}  \\
  \text{94} & \text{332} & \text{899} & \text{1013} & \text{1512}  \\
  \text{196} & \text{347} & \text{945} & \text{1152} & \text{1525}  \\
  \text{268} & \text{544} & \text{950} & \text{1362} & \text{1539}  \\
  \text{290} & \text{732} & \text{955} & \text{1459} & {}  \\
\end{matrix}</math>
 
 
1. Estimate the GRP model parameters using the Type I virtual age option.
 
2. Plot the failure number and instantaneous failure intensity vs. time with 90% two-sided confidence bounds.
 
3. Plot the conditional reliability vs. time with 90% two-sided confidence bounds. The mission start time is 40 and mission time is varying.
 
4. Using the QCP, calculate the expected failure number and expected instantaneous failure intensity by time 1800.
 
====Solution to Example 4====
Enter the data into a Parametric RDA Specialized Folio in Weibull++. Choose 3 under Parameters and Type I under Settings. Keep the default simulation settings.
 
1. The estimated parameters are  <math>\hat{\beta }=1.1976,</math>  <math>\hat{\lambda }=4.94E-03,</math>  <math>\hat{q}=0.1344</math> .
 
2. The failure number and instantaneous failure intensity are given in the following plots.
 
3. The conditional reliability is plotted below.
 
4. Using QCP, the failure number and instantaneous failure intensity are:

Latest revision as of 10:53, 25 June 2012

Redirect to:

Retrieved from "https://www.reliawiki.com/index.php?title=Non-Parametric_Recurrent_Events_Data_Analysis&oldid=25714"