Template:Arrhenius confidence bounds: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(Created page with '=Appendix 6.A: Arrhenius Confidence Bounds= <br> ==Approximate Confidence Bounds for the Arrhenius-Exponential== <br> There are different methods for computing confidence bounds.…')
 
 
(4 intermediate revisions by 2 users not shown)
Line 1: Line 1:
=Appendix 6.A: Arrhenius Confidence Bounds=
#REDIRECT [[Arrhenius_Relationship#Arrhenius_Confidence_Bounds]]
<br>
==Approximate Confidence Bounds for the Arrhenius-Exponential==
<br>
There are different methods for computing confidence bounds. ALTA utilizes confidence bounds that are based on the asymptotic theory for maximum likelihood estimates, most commonly referred to as the Fisher matrix bounds.
<br>
===Confidence Bounds on the Mean Life===
<br>
The Arrhenius-exponential distribution is given by Eqn. (arrhenius) by setting  <math>m=L(V)</math>  as shown in Eqn. (Arrean). The upper  <math>({{m}_{U}})</math>  and lower  <math>({{m}_{L}})</math>  bounds on the mean life are then estimated by:
 
<br>
::<math>\begin{align}
  & {{m}_{U}}= & \widehat{m}\cdot {{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{m})}}{\widehat{m}}}} \\
& {{m}_{L}}= & \widehat{m}\cdot {{e}^{-\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{m})}}{\widehat{m}}}} 
\end{align}</math>
 
<br>
where  <math>{{K}_{\alpha }}</math>  is defined by:
 
<br>
::<math>\alpha =\frac{1}{\sqrt{2\pi }}\mathop{}_{{{K}_{\alpha }}}^{\infty }{{e}^{-\tfrac{{{t}^{2}}}{2}}}dt=1-\Phi ({{K}_{\alpha }})</math>
 
<br>
If  <math>\delta </math>  is the confidence level (i.e., 95%=0.95), then  <math>\alpha =\tfrac{1-\delta }{2}</math>  for the two-sided bounds, and  <math>\alpha =1-\delta </math>  for the one-sided bounds. The variance of  <math>\widehat{m}</math>  is given by:
 
<br>
::<math>\begin{align}
  & Var(\widehat{m})= & {{\left( \frac{\partial m}{\partial C} \right)}^{2}}Var(\widehat{C})+{{\left( \frac{\partial m}{\partial B} \right)}^{2}}Var(\widehat{B}) \\
&  & +2\left( \frac{\partial m}{\partial C} \right)\left( \frac{\partial m}{\partial B} \right)Cov(\widehat{B},\widehat{C}) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>Var(\widehat{m})={{e}^{\tfrac{2\widehat{B}}{V}}}\left[ Var(\widehat{C})+\frac{{{\widehat{C}}^{2}}}{{{V}^{2}}}Var(\widehat{B})+\frac{2\widehat{C}}{V}Cov(\widehat{B},\widehat{C}) \right]</math>
 
<br>
The variances and covariance of  <math>B</math>  and  <math>C</math>  are estimated from the local Fisher matrix (evaluated at  <math>\widehat{B}</math> ,  <math>\widehat{C})</math>  as follows:
 
<br>
::<math>\left[ \begin{matrix}
  Var(\widehat{B}) & Cov(\widehat{B},\widehat{C})  \\
  Cov(\widehat{C},\widehat{B}) & Var(\widehat{C})  \\
\end{matrix} \right]={{\left[ \begin{matrix}
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{B}^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial C}  \\
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial C\partial B} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{C}^{2}}}  \\
\end{matrix} \right]}^{-1}}</math>
 
===Confidence Bounds on Reliability===
<br>
The bounds on reliability for any given time,  <math>T</math> , are estimated by:
<br>
::<math>\begin{align}
  & {{R}_{U}}(T)= & {{e}^{-\tfrac{T}{{{m}_{U}}}}} \\
