Template:Weibull bayesian confidence bounds

Bounds on Parameters
Bayesian Bounds use non-informative prior distributions for both parameters. From Chapter 5, we know that if the prior distribution of η and β are independent, the posterior joint distribution of η and β can be written as:


 * $$ f(\eta ,\beta |Data)= \dfrac{L(Data|\eta ,\beta )\varphi (\eta )\varphi (\beta )}{\int_{0}^{\infty }\int_{0}^{\infty }L(Data|\eta ,\beta )\varphi (\eta )\varphi (\beta )d\eta d\beta } $$

The marginal distribution of η is:


 * $$ f(\eta |Data) =\int_{0}^{\infty }f(\eta ,\beta |Data)d\beta =

\dfrac{\int_{0}^{\infty }L(Data|\eta ,\beta )\varphi (\eta )\varphi (\beta )d\beta }{\int_{0}^{\infty }\int_{-\infty }^{\infty }L(Data|\eta ,\beta )\varphi (\eta )\varphi (\beta )d\eta d\beta } $$

where: $$ \varphi (\beta )=\frac{1}{\beta } $$ is the non-informative prior of β. $$ \varphi (\eta )=\frac{1}{\eta } $$ is the non-informative prior of η. Using these non-informative prior distributions, $$f(\eta|Data)$$ can be rewritten as:


 * $$ f(\eta |Data)=\dfrac{\int_{0}^{\infty }L(Data|\eta ,\beta )\frac{1}{\beta } \frac{1}{\eta }d\beta }{\int_{0}^{\infty }\int_{0}^{\infty }L(Data|\eta ,\beta )\frac{1}{\beta }\frac{1}{\eta }d\eta d\beta } $$

The one-sided upper bounds of η is:


 * $$ CL=P(\eta \leq \eta _{U})=\int_{0}^{\eta _{U}}f(\eta |Data)d\eta $$

The one-sided lower bounds of η is:


 * $$ 1-CL=P(\eta \leq \eta _{L})=\int_{0}^{\eta _{L}}f(\eta |Data)d\eta $$ The two-sided bounds of η is:


 * $$ CL=P(\eta _{L}\leq \eta \leq \eta _{U})=\int_{\eta _{L}}^{\eta _{U}}f(\eta |Data)d\eta $$ Same method is used to obtain the bounds of β.

Bounds on Time (Type 1)
From Chapter 5, we know that:


 * $$ CL=\Pr (T\leq T_{U})=\Pr (\eta \leq T_{U}\exp (-\frac{\ln (-\ln R)}{\beta })) $$ From the posterior distribution of η, we have:


 * $$ CL=\dfrac{\int\nolimits_{0}^{\infty }\int\nolimits_{0}^{T_{U}\exp (-\dfrac{ \ln (-\ln R)}{\beta })}L(\beta ,\eta )\frac{1}{\beta }\frac{1}{\eta }d\eta d\beta }{\int\nolimits_{0}^{\infty }\int\nolimits_{0}^{\infty }L(\beta ,\eta )\frac{1}{\beta }\frac{1}{\eta }d\eta d\beta } EQNREF BayesCLT $$

Eqn. (EQNREF BayesCLT ) is solved numerically for TU. The same method can be applied to calculate one sided lower bounds and two-sided bounds on time.

Bounds on Reliability (Type 2)

 * $$ CL=\Pr (R\leq R_{U})=\Pr (\eta \leq T\exp (-\frac{\ln (-\ln R_{U})}{\beta })) $$ From the posterior distribution of η, we have:


 * $$ CL=\dfrac{\int\nolimits_{0}^{\infty }\int\nolimits_{0}^{T\exp (-\dfrac{\ln (-\ln R_{U})}{\beta })}L(\beta ,\eta )\frac{1}{\beta }\frac{1}{\eta }d\eta d\beta }{\int\nolimits_{0}^{\infty }\int\nolimits_{0}^{\infty }L(\beta ,\eta )\frac{1}{\beta }\frac{1}{\eta }d\eta d\beta } EQNREF BayesCLR $$

Eqn. (EQNREF BayesCLR ) is solved numerically for RU. The same method can be used to calculate the one sided lower bounds and two-sided bounds on reliability.