09.12.2012 Views

The Weibull Distribution: A Handbook - Index of

The Weibull Distribution: A Handbook - Index of

The Weibull Distribution: A Handbook - Index of

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

17.1 Classical estimation approaches 559<br />

RIDER (1961) gave the following solution <strong>of</strong> the system (17.3a–c):<br />

�bj = d1 d2 − d3 ± � d2 3 − 6d1 d2 d3 − 4d2 1 d22 + 4d31 d3 + 4d3 3<br />

d2 1 − d2<br />

where for simplification he has set<br />

<strong>The</strong> estimator <strong>of</strong> the mixing proportion is<br />

© 2009 by Taylor & Francis Group, LLC<br />

dj = m ′ j Γ<br />

�<br />

1 + j<br />

�<br />

; j = 1,2,3.<br />

c<br />

, (17.4a)<br />

�p = d1 −�b2 �b1 −� . (17.4b)<br />

b2<br />

We now turn to the general case where all five parameters p,b1,b2,c1,c2 are unknown. In<br />

order to avoid the solution <strong>of</strong> a system <strong>of</strong> five dependent non–linear equations<br />

m ′ r = p b r �<br />

1 Γ 1 + r<br />

�<br />

+ (1 − p)b<br />

c1<br />

r �<br />

2 Γ 1 + r<br />

�<br />

; r = 1,... ,5.<br />

c2<br />

FALLS (1970) suggested to estimate p by the graphical approach <strong>of</strong> KAO (1959):<br />

1. Plot the sample CDF for the mixed data on WEIBULL–probability–paper (see<br />

Fig. 9/4) and visually fit a curve among these points (= WEIBULL plot).<br />

2. Starting at each end <strong>of</strong> the WEIBULL plot, draw two tangent lines and denote them<br />

by � p F1(x) and � (1 − p)F2(x), which are estimates <strong>of</strong> p F1(x) and (1 − p)F2(x),<br />

respectively.<br />

3. At the intersection <strong>of</strong> both tangent lines drop a vertical line on the percent scale which<br />

gives the estimate �p <strong>of</strong> p.<br />

<strong>The</strong> solution <strong>of</strong> m ′ r = µ ′ r (r = 1,... ,4) with p substituted by �p will not be unique. When<br />

confronted with more than one set <strong>of</strong> acceptable estimates �b1, �b2,�c2,�c2, FALLS (1970) suggested<br />

to adopt PEARSON’s procedure and chose the set which produces the closest agreement<br />

between m ′ 5 and the theoretical moment µ′ 5 when evaluated by the estimates.<br />

17.1.2 Estimation by maximum likelihood<br />

We will first present the solution <strong>of</strong> the special case <strong>of</strong> a two–fold mixture (Sect. 17.1.2.1)<br />

and then turn to the general case with more than two subpopulations (Sect. 17.1.2.2).<br />

17.1.2.1 <strong>The</strong> case <strong>of</strong> two subpopulations<br />

When we have postmortem data that are type–I censored, T being the censoring time,<br />

there are r1 failures in subpopulation 1 and r2 in subpopulation 2, found in a sample <strong>of</strong> n

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!