Template:Normal distribution estimation of the parameters: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
(Redirected page to The Normal Distribution)
 
(One intermediate revision by one other user not shown)
Line 1: Line 1:
==Estimation of the Parameters==
#REDIRECT [[The_Normal_Distribution]]
 
{{normal distribution probability plotting}}
 
{{normal distribution rank regression on Y}}
 
{{normal distribution rank regression on X}}
 
===Maximum Likelihood Estimation===
As it was 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. 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. The log-likelihood functions and associated partial derivatives used to determine maximum likelihood estimates for the normal distribution are covered in Appendix C.
 
===Special Note About Bias===
Estimators (i.e. parameter estimates) have properties such as unbiasedness, minimum variance, sufficiency, consistency, squared error constancy, efficiency and completeness [7][5]. Numerous books and papers deal with these properties and this coverage is beyond the scope of this reference.
 
However, we would like to briefly address one of these properties, unbiasedness. An estimator is said to be unbiased if the estimator  <math>\widehat{\theta }=d({{X}_{1,}}{{X}_{2,}}...,{{X}_{n)}}</math>  satisfies the condition  <math>E\left[ \widehat{\theta } \right]</math>  <math>=\theta </math>  for all  <math>\theta \in \Omega .</math>
Note that  <math>E\left[ X \right]</math>  denotes the expected value of X and is defined (for continuous distributions) by:
 
::<math>\begin{align}
  E\left[ X \right]=  \int_{\varpi }x\cdot f(x)dx \\
  X\in  & \varpi . 
\end{align}</math>
 
 
It can be shown [7][5] that the MLE estimator for the mean of the normal (and lognormal) distribution does satisfy the unbiasedness criteria, or  <math>E\left[ \widehat{\mu } \right]</math>  <math>=\mu .</math>  The same is not true for the estimate of the variance  <math>\hat{\sigma }_{T}^{2}</math> . The maximum likelihood estimate for the variance for the normal distribution is given by:
 
::<math>\hat{\sigma }_{T}^{2}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}}</math>
 
with a standard deviation of:
 
::<math>{{\hat{\sigma }}_{T}}=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}}}</math>
 
 
These estimates, however, have been shown to be biased. It can be shown [7][5] that the unbiased estimate of the variance and standard deviation for complete data is given by:
 
 
::<math>\begin{align}
  \hat{\sigma }_{T}^{2}= & \left[ \frac{N}{N-1} \right]\cdot \left[ \frac{1}{N}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}} \right]=\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}} \\
  {{{\hat{\sigma }}}_{T}}= & \sqrt{\left[ \frac{N}{N-1} \right]\cdot \left[ \frac{1}{N}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}} \right]} \\
  = & \sqrt{\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathop \sum }}\,{{({{T}_{i}}-\bar{T})}^{2}}} 
\end{align}</math>
 
 
Note that for larger values of  <math>N</math> ,  <math>\sqrt{\left[ N/(N-1) \right]}</math>  tends to 1.
 
Weibull++ by default returns the standard deviation as defined by Eqn. (NormSt2).
The Use Unbiased Std on Normal Data option in the User Setup under the Calculations tab allows biasing to be considered when estimating the parameters.
 
When this option is selected, Weibull++ returns the standard deviation as defined by Eqn. (NormSt2). This is only true for complete data sets. For all other data types, Weibull++ by default returns the standard deviation as defined by Eqn. (normbias2) regardless of the selection status of this option.  The next figure shows this setting in Weibull++.
 
 
<math></math>
[[Image:ldachp8fig3.gif|thumb|center|300px| ]]

Latest revision as of 09:06, 3 August 2012

Redirect to:

Retrieved from "https://www.reliawiki.com/index.php?title=Template:Normal_distribution_estimation_of_the_parameters&oldid=29209"