& {{R}_{L}}(T)= & {{e}^{-\tfrac{T}{{{m}_{L}}}}} 
\end{align}</math>
 
<br>
where  <math>{{m}_{U}}</math>  and  <math>{{m}_{L}}</math>  are estimated using Eqns. (ArrhuUpper) and (ArrhuLower).
<br>
 
===Confidence Bounds on Time===
<br>
The bounds on time (ML estimate of time) for a given reliability are estimated by first solving the reliability function with respect to time:
 
<br>
::<math>\widehat{T}=-\widehat{m}\cdot \ln (R)</math>
 
<br>
The corresponding confidence bounds are then estimated from:
 
<br>
::<math>\begin{align}
  & {{T}_{U}}= & -{{m}_{U}}\cdot \ln (R) \\
& {{T}_{L}}= & -{{m}_{L}}\cdot \ln (R) 
\end{align}</math>
 
 
where  <math>{{m}_{U}}</math>  and  <math>{{m}_{L}}</math>  are estimated using Eqns. (ArrhuUpper) and (ArrhuLower).
 
==Approximate Confidence Bounds for the Arrhenius-Weibull:==
<br>
===Bounds on the Parameters===
<br>
From the asymptotically normal property of the maximum likelihood estimators, and since  <math>\widehat{\beta },</math>  and  <math>\widehat{C}</math>  are positive parameters,  <math>\ln (\widehat{\beta }),</math>  and  <math>\ln (\widehat{C})</math>  can then be treated as normally distributed. After performing this transformation, the bounds on the parameters can be estimated from:
 
<br>
::<math>\begin{align}
  & {{\beta }_{U}}= & \widehat{\beta }\cdot {{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\beta })}}{\widehat{\beta }}}} \\
& {{\beta }_{L}}= & \widehat{\beta }\cdot {{e}^{-\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{\beta })}}{\widehat{\beta }}}} 
\end{align}</math>
 
<br>
also:
 
<br>
::<math>\begin{align}
  & {{B}_{U}}= & \widehat{B}+{{K}_{\alpha }}\sqrt{Var(\widehat{B})} \\
& {{B}_{L}}= & \widehat{B}-{{K}_{\alpha }}\sqrt{Var(\widehat{B})} 
\end{align}</math>
 
<br>
and:
 
<br>
::<math>\begin{align}
  & {{C}_{U}}= & \widehat{C}\cdot {{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{C})}}{\widehat{C}}}} \\
& {{C}_{L}}= & \widehat{C}\cdot {{e}^{-\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{C})}}{\widehat{C}}}} 
\end{align}</math>
 
<br>
The variances and covariances of  <math>\beta ,</math>  <math>B,</math>  and  <math>C</math>  are estimated from the local Fisher matrix (evaluated at  <math>\widehat{\beta },</math>  <math>\widehat{B},</math>  <math>\widehat{C})</math> , as follows:
 
<br>
::<math>\left[ \begin{matrix}
  Var(\widehat{\beta }) & Cov(\widehat{\beta },\widehat{B}) & Cov(\widehat{\beta },\widehat{C})  \\
  Cov(\widehat{B},\widehat{\beta }) & Var(\widehat{B}) & Cov(\widehat{B},\widehat{C})  \\
  Cov(\widehat{C},\widehat{\beta }) & Cov(\widehat{C},\widehat{B}) & Var(\widehat{C})  \\
\end{matrix} \right]={{\left[ \begin{matrix}
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\beta }^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial \beta \partial B} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial \beta \partial C}  \\
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial \beta } & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{B}^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial C}  \\
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial C\partial \beta } & -\tfrac{{{\partial }^{2}}\Lambda }{\partial C\partial B} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{C}^{2}}}  \\
\end{matrix} \right]}^{-1}}</math>
 
===Confidence Bounds on Reliability===
 
