Template:Weibull parameters maximum likelihood estimation: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
(Redirected page to The Weibull Distribution)
 
(5 intermediate revisions by one other user not shown)
Line 1: Line 1:
=== Maximum Likelihood Estimation ===
#REDIRECT [[The Weibull Distribution]]
 
As outlined in Chapter 3, maximum likelihood estimation works by developing a likelihood function based on the available data and finding the values of the parameter estimates that maximize the likelihood function.  This can be achieved by using iterative methods to determine the parameter estimate values that maximize the likelihood function, but this can be rather difficult and time-consuming, particularly when dealing with the three-parameter distribution.  Another method of finding the parameter estimates involves taking the partial derivatives of the likelihood function with respect to the parameters, setting the resulting equations equal to zero and solving simultaneously to determine the values of the parameter estimates. ( Note that MLE asymptotic properties do not hold when estimating <span class="texhtml">γ</span> using MLE [27].) The log-likelihood functions and associated partial derivatives used to determine maximum likelihood estimates for the Weibull distribution are covered in Appendix C.
 
 
'''Example 5:'''
 
{{Example: Weibull MLE}}
 
'''Maximum Likelihood Estimation Example'''
Repeat Example 1 using maximum likelihood estimation.
 
 
'''Solution'''
 
In this case, we have non-grouped data with no suspensions or intervals, i.e. complete data. The equations for the partial derivatives of the log-likelihood function are derived in Appendix C and given next:
::<math> \frac{\partial \Lambda }{\partial \beta }=\frac{6}{\beta } +\sum_{i=1}^{6}\ln \left( \frac{T_{i}}{\eta }\right) -\sum_{i=1}^{6}\left( \frac{T_{i}}{\eta }\right) ^{\beta }\ln \left( \frac{T_{i}}{\eta }\right) =0
</math>
 
:and:
 
::<math> \frac{\partial \Lambda }{\partial \eta }=\frac{-\beta }{\eta }\cdot 6+\frac{ \beta }{\eta }\sum\limits_{i=1}^{6}\left( \frac{T_{i}}{\eta }\right) ^{\beta }=0 </math>
 
Solving the above equations simultaneously we get:
 
::<math> \hat{\beta }=1.933,</math> <math>\hat{\eta }=73.526 </math>
 
<br>
The variance/covariance matrix is found to be,
 
::<math> \left[ \begin{array}{ccc} \hat{Var}\left( \hat{\beta }\right) =0.4211 & \hat{Cov}( \hat{\beta },\hat{\eta })=3.272  \\
 
\hat{Cov}(\hat{\beta },\hat{\eta })=3.272 & \hat{Var} \left( \hat{\eta }\right) =266.646 \end{array} \right] </math>
 
The results and the associated graph using Weibull++ (MLE) are shown next.
 
[[Image:weibullfolio16plot.png|thumb|center|400px| ]]
 
You can view the variance/covariance matrix directly by clicking the ''Quick Calculation Pad ''(QCP) icon
 
[[Image:qcpicon.gif|thumb|center|400px| ]]
 
 
[[Image:qcpfolio16.png|thumb|center|400px| ]]
 
 
<br> Note that the decimal accuracy displayed and used is based on your individual User Setup.

Latest revision as of 09:08, 3 August 2012