arxiv:math/ v1 [math.st] 12 Nov 2004 ON REPRESENTING THE MEAN RESIDUAL LIFE IN TERMS OF THE FAILURE RATE

Transcription

1 arxiv:mah/411297v1 [mah.st] 12 Nov 24 ON REPRESENTING THE MEAN RESIDUAL LIFE IN TERMS OF THE FAILURE RATE RAMESH C. GUPTA AND DAVID M. BRADLEY Absrac. In survival or reliabiliy sudies, he mean residual life or life expecancy is an imporan characerisic of he model. Whereas he failure rae can be expressed quie simply in erms of he mean residual life and is derivaive, he inverse problem namely ha of expressing he mean residual life in erms of he failure rae ypically involves an inegral of a complicaed expression. In his paper, we obain simple expressions for he mean residual life in erms of he failure rae for cerain classes of disribuions which subsume many of he sandard cases. Several resuls in he lieraure can be obained using our approach. Addiionally, we develop an expansion for he mean residual life in erms of Gaussian probabiliy funcions for a broad class of ulimaely increasing failure rae disribuions. Some examples are provided o illusrae he procedure. 1. Inroducion In life esing siuaions, he expeced addiional lifeime given ha a componen has survived unil ime is a funcion of, called he mean residual life. More specifically, if he random variable X represens he life of a componen, hen he mean residual life is given by m() = E(X X > ). The mean residual life has been employed in life lengh sudies by various auhors, e.g. Bryson and Siddiqui (1969), Hollander and Proschan (1975), and Muh (1977). Limiing properies of he mean residual life have been sudied by Meilijson (1972), Balkema and de Hann (1974), and more recenly by Bradley and Gupa (22). A smooh esimaor of he mean residual life is given by Chaubey and Sen (1999). I is well known ha he failure rae can be expressed quie simply in erms of he mean residual life and is derivaive: see (2.4) below. However, he inverse problem namely ha of expressing he mean residual life in erms of he failure rae ypically involves an inegral of a complicaed expression. In his paper, we obain a simple expression for he mean residual life in erms of he failure rae for cerain classes of disribuions. Many of he sandard cases are subsumed, and several resuls in he lieraure can be obained Dae: February 1, Mahemaics Subjec Classificaion. Primary: 62N2; Secondary: 62N5. Key words and phrases. Pearson family of disribuions, Gaussian probabiliy funcion, ulimaely increasing failure rae. 1

2 2 RAMESH C. GUPTA AND DAVID M. BRADLEY using our approach. However, he emphasis here is o express he mean residual life in erms of he failure rae. For he class of ulimaely increasing failure rae disribuions, we also provide sufficien condiions under which he mean residual life can be expanded in erms of Gaussian probabiliy funcions. Finally, some examples are presened o illusrae he procedure. 2. Background and Definiions Le F : [, ) [, ) be a nondecreasing, righ coninuous funcion wih F() =, lim x F(x) = 1, and le ν denoe he induced Lebesgue-Sieljes measure. (Equivalenly, le ν be a probabiliy measure on [, ) and le F be he cumulaive disribuion funcion of ν.) If X is a nonnegaive random variable represening he life of a componen having disribuion funcion F, he mean residual life is defined by m() = E(X X > ) = 1 F() (x ) dν(x),, where F = 1 F is he so-called survival funcion. Wriing x = x du and employing Tonelli s heorem yields he equivalen formula m() = 1 F() x du dν(x) = 1 F() u dν(x) du = 1 F() F(u) du, (2.1) which is someimes also used as a definiion. The cumulaive hazard funcion may be defined by R = log F. Then (2.1) implies ha m() = exp { R() R( + x) } dx. (2.2) If F (equivalenly, ν) is also absoluely coninuous, hen he probabiliy densiy funcion f and he failure rae (hazard funcion) r are defined almos everywhere by f = F and r = f/ F = R, respecively, and hen R() = log F() = df(x) F(x) = r(x) dx. (2.3) In view of (2.2) and (2.3), we have expressed m in erms of r, albei somewha indirecly. Ideally, we d like o express he mean residual life in erms of known funcions of he failure rae and is derivaives wihou he use of inegrals. In any case, i is useful o have alernaive represenaions of he mean residual life. We noe ha he converse problem, ha of expressing he failure rae in erms of he mean residual life and is derivaives is rivial, for (2.1) and (2.3) imply ha m () = r()m() 1. (2.4)

