Template:Example: Wei-Bayesian Log-normal Prior: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
 
(6 intermediate revisions by 3 users not shown)
Line 1: Line 1:
'''Weibull-Bayesian with Lognormal Prior Example'''
#REDIRECT [[Bayesian-Weibull_Lognormal_Prior_Example]]
 
A manufacturer has tested prototypes of a modified product. The test was terminated at 2000 hours, with only two failures observed from a sample size of eighteen.
 
{|  border=1 cellspacing=1 align="center"
|-
!Number of State
!State of F or S
!State End Time
|-
| 1 || F || 1180
|-
| 1 || F || 1842 
|-
| 16 || S || 2000 
|}
 
Because of the lack of failure data in the prototype testing, the manufacturer decided to use information gathered from prior tests on this product to increase the confidence in the results of the prototype testing. This decision was made because failure analysis indicated that the failure mode of these two failures is the same as the one observed in previous tests. In other words, it is expected that the shape of the distribution hasn't changed, but hopefully the scale has, indicating longer life. The two-parameter Weibull distribution have been used to model all prior tests results. The list of the estimated <span class="texhtml">β</span> parameter is as follows:
 
{| border=1 cellspacing=1 align="center"
|-
!Betas Obtained for Similar Mode
|-
| 1.7
|-
| 2.1
|-
| 2.4
|-
|3.1
|-
|3.5
|}
 
'''Solution'''
 
First, in order to fit the data to a Weibull-Bayesian model, a prior distribution for <span class="texhtml">β</span> needs to be determined. Based on the prior tests' <span class="texhtml">β</span> values, the prior distribution for <span class="texhtml">β</span> was found to be a lognormal distribution with <span class="texhtml">μ = 0.9064</span>, <span class="texhtml">σ = 0.3325</span> (obtained by entering the <span class="texhtml">β</span> values into a Weibull++ ''Standard Folio'' and analyzing it based on the RRX analysis method.)
 
the test data is entered into a ''Standard Folio'', the '''B-W Lognormal Prior''' is selected under '' Distribution'' and the <span class="texhtml">β</span> prior distribution is entered after clicking the ''Calculate'' button.
 
[[Image:Weibull Distribution Example 6 Data and Prior.png|thumb|center|500px| ]]
 
Suppose that the reliability at 3000hr is the metric of interest in this example. This reliability can be obtained using the equation for median reliability, resulting in the median value of the posterior  of the reliability at 3000hr. Using the ''QCP'', this value is calculated to be 76.97%. ( By default Weibull++ returns the median values of the posterior distribution. )
 
The posterior <math>pdf</math> of the reliability function at 3000hrs can be obtained using the equation at the begining of this section. In Figure 6-10 the posterior <math>pdf</math> of the reliability at 3000hrs is plotted, with the corresponding median value as well as the 10th percentile value shown. The 10th percentile constitutes the 90 Lower 1-Sided bound on the reliability at 3000hrs, which is calculated to be 50.77%.
 
[[Image:lda6.10.gif|thumb|center|300px| ]]
 
Notice that the <math>pdf</math> plotted in the above Figure is of the reliability at 3000hrs, and not the <math>pdf</math> of the times-to-failure data. The <math>pdf</math> of the times-to-failure data can be obtained using the equation at the begining of this section and plotted using Weibull++, as shown next:
 
[[Image:Weibull Distribution Example 6 pdf.png|thumb|center|400px| ]]
 
<br>
{{RS Copyright}} 
 
[[Category:Life_Data_Analysis_Reference]]

Latest revision as of 01:33, 13 August 2012