Template:Goodness-of-fit tests for rsa - Revision history

Lisa Hacker: Redirected page to RGA Models for Repairable Systems Analysis

2012-08-23T06:18:26Z

Redirected page to RGA Models for Repairable Systems Analysis

← Older revision		Revision as of 06:18, 23 August 2012
Line 1:		Line 1:
	#REDIRECT [[Repairable Systems Analysis]]		#REDIRECT [[RGA Models for Repairable Systems Analysis]]

Lisa Hacker: Redirected page to Repairable Systems Analysis

2012-08-23T06:17:49Z

Redirected page to Repairable Systems Analysis

← Older revision		Revision as of 06:17, 23 August 2012
Line 1:		Line 1:
	~~===Goodness-of-Fit Tests for~~ Repairable ~~System~~ Analysis~~===~~		#REDIRECT [[Repairable Systems Analysis]]
	~~<br>~~
	It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A parametric Cramér-von Mises goodness-of-fit test is used for the multiple system and repairable system Power Law model, as proposed by Crow in [17]~~. This goodness-of-fit test is appropriate whenever the start time for each system is 0 and the failure data is complete over the continuous interval <math>[0,{{T}_{q}}~~]</math> with no gaps in the data. The Chi-Squared test is a goodness-of-fit test that can be applied under more general circumstances. In addition, the Common Beta Hypothesis test also can be used to compare the intensity functions of the individual systems by comparing the <math>{{\beta }_{q}}</math> values of each system. Lastly, the Laplace Trend test checks for trends within the data. Due to their general applicatoin, the Common Beta Hypothesis test and the Laplace Trend test are both presented in Appendix B. The Cramér-von Mises and Chi-Squared goodness-of-fit tests are illustrated next.
	~~<br>~~
	~~<br>~~
	~~{{cramer-con mises test rsa}}~~

	~~{{chi-squared test rsa}}~~

Nicolette Young: /* Chi-Squared Test */

2012-01-10T22:03:21Z

Chi-Squared Test

← Older revision		Revision as of 22:03, 10 January 2012
Line 6:		Line 6:
	{{cramer-con mises test rsa}}		{{cramer-con mises test rsa}}

	~~====Chi-Squared Test====~~		{{chi-squared test rsa}}
	~~<br>~~
	~~The parametric Cramér-von Mises test described above requires that the starting time, <math>~~{{S}_{q}}</math> , be equal to 0 for each of the <math>K</math> systems. Although not as powerful as the Cramér-von Mises test, the Chi-Squared test can be applied regardless of the starting times. The expected number of failures for a system over its age <math>(a,b)</math> for the Chi-~~Squared~~ test is estimated by <math>\widehat{\lambda }{{b}^{\widehat{\beta }}}-\widehat{\lambda }{{a}^{\widehat{\beta }}}=\widehat{\theta }</math> , where <math>\widehat{\lambda }</math> and <math>\widehat{\beta }</math> are the maximum likelihood estimates.
	~~The computed <math>{{\chi }^{2}}</math> statistic is:~~

	~~::<math>{{\chi }^{2}}=\underset{j=1}{\overset{d}{\mathop \sum }}\,{{\frac{\left[ N(j)-\theta (j) \right]}{\widehat{\theta }(j)}}^{2}}</math>~~

	where <math>d</math> is the total number of intervals. The random variable <math>{{\chi }^{2}}</math> is approximately Chi-Square distributed with <math>df=d-2</math> degrees of freedom. There must be at least three intervals and the length of the intervals do not have to be equal. It is common practice to require that the expected number of failures for each interval, <math>\theta (j)</math> , be at least five. If <math>\chi _{0}^{2}>\chi _{\alpha /2,d-2}^{2}~~</math> or if <math>\chi _{0}^{2}<\chi _{1-(\alpha /2),d-2}^{2~~}~~</math> , reject the null hypothesis.~~

Nicolette Young: /* Cramér-von Mises Test */

2012-01-10T22:02:35Z

Cramér-von Mises Test

@@ Line 4: / Line 4: @@
 <br>
 <br>
-====Cramér-von Mises Test====
+{{cramer-con mises test rsa}}
 ====Chi-Squared Test====

