Simulation Based Bounds Example: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 3: Line 3:
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 10,000. The SimuMatic inputs are shown next.  
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 10,000. The SimuMatic inputs are shown next.  


[[Image:Weibulltoolssimumatic.png|thumb|center|300px]]  
[[Image: SimuMatic Parameters.png|thumb|center|400px]]
 
[[Image: SimuMatic Censoring.png|thumb|center|400px]]
 
[[Image: SimuMatic Number of Data Sets.png|thumb|center|400px]]


{| align="center" border="1"
|-
| [[Image:Weibullsimumatic.png|thumb|center|300px]]
| [[Image:Weibullsimumatic2.png|thumb|center|300px]]
|-
| [[Image:Weibullsimumatic3.png|thumb|center|300px]]
| [[Image:Weibullsimumatic4.png|thumb|center|300px]]
|-
| [[Image:Weibullsimumatic5.png|thumb|center|300px]]
| [[Image:Weibullsimumatic6.png|thumb|center|300px]]
|}


The parameters are estimated using RRX, RRY and MLE. The plotted results generated by SimuMatic are shown next.  
The parameters are estimated using RRX, RRY and MLE. The plotted results generated by SimuMatic are shown next.  
Line 21: Line 14:
:Using RRX:
:Using RRX:


[[Image:Probabilityweibull.gif|thumb|center|300px]]  
[[Image:SimuMatic RRX.png|thumb|center|300px]]  


<br>  
<br>  

Revision as of 21:02, 1 March 2012

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 following a Weibull distribution with [math]\displaystyle{ \beta =2 }[/math] and [math]\displaystyle{ \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 10,000. The SimuMatic inputs are shown next.

SimuMatic Parameters.png
SimuMatic Censoring.png
SimuMatic Number of Data Sets.png


The parameters are estimated using RRX, RRY and MLE. The plotted results generated by SimuMatic are shown next.

Using RRX:
SimuMatic RRX.png


Using RRY:


Probabilityweibull2.gif


Using MLE:


Probabilityweibull3.gif


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.


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.