3 REPRESENTING THE MEAN RESIDUAL LIFE 3 3. A General Family of Disribuions, Including he Pearson Family Consider he family of disribuions whose probabiliy densiy funcion f is differeniable. Wrie f(x) = µ x g(x) g (x) g(x), (3.1) where µ is a consan, and g saisfies he firs order linear differenial equaion g (x) + g(x) = µ x. f(x) Alernaively, we may view g as given; hen f is uniquely deermined by (3.1). Clearing he fracions in (3.1) and inegraing yields In oher words, or equivalenly, xf(x) dx = µf() + g()f(),. E(X X > ) = µ + g()r(),, (3.2) m() = µ + g()r(),. (3.3) Thus, he mean residual life has been expressed in erms of he failure rae r, he given funcion g, and he consan µ. By appropriaely specializing g in (3.3), one can obain many of he imporan cases ha have appeared in he lieraure. We noe ha a resul similar o (3.3) was obained by Ruiz and Navarro (1994), bu heir emphasis was differen The Pearson Family. The quadraic funcion g(x) = a + a 1 x + a 2 x 2 wih a 2 1/2 yields he Pearson family, special cases of which include he bea disribuions, he gamma disribuion, and he normal disribuion. From (3.1), we have f(x) = µ x a + a 1 x + a 2 x 2 a 1 + 2a 2 x a + a 1 x + a 2 x = x + d 2 A + A 1 x + A 2 x2, (3.4) say, where A j = a j /(1 + 2a 2 ) for j =, 1, 2 and d = (a 1 µ)/(1 + 2a 2 ). For he Pearson family, (3.2) gives E(X X > ) = A 1 d 1 2A 2 + A + A 1 + A A 2 r(). (3.5) See Nair and Sankaran (1991) and Ruiz and Navarro (1994). Noe ha since A 2 = a 2 /(1 + 2a 2 ), i is impossible for 1 2A 2 o vanish in (3.5).

4 4 RAMESH C. GUPTA AND DAVID M. BRADLEY The Bea Disribuions. Wih g(x) = x(1 x)/(a+b), we have he bea disribuion of he firs kind: I saisfies f(x) = Γ(a + b) Γ(a)Γ(b) xa 1 (1 x) b 1, < x < 1, a >, b >. f(x) = x + (1 a)/(a + b 2) x(1 x)/(a + b 2), a + b 2, and hence belongs o he Pearson family given by (3.4) wih d = (1 a)/(a + b 2), A =, A 1 = A 2 = 1/(a + b 2), and a + b 2. Thus, E(X X > ) = µ + where µ = a/(a + b). If a + b = 2, we have f(x) = 1 Γ(a)Γ(2 a) (1 ) r(), < < 1, a + b ( x 1 x and g(x) = x(1 x)/2, µ = a/2, from which (3.2) gives ) a 1, < a < 2, E(X X > ) = 1a + 1 (1 )r(), < < The bea disribuion of he second kind has he form Since hen f(x) = f(x) cxβ 1, c, α, β, γ, x >. (γ + x) α+β = x + (1 β)γ/(α + 1) x 2 /(α + 1) + γx/(α + 1), f belongs o he Pearson family (3.4) wih d = (1 β)γ/(1+α), A =, A 1 = γ/(1+α), and A 2 = 1/(1 + α). Hence, ( ) 2 + γ E(X X > ) = µ + r(), α 1, α 1 where µ = βγ/(α 1) is he mean. A similar resul was obained by Ahmed (1991) using a compleely differen approach The Gamma Disribuion. Le B >. The linear funcion g(x) = Bx for x > yields he gamma disribuion. In his case, (3.1) gives f(x) = x + B µ, Bx

