Simulation Based Bounds Example: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
(updated example)
Line 1: Line 1:
'''Comparing Parameter Estimation Methods Using Simulation Based Bounds'''
<noinclude>{{Banner Weibull Examples}}
''This example appears in the [[Introduction to Confidence Bounds]] article.''


The purpose of this example is to determine the best parameter estimation method for a sample of ten units following a Weibull distribution with <math>\beta =2</math> and <math>\eta =100</math> and with complete time-to-failure data for each unit (i.e. no censoring). The number of generated data sets is set to 1,000. The SimuMatic inputs are shown next.


[[Image: SimuMatic Parameters.png|thumb|center|500px]]
</noinclude>'''Comparing Parameter Estimation Methods Using Simulation Based Bounds'''


[[Image: SimuMatic Censoring.png|thumb|center|500px]]
The purpose of this example is to determine the best parameter estimation method for a sample of ten units with complete time-to-failure data for each unit (i.e., no censoring). The data set follows a Weibull distribution with <math>\beta =2</math> and <math>\eta =100</math> hr.


[[Image: SimuMatic Number of Data Sets.png|thumb|center|500px]]
The confidence bounds for the data set could be obtained by using '''Weibull++''''s SimuMatic utility. To obtain the results, use the following settings in SimuMatic.


:# On the Main tab, choose the '''2P-Weibull''' distribution and enter the given parameters (i.e., <math>\beta =2</math> and <math>\eta =100</math> hr)
:# On the Censoring tab, select the '''No censoring''' option.
:# On the Settings tab, set the number of data sets to '''1,000''' and the number of data points to '''10'''.
:# On the Analysis tab, choose the '''RRX''' analysis method and set the confidence bounds to '''90'''.


The parameters are estimated using RRX, RRY and MLE. The plotted results generated by SimuMatic are shown next.  
The following plot shows the simulation-based confidence bounds for the RRX parameter estimation method, as well as the expected variation due to sampling error.


:Using RRX:


[[Image:SimuMatic RRX.png|thumb|center|500px|'''RRX''']]  
[[Image:SimuMatic RRX.png|center|500px|'''RRX''']]  


<br>


:Using RRY:
Create another SimuMatic folio and generate a second data using the same settings, but this time, select the RRY analysis method on the Analysis tab. The following plot shows the result.


[[Image:SimuMatic RRY.png|thumb|center|500px|'''RRY''']]


<br>
[[Image:SimuMatic RRY.png|center|500px|'''RRY''']]


:Using MLE:
The following plot shows the results using the MLE analysis method.


<br> [[Image:SimuMatic MLE.png|thumb|center|500px|'''MLE''']]


<br> The results clearly demonstrate that the median RRX estimate provides the least deviation from the truth for this sample size and data type. However, the MLE outputs are grouped more closely together, as evidenced by the bounds. The previous figures also show the simulation-based bounds, as well as the expected variation due to sampling error.  
[[Image:SimuMatic MLE.png|center|500px|'''MLE''']]
 
 
The results clearly demonstrate that the median RRX estimate provides the least deviation from the truth for this sample size and data type. However, the MLE outputs are grouped more closely together, as evidenced by the bounds.


This experiment can be repeated in SimuMatic using multiple censoring schemes (including Type I and Type II right censoring as well as random censoring) with various distributions. Multiple experiments can be performed with this utility to evaluate assumptions about the appropriate parameter estimation method to use for data sets.
This experiment can be repeated in SimuMatic using multiple censoring schemes (including Type I and Type II right censoring as well as random censoring) with various distributions. Multiple experiments can be performed with this utility to evaluate assumptions about the appropriate parameter estimation method to use for data sets.

Revision as of 04:44, 6 August 2012

Weibull Examples Banner.png


New format available! This reference is now available in a new format that offers faster page load, improved display for calculations and images and more targeted search.

As of January 2024, this Reliawiki page will not continue to be updated. Please update all links and bookmarks to the latest references at Weibull examples and Weibull reference examples.




This example appears in the Introduction to Confidence Bounds article.


Comparing Parameter Estimation Methods Using Simulation Based Bounds

The purpose of this example is to determine the best parameter estimation method for a sample of ten units with complete time-to-failure data for each unit (i.e., no censoring). The data set follows a Weibull distribution with [math]\displaystyle{ \beta =2 }[/math] and [math]\displaystyle{ \eta =100 }[/math] hr.

The confidence bounds for the data set could be obtained by using Weibull++'s SimuMatic utility. To obtain the results, use the following settings in SimuMatic.

  1. On the Main tab, choose the 2P-Weibull distribution and enter the given parameters (i.e., [math]\displaystyle{ \beta =2 }[/math] and [math]\displaystyle{ \eta =100 }[/math] hr)
  2. On the Censoring tab, select the No censoring option.
  3. On the Settings tab, set the number of data sets to 1,000 and the number of data points to 10.
  4. On the Analysis tab, choose the RRX analysis method and set the confidence bounds to 90.

The following plot shows the simulation-based confidence bounds for the RRX parameter estimation method, as well as the expected variation due to sampling error.


RRX


Create another SimuMatic folio and generate a second data using the same settings, but this time, select the RRY analysis method on the Analysis tab. The following plot shows the result.


RRY


The following plot shows the results using the MLE analysis method.


MLE


The results clearly demonstrate that the median RRX estimate provides the least deviation from the truth for this sample size and data type. However, the MLE outputs are grouped more closely together, as evidenced by the bounds.

This experiment can be repeated in SimuMatic using multiple censoring schemes (including Type I and Type II right censoring as well as random censoring) with various distributions. Multiple experiments can be performed with this utility to evaluate assumptions about the appropriate parameter estimation method to use for data sets.