<br>
The reliability function for the Arrhenius-Weibull model (ML estimate) is given by:
 
<br>
::<math>\widehat{R}(T,V)={{e}^{-{{\left( \tfrac{T}{\widehat{C}\cdot {{e}^{\tfrac{\widehat{B}}{V}}}} \right)}^{\widehat{\beta }}}}}</math>
 
<br>
or:
 
<br>
::<math>\widehat{R}(T)={{e}^{-{{e}^{\ln \left[ {{\left( \tfrac{T}{\widehat{C}\cdot {{e}^{\tfrac{\widehat{B}}{V}}}} \right)}^{\widehat{\beta }}} \right]}}}}</math>
 
<br>
Setting:
 
<br>
::<math>\widehat{u}=\ln \left[ {{\left( \frac{T}{\widehat{C}\cdot {{e}^{\tfrac{\widehat{B}}{V}}}} \right)}^{\widehat{\beta }}} \right]</math>
 
<br>
or:
 
<br>
::<math>\widehat{u}=\widehat{\beta }\left[ \ln (T)-\ln (\widehat{C})-\frac{\widehat{B}}{V} \right]</math>
 
<br>
The reliability function now becomes:
 
<br>
::<math>\widehat{R}(T,V)={{e}^{-{{e}^{\widehat{u}}}}}</math>
 
<br>
The next step is to find the upper and lower bounds on  <math>\widehat{u}\ \ :</math>
 
<br>
::<math>{{u}_{U}}=\widehat{u}+{{K}_{\alpha }}\sqrt{Var(\widehat{u})}</math>
 
<br>
::<math>{{u}_{L}}=\widehat{u}-{{K}_{\alpha }}\sqrt{Var(\widehat{u})}</math>
 
<br>
where:
 
<br>
::<math>\begin{align}
  & Var(\widehat{u})= & {{\left( \frac{\partial \widehat{u}}{\partial \beta } \right)}^{2}}Var(\widehat{\beta })+{{\left( \frac{\partial \widehat{u}}{\partial B} \right)}^{2}}Var(\widehat{B})+{{\left( \frac{\partial \widehat{u}}{\partial C} \right)}^{2}}Var(\widehat{C}) \\
&  & +2\left( \frac{\partial \widehat{u}}{\partial \beta } \right)\left( \frac{\partial \widehat{u}}{\partial B} \right)Cov(\widehat{\beta },\widehat{B})+2\left( \frac{\partial \widehat{u}}{\partial \beta } \right)\left( \frac{\partial \widehat{u}}{\partial C} \right)Cov(\widehat{\beta },\widehat{C}) \\
&  & +2\left( \frac{\partial \widehat{u}}{\partial B} \right)\left( \frac{\partial \widehat{u}}{\partial C} \right)Cov(\widehat{B},\widehat{C}) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>\begin{align}
  & Var(\widehat{u})= & {{\left( \frac{\widehat{u}}{\widehat{\beta }} \right)}^{2}}Var(\widehat{\beta })+{{\left( \frac{\widehat{\beta }}{V} \right)}^{2}}Var(\widehat{B})+{{\left( \frac{\widehat{\beta }}{\widehat{C}} \right)}^{2}}Var(\widehat{C}) \\
&  & -\frac{2\widehat{u}}{V}Cov(\widehat{\beta },\widehat{B})-\frac{2\widehat{u}}{\widehat{C}}Cov(\widehat{\beta },\widehat{C})+\frac{2{{\widehat{\beta }}^{2}}}{V\widehat{C}}Cov(\widehat{B},\widehat{C}) 
\end{align}</math>
 
<br>
The upper and lower bounds on reliability are:
 
<br>
::<math>\begin{align}
  & {{R}_{U}}(T,V)= & {{e}^{-{{e}^{\left( {{u}_{L}} \right)}}}} \\
