Image Anal Stereol 2006;25:145-154 Original Research Paper SIZE AND SHAPE FACTOR EXTREMES OF SPHEROIDS Daniel Hlubinka Charles University in Prague, Faculty of Mathematics and Physics, Department of Probability and Statistics, Sokolovska´ 83, 186 75 Praha 8, Czech Republic e-mail: daniel.hlubinka@mff.cuni.cz (Accepted October 31, 2006) ABSTRACT In the paper we consider random prolate (oblate) spheroids and their random profiles. The limiting distribution of the extremal characteristics of the spheroids is related to the limiting distribution of the corresponding extremal characteristics of the profiles. The difference between the analysis of the prolate and oblate spheroids is discussed. We propose the possible application of the theoretical results. Keywords: maximum domain of attraction, shape factor of spheroids, size of spheroids, stereology of extremes. INTRODUCTION Microscopic particles (inclusions) in a material may cause a serious damage to that material. The study of the extremes of certain characteristics of such particles is important since material damage is related to the extremal rather than to the mean characteristics of the microstructure. The principal problem is typically that the particles are not observed directly. Only the profiles of the particles appearing on a random planar section of the material can be observed. The basic model of the particles is random balls characterized by their size. The problem of estimating the size of random balls observed through their profiles has been well known for a long time as Wicksell’s corpuscle problem (see Wicksell, 1925; 1926). Recently the problem of predicting extremal size based on profiles has been studied extensively; (see e.g., Drees and Reiss, 1992; Takahashi and Sibuya, 1996; 2002; Ko¨tzer and Molchanov, 2006), where both the theoretical aspects and the applications to materials science are discussed. In a series of papers, Hlubinka (2003a;b; 2006); Hlubinka and Kotz (2006), we have discussed the extremes of size or shape factor of oblate spheroids rather than balls. We present a new development of the topic in this paper. In particular we complete the study of the extremal characteristics of prolate spheroids; some new results are also proved for oblate spheroids. The analyses of the extremal characteristics of the two classes of spheroids are compared. PROLATE AND OBLATE SPHEROIDS Random microscopic particles of spheroidal shape in a volume of an opaque material are considered. We shall further restrict our study either to prolate or to oblate spheroids (see e.g., Cruz-Orive (1976) for arguments why the restriction is used). Let us suppose that the particle arrangement is stationary and isotropic, which means that the underlying point process of their centers is stationary and their orientation is a uniformly distributed random variable. Recall that a prolate spheroid is a spheroid with one major semi-axis whose length is denoted by X and two equal minor semi-axes with length V. An oblate spheroid is characterized by two equal major semi-axes (denoted by X) and one minor semi-axis (denoted by V, again). In the paper we use yet another set of characteristics of such spheroids. Define the shape factor of the spheroid by S = X2/V2 - 1. In Cruz-Orive (1976) it is proposed to use the characteristics (X,S) for oblate spheroids, while the pair (V,S) is used for prolate spheroids. However, we use the pair (X,S) beside the usual characteristics (V,S) for prolate spheroids later in the study when it is more appropriate. STEREOLOGY OF SPHEROIDS It is a usual complication in materials science that spheroids cannot be observed directly. It is, however, possible to observe in a planar section of the material, the profiles of spheroids. The random profiles of both oblate and prolate spheroids are ellipses. These ellipses are usually characterized by the lengths Y and W of their major and minor semi-axes respectively. The characteristics used further in this paper are an analogy to the characteristics of spheroids, namely the shape factor T =Y2/W2 -1 and the size (Y for oblate andW for prolate spheroids). Remark 1. Note that for a prolate spheroid the length of the minor semi-axis is used as a size 145 Hlubinka D: Size and shape factor extremes characteristic rather than the length of the major semi-axis, which seems to be the natural choice. Cruz-Orive (1976) explains why these characteristics are used in stereology. The main reason is that “the manipulation becomes much simpler” (Cruz-Orive, 1976). We will discuss the question of the size of prolate spheroid, in more detail later in the paper. As used throughout the paper, the spheroid population consists either of prolate, or of oblate spheroids exclusively. Further, we consider the characteristics (V or X, and S) to be a random vector independent of the position and orientation of the particle. The bivariate distribution of the vector is considered to be absolutely continuous with the probability density function (density) denoted further by g(v,s) and g(x,s), respectively. The density of the profile characteristics (W,T) or (Y,T) is denoted by f (w,t) and f (y,t), respectively. Let M be the population mean size of the particles (half of the mean caliper diameter) and xF and sF the upper endpoint of the marginal distribution of the size X and of the shape factor S respectively. It can be determined e.g., in Hlubinka (2003a), that the density functions for the profile characteristics are f(w,t) f (y,t) w xF ^(1 + sf/2g(v,s)dsdv yV1Tt xF sF y t 2M xp w t xf y t s{1+s){s-t)^ (P) Vs g(x,s)dsdx y2 (O) Remark 2. Note that the size V of a prolate spheroid with a given shape factor S = s cannot exceed the value vF(s) = xF/y/1+s. Similarly the shape factor S of a prolate spheroid with a given size V = v cannot exceed the value sF(v) = x2F/v2 - 1 as follows from the definition. Hence the transformation (P) should be rewritten as f (w,t) w 2M(1+t)2 sf ç t w sF vF (s) (1+s)3/2g(v,s)dvds tw w 2M{1 + tf fxF w t (P1) s v)(1+s^ßg(v s)dsdv (P2) This fact is not important when xF = 8 but it may otherwise cause a serious complication, as we shall see in the section “Application of the results”. STABILITY OF MAXIMUM DOMAIN OF ATTRACTION Recall that if for a univariate distribution function H there exist normalizing constants an and bn such that, as n ^ oo it holds that A(u) = exp(-e-u), ueR, Hn{anu + bn) —? 0, ya(u) = exp(-(-u)a), u<0, (1) for some a > 0, then H is said to belong to a maximum domain of attraction (H G MDA(-)) of the limiting distribution function A (Gumbel distribution), O (Fre´chet distribution) or ¥ (Weibull distribution), respectively. There are no other possible limiting distributions. There are well known sufficient conditions for H G MDA(-) if H is an absolutely continuous distribution function with a density h and an upper endpoint zf. Consider the conditions h(z + ub(z)) _ z/zf h(z) ueR, (2) z l im h uz) =u-(«+i); u>0,zF = +oo, (3) h(zF-uz) _ ! z\0h(zF-z) u > 0,z > 0, zF < +8, (4) where a > 0. The detailed description of the auxiliary function b(-) in condition Eq. 2 can be found e.g., in Embrechts et al. (1997). In particular b(-) can be chosen such that it is differentiable for zH£MDA(A), Eq.3^HeMDA(a), and Eq.4=>H£MDA(¥a) Remark 3. Throughout the paper we shall use the following notation for probability density functions: – gS(s), fT(t) are the marginal density functions of the shape factor of the spheroids and of the profiles, respectively; _ _ _ _ 146 Image Anal Stereol 2006;25:145-154 – g·(s) and f·(t) are the conditional densities of the spheroid’s shape factor and of their profiles respectively, given an exact value of the size; and – g>·(s) and f>·(t) are the conditional densities of the shape factor of the oblate (prolate) spheroids and of their profiles respectively, given the size exceeds a given threshold. This case may be of particular interest for an application in the material sciences as the extremal shape factor is restricted to large particles only; – We use the superscript f o(·) for the density of oblate spheroids and f p(·) for the density of prolate spheroids if it is not quite clear from the context which spheroid class is being considered. We are now ready to prove a series of results concerning the stability of MDAs. Theorem 4 (Shape factor of prolate spheroids I). Assume that the upper endpoint of the size distribution is xF = 8. Suppose further that the conditional density function gv(s) of the shape factor given the size V = v satisfies one of the conditions Eq. 2, Eq. 3 or Eq. 4 uniformly in size. Then for the conditional density fw(t) of the profile shape factor given the size the same condition holds with a parameterß replacing a. The value of ß is then ß = a+ 1/2 for both the Fre´ chet MDA and the Weibull MDA. Remark 5. In fact, we will prove the “stability of the sufficient conditions” (2–4) up to a parameter rather than “the stability of the MDA”. On the other hand it is possible but a bit more difficult to prove the stability of the equivalent conditions for the distribution functions; for these conditions see, e.g., Drees and Reiss (1992). Since the idea of such a proof is quite similar to the one presented in our proof, we prefer to use the simpler conditions (2–4). Since the MDA can be estimated from the observed data (profiles) and since for the Weibull and Fre´chet distributions the parameter can be also estimated based on the profiles, one can conclude that under the uniformity condition we also know the limiting distribution, and its parameter respectively, for the particle shape factor extremes. Remark 6. The uniformity assumption of Theorem 4 (and of the similar theorems below) implies in particular that the parameter a of Eq. 3 and Eq. 4 is the same for all values of the size (x or v respectively) and also that the auxiliary function b(·) of Eq. 2 is the same for all sizes. It also means that the upper endpoint xF does not depend on the size. There are classes of bivariate distributions constructed by generalized Farlie-Gumbel-Morgenstern copulas for which the tail uniformity holds. The reader is referred to Hlubinka and Kotz (2006) for more details. The tail uniformity condition in a similar sense may also be found for the extremal quantile regression, e.g., Chernozhukov (2005). Remark 7. There is a natural question of what can happen without the tail uniformity assumption. In such a case it is possible that the spheroid’s shape factor may belong to a different maximum domain of attractions for different sizes. As for the profile shape factor (when the profile size is fixed), the respective domain of attraction should reflect the heaviest tail of the spheroid’s shape factor for which the spheroid size exceeds the observed profile size. Proof of Theorem 4. We shall treat the three limiting cases separately. We use the fact that since xF = °° the upper bound of the inner integral in (P1) is also infinite. Denote {v)ds dv iw it v s(s-t) lim a-1/2 (1 +t)2 8 8 (a-1+s)3/2gv(as)f(v)dsdv OO t w -((a+l/2)+l (a-1 +t)2 v s(s-t) J'OO r w i w t (l+s)3/2gv(s)(j>(v)dsdv ys(s - t) (5) The last equality follows from the fact that (a-1 + t)/(1+t) ? 1 as t ? 8 and gv(as)/gv(s) ? a-(a+1) as s = t ? 8 uniformly in v and using Lemma 1.2.1 of de Haan (1970) for the limit of the ratio of integrals. _ _ 147 Hlubinka D: Size and shape factor extremes Weibull limit: Using similar arguments as for Eq. 5 we obtain limfw(sF-at)/fw(sF-t) -O2 lim~----------^ ft» rsF w sF-at w (\+s)3l2g(v,s)dsdv ^s -y/s— (sF —at) y v2 — i (1+s)3/2g(v,s)dsdv a(1+sF limal/2^ t\0 °° t (l+sF-s)3/2gv(sF-s)^(v)dsdv t ( 1 +sF —as) ' gvjsF —as)(j>(v)dsdv ^sF^as^/{sF-as)-{sF-at) t {ì+sF-s)3/2gv(sF-s)(j>(v)dsd |3/ limf1+ F -at^ r ft a+sF-sgv sF—s)(j>(v)dsdv h oo pt (1+sF—as) / gv(sF—as)(j>(v)dsdv 8t w0 vv sF-s t-s aa-\l2 X a (a+l/2)-l (6) The second equality follows from the substitutions s <-> sF - as and s <-> sF - s in the numerator and denominator, respectively. Lemma 1.2.1 of de Haan (1970) and the uniformity (in v) of the limit gv(sF -as)/gv(sF -s)-»¦ a«-1 as s \ 0 complete the proof. Gumbel limit: The situation is now complicated by the auxiliary function b(-). Recall (see e.g., Chapter 3 in Embrechts et al, 1997) that b(s) can be chosen such that it is differentiable for s < sF and s limb'(s) = 0, limsr1b(s = 0 if sF 00; or limb'{s) = 0, lim(sF-s)-1b{s) = 0 ifsF<°°. s^sF s ^sF Hence, using the second equality in Eq. 7 below the substitution s^s + ab(s) and the limit lim t/sF 1 + t + ab(t) 1+t 1 following from the properties of b(-), we obtain lim Ut+ ab(t))/fw(t) t/sF (1 +t)2 r°° sF (l+s)3/2gv(s)(l>(v)dsdv Jw Jt+ab(t) ss-{t+ab{t)) lim sF(1+t + ab(t)) 2' r rsF (l+s)3/2gv(s)^v)dsdv Jw t r°o rsF (l+s+ab(s))3'2(l+ab'(s))gv(s+ab(s))(l>(v)dsdv oo s w t lim------ t/sF exp{-4 v/s+ab(s)v/(s+ab(s))-(t+ab(t)) o, sF (l+s)3/2gv(s)t. But this follows from the properties of b(-): It holds that (s + ab(s))-(t+ab(t)) s-t 1+ab'^), for some t < £, < s and since b'{t) -^0 ast^sF, the proof is complete. D Remark 8. The case of oblate spheroids is simpler. There is no need to assume xF = °° in analogy to Theorem 4; see Hlubinka (2003a). There is, however, a difference in the parameter change. In particular, the adjusted parameter is ß = a for the Fre´chet MDA and ß = a + 1/2 for the Weibull MDA. Very similar arguments can be used to prove the stability of MDA for the marginal density functions gS(s),Mt) and for the densities g>.{t) of the shape factor conditioned by the size exceeding a given threshold. Theorem 9 (Shape factor of prolate spheroids II). Suppose that the conditional density function gv(s) of the shape factor given the size V = v satisfies one of the conditions Eq. 2, Eq. 3 or Eq. 4 uniformly in v and for some parameter a > 0, respectively. Suppose moreover that the upper end point of the size distribution is xF = °°. Then i. the marginal density gS (s) of the shape factor and the conditional densities g>v(s) of the shape factor given that the size exceeds a threshold, satisfy the same conditions (2-4) as gv(s) does. The respective parameter a of the limiting distribution is unchanged. ii. the marginal density f(t) of the profile shape factor and the conditional densities f>w(t) of the profile shape factor given the profile size exceeds a threshold satisfy the same conditions (2-4) as gv(s) does. The parameter ß of the limiting distribution isß = a+ 1 /2for the Weibull and Fre´chetMDAs. Proof. The proof is an analogy to the proof of Theorem 4. We shall briefly show the main idea for the Fre´chet limiting case. The marginal density of the profile shape factor is fT (t) 1 2M(1+t)2 io t (1 + s)3/2vg(v,s)dsdv (8) _ _ _ _ _ _ 148 Image Anal Stereol 2006;25:145-154 The arguments of Eq. 5 can be repeated again and the proof for the marginal density functions can be concluded. For the conditional density functions let us first note that g>v(s) f>w(t) S? gu(s)gV(u)du vgV(u)du ' w fu(t)fW(u)du IwfW(u)du and (9) and that the denominators do not depend on the shape factor. Hence it is easy to find g>v(as) = Çgu(as)gV(u)du ^»g^v ( s ) vgu(s)gV(u)du lim a-a+1 again as a simple consequence of the uniformity assumption. The same argument may be used for the other two limiting cases. For the last part of the proof it is sufficient to repeat the arguments in Eq. 5 line by line and to use the fact thatf>w(t ) = w fW(u)du D Quite similar arguments hold for the size characteristics. There is no problem for the size of prolate spheroids if xF = 8 (hence the Weibull MDA is excluded). Theorem 12 (Size of prolate spheroids I). Suppose that the upper endpoint xF = 8. Suppose that the conditional density function gs(v) of the size given the shape factor S = s satisfies one of the conditions Eq. 2 or Eq. 3 uniformly in s and for some parameter a> 1 for Eq. 3. Then i. the marginal density gV (v) of the size and the conditional densities g>s(v) of the size given that the shape factor exceeds the threshold s, satisfy the same of the conditions Eq. 2 or Eq. 3 as gs(v) does. The respective parameter a of the limiting distribution is unchanged. ii. all the conditional densities ft (w) of the profile size given that the shape factor T = t, the marginal density fW(w) of the profile size and the conditional densities f>t(w) of the profile size given the profile shape factor exceeds the threshold t, satisfy the same of the conditions Eq. 2 or Eq. 3 as gs(v) does. The parameter ß of the limiting distribution is ß = a-1 for the Fre´ chet MDA. Remark 10. There is an analogy to Theorem 9 for oblate spheroids again. The assumption xF = 8 may be omitted, as in Remark 8. The parameter change is ß = a for the Fre´chet MDA and ß = a+1/2 for the Weibull MDA exactly as in Remark 8. Remark 11. The reason for the uniformity assumption of course can be seen from the proofs of Theorems 4 and 9. There is, however, also an intuitive reason for the uniformity assumption. Assume that we observe a profile with size W = w. Since it holds that 0 = w =V, the observed spheroid may be of any size which is greater than or equal to the observed one. Hence the limiting behavior of the profile shape factor extreme is influenced by the limiting behavior of the shape factor extreme of any spheroid whose size is larger than the size of the observed profile. The uniformity is assumed in order to obtain just one MDA for the profile shape factor regardless of the actually observed size. If the uniformity assumption is not valid then we may expect that there can be some dominating extreme value behavior for the shape factor extremes and that this dominating behavior is connected to some values of the size. But then it would be impossible to recover the non-dominating limiting behavior of the shape factor for the other (conditioning) values of the size. Proof. The proof follows similar arguments as the proofs of Theorems 4 and 9, so it is sufficient to show the main steps only. We shall again distinguish the two limiting cases. For the Fre´chet MDA we can write (1+s)3/2g(v,s)dvds ft (aw) lim lim a sF t aw ¦/s^/s-tyv2 -a2w2 t~ t sF wc (1+s)3/2g(v,s)dvds rsp p (1+s) / gs(av)gs(s)advds Jt Jw vsvs - tVa2 v 2-a2w2 a lim w" CsF roo (1+s)3/2gs(v)gS(s)dsdv sF tw -s s - tv- v2-w2 sf (1+s)3/2gS(s) ts s - t wjv 2- w 2 dvds a lim — w^oo ,sF (1+s)3/2gs(s) woo gs(v) dvds -((a-1)+1 as follows from the assumptions of the theorem. (10) _ _ 149 Hlubinka D: Size and shape factor extremes We can continue with the Gumbel MDA. Now due to the properties of b(-), we obtain lim ft (w+ab(w))/ft(w) w" sF tw (1+s)3/2g(v,s)dvds w li m rsp r°° (1+s)3'2 g(v,s)dsdv w t w = lim t—sp (l+s)3i2(l+ab'(s))gs(v+ab(v))gs(s)dvds ^Js\/s—t ^(v+ab(vyf-(w+ab(wyf rsp roo (l-\-s)3/2gs(v)gs(s)dvds exp{-4 (11) by the Lemma 1.2.1 of de Haan (1970), again. In particular we need to check that w li m [v + ab(v))2 - (w + ab(w)f v2w2 1 when v > w ^ oo. But this follows from the properties ofb(-). It holds that (v + ab(v))2-(w + ab(w))2 v2w2 vh( 1+2a vb(v) —wb(w) (v — w)(v + w) 2(b(v)-b(w))(b(v) + b(w)) + a (v-w)(v + w) 1 as v > w ->¦ o« and since there exists £ G (w,v) such that b(v)-b(w) v — w b'&^0 as w —> °° and also b(w)/w -? 0, as w -^ oo. The proof for the conditional density f>t(w) is easy and follows from the fact that f>t{w) = where the term 1 1-FT(t)' ?(s,t) 1 2^ 2^-0(1 + s) 1 + t +log 1 + tts does not depend on w. The rest of the proof is the same as the proof of Theorem 9 for the Fre´chet MDA. Since the arguments for the Gumbel MDA are similar, they can be omitted and the proof is complete. D Remark 13 (Size of oblate spheroids). There is an analogy to Theorem 12 for oblate spheroids again. For the marginal and conditional densities gS(s),fY(y) and ft{y) it is Hlubinka (2003b, Theorem 2 in). The same result can be proved for the conditional densities g>s{x) and f>t{y) using similar arguments as in the proof of Theorem 12. Recall that the parameter change is ß = a - 1 for the Fre´chet MDA (compare with the prolate case) and ß = a + 1/2 for the Weibull MDA. PROLATE SPHEROIDS WITH BOUNDED SIZE We turn our attention to prolate spheroids with the size bounded from above. The situation is now much more difficult than the unbounded case. We summarize the essential differences in a remark. Remark 14. The pairs (X,V) and (X,S) are in a one-to-one correspondence given by (V1TS,S) = (X,S) = (X,X2/V2-1), and so are (W,Y) and (T,Y). In particular, having a spheroid with a given shape factor S = s there is an upper bound for the minor semi-axis, namely V < xF/^1+s. In this respect, there is no difference between oblate and prolate spheroids. But, and this is the main difference, the value of V is also bounded from below for prolate spheroids by the value of the profile size W. It follows that observing a profile of a prolate spheroid with (profile) size W = w the shape factor of the spheroid is limited from above by S < x2F/w2 - 1. There is no similar bound for oblate spheroids. On the other hand, observing the profile of a prolate spheroid with (profile) shape factor T = t it is clear that the size of the spheroid is limited from above by V < xF/1Tt. Again there is no similar bound for < oblate spheroids. We shall see later that the fact V < xF/1+t causes particular technical problems for the method of proof we have applied in section “Stability of maximum domain of attraction” when we were studying the profile size (or shape factor) given the profile shape factor (or size, respectively). There is, however, no evidence that “MDA stability” does not hold in such a case and in particular we suspect that, possibly under additional conditions, the analogies to Theorems 4, 9 and 12 may be proved using a different method of proof. _ _ _ _ _ 150 Image Anal Stereol 2006;25:145-154 SIZE where Let us now consider the length of the major semi-axis X and the shape factor S to be the characteristics of a prolate spheroid with joint density g?(x,s). The joint density of the analogous characteristics of the profiles then becomes (we use the superscript ? to emphasize that x and y are used as the size) f?(y,t) y 2M(1+t)5/2 sF t 1 xF 1 [^/ï+s<^-^/ï+t] ](1 + s)2g*(x,s)dxds y/& s ^ - t) V(1+0x2-(1+s y2 ' (12) for 0 < t < sF, and 0 otherwise. We have denoted by 1 [A] the indicator function of a set A. The main problem lies in the fact that for a given prolate spheroid the upper bound of the profile size Y depends on the unknown shape factor S. In particular assuming a spheroid whose shape factor is S = s then Y 0 the condition Eq. 4 uniformly in s. Then i. the marginal density g*X (x) of the size also satisfies the condition Eq. 4. The parameter a of the limiting distribution is unchanged. ii. the marginal density fY£ (y) of the profile size also satisfies the conditions Eq. 4 with parameter ß. The parameter ß of the limiting distribution is ß = a+1. Proof. We show only the principal argument for (ii) since the idea of the proof is essentialy the same as for the other theorems. We need to study ff (xF - az) /fY (xF - z) ~ (xF-zf £F V(s) %ag*(xF-ax)Çaz(ax)dxds (xF-az)3 JoF W(s) Jog*s(xF -x)^z{x)dxds -* aaa-1 = aa, as z \ 0, where “A(z) ~ B(z)” is used in the sense |A{z)/B{z) -1| -^0, asz^0, Çz(x)=(2(xF-x) 2 -(xF-z) 2) and ijf(s) =g*S(s)[s(1 + s)}-1/2 does not depend on z. n Remark 16. Comparing Theorem 15 and Remark 13 we can see that the parameter ß is different for oblate and prolate spheroids. In particular the tail of the profile size is lighter in the prolate case. The reason is in the fact that there are two major semiaxes in the oblate case and hence the profile size can attain its maximum xF regardless of the orientation of the sectioning plane respective to the orientation of the spheroid. On the other hand, in the prolate case the size Y can reach the upper endpoint xF if and only if the sectioning plane is parallel to the major semiaxis of the spheroid. _ y 151 Hlubinka D: Size and shape factor extremes SHAPE FACTOR The situation of the shape factor is quite similar. Hence we have to restrict our study to the MDA of the marginal distribution of the shape factor. For the conditional distribution of the shape factor given the size we encounter similar difficulties as for the size. We shall use the length of the major semi-axis X as the more natural choice of the size characteristic of the spheroid. The marginal distribution of the profile shape factor T is obtained by integrating of f ?(y,t), see Eq. 12, over y. Hence fT?(t) 1 sF 1 + s xf xg*(x,s)dxds, (14) where g* (x, s) is the joint density of the random vector (X,S). We can see that the limits of the inner integral do not depend on the unknown value of S. Now we can prove the following theorem. Theorem 17 (Shape factor of prolate spheroids III). Assume that the upper endpoint of the size distribution satisfies xF < °°. Suppose further that the conditional density function g*(s) of the shape factor given the size satisfies one of the conditions Eq. 2, Eq. 3 or Eq. 4 uniformly in size. Then, for the marginal density fT (t) of the profile shape factor, the same conditions hold with a respective parameter ß. The parameter satisfies ß = a+1for the Fre´chet MDA andß = a+ 1/2for the WeibullMDA. Remark 18. Comparing Theorems 4 and 17 we can see that only the result for the Fre´chet MDA is different. In particular, the tail of the profile shape factor is lighter than the tail of the spheroid shape factor. The reason is probably in the fact that having an unbounded shape factor and bounded size the shape factor is extremal when V is negligible. Note that a large profile shape factor is observed if the sectioning plane is almost parallel to the major semi-axis. But such an event is quite rare since V is very small. Proof We show only the essential part, as the proof is quite similar to the proof of Theorem 4. For the Fre´chet MDA we get t l im fT?(at)/f?(0 = (1+tf xoF x*g(x) at (1 s s g s ds dx lim lim (1TatF ~r x* (x) r (1+s)gi(s)dsdx O g X t \fs\/s—t (1+t)2 SoFxg*X(x)Jt lim (l-\-as)g*(as)ads ^/as^Jas—at (15) dx lim~ (1 + a)! fxF {x r (i+s)gi(s)dsdx o xgX\ > t *ys^s—t a-(a+2) • For the Weibull MDA it holds that limfî(sF-at)/fî(sF-t) = t\0 (1+sF -tf [xF rsF xg*x(x)(l+s)g*(s)dsdx O Jsf-at s /s-(s F -at) t l i\p(1+s F -at) 2 (xF fsF xgj(x)(l+s)gis dsdx O - ' fxF ft (1+sF-t)2 Jo Jo lim--------rat aa-1/2 aa+\l2)-\ s F -t ^sy/s-(s F -t) xF ft xg^(x)(l+sF-as)gì(sF-as)adsdx s=asy/(sF-as)-(sF-at) *v(x)(l+sFs)g1c(sF-s)dsdx sFSy/(sF-s)-(sF-t) t lim 0 (1+sF-at)2 fxF ft xgi(x)(l+sF-s)g*(sF-s)dsdx (16) Finally for the Gumbel case it holds that lim fr(t + ab(t))/fì(t) t/s F (1 +t)2 rxF rsF xgx(x)(\+s)g*(s)dsdx O t+ab(t) s^s-(t+ab(tj) t l im s{1+t + ab{t)Y' (xF rsF xgi{x)(l+s)gî(s)dsdx Jo Jt js^t [xp rsF xgx(x)(\+s+ab(s))(\+ab'(s))g*(s+ab(s))dsdx h t lim t/sF sTab s) ^(s+ab(s)) - (t+ab(t)) fxF rsF xg*x(x)(l+s)g*(s)dsdx Jo t vv s s-t = exp{-s}. The proof is complete. (17) D APPLICATION OF THE RESULTS In this section we discuss the possible application of the above theorems. We are interested in the distribution of the extreme spheroid characteristics. Let us denote by C the spheroid characteristics of interest and by Cn:n the sample maximum. The corresponding profile characteristics are denoted by D and the sample maximum is Dm:m. What can be said about P[Cn:n < x] for large n, based on the profile sample and using the results of the above theorems? First of all, the extreme value index k for the profile characteristics D can be estimated from the observed profiles. Recall that the sign of the extreme value index corresponds to the maximum domain of attraction. In particular for the Weibull MDA it holds that k < 0, k> 0 for the Fre´chet MDA and k = 0 for the Gumbel MDA. The estimate ß of the respective parameter of the Fre´chet or Weibull limiting distribution can also be based on the observed profiles. There are classical 152 Image Anal Stereol 2006;25:145-154 methods for these estimates; see e.g., Coles (2001) or Beirlant et al. (2004) for recent references. Then, under the uniformity conditions of the theorems, we can estimate the maximum domain of attraction for C, and its respective parameter is estimated by a simple transformation of ß as proved in the Theorems above. Assume that the characteristic C (and also D) belongs to the MDA T, where T is one of the three possible limiting distributions. Hence, there exist normalizing constants (an,bn) suchthat P[Cn:n0,j3,7>0,<5>0 such that 1-K(v) 1, then the normalizing constants can be chosen as an bn logn 1/5-1 1 1/5 g + ¦ logn1/ü f (log log n-log y) + log a Ï logn N 1-1/5 y8 This lemma and its proof can be found in Takahashi (1987). Similar results are also known for the other two MDA; see e.g., Embrechts et al. (1997) for a possible choice of the normalizing constants based on the quantile function. Lemma 19 is generalized by Lemma 3.8 in Hlubinka and Kotz (2006) which is of particular interest if the generalized Farlie-Gumbel-Morgenstern bivariate distribution is used for the spheroid’s characteristics. The procedure for the case of the Gumbel MDA can be the following. For some parametric model describing the asymptotic tail behavior of the spheroid characteristics we shall derive the formula for the normalizing constants (asn,bsn) using Lemma 19, where n is considered to be a variable rather than a constant. The same Lemma 19 applied to the appropriate density function based on transformation (O) or (P), respectively, can be used to derive the formula for the normalizing constants (amp ,bmp ) for the profile characteristics. Comparison of these two pairs of normalizing constants will suggest how to adjust the estimated normalizing constants of the profile characteristics to the estimate of the normalizing constants of the spheroid characteristics. For a more detailed discussion of the recalculation of the normalizing constants, see Hlubinka (2006) and Benes? et al. (2003). CONCLUSION We have shown that under a tail equivalence assumption the maximum domain of attraction of spheroid characteristics is preserved in their profiles up to a simple change of the respective parameter. These results can be used for the estimation of the distribution of spheroid characteristics extremes based on the observed profile characteristics extremes. The exception is the class of prolate profiles with a bounded size. In this case it seems that the proposed analysis can only be used for the marginal distribution of the size or the shape factor. The study of the conditional distribution of the size given the shape factor and vice versa must probably be based on a different approach in this situation. _ _ _ 153 Hlubinka D: Size and shape factor extremes ACKNOWLEDGMENT The work is a part of the research project MSM 0021620839 financed by MSMT The author thanks the referees for their very careful reviews and valuable comments. REFERENCES Benes? V, Bodla´k K, Hlubinka D (2003). Stereology of extremes; FGM bivariate distributions. Method Comp Appl Prob 5:289-308. Beirlant J, Goegebeur Y, Teugels J, Segers J, De Waal D, Ferro C (2004). Statistics of extremes. Theory and applications. Chichester: Wiley. Chernozhukov V (2005). Extremal quantile regression. Ann Stat 33:806-39. Coles S (2001). An introduction to statistical modeling of extreme values. London: Springer Verlag. Cruz-Orive L-M (1976). Particle size-shape distributions; the general spheroid problem. J Microsc 107:235-53. Drees H, Reiss R-D (1992). Tail behavior in Wicksell’s corpuscle problem. In: Galambos J, Ka´tai J, eds. Probability Theory and Applications. Dordrecht: Kluwer, 205-20. Embrechts P, Klu¨ppelberg C, Mikosch T (1997). Modelling extremal events for insurance and finance. Berlin: Springer Verlag. de Haan L (1970). On regular variation and its application to the weak convergence of sample extremes. Math Centre Tracts 32. Amsterdam: Mathematisch Centrum Amsterdam. Hlubinka D (2003a). Stereology of extremes; shape factor of spheroids. Extremes 5:5-24. Hlubinka D (2003b). Stereology of extremes; size of spheroids. Math Bohem 128:419-38. Hlubinka D (2006). Extremes of spheroid shape factor based on two dimensional profiles. Kybernetika 42:77-94. Hlubinka D, Kotz S (2006) The generalized FGM distribution and its application to stereology of extremes. To appear. Ko¨tzer S, Molchanov I (2006). On the domain of attraction for the lower tail in Wicksell’s corpuscle problem. In: Benes? V, Lechnerova´ R, Saxl I, eds. Proceedings S4G. ? Praha: JCMF, 91-6. Takahashi R (1987). Normalizing constants of a distribution which belongs to the domain of attraction of the Gumbel distribution. Stat Probabil Lett 5:197-200. TakahashiR and Sibuya M (1996). The maximum size of the planar sections of random spheres and its application to metalurgy. Ann Inst Stat Math 48:127-44. Takahashi R, Sibuya M (2002). Maximum size prediction in Wicksell’s corpuscle problem for the exponential tail data. Extremes 5:55-70. Wicksell SD (1925). The corpuscle problem I. Biometrika 17:84-99. Wicksell SD (1926). The corpuscle problem II. Biometrika 18:152-72. 154