5 REPRESENTING THE MEAN RESIDUAL LIFE 5 so ha f belongs o he Pearson family (3.4) wih d = B µ, A =, A 1 = B, and A 2 =. In fac, f(x) = B µ/b Γ(µ/B) xµ/b 1 e x/b, a gamma disribuion wih mean µ. Hence (3.2) akes he form E(X X > ) = µ + Br() in his case. This is Theorem 1 of Osaki and Li (1988). See also El-Arishy (1995) The Normal Disribuion. Le σ >. The consan funcion g(x) = σ 2 for < x < yields he normal disribuion. In his case, (3.1) gives f(x) = µ x, σ 2 so ha f belongs o he Pearson family (3.4) wih d = µ, A = σ 2, and A 1 = A 2 =. I follows ha f(x) = σ 1 (2π) 1/2 exp { (x µ) 2 /2σ 2}, a normal disribuion wih mean µ and variance σ 2. Hence (3.2) can be wrien as E(X X > ) = µ + σ 2 r(). This is Theorem 3.1 of Ahmed and Abdul-Rahman (1993). See also McGill (1992) and Koz and Shanbhag (198) The Maxwell Disribuion. The Maxwell disribuion has he form f(x) = 4b 3 π 1/2 x 2 exp{ x 2 /b 2 }, x >, b >. See El-Arishy (1993). We hen have f(x) = 2 x 2x b 2, which can be wrien in he form (3.1) wih µ = and g(x) = (1 + b 2 /x 2 )b 2 /2. Hence, E(X X > ) = (2 + b 2 )b 2 r() 2 2. The fac ha he corresponding failure rae is increasing can be seen by using Glaser s (198) resul.

6 6 RAMESH C. GUPTA AND DAVID M. BRADLEY 4. Ulimaely Increasing Failure Rae Disribuions Consider he class of disribuions whose failure rae is ulimaely increasing. More specifically, he failure rae should be sricly increasing from some poin onward. Obviously, he imporan class of lifeime disribuions having a bahub-shaped failure rae wih change poins 1 2 < (i.e. for which he failure rae is sricly decreasing on he inerval [, 1 ], consan on [ 1, 2 ] and sricly increasing on [ 2, )) consiues a proper subclass of he disribuions we consider here. We ll see ha if he failure rae is sricly increasing from some poin onward, hen under cerain addiional condiions he mean residual life can be expanded in erms of Gaussian probabiliy funcions. Noaion. Our convenions regarding he Bachmann-Landau O-noaion, he Vinogradov -noaion, and o-noaion are fairly sandard. Thus, if h is a funcion of a posiive real variable, he symbol O(h()),, denoes an unspecified funcion g for which here exis posiive real numbers and B such ha g() B h() for all real >. For such g we wrie g() h() or g() = O(h()). The noaion g() = o(h()),, means ha for every real ε >, no maer how small, here exiss a posiive real number such ha g() ε h() whenever >. Theorem 1. Suppose ha from some poin onward, he failure rae r increases (sricly) wihou bound. Suppose furher ha for some posiive ineger n, he n 1 derivaive is coninuous and saisfies r (n 1) ( + x) r (n 1) (),, (4.1) uniformly in x for x min ( 1, r () 1/3), and r (j) () max ( 1, r () (j+1)/3),, 3 j n 1. Finally, suppose here exiss a posiive real number ε such ha for each ineger j in he range 3 j n, r (j 1) () = o ( (r()) j jε),. (4.2) Then we have he following expansion for he mean residual life: m() = b k ()ϕ k () + o ( (r()) 1 nε),, (4.3) k= where he coefficiens b k () are given by he formal power series ideniy k= { b k x k = exp k=3 } r (k 1) () xk, (4.4) k!

7 REPRESENTING THE MEAN RESIDUAL LIFE 7 and ϕ k () = x k exp{ xr() 1 2 x2 r ()} dx [ ( ( )) { ] = ( 1) k 2π k p p 2 1 Φ exp. r () p k r () 2r ()} p=r() (4.5) Here, Φ(x) = 1 2π x e v2 /2 dv (4.6) is he Gaussian probabiliy funcion, i.e. he cumulaive disribuion funcion of he sandard normal disribuion. Before proving Theorem 1, we make some preliminary remarks and give wo illusraive examples. Firs, if r () =, hen he uniformiy condiion on x in (4.1) should be inerpreed as x 1. Nex, observe ha he hypohesis (4.2) becomes more resricive as ε increases. In paricular ε > 1 implies lim r (j 1) () =. Of course, as ε increases, he error erm in (4.3) decreases. On he oher hand, if < ε < 1, hen r (j 1) () does no necessarily approach zero, bu he correspondingly weaker hypohesis implies a weaker conclusion (larger error erm). In any case, since r increases wihou bound, he error erm ends o zero as. Addiionally, if r is infiniely differeniable we may le n in (4.3) o obain he convergen infinie series expansion m() = b k ()ϕ k (), k= valid for all sufficienly large values of. (More specifically, for hose for which r() > 1.) In general, however, we do no assume he failure rae has infiniely many derivaives; n is fixed and he generaing funcion (4.4) is a formal power series. Expanding (4.4) o compue b k () in erms of r() and is derivaives shows ha if r has only n 1 derivaives, hen b k () is undefined if k n. Differeniaing (4.4) leads o he recurrence b k+1 () = 1 k + 1 k j=2 r (j) () b k j (), k 2, (4.7) j! from which he coefficiens b k () may be successively deermined, saring wih he iniial values b () = 1, b 1 () = b 2 () =. On he oher hand, an applicaion of he mulinomial heorem yields he explici represenaion k/3 b k () = ( 1) ( ) p 1 r (j) αj (), (4.8) α p= j 2 j! (j + 1)!

