Simulation Based Bounds Example

From ReliaWiki
Revision as of 18:33, 28 February 2012 by Harry Guo (talk | contribs)
Jump to navigation Jump to search

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.

Weibulltoolssimumatic.png
Weibullsimumatic.png
Weibullsimumatic2.png
Weibullsimumatic3.png
Weibullsimumatic4.png
Weibullsimumatic5.png
Weibullsimumatic6.png

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

Using RRX:
Probabilityweibull.gif


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.