Template:Example: Bayesian Test Design with Prior Information from Expert Opinion

From ReliaWiki
Revision as of 06:33, 12 July 2012 by Harry Guo (talk | contribs)
Jump to navigation Jump to search

Bayesian Test Design with Prior Information from Expert Opinion

Suppose you wanted to know the reliability of a system and you had the following prior knowledge of the system:

  • Lowest possible reliability: a = 0.8
  • Most likely reliability: b = 0.85
  • Highest possible reliability: c = 0.97


This information can be used to approximate the expected value and the variance of the prior system reliability.

[math]\displaystyle{ E\left(R_{0}\right)=\frac{a+4b+c}{6}=0.861667 }[/math]


[math]\displaystyle{ \lt math\gt Var({{R}_{0}})={{\left( \frac{c-a}{6} \right)}^{2}} }[/math]=0.000803 </math>


These approximations of the expected value and variance of the prior system reliability can then be used to estimate [math]\displaystyle{ \alpha_{0} }[/math] and [math]\displaystyle{ \beta_{0} }[/math] used in the beta distribution for the system reliability, as given next:

[math]\displaystyle{ \alpha\,\!_{0}=E\left(R_{0}\right)\left[\frac{E\left(R_{0}\right)-E^{2}\left(R_{0}\right)}{Var\left(R_{0}\right)}-1\right]=127.0794 }[/math]


[math]\displaystyle{ \beta\,\!_{0}=\left(1-E\left(R_{0}\right)\right)\left[\frac{E\left(R_{0}\right)-E^{2}\left(R_{0}\right)}{Var\left(R_{0}\right)}-1\right]=20.40153 }[/math]

With [math]\displaystyle{ \alpha_{0} }[/math] and [math]\displaystyle{ \beta_{0} }[/math] known, any single value of the four quantities system reliability R, confidence level CL, number of units n, or number of failures r can be calculated from the other three using the beta distribution function:

[math]\displaystyle{ 1-CL=\text{Beta}\left(R,\alpha,\beta\right)=\text{Beta}\left(R,n-r+\alpha_{0},r+\beta_{0}\right) }[/math]


Solve for System Reliability R

Given CL = 0.8, n = 20, and r = 1, using the above prior information to solve R.

First, we get the number of successes: [math]\displaystyle{ s = n – r = 19 }[/math]. Then the parameters in the posterior beta distribution for R are calculated as:

[math]\displaystyle{ \alpha\,\!=\alpha\,\!_{0}+s=146.0794 }[/math]


[math]\displaystyle{ \beta\,\!=\beta\,\!_{0}+r=21.40153 }[/math]

Finally, from this posterior distribution, the system reliability R at a confidence level of CL=0.8 is solved as:

[math]\displaystyle{ R=\text{BetaINV}\left(1-CL,\alpha\,\!,\beta\,\!\right)=0.838374 }[/math]


Solve for Confidence Level CL

Given R = 0.85, n = 20, and r = 1, using the above prior information on system reliability to solve for CL.

First, we get the number of successes: [math]\displaystyle{ s = n – r = 19 }[/math]. Then the parameters in the posterior beta distribution for R are calculated as:

[math]\displaystyle{ \alpha\,\!=\alpha\,\!_{0}+s=146.07943 }[/math]


[math]\displaystyle{ \beta\,\!=\beta\,\!_{0}+r=21.40153 }[/math]

Finally, from this posterior distribution, the corresponding confidence level for reliability R=0.85 is:

[math]\displaystyle{ CL=\text{Beta}\left(R,\alpha,\beta\right)=0.81011 }[/math]


Solve for Sample Size n

Given R = 0.9, CL = 0.8, and r = 1, using the above prior information on system reliability to solve the required sample size in the demonstration test.

Again, the above beta distribution equation for the system reliability can be utilized. The figure below shows the result from Weibull++.

Bayesian Prior Expert Solve for n.png

The results show that the required sample size is 103. Weibull++ always displays the sample size as an integer.