Template:Cb for the eyring-log

Approximate Confidence Bounds for the Eyring-Lognormal

Bounds on the Parameters

The lower and upper bounds on [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math] are estimated from:

[math]\displaystyle{ \begin{align} & {{A}_{U}}= & \widehat{A}+{{K}_{\alpha }}\sqrt{Var(\widehat{A})}\text{ (Upper bound)} \\ & {{A}_{L}}= & \widehat{A}-{{K}_{\alpha }}\sqrt{Var(\widehat{A})}\text{ (Lower bound)} \end{align} }[/math]

and:

[math]\displaystyle{ \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]

Since the standard deviation, [math]\displaystyle{ {{\widehat{\sigma }}_{{T}',}} }[/math] is a positive parameter, [math]\displaystyle{ \ln ({{\widehat{\sigma }}_{{{T}'}}}) }[/math] is treated as normally distributed, and the bounds are estimated from:

[math]\displaystyle{ \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]

The variances and covariances of [math]\displaystyle{ A, }[/math] [math]\displaystyle{ B, }[/math] and [math]\displaystyle{ {{\sigma }_{{{T}'}}} }[/math] are estimated from the local Fisher matrix (evaluated at [math]\displaystyle{ \widehat{A}, }[/math] [math]\displaystyle{ \widehat{B} }[/math] , [math]\displaystyle{ {{\widehat{\sigma }}_{{{T}'}}}) }[/math] as follows:

[math]\displaystyle{ \left( \begin{matrix} Var\left( {{\widehat{\sigma }}_{{{T}'}}} \right) & Cov\left( \widehat{A},{{\widehat{\sigma }}_{{{T}'}}} \right) & Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right) \\ Cov\left( {{\widehat{\sigma }}_{{{T}'}}},\widehat{A} \right) & Var\left( \widehat{A} \right) & Cov\left( \widehat{A},\widehat{B} \right) \\ Cov\left( {{\widehat{\sigma }}_{{{T}'}}},\widehat{B} \right) & Cov\left( \widehat{B},\widehat{A} \right) & Var\left( \widehat{B} \right) \\ \end{matrix} \right)={{[F]}^{-1}} }[/math]

where:

[math]\displaystyle{ F=\left( \begin{matrix} -\tfrac{{{\partial }^{2}}\Lambda }{\partial \sigma _{{{T}'}}^{2}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\sigma }_{{{T}'}}}\partial A} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{\sigma }_{{{T}'}}}\partial B} \\ -\tfrac{{{\partial }^{2}}\Lambda }{\partial A\partial {{\sigma }_{{{T}'}}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{A}^{2}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial A\partial B} \\ -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial {{\sigma }_{{{T}'}}}} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial B\partial A} & -\tfrac{{{\partial }^{2}}\Lambda }{\partial {{B}^{2}}} \\ \end{matrix} \right) }[/math]

Bounds on Reliability

The reliability of the lognormal distribution is given by:

[math]\displaystyle{ R({T}',V;A,B,{{\sigma }_{{{T}'}}})=\mathop{}_{{{T}'}}^{\infty }\frac{1}{{{\widehat{\sigma }}_{{{T}'}}}\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{\left( \tfrac{t+\ln (V)+\widehat{A}-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}} \right)}^{2}}}}dt }[/math]

