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

From ReliaWiki
Jump to navigation Jump to search
 
(13 intermediate revisions by 4 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 Weibull-Bayesian is selected under '' Distribution'' and the <span class="texhtml">β</span> prior distribution is entered after clicking the ''Calculate'' button.
 
[[Image:chp6folio2data1.gif|thumb|center|300px| ]]
 
Suppose that the reliability at 3000hr is the metric of interest in this example. This reliability can be obtained using Eqn. (EQNREF MedRel ), 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 Eqn. (EQNREF Rpdf ). 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 Fig. 6-10 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 Eqn. (EQNREF WeibBayesPDF ) and plotted using Weibull++, as shown next:
 
[[Image:chp6plotofdata1.gif|thumb|center|400px| ]]
 
<br>
{{RS Copyright}} 
 
[[Category:Life_Data_Analysis_Reference]]

Latest revision as of 01:33, 13 August 2012