For large, complete data sets, both the LSE method and the MLE method provide consistent results. In reliability applications, data sets are typically small or moderate in size. Extensive simulation studies show that in small sample designs where there are only a few failures, the MLE method is better than the LSE method.1 Thus, the default estimation method in Minitab is MLE.
When there are only a few failures because the data are heavily censored, the MLE method uses the information in the entire data set, including the censored values. The LSE method ignores the information in the censored observations.1
Usually, the advantages of the MLE method outweigh the advantages of the LSE method. The LSE method is easier to calculate by hand and easier to program. The LSE method is also traditionally associated with the use of probability plots to assess goodness-of-fit. However, the LSE method can provide misleading results on a probability plot. Examples exist where the points on a Weibull probability plot that uses the LSE method fall along a line when the Weibull model is actually inappropriate.1
1. Genschel, U. and Meeker, W.Q. (2010). A Comparison of Maximum Likelihood and Median-Rank Regression for Weibull Estimation. Quality Engineering, 22(4): 236–255.
In earlier releases, Minitab provided calculated results for standard errors, confidence intervals, and tests for model parameters when using the LSE method. These calculated results were based on an ad-hoc method. However, there is no established, accepted statistical method for calculating standard errors for model parameters using the LSE method. Therefore, if you change the default method of estimation and select Least Squares (failure time(X) on rank(Y)), your output will not include calculated results for standard errors, confidence intervals, and tests for the model parameters. If you want to include confidence intervals and tests for model parameters in your results, you must use the MLE (default) method.
To change the parameter estimation method from MLE to LSE when using a parametric distribution analysis, a distribution ID plot, or a distribution overview plot, do the following:
Analysis | Button |
---|---|
Parametric Distribution Analysis | Estimate |
Distribution ID Plot | Options |
Distribution Overview Plot | Options |
If you use the least squares estimation method, estimates are calculated by fitting a regression line to the points in a probability plot. The line is formed by regressing time to failure or log (time to failure) (X) on the transformed percent (Y).
Because the percentiles of the distribution are based on the estimated distribution parameters, differences in the estimated parameters will cause differences in estimated percentiles.
When you estimate the parameters using the maximum likelihood estimation method, you can specify starting values for the algorithm and specify the maximum number of iterations.
Distribution | Parameters |
---|---|
Weibull | Enter shape and scale |
Exponential | Enter mean |
Other 2-parameter distributions | Enter location and scale |
2-parameter exponential | Enter scale and threshold |
3-parameter Weibull | Enter shape, scale, and threshold |
Other 3-parameter distributions | Enter location, scale, and threshold |
Minitab obtains maximum likelihood estimates through an iterative process. If the maximum number of iterations is obtained before convergence, the algorithm stops.
You can use either of the estimation methods in Parametric Distribution Analysis (Right Censoring) and Parametric Distribution Analysis (Arbitrary Censoring). However, instead of having Minitab estimate the parameters using one of these methods, you can also specify some parameters or all the parameters. If you choose to specify parameters, the calculated results—such as the percentiles—are based on the values of the parameters that you entered for the analysis.
You can specify some of the parameters for your distribution and have Minitab estimate the others from the data. Usually, you estimate some parameters to perform a Bayes Analysis when the data have few or no failures. See How to perform a reliability analysis with few or no failures for more details.
Distribution | Parameters that you can specify |
---|---|
Weibull | Shape |
3-parameter Weibull | Shape, threshold, or both |
Exponential | None |
2-parameter Exponential | Threshold |
Other distributions without a threshold | Scale |
Other distributions with a threshold | Scale, threshold, or both |
You always estimate the scale parameter for the Weibull distribution. For distributions that have a location parameter, you always estimate the location parameter.
You can specify all of the parameters instead of estimating them from the data. You can specify historical parameters to do things like compare the estimates you based on historical parameters to estimates based on the current data or see how the current data fit a probability plot based on the historical parameters.
Distribution | Parameters |
---|---|
Weibull | Enter shape and scale |
Exponential | Enter mean |
Other 2-parameter distributions | Enter location and scale |
2-parameter exponential | Enter scale and threshold |
3-parameter Weibull | Enter shape, scale, and threshold |
Other 3-parameter distributions | Enter location, scale, and threshold |
When you perform parametric distribution analysis, you can have Minitab assume common shape or scale parameters for the estimates.
Minitab then assumes common shape or scale parameters when calculating the estimates. For example, suppose that you have 2 (or more generally k>2) independent normally distributed samples with different means but the same variance. To estimate the mean for each sample Minitab uses a pooled estimate of the variance. This approach is generalized to other distributions as well. The specific result, however, depends on the estimation method that you have selected for the analysis.
For the maximum likelihood method, Minitab uses the log likelihood function. In this case, the log likelihood function of the model is the sum of the individual log likelihood functions, with the same shape parameter assumed in each individual log likelihood function. The resulting overall log likelihood function is maximized to obtain the scale parameters associated with each group and the common shape parameter. For more information, see the following reference: W. Nelson (1982). Applied Life Data Analysis, Chapter 12. John Wiley & Sons.
Minitab first calculates the y-coordinate and x-coordinate for each group (for details, see the "Plot points" and "Fitted line" topics in Methods and formulas for probability plot in Parametric Distribution Analysis (Right Censoring). Then, to obtain the LSE estimates, Minitab performs the following steps:
For log-location-scale distributions (for example, Weibull), the x-coordinates must be log-transformed. The groups should have the same slope, which is the inverse of the common shape parameter. The scale parameter of each group is obtained by exponentiation of the intercept for each group.