8 8 RAMESH C. GUPTA AND DAVID M. BRADLEY in which k/3 is he greaes ineger no exceeding k/3 and he inner sum is over all non-negaive inegers α 2, α 3,... such ha j 2 α j = p and j 2 (j + 1)α j = k. Finally, we noe ha he funcions ϕ k of (4.5) may also be given more explicily. By seing a = r() and b = 2r () in Lemma 1 below, we find ha ) (k+1)/2 ( 2 ϕ k () = ( 1) k r () { k/2 ( ) ( k (1 ( )) λ k/2 h Γ(h + 1/2) Φ 2λ e λ + 1 2h 2 h= 1 2 k/2 h= h 1 j= ( k h! 2h + 1 ) ) λ j+1/2 Γ(j + 3/2) λ (k 1)/2 h h j= } λ j, j! where λ = (r()) 2 /2r (), Γ(h + 1/2) = π 1/2 h j=1 (j 1/2), and Φ denoes he Gaussian probabiliy funcion (4.6). Lemma 1. Le a be a real number, le b be a posiive real number, and le k be a nonnegaive ineger. Then x k exp{ ax bx 2 } dx { k/2 = ( 1) k b (k+1)/2 where λ = a 2 /4b. h= ( ) ( k (1 ( )) λ k/2 h Γ(h + 1/2) Φ 2λ e λ + 1 h 1 2h 2 j= 1 2 k/2 h= ( ) k h! λ (k 1)/2 h 2h + 1 Proofs of Theorem 1 and Lemma 1 are relegaed o 6, he final secion. ) λ j+1/2 Γ(j + 3/2) h j= } λ j, j! 5. Applicaions We provide wo examples indicaing how Theorem 1 may be applied. Example 1. Consider a linear failure rae of he form r() = α + β, β >. (5.1) The moivaion and applicaion of (5.1) o analyzing various daa ses has been demonsraed by Kodlin (1967) and Carbone e al. (1967). Saisical inference relaed o he

9 REPRESENTING THE MEAN RESIDUAL LIFE 9 linear failure rae model has been sudied by Bain (1974), Shaked (1974) and more recenly by Sen and Bhaacharya (1995). For his model, he hypoheses of Theorem 1 are rivially saisfied for any posiive ineger n and any posiive real number ε. Since r vanishes idenically in his case, we see ha b k () = for k > in (4.3) and in fac we have he exac resul { }( ( )) (α + β) m() = exp{ (α + β)x βx 2 2 α + β 2π /2} dx = exp 1 Φ 2β β β. Example 2. Chen (2) proposes he wo-parameer disribuion wih cumulaive disribuion funcion given by F() = 1 exp {( 1 exp( β ) ) λ }, >, where λ > and β > are parameers. The corresponding hazard funcion is he ulimaely sricly increasing funcion of given by r() = λβ β 1 exp( β ), >. (5.2) I is sraighforward, albei somewha edious, o verify ha Chen s failure rae (5.2) saisfies he hypoheses of Theorem 1 wih n > 2 and < ε 2/3. Clearly ε = 2/3 is opimal here. Thus, wih derivaives of r in (4.4) and (4.5) now coming from (5.2), we see ha he asympoic formula m() = b k ϕ k () + o ( (r()) 1 2n/3),, k= holds for all inegers n > 2. In paricular, as he error erm in (4.3) ends o zero in he limi as n, we obain he convergen infinie series represenaion m() = b k ()ϕ k (), k= valid for all sufficienly large values of. There is no need o work ou he coefficiens b k () explicily in his case. One can simply use he recurrence (4.7) o generae hem. 6. Proofs 6.1. Proof of Theorem 1. Since he failure rae is sricly increasing from some poin onward, here exiss such ha r () > for all. Also, since lim r() =, here exiss 1 such ha r() 1 for 1. Now le max(, 1 ), δ = δ() = min ( 1, 1/ 3 r () ), and se I() := δ exp{r() R( + x)} dx, J() := δ exp{r() R( + x)} dx,

10 1 RAMESH C. GUPTA AND DAVID M. BRADLEY so ha m() = I() + J(). We have J() δ r( + x) exp {R() R( + x)} dx = exp {R() R( + δ)} { δ } = exp δr() xr ( + δ x) dx exp { δr()}. Bu r = o(r 3 3ε ). By definiion of δ, i follows ha from some poin onward we mus have δr min(r, r ε ). Therefore, if we se ν = min(1, ε), hen ν > and J() exp { (r()) ν}, (6.1) for all sufficienly large values of. Nex, we wrie where I() = = δ δ { exp r (k 1) () xk k! 1 (n 1)! k=1 ( ) b k ()x k + E n (x, ) k= { E n (x, ) = exp r (k 1) () xk k! 1 (n 1)! k=3 ( n 1 = O x n max r (j) () ) α j, j=2 x } u n 1 r (n 1) ( + x u) du dx exp { xr() 1 2 x2 r () } dx, x } u n 1 r (n 1) ( + x u) du b k ()x k k= and he maximum is aken over all non-negaive inegers α j saisfying n 1 j=2 (j+1)α j = n. In view of he fac ha r (j 1) = o(r j jε ) for 3 j n, i follows ha E n (x, ) = o(x n (r()) n nε ), x δ.

11 REPRESENTING THE MEAN RESIDUAL LIFE 11 If we now wrie I() = b k () k= hen we find ha x k exp { xr() 1 2 x2 r () } dx k= I = b k ϕ k + k= b k () x k exp { xr() 1 2 x2 r () } dx δ k= + δ E n (x, ) exp { xr() 1 2 x2 r () } dx, O ( b k r k e δr) δ ) + o (r n nε x n e xr dx. The hypoheses on r and he definiion of he coefficiens b k imply ha b k = O(r k kε ). From he derivaion of he esimae (6.1) for J we recall ha exp( δr) exp( r ν ), a leas from some poin onward. Finally, as i follows ha δ x n e xr dx x n e xr dx = n! r n 1, I = b k ϕ k + o ( r 1 nε). (6.2) k= Since m = I + J, combining (6.1) and (6.2) gives he saed resul for m. To complee he proof, i remains only o esablish he assered evaluaion of he inegrals ϕ k. Bu his is readily obained by compleing he square in he exponenial and differeniaing under he inegral Proof of Lemma 1. By a sraighforward change of variables, we find ha where ( 2b (k+1)/2 x k exp{ ax bx 2 } dx = e λ 1/2 λ 1/2) k e 1/2 d Γ(α, λ) = λ = e λ k λ h= ( ) k ( 1) k h λ (k h)/2 Γ ( h+1 2 h, λ), (6.3) α 1 e d (6.4)

12 12 RAMESH C. GUPTA AND DAVID M. BRADLEY is he incomplee gamma funcion. If we inegrae (6.4) by pars and hen divide boh sides of he resul by Γ(α + 1) = αγ(α), we obain he recurrence formula which can be ieraed o give Γ(α + 1, λ) Γ(α + 1) = λα e λ Γ(α) + Γ(α, λ) Γ(α), Γ(α + k, λ) Γ(α + k) k 1 = e λ λ α+h Γ(α, λ) + Γ(α + h + 1) Γ(α), (6.5) valid for any non-negaive ineger k. In paricular, when α =, Γ(k, λ) Γ(k) h= k 1 = e λ h= λ h h!. (6.6) Equaion (6.6) is valid for all posiive inegers k if λ ; i is also valid when k = if λ >. Subsiuing α = 1/2 in (6.5) yields Γ(k + 1/2, α) Γ(k + 1/2) k 1 = e λ λ k+1/2 Γ(1/2, λ) + Γ(k + 3/2) Γ(1/2) h= k 1 = e λ λ k+1/2 ( Γ(k + 3/2) Φ ( ) ) 2λ. (6.7) h= Using (6.6) and (6.7), we ge he saed resul from (6.3). References [1] Ahmed, A. N. (1991). Characerizaion of bea, binomial and Poisson disribuions. IEEE Transacions on Reliabiliy, 4 (3), [2] Ahmed, A. N. and Abdul-Rahman, A. A. (1993). On characerizaion of he normal, one-sided normal, lognormal, and logisic disribuions via condiional expecaions. Pakisan Journal of Saisics, 9 (3)B, [3] Bain, L. J. (1974). Analysis for he linear failure rae life esing disribuion. Technomerics, 16 (4), [4] Balkema, A. A. and de Hann, L. (1974). Residual life ime a grea age. Annals of Probabiliy, 2 (5), [5] Bradley, D. M. and Gupa, R. C. (22). Limiing behaviour of he mean residual life. Annals of he Insiue of Saisical Mahemaics, o appear. [6] Bryson, C. and Siddiqui, M. M. (1969). Some crieria for aging. Journal of he American Saisical Associaion, 64, [7] Carbone, P. O., Kellerhouse, L. E. and Gehan, E. A. (1967). Plasmacyic myeloma: A sudy of relaionship of survival o various clinical manifesaions and anomalous proein ype in 112 paiens. American Journal of Medicine, 42,

13 REPRESENTING THE MEAN RESIDUAL LIFE 13 [8] Chaubey, Y. P. and Sen, P. K. (1999). On smooh esimaion of mean residual life. Journal of Saisical Planning and Inference, 75, [9] Chen, Z. (2). A new wo-parameer lifeime disribuion wih bahub shape or increasing failure rae funcion. Saisics & Probabiliy Leers, 49, [1] El-Arishy, S. M. (1993). Characerizaions of Weibull, Rayleigh and Maxwell disribuions using runcaed momens and order saisics. Egypian Journal of Saisics, 37 (2), [11] El-Arishy, S. M. (1995). Useful relaionships beween he log-survival funcion and runcaed momens, wih applicaions. Saisical Papers, 36, [12] Glaser, R. E. (198). Bahub and relaed failure rae characerizaions. Journal of he American Saisical Associaion, 75, [13] Hollander, W. and Proschan, F. (1975). Tess for mean residual life. Biomerika, 62, [14] Kodlin, D. (1967). A new response ime disribuion. Biomerics, 2, [15] Koz, S. and Shanbhag, D. N. (198). Some new approaches o probabiliy disribuions. Advances in Applied Probabiliy, 12, [16] McGill, J. I. (1992). The mulivariae hazard gradien and momens of he runcaed mulinormal disribuion. Communicaions in Saisics Theory and Mehods, 21 (11), [17] Meilijson, I. (1972). Limiing properies of he mean residual lifeime funcion. Annals of Mahemaical Saisics, 43 (1), [18] Muh, E. J. (1977). Reliabiliy models wih posiive memory derived from he mean residual life funcion, in Theory and Applicaions of Reliabiliy, (C. P. Tsokos and I. N. Shimi eds.), Academic Press, [19] Nair, U. N. and Sankaran, P. G. (1991). Characerizaion of he Pearson family of disribuions. IEEE Transacions on Reliabiliy, 4 (1), [2] Osaki, S. and Li, X. (1988). Characerizaions of gamma and negaive binomial disribuions. IEEE Transacions on Reliabiliy, 37 (4), [21] Ruiz, J. M. and Navarro, J. (1994). Characerizaion of disribuions by relaionships beween failure rae and mean residual life. IEEE Transacions on Reliabiliy, 43 (4), [22] Sen, A. and Bhaacharya, G. K. (1995). Inference procedures for he linear failure rae model. Journal of Saisical Planning and Inference, 46, [23] Shaked, M. (1978). Acceleraed life esing for a class of linear hazard rae ype disribuions. Technomerics, 2 (4), Deparmen of Mahemaics & Saisics, Universiy of Maine, 5752 Neville Hall Orono, Maine , U.S.A. address, Ramesh C. Gupa: rcgupa@maine.maine.edu address, David M. Bradley: bradley@mah.umaine.edu, dbradley@member.ams.org