& {{R}_{L}}(T,V)= & {{e}^{-{{e}^{\left( {{u}_{U}} \right)}}}} 
\end{align}</math>
 
 
where  <math>{{u}_{U}}</math>  and  <math>{{u}_{L}}</math>  are estimated from Eqns. (ArreibRupper) and (ArreibRlower).
<br>
 
===Confidence Bounds on Time===
<br>
The bounds on time for a given reliability are estimated by first solving the reliability function with respect to time:
 
<br>
::<math>\begin{align}
  & \ln (R)= & -{{\left( \frac{\widehat{T}}{\widehat{C}\cdot {{e}^{\tfrac{\widehat{B}}{V}}}} \right)}^{\widehat{\beta }}} \\
& \ln (-\ln (R))= & \widehat{\beta }\left( \ln \widehat{T}-\ln \widehat{C}-\frac{\widehat{B}}{V} \right) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>\widehat{u}=\frac{1}{\widehat{\beta }}\ln (-\ln (R))+\ln \widehat{C}+\frac{\widehat{B}}{V}</math>
<br>
 
where  <math>\widehat{u}=\ln \widehat{T}</math> .
 
<br>
The upper and lower bounds on  <math>u</math>  are estimated from:
 
 
::<math>{{u}_{U}}=\widehat{u}+{{K}_{\alpha }}\sqrt{Var(\widehat{u})}</math>
 
<br>
::<math>{{u}_{L}}=\widehat{u}-{{K}_{\alpha }}\sqrt{Var(\widehat{u})}</math>
 
<br>
where:
 
<br>
::<math>\begin{align}
  & Var(\widehat{u})= & {{\left( \frac{\partial \widehat{u}}{\partial \beta } \right)}^{2}}Var(\widehat{\beta })+{{\left( \frac{\partial \widehat{u}}{\partial B} \right)}^{2}}Var(\widehat{B})+{{\left( \frac{\partial \widehat{u}}{\partial C} \right)}^{2}}Var(\widehat{C}) \\
&  & +2\left( \frac{\partial \widehat{u}}{\partial \beta } \right)\left( \frac{\partial \widehat{u}}{\partial B} \right)Cov(\widehat{\beta },\widehat{B})+2\left( \frac{\partial \widehat{u}}{\partial \beta } \right)\left( \frac{\partial \widehat{u}}{\partial C} \right)Cov(\widehat{\beta },\widehat{C}) \\
&  & +2\left( \frac{\partial \widehat{u}}{\partial B} \right)\left( \frac{\partial \widehat{u}}{\partial C} \right)Cov(\widehat{B},\widehat{C}) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>\begin{align}
  & Var(\widehat{u})= & \frac{1}{{{\widehat{\beta }}^{4}}}{{\left[ \ln (-\ln (R)) \right]}^{2}}Var(\widehat{\beta })+\frac{1}{{{V}^{2}}}Var(\widehat{B})+\frac{1}{{{\widehat{C}}^{2}}}Var(\widehat{C}) \\
&  & -\frac{2\ln (-\ln (R))}{{{\widehat{\beta }}^{2}}V}Cov(\widehat{\beta },\widehat{B})-\frac{2\ln (-\ln (R))}{{{\widehat{\beta }}^{2}}\widehat{C}}Cov(\widehat{\beta },\widehat{C}) \\
&  & +\frac{2}{V\widehat{C}}Cov(\widehat{B},\widehat{C}) 
\end{align}</math>
 
<br>
The upper and lower bounds on time can then found by:
 
<br>
::<math>\begin{align}
  & {{T}_{U}}= & {{e}^{{{u}_{U}}}} \\
& {{T}_{L}}= & {{e}^{{{u}_{L}}}} 
\end{align}</math>
 
<br>
where  <math>{{u}_{U}}</math>  and  <math>{{u}_{L}}</math>  are estimated using Eqns. (ArreibTupper) and (ArreibTlower).
 
