Template:Beta Binomial Confidence Bounds: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
 
Line 1: Line 1:
== Beta Binomial Confidence Bounds  ==
#REDIRECT [[Confidence_Bounds#Beta_Binomial_Confidence_Bounds]]
 
Another less mathematically intensive method of calculating confidence bounds involves a procedure similar to that used in calculating median ranks (see Chapter [[Parameter Estimation]]). This is a non-parametric approach to confidence interval calculations that involves the use of rank tables and is commonly known as beta-binomial bounds (BB). By non-parametric, we mean that no underlying distribution is assumed. (Parametric implies that an underlying distribution, with parameters, is assumed.) In other words, this method can be used for any distribution, without having to make adjustments in the underlying equations based on the assumed distribution. Recall from the discussion on the median ranks that we used the binomial equation to compute the ranks at the 50% confidence level (or median ranks) by solving the cumulative binomial distribution for <span class="texhtml">''Z''</span> (rank for the <span class="texhtml">''j''<sup>''t''''h'''</sup></span> failure):
 
::<math>P=\underset{k=j}{\overset{N}{\mathop \sum }}\,\left( \begin{matrix}
  N  \\
  k  \\
\end{matrix} \right){{Z}^{k}}{{\left( 1-Z \right)}^{N-k}}</math>
 
<br>where <span class="texhtml">''N''</span> is the sample size and <span class="texhtml">''j''</span> is the order number. <br>The median rank was obtained by solving the following equation for <span class="texhtml">''Z''</span>:
 
::<math>0.50=\underset{k=j}{\overset{N}{\mathop \sum }}\,\left( \begin{matrix}
  N  \\
  k  \\
\end{matrix} \right){{Z}^{k}}{{\left( 1-Z \right)}^{N-k}}</math>
 
<br>The same methodology can then be repeated by changing <span class="texhtml">''P''</span> for <span class="texhtml">0.50</span> <span class="texhtml">(50%)</span> to our desired confidence level. For <span class="texhtml">''P'' = 90%</span> one would formulate the equation as
 
::<math>0.90=\underset{k=j}{\overset{N}{\mathop \sum }}\,\left( \begin{matrix}
  N  \\
  k  \\
\end{matrix} \right){{Z}^{k}}{{\left( 1-Z \right)}^{N-k}}</math>
 
Keep in mind that one must be careful to select the appropriate values for <span class="texhtml">''P''</span> based on the type of confidence bounds desired. For example, if two-sided 80% confidence bounds are to be calculated, one must solve the equation twice (once with <span class="texhtml">''P'' = 0.1</span> and once with <span class="texhtml">''P'' = 0.9</span>) in order to place the bounds around 80% of the population. <br><br>Using this methodology, the appropriate ranks are obtained and plotted based on the desired confidence level. These points are then joined by a smooth curve to obtain the corresponding confidence bound.
 
In Weibull++, this non-parametric methodology is&nbsp;used only&nbsp;when plotting bounds on the mixed Weibull distribution. Full details on this methodology can be found in [[Appendix D: Weibull References|Kececioglu [20]]]. These binomial equations can again be transformed using the beta and F distributions, thus the name beta binomial confidence bounds.

Latest revision as of 23:34, 12 August 2012