Nicolette Young: Created page with '===Goodness-of-Fit Tests for Repairable System Analysis===
It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A p…'

2012-01-10T22:01:30Z

Created page with '===Goodness-of-Fit Tests for Repairable System Analysis=== It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A p…'

New page

===Goodness-of-Fit Tests for Repairable System Analysis===
 
It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A parametric Cramér-von Mises goodness-of-fit test is used for the multiple system and repairable system Power Law model, as proposed by Crow in [17]. This goodness-of-fit test is appropriate whenever the start time for each system is 0 and the failure data is complete over the continuous interval <math>[0,{{T}_{q}}]</math> with no gaps in the data. The Chi-Squared test is a goodness-of-fit test that can be applied under more general circumstances. In addition, the Common Beta Hypothesis test also can be used to compare the intensity functions of the individual systems by comparing the <math>{{\beta }_{q}}</math> values of each system. Lastly, the Laplace Trend test checks for trends within the data. Due to their general applicatoin, the Common Beta Hypothesis test and the Laplace Trend test are both presented in Appendix B. The Cramér-von Mises and Chi-Squared goodness-of-fit tests are illustrated next.
 
 
====Cramér-von Mises Test====
 
To illustrate the application of the Cramér-von Mises statistic for multiple system data, suppose that <math>K</math> like systems are under study and you wish to test the hypothesis <math>{{H}_{1}}</math> that their failure times follow a non-homogeneous Poisson process. Suppose information is available for the <math>{{q}^{th}}</math> system over the interval <math>[0,{{T}_{q}}]</math> , with successive failure times , <math>(q=1,2,\ldots ,\,K)</math> . The Cramér-von Mises test can be performed with the following steps:
 
 
Step 1: If <math>{{x}_{{{N}_{q}}q}}={{T}_{q}}</math> (failure terminated) let <math>{{M}_{q}}={{N}_{q}}-1</math> , and if <math>{{x}_{{{N}_{q}}q}}<T</math> (time terminated) let <math>{{M}_{q}}={{N}_{q}}</math> . Then:

::<math>M=\underset{q=1}{\overset{K}{\mathop \sum }}\,{{M}_{q}}</math>

Step 2: For each system divide each successive failure time by the corresponding end time <math>{{T}_{q}}</math> , <math>\,i=1,2,...,{{M}_{q}}.</math> Calculate the <math>M</math> values:

::<math>{{Y}_{iq}}=\frac{{{X}_{iq}}}{{{T}_{q}}},i=1,2,\ldots ,{{M}_{q}},\text{ }q=1,2,\ldots ,K</math>

Step 3: Next calculate <math>\overline{\beta }</math> , the unbiased estimate of <math>\beta </math> , from:

::<math>\overline{\beta }=\frac{M-1}{\underset{q=1}{\overset{K}{\mathop{\sum }}}\,\underset{i=1}{\overset{Mq}{\mathop{\sum }}}\,\ln \left( \tfrac{{{T}_{q}}}{{{X}_{i}}{{}_{q}}} \right)}</math>

Step 4: Treat the <math>{{Y}_{iq}}</math> values as one group and order them from smallest to largest. Name these ordered values <math>{{z}_{1}},\,{{z}_{2}},\ldots ,{{z}_{M}}</math> , such that <math>{{z}_{1}}<\ \ {{z}_{2}}<\ldots <{{z}_{M}}</math> .
 
 
Step 5: Calculate the parametric Cramér-von Mises statistic.

::<math>C_{M}^{2}=\frac{1}{12M}+\underset{j=1}{\overset{M}{\mathop \sum }}\,{{(Z_{j}^{\overline{\beta }}-\frac{2j-1}{2M})}^{2}}</math>

Critical values for the Cramér-von Mises test are presented in Table B.2 of Appendix B.
 
 
Step 6: If the calculated <math>C_{M}^{2}</math> is less than the critical value then accept the hypothesis that the failure times for the <math>K</math> systems follow the non-homogeneous Poisson process with intensity function <math>u(t)=\lambda \beta {{t}^{\beta -1}}</math> .
 
 
=====Example 2=====
 
For the data from Example 1, use the Cramér-von Mises test to examine the compatibility of the model at a significance level <math>\alpha =0.10</math>
 
 
''Solution''
 