==Approximate Confidence Bounds for the Arrhenius-Lognormal==
<br>
 
===Bounds on the Parameters===
 
<br>
The lower and upper bounds on  <math>B</math>  are estimated from:
 
<br>
::<math>\begin{align}
  & {{B}_{U}}= & \widehat{B}+{{K}_{\alpha }}\sqrt{Var(\widehat{B})}\text{ (Upper bound)} \\
& {{B}_{L}}= & \widehat{B}-{{K}_{\alpha }}\sqrt{Var(\widehat{B})}\text{ (Lower bound)} 
\end{align}</math>
 
<br>
Since the standard deviation,  <math>{{\widehat{\sigma }}_{{{T}'}}}</math> , and the parameter  <math>C</math>  are positive parameters,  <math>\ln ({{\widehat{\sigma }}_{{{T}'}}})</math>  and  <math>\ln (C)</math>  are treated as normally distributed. The bounds are estimated from:
 
<br>
::<math>\begin{align}
  & {{C}_{U}}= & \widehat{C}\cdot {{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{C})}}{\widehat{C}}}}\text{ (Upper bound)} \\
& {{C}_{L}}= & \frac{\widehat{C}}{{{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var(\widehat{C})}}{\widehat{C}}}}}\text{ (Lower bound)} 
\end{align}</math>
 
<br>
and:
 