Let [math]\displaystyle{ \widehat{z}(t,V;A,B,{{\sigma }_{T}})=\tfrac{t+\ln (V)+\widehat{A}-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}}, }[/math] then [math]\displaystyle{ \tfrac{d\widehat{z}}{dt}=\tfrac{1}{{{\widehat{\sigma }}_{{{T}'}}}}. }[/math]
For [math]\displaystyle{ t={T}' }[/math] , [math]\displaystyle{ \widehat{z}=\tfrac{{T}'+\ln (V)+\widehat{A}-\tfrac{\widehat{B}}{V}}{{{\widehat{\sigma }}_{{{T}'}}}} }[/math] , and for [math]\displaystyle{ t=\infty , }[/math] [math]\displaystyle{ \widehat{z}=\infty . }[/math] The above equation then becomes:

[math]\displaystyle{ R(\widehat{z})=\mathop{}_{\widehat{z}({T}',V)}^{\infty }\frac{1}{\sqrt{2\pi }}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz }[/math]

The bounds on [math]\displaystyle{ z }[/math] are estimated from:

[math]\displaystyle{ \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]

where:

[math]\displaystyle{ \begin{align} & Var(\widehat{z})= & \left( \frac{\partial \widehat{z}}{\partial A} \right)_{\widehat{A}}^{2}Var(\widehat{A})+\left( \frac{\partial \widehat{z}}{\partial B} \right)_{\widehat{B}}^{2}Var(\widehat{B})+\left( \frac{\partial \widehat{z}}{\partial {{\sigma }_{{{T}'}}}} \right)_{{{\widehat{\sigma }}_{{{T}'}}}}^{2}Var({{\widehat{\sigma }}_{T}}) \\ & & +2{{\left( \frac{\partial \widehat{z}}{\partial A} \right)}_{\widehat{A}}}{{\left( \frac{\partial \widehat{z}}{\partial B} \right)}_{\widehat{B}}}Cov\left( \widehat{A},\widehat{B} \right) \\ & & +2{{\left( \frac{\partial \widehat{z}}{\partial A} \right)}_{\widehat{A}}}{{\left( \frac{\partial \widehat{z}}{\partial {{\sigma }_{{{T}'}}}} \right)}_{{{\widehat{\sigma }}_{{{T}'}}}}}Cov\left( \widehat{A},{{\widehat{\sigma }}_{T}} \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) \end{align} }[/math]

or:

[math]\displaystyle{ \begin{align} & Var(\widehat{z})= & \frac{1}{\widehat{\sigma }_{{{T}'}}^{2}}[Var(\widehat{A})+\frac{1}{{{V}^{2}}}Var(\widehat{B})+{{\widehat{z}}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\ & & -\frac{2}{V}Cov\left( \widehat{A},\widehat{B} \right)-2\widehat{z}Cov\left( \widehat{A},{{\widehat{\sigma }}_{{{T}'}}} \right)+\frac{2\widehat{z}}{V}Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right)] \end{align} }[/math]

The upper and lower bounds on reliability are:

[math]\displaystyle{ \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

The bounds around time for a given lognormal percentile (unreliability) are estimated by first solving the reliability equation with respect to time as follows:

[math]\displaystyle{ {T}'(V;\widehat{A},\widehat{B},{{\widehat{\sigma }}_{{{T}'}}})=-\ln (V)-\widehat{A}+\frac{\widehat{B}}{V}+z\cdot {{\widehat{\sigma }}_{{{T}'}}} }[/math]

where:

[math]\displaystyle{ \begin{align} & {T}'(V;\widehat{A},\widehat{B},{{\widehat{\sigma }}_{{{T}'}}})= & \ln (T) \\ & z= & {{\Phi }^{-1}}\left[ F({T}') \right] \end{align} }[/math]

and:

[math]\displaystyle{ \Phi (z)=\frac{1}{\sqrt{2\pi }}\mathop{}_{-\infty }^{z({T}')}{{e}^{-\tfrac{1}{2}{{z}^{2}}}}dz }[/math]

The next step is to calculate the variance of [math]\displaystyle{ {T}'(V;\widehat{A},\widehat{B},{{\widehat{\sigma }}_{{{T}'}}}): }[/math]

[math]\displaystyle{ \begin{align} & Var({T}')= & {{\left( \frac{\partial {T}'}{\partial A} \right)}^{2}}Var(\widehat{A})+{{\left( \frac{\partial {T}'}{\partial B} \right)}^{2}}Var(\widehat{B})+{{\left( \frac{\partial {T}'}{\partial {{\sigma }_{{{T}'}}}} \right)}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\ & & +2\left( \frac{\partial {T}'}{\partial A} \right)\left( \frac{\partial {T}'}{\partial B} \right)Cov\left( \widehat{A},\widehat{B} \right) \\ & & +2\left( \frac{\partial {T}'}{\partial A} \right)\left( \frac{\partial {T}'}{\partial {{\sigma }_{{{T}'}}}} \right)Cov\left( \widehat{A},{{\widehat{\sigma }}_{{{T}'}}} \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) \end{align} }[/math]

or:

[math]\displaystyle{ \begin{align} & Var({T}')= & Var(\widehat{A})+\frac{1}{V}Var(\widehat{B})+{{\widehat{z}}^{2}}Var({{\widehat{\sigma }}_{{{T}'}}}) \\ & & -\frac{2}{V}Cov\left( \widehat{A},\widehat{B} \right) \\ & & -2\widehat{z}Cov\left( \widehat{A},{{\widehat{\sigma }}_{{{T}'}}} \right) \\ & & +\frac{2\widehat{z}}{V}Cov\left( \widehat{B},{{\widehat{\sigma }}_{{{T}'}}} \right) \end{align} }[/math]

The upper and lower bounds are then found by:

[math]\displaystyle{ \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]

Solving for [math]\displaystyle{ {{T}_{U}} }[/math] and [math]\displaystyle{ {{T}_{L}} }[/math] yields:

[math]\displaystyle{ \begin{align} & {{T}_{U}}= & {{e}^{T_{U}^{\prime }}}\text{ (Upper bound)} \\ & {{T}_{L}}= & {{e}^{T_{L}^{\prime }}}\text{ (Lower bound)} \end{align} }[/math]