Step 1:

::<math>\begin{align}
& {{X}_{9,1}}= & 1913.5<2000,\,\ {{M}_{1}}=9 \\
& {{X}_{11,2}}= & 1867<2000,\,\ {{M}_{2}}=11 \\
& {{X}_{14,3}}= & 1604.8<2000,\,\ {{M}_{3}}=14 \\
& M= & \underset{q=1}{\overset{3}{\mathop \sum }}\,{{M}_{q}}=34
\end{align}</math>

Step 2: Calculate <math>{{Y}_{iq}},</math> treat the <math>{{Y}_{iq}}</math> values as one group and order them from smallest to largest. Name these ordered values <math>{{z}_{1}},\,{{z}_{2}},\ldots ,{{z}_{M}}</math> .
 
 
Step 3: Calculate <math>\overline{\beta }=\tfrac{M-1}{\underset{q=1}{\overset{K}{\mathop{\sum }}}\,\underset{i=1}{\overset{Mq}{\mathop{\sum }}}\,\ln \left( \tfrac{{{T}_{q}}}{{{X}_{i}}{{}_{q}}} \right)}=0.4397</math>
 
 
Step 4: Calculate <math>C_{M}^{2}=\tfrac{1}{12M}+\underset{j=1}{\overset{M}{\mathop{\sum }}}\,{{(Z_{j}^{\overline{\beta }}-\tfrac{2j-1}{2M})}^{2}}=0.0611</math>
 
 
Step 5: Find the critical value (CV) from Table B.2 for <math>M=34</math> at a significance level <math>\alpha =0.10</math> . <math>CV=0.172</math> .
 
 
Step 6: Since <math>C_{M}^{2}<CV</math> , accept the hypothesis that the failure times for the <math>K=3</math> repairable systems follow the non-homogeneous Poisson process with intensity function <math>u(t)=\lambda \beta {{t}^{\beta -1}}</math> .
 
 

====Chi-Squared Test====
 
The parametric Cramér-von Mises test described above requires that the starting time, <math>{{S}_{q}}</math> , be equal to 0 for each of the <math>K</math> systems. Although not as powerful as the Cramér-von Mises test, the Chi-Squared test can be applied regardless of the starting times. The expected number of failures for a system over its age <math>(a,b)</math> for the Chi-Squared test is estimated by <math>\widehat{\lambda }{{b}^{\widehat{\beta }}}-\widehat{\lambda }{{a}^{\widehat{\beta }}}=\widehat{\theta }</math> , where <math>\widehat{\lambda }</math> and <math>\widehat{\beta }</math> are the maximum likelihood estimates.
The computed <math>{{\chi }^{2}}</math> statistic is:

::<math>{{\chi }^{2}}=\underset{j=1}{\overset{d}{\mathop \sum }}\,{{\frac{\left[ N(j)-\theta (j) \right]}{\widehat{\theta }(j)}}^{2}}</math>

where <math>d</math> is the total number of intervals. The random variable <math>{{\chi }^{2}}</math> is approximately Chi-Square distributed with <math>df=d-2</math> degrees of freedom. There must be at least three intervals and the length of the intervals do not have to be equal. It is common practice to require that the expected number of failures for each interval, <math>\theta (j)</math> , be at least five. If <math>\chi _{0}^{2}>\chi _{\alpha /2,d-2}^{2}</math> or if <math>\chi _{0}^{2}<\chi _{1-(\alpha /2),d-2}^{2}</math> , reject the null hypothesis.

Template:Goodness-of-fit tests for rsa - Revision history

Lisa Hacker: Redirected page to RGA Models for Repairable Systems Analysis

Lisa Hacker: Redirected page to Repairable Systems Analysis

Nicolette Young: /* Chi-Squared Test */

Nicolette Young: /* Cramér-von Mises Test */

Nicolette Young: Created page with '===Goodness-of-Fit Tests for Repairable System Analysis=== It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A p…'

Nicolette Young: Created page with '===Goodness-of-Fit Tests for Repairable System Analysis===
It is generally desirable to test the compatibility of a model and data by a statistical goodness-of-fit test. A p…'