<br>
::<math>\begin{align}
  & {{\sigma }_{U}}= & {{\widehat{\sigma }}_{{{T}'}}}\cdot {{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var({{\widehat{\sigma }}_{{{T}'}}})}}{{{\widehat{\sigma }}_{{{T}'}}}}}}\text{ (Upper bound)} \\
& {{\sigma }_{L}}= & \frac{{{\widehat{\sigma }}_{{{T}'}}}}{{{e}^{\tfrac{{{K}_{\alpha }}\sqrt{Var({{\widehat{\sigma }}_{{{T}'}}})}}{{{\widehat{\sigma }}_{{{T}'}}}}}}}\text{ (Lower bound)} 
\end{align}</math>
 
<br>
The variances and covariances of  <math>B,</math>  <math>C,</math>  and  <math>{{\sigma }_{{{T}'}}}</math>  are estimated from the local Fisher matrix (evaluated at  <math>\widehat{B},</math>  <math>\widehat{C}</math> ,  <math>{{\widehat{\sigma }}_{{{T}'}}}),</math>  as follows:
<br>
 
::<math>\left[ \begin{matrix}
  Var\left( {{\widehat{\sigma }}_{{{T}'}}} \right) & Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right) & Cov\left( \widehat{C},{{\widehat{\sigma }}_{{{T}'}}} \right)  \\
  Cov\left( {{\widehat{\sigma }}_{{{T}'}}},\widehat{B} \right) & Var\left( \widehat{B} \right) & Cov\left( \widehat{B},\widehat{C} \right)  \\
  Cov\left( {{\widehat{\sigma }}_{{{T}'}}},\widehat{C} \right) & Cov\left( \widehat{C},\widehat{B} \right) & Var\left( \widehat{C} \right)  \\
\end{matrix} \right]=</math>
<br>
::<math>={{\left[ \begin{matrix}
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial \sigma _{{{T}'}}^{2}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\sigma }_{{{T}'}}}\partial B} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\sigma }_{{{T}'}}}\partial C}  \\
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial {{\sigma }_{{{T}'}}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{B}^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial C}  \\
  -\tfrac{{{\partial }^{2}}\Lambda }{\partial C\partial {{\sigma }_{{{T}'}}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial C\partial B} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{C}^{2}}}  \\
\end{matrix} \right]}^{-1}}</math>
<br>
 
===Bounds on Reliability===
<br>
The reliability of the lognormal distribution is:
 
<br>
::<math>R({T}',V;B,C,{{\sigma }_{{{T}'}}})=\mathop{}_{{{T}'}}^{\infty }\frac{1}{{{\widehat{\sigma }}_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{t-\ln (\widehat{C})-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}} \right)}^{2}}}}dt</math>
 
<br>
Let  <math>\widehat{z}(t,V;B,C,{{\sigma }_{T}})=\tfrac{t-\ln (\widehat{C})-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}},</math>  then  ..
For  <math>t={T}'</math> ,  <math>\widehat{z}=\tfrac{{T}'-\ln (\widehat{C})-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}}</math> , and for  <math>t=\infty ,</math>  <math>\widehat{z}=\infty .</math>  The above equation then becomes:
 
<br>
::<math>R(\widehat{z})=\mathop{}_{\widehat{z}({T}')}^{\infty }\frac{1}{\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz</math>
 
<br>
The bounds on  <math>z</math>  are estimated from:
 
<br>
::<math>\begin{align}
  & {{z}_{U}}= & \widehat{z}+{{K}_{\alpha }}\sqrt{Var(\widehat{z})} \\
& {{z}_{L}}= & \widehat{z}-{{K}_{\alpha }}\sqrt{Var(\widehat{z})} 
\end{align}</math>
 
<br>
where:
 
<br>
::<math>\begin{align}
  & Var(\widehat{z})= & \left( \frac{\partial \widehat{z}}{\partial B} \right)_{\widehat{B}}^{2}Var(\widehat{B})+\left( \frac{\partial \widehat{z}}{\partial C} \right)_{\widehat{C}}^{2}Var(\widehat{C})+\left( \frac{\partial \widehat{z}}{\partial {{\sigma }_{{{T}'}}}} \right)_{{{\widehat{\sigma }}_{{{T}'}}}}^{2}Var({{\widehat{\sigma }}_{T}}) \\
&  & +2{{\left( \frac{\partial \widehat{z}}{\partial B} \right)}_{\widehat{B}}}{{\left( \frac{\partial \widehat{z}}{\partial C} \right)}_{\widehat{C}}}Cov\left( \widehat{B},\widehat{C} \right) \\
&  & +2{{\left( \frac{\partial \widehat{z}}{\partial B} \right)}_{\widehat{B}}}{{\left( \frac{\partial \widehat{z}}{\partial {{\sigma }_{{{T}'}}}} \right)}_{{{\widehat{\sigma }}_{{{T}'}}}}}Cov\left( \widehat{B},{{\widehat{\sigma }}_{T}} \right) \\
&  & +2{{\left( \frac{\partial \widehat{z}}{\partial C} \right)}_{\widehat{C}}}{{\left( \frac{\partial \widehat{z}}{\partial {{\sigma }_{{{T}'}}}} \right)}_{{{\widehat{\sigma }}_{{{T}'}}}}}Cov\left( \widehat{C},{{\widehat{\sigma }}_{T}} \right) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>\begin{align}
  & Var(\widehat{z})= & \frac{1}{\widehat{\sigma }_{{{T}'}}^{2}}[\frac{1}{{{V}^{2}}}Var(\widehat{B})+\frac{1}{{{C}^{2}}}Var(\widehat{C})+{{\widehat{z}}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\
&  & +\frac{2}{C\cdot V}Cov\left( \widehat{B},\widehat{C} \right)+\frac{2\widehat{z}}{V}Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right)+\frac{2\widehat{z}}{C}Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right)] 
\end{align}</math>
 
<br>
The upper and lower bounds on reliability are:
 
<br>
::<math>\begin{align}
  & {{R}_{U}}= & \mathop{}_{{{z}_{L}}}^{\infty }\frac{1}{\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz\text{ (Upper bound)} \\
& {{R}_{L}}= & \mathop{}_{{{z}_{U}}}^{\infty }\frac{1}{\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz\text{ (Lower bound)} 
\end{align}</math>
 
 
===Confidence Bounds on Time===
<br>
The bounds around time, for a given lognormal percentile (unreliability), are estimated by first solving the reliability equation with respect to time, as follows:
 
<br>
::<math>{T}'(V;\widehat{B},\widehat{C},{{\widehat{\sigma }}_{{{T}'}}})=\ln (\widehat{C})+\frac{\widehat{B}}{V}+z\cdot {{\widehat{\sigma }}_{{{T}'}}}</math>
 
<br>
where:
 
<br>
::<math>\begin{align}
  & {T}'(V;\widehat{B},\widehat{C},{{\widehat{\sigma }}_{{{T}'}}})= & \ln (T) \\
& z= & {{\Phi }^{-1}}\left[ F({T}') \right]
\end{align}</math>
 
<br>
and:
 
<br>
::<math>\Phi (z)=\frac{1}{\sqrt{2\pi }}\mathop{}_{-\infty }^{z({T}')}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz</math>
 
<br>
The next step is to calculate the variance of  <math>{T}'(V;\widehat{B},\widehat{C},{{\widehat{\sigma }}_{{{T}'}}}):</math>
 
<br>
::<math>\begin{align}
  & Var({T}')= & {{\left( \frac{\partial {T}'}{\partial B} \right)}^{2}}Var(\widehat{B})+{{\left( \frac{\partial {T}'}{\partial C} \right)}^{2}}Var(\widehat{C})+{{\left( \frac{\partial {T}'}{\partial {{\sigma }_{{{T}'}}}} \right)}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\
&  & +2\left( \frac{\partial {T}'}{\partial B} \right)\left( \frac{\partial {T}'}{\partial C} \right)Cov\left( \widehat{B},\widehat{C} \right) \\
&  & +2\left( \frac{\partial {T}'}{\partial B} \right)\left( \frac{\partial {T}'}{\partial {{\sigma }_{{{T}'}}}} \right)Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right) \\
&  & +2\left( \frac{\partial {T}'}{\partial C} \right)\left( \frac{\partial {T}'}{\partial {{\sigma }_{{{T}'}}}} \right)Cov\left( \widehat{C},{{\widehat{\sigma }}_{{{T}'}}} \right) 
\end{align}</math>
 
<br>
or:
 
<br>
::<math>\begin{align}
  & Var({T}')= & \frac{1}{{{V}^{2}}}Var(\widehat{B})+\frac{1}{{{C}^{2}}}Var(\widehat{C})+{{\widehat{z}}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\
&  & +\frac{2}{B\cdot C}Cov\left( \widehat{B},\widehat{C} \right) \\
&  & +\frac{2\widehat{z}}{V}Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right) \\
&  & +\frac{2\widehat{z}}{C}Cov\left( \widehat{C},{{\widehat{\sigma }}_{{{T}'}}} \right) 
\end{align}</math>
 
 
<br>
The upper and lower bounds are then found by:
 
 
<br>
::<math>\begin{align}
  & T_{U}^{\prime }= & \ln {{T}_{U}}={T}'+{{K}_{\alpha }}\sqrt{Var({T}')} \\
& T_{L}^{\prime }= & \ln {{T}_{L}}={T}'-{{K}_{\alpha }}\sqrt{Var({T}')} 
\end{align}</math>
 
<br>
Solving for  <math>{{T}_{U}}</math>  and  <math>{{T}_{L}}</math>  yields:
 
<br>
::<math>\begin{align}
  & {{T}_{U}}= & {{e}^{T_{U}^{\prime }}}\text{ (Upper bound)} \\
& {{T}_{L}}= & {{e}^{T_{L}^{\prime }}}\text{ (Lower bound)} 
\end{align}</math>
 
 
 
 
 
{{RS Copyright}}

Latest revision as of 05:45, 16 August 2012

Redirect to:

Retrieved from "https://www.reliawiki.com/index.php?title=Template:Arrhenius_confidence_bounds&oldid=32278"