Parameter Estimation
The estimate of the parameters of the lognormal distribution can be found graphically on probability plotting paper or analytically using either least squares or maximum likelihood. Parameter estimation methods are presented in detail in Appendix B: Parameter Estimations.
Probability Plotting
One method of calculating the parameter of the lognormal distribution is by using probability plotting. To better illustrate this procedure, consider the following example.
Example 5
Let's assume six identical units are being reliability tested at the same application and operation stress levels. All of these units fail during the test after operating the following times (in hours), Ti: 144, 385, 747, 1,144, 1,576 and 2,616.
The steps for determining the parameters of the lognormal pdf representing the data, using probability plotting, are as follows:
Rank the times-to-failure in ascending order as shown next.
Obtain their median rank plotting positions. The times-to-failure, with their corresponding median ranks, are shown next:
On a lognormal probability paper, plot the times and their corresponding rank value. Figure 10 displays an example of a lognormal probability paper. The paper is simply a log-log paper. (The solution is given in Figure 11.)
Fig. 10: Sample Lognormal Probability Plotting Paper
Draw the best possible straight line that goes through the t = 0 and R(t) = 100% point and through these points (as shown in Figure 11).
At the Q(t) = 50% ordinate point, draw
a straight horizontal line until this line intersects the fitted straight
line. Draw a vertical line through this intersection until it crosses
the abscissa. The value at the intersection of the abscissa is the
estimate of the median. For this case, 
= 760 hr which means that 
= ln () = 6.633 (see Eqn.
( 15)).
Fig. 11: Probability plot for example 5
The standard deviation, 
,
can be found using the following equation:
Now any reliability value for any mission time t can be obtained. For example, the reliability for a mission of 200 hr, or any other time, can now be obtained either from the plot or analytically.
To obtain the value from the plot, draw a vertical line from the abscissa, at t = 200 hr, to the fitted line. Draw a horizontal line from this intersection to the ordinate and read Q(t). In this case, R(t = 200) = 1 - Q(t = 200) = 92%.This can also be obtained analytically, from the lognormal reliability function. However, standard normal tables (or the Quick Statistical Reference in ALTA) must be used.
MLE Parameter Estimation
The parameters of the lognormal distribution can also be estimated using maximum likelihood estimation (MLE). This general log-likelihood function is:
where:
and:
Fe is the number of groups of times-to-failure data points.
Ni is the number of times-to-failure in the ith time-to-failure data group.
is the mean of the natural logarithms of the times-to-failure (unknown
a priori, the first of two
parameters to be found).
is the standard deviation of the natural logarithms of the times-to-failure
(unknown a priori, the second
of two parameters to be found).
Ti is the time of the ith group of time-to-failure data.
S is the number of groups of suspension data points.
is the number of suspensions in ith
group of suspension data points.
is the time of the ith
suspension data group.
FIis the number of interval data groups.
is the number of intervals in the ith
group of data intervals.
is the beginning of the ith
interval.
is the ending of the ith
interval.
The solution will be found by solving for a pair of parameters (,
) so that 
= 0 and 
= 0, where:
and
Example 6
Using the same data as in the probability plotting example (Example 5) and assuming a lognormal distribution, estimate the parameters using the MLE method.
Solution
In this example we have non-grouped data, without suspensions. Thus, the partials reduce to,
Substituting the values of Ti and solving the above system simultaneously, we get:
The mean and standard deviation of the times-to-failure can be estimated using Eqns. ( 12) and ( 14),
and
See Also:
Lognormal Distribution
Go
to Weibull.com
Go to ReliaSoft.com
©1998-2010. ReliaSoft Corporation. ALL RIGHTS RESERVED.