Bias-corrected Pearson estimating functions for Taylor’s power law applied to benthic macrofauna data Bent Jørgensen, Clarice G.B. Demétrio, Erik Kristensen, Gary T. Banta, Hans Christian Petersen, Matthieu Delefosse To cite this version: Bent Jørgensen, Clarice G.B. Demétrio, Erik Kristensen, Gary T. Banta, Hans Christian Petersen, et al.. Bias-corrected Pearson estimating functions for Taylor’s power law applied to benthic macrofauna data. Statistics and Probability Letters, Elsevier, 2011, ฀10.1016/j.spl.2011.01.005฀. ฀hal-00746104฀ HAL Id: hal-00746104 https://hal.archives-ouvertes.fr/hal-00746104 Submitted on 27 Oct 2012 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L’archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d’enseignement et de recherche français ou étrangers, des laboratoires publics ou privés. Accepted Manuscript Bias-corrected Pearson estimating functions for Taylor’s power law applied to benthic macrofauna data Bent Jørgensen, Clarice G.B. Demétrio, Erik Kristensen, Gary T. Banta, Hans Christian Petersen, Matthieu Delefosse PII: DOI: Reference: S0167-7152(11)00012-5 10.1016/j.spl.2011.01.005 STAPRO 5880 To appear in: Statistics and Probability Letters Please cite this article as: Jørgensen, B., Demétrio, C.G.B., Kristensen, E., Banta, G.T., Petersen, H.C., Delefosse, M., Bias-corrected Pearson estimating functions for Taylor’s power law applied to benthic macrofauna data. Statistics and Probability Letters (2011), doi:10.1016/j.spl.2011.01.005 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. a∗ c ! a % a c d ' b ' # + ' d b & & ) * ", ' " * b " #$ "'-./ ) * ' ( # '0 2 # 1 ) * # 3 4566 " * 7 8 9 : * * ; < * , - * * j . , . . . , ni i Yij i ,...,k Yij * * a> b Yij # * i $ a bi , = 6?@6> 9 * CCC 7 =455C> 1 # =6 6> , =6?AB> , ∗ DE #= > 6 * =6?AC> =" =6??F C> 455AG 9 1 #H9 = 1 # 4565> b b : # =6 6> =6?@6> , * * 6?A@G $ # H (* b * = 6?ABG 455B> ( b > b J CA> I 6??C> * * * b = H * * * 6??C * 4 b =6?AB> * ( * 6?A6> = = b 4> ' 7 ( ni =6?A6> * * 9 * # #H J =455C> : , * 1 a b Yij J * : # # # # , 7 * 8 , H $ =455C> % H a =4565> BKC> b =' * ( a b * b =' * 3> 8 =' 8 , @> * , : * =6 6> . ( : n Yi ni n Yij Si ni − j 4 Yij − Y i j =4 6> * i Y i, Si a b i * * i ,...,k , a Si ǫi , b ǫi i * =4 4> * ni − Yi =4 4> % Si ni − * , * =6?A6> ǫi Yi 1 * b Si , 1 ni − * =4 4> i 1 * * , * ' * * # J @ * 8 * fi β , i ℓ* β * * fi * 1 * * * : # ' * : * 8 b a * i, a b i , =6 6> = 6??F * C> = Yij ∼ b 6??C> " i, a # , b * #H , * β 8 ' β ,..., k # b B a b a a β # =455C> . * J a b , =455?> # R : * b k ⊤ ' B b ∈ 1 ' , @ ( = * 6??F C> Y * N Y XN X B =4 B> 9 6L * I Y N 1 % * b N −b / * −b Xm b− −b ; * b− −b I Xm 8 b= * I Xm 9 6> Y - Y 1 Y > * 1 * * 9 J 4 b , 1 b Y I . 8 b * =4 B> 9 bD J b * 6 - b : I b 4 * , * * b b > I = . = = 9 9 6> I 4> - b * : > 7 * = 6??F b C> 8 = C 6??F * C> 9 = 4L > b 7 * a= 1 > . 7 1 * , * b * * * 6 * 8 H$ : Yij * =455C> =6 6> % H 9 Yi . i λ γ λ, b a ⊤ * 3 * =4565> 8 fi β i J β n k Yij − i xi eλ bi ψβ β, γ i ℓ* * * xi β ψβ β, γ , * k i * Ci fi , ni xi x⊤ i λ e bi Jβ . j ni b i * γ J : Yij , ∂Ci ∂γ Ci− Wiγ Ci− ∂ ∂ Ci H$ ∂ Ci Ci i γ n Yij − Wiγ i i n ∂ Ci ∂γ j Yij − Ci γ Ci J * * J : * i =B B> − Yij − #H # Yij − # * ∇γ ψ λ β, γ i =455C> * i J β 455C> ( i J * γ − Ci i j k J = * =455C> ψγ β, γ β, γ ψ γ β, γ ψβ β, γ ψ ∇ γ β β, γ β,γ # H$ Yi i =B 4> ∂γ ∂b k . β Wiγ * ∂λ 9 # J γ =B 6> b eλ Y i * eλ β ψβ β, γ ψb β, γ 8 H$ ψ γ β, γ J =455C> * γ J # γ # fi β ψ γ β, γ γ =B B> * 7 γ ψγ β, γ @ 9 ψγ β, γ H$ =455C> * β, γ γ ψ β β, γ J . ψγ β, γ * % =B C> H =4565> γ * 9 * : γ ψγ β, γ γ γ Jβ * * λ Jβ J− β Jβ ' =B B> Jβ =B 3> λ n Yij − i eλ bi k i j n k i Yij − i eλ bi j k ni − eλ n fi β n−ℓ i : ℓ − − n−ℓ =B @> − n−k =B F> λ * k Yij − i n−k b i j =B 3> −J β 1 ni a * k i Si b eλ Y i i a γ ψ γ β, γ − ψ λ β, γ n ni − β =B B> n k Yij − i eλ bi ψ b β, γ i j * − * * i ni eλ i b i * xi x⊤ i F =B ?> =B B> b Jβ − * i b − =B A> b b k Si Yi i * . b Jβ 8 =B 3> * ni Yi b eλ Y i =B 65> % k b Jβ J− β − i ni eλ b i ni eλ b i k − i k ni eλ − i i i − xi x⊤ i Jβ i − x⊤ i J β xi − x⊤ i J β xi b i =B 66> k Yi − i b n k i j n k i Yij − − − x⊤ i J β xi i eλ Si ni − i b i − x⊤ i J β xi b : =B 64> * λ − x⊤ i J β xi =B 64> i =B 6B> /ni b =B F> * =B 6B> * b ,...,k Yi J i J k ni − ∂ i γ i b i Yi ) i − − ni − b eλ Y i ni eλ i i j k # b k Yij − i eλ bi ψb β, γ λ b Ci Si − ∂γ Ci =B 6C> ni → ∞ i =:, β ,..., k ⊤ M< ' M < ' =B 6C> =B B> Si * b , Ci * k = * Ci / ni − , =B 6C> * ' b k> 3 J > k : = . =455B> * * A ' 4> , + * a =B A> , * ( * * * ! 8 * < * < J =B C> , # * * H$ * γ β, γ =455C> ψβ β, γ = J * *> β L γ * Sγ 3 ( ' β∗ β − J− β ψβ γ∗ γ − S γ ψγ =C 6> − , = , =C 4> * * > =C 4> =C 6> ψγ * β γ J * γ * * * : J λ * 8 Sγ =B 3> β,γ k ni Ci− − γ i γ γ =B A> , λ Sγ b − : ψ λ β, γ ni − ∂ =C B> * Ci ∂ ∂γ i Ci ∂γ * k i −n ni i k i k i − − ni ni i i * 1 * ' # * k Sγ ψλ - Sγ − T γ , , ψγ β, γ ∂Ci ∂Ci ∂γ ∂γ β : =C 4> ∇γ ψ γ β, γ * Sγ Sγ γ * * λ 1 , ψb b ) J 3 6A @ Tγ n−ℓ =B @> + ? b, b 1 Tγ 1 b b b,b − Jβ J− β Jβ Jβ − Jβ =B 6> =B 65> J− β = % k ni b,b i eλ i * 4565> b Jβ Jβ * Jβ * H xi x⊤ i b i b ( Jβ 1 =B 66> " # 8 * * −⊤ S− J− θ θ V θ Sθ Sθ θ ∇θ ψ θ θ * θ ψθ θ , θ * * −⊤ θ % Vθ Vθ * −S θ * γ 8 % H ( Vθ = Sθ Sθ − T θ * γ 7 =455C> J− β * =4565> * , Jβ, ( V βγ * Vγ * V β V βγ V γβ V γ V⊤ γβ Yij Vγ * V γβ k ⊤ ψ γij β, γ ψγij β, γ ' , J γθ J ψγij β, γ ψ⊤ βij β, γ i # γ − Sγ * Vγ * * J− θ n V γβ j : * θ n k Vγ Sγ * =B 3>> =C B> γ θ γ 7 i β γ * * =4565> =B 6> θ ψβ * Vθ Vβ * , , ( , ( H$ β * * J ψθ Jθ * * θ ) j , γ % H * ⊤ S γβ J − β S γβ S γβ J − β V βγ ⊤ * Vγ V γβ # S γβ 65 − V γβ J − β S γβ S γ =3 6> + =3 6> γ Jθ =3 6> ψγ β, γ * * , =B B> n k ∂ ∇β ψ γ β, γ θ γ γ ) : i Ci ∂γ j k ni − i k ni − ∂ k S γβ Ci Ci b ∂γ i ( Ci ∂γ i x⊤ i x⊤ i i i i * T γβ β γ γ,β Jβ J− β Jβ βm J− β − Jβ J− β γ β λ β λ,β Jβ J− β Jβ b k β Jβ J− β − Jβ β i k i ni λ e bi xi xi k b,β Jβ − βm ni xim −b eλ bi i − β nm −b eλ Ymb Ym xi β x⊤ i * xi xi 1 β J− β Jβ ⊤ − β m xi * k i λ J− β J− β 6 xm m −J β J − β Jβ J− β − Jβ ni bxim xi x⊤ i eλ bi − β J− β b,β Jβ J− β Jβ b − Ci i ∇β Ci ni b − Yij − ∇β θ Ci ∂ : , S γβ xi x⊤ i − i ni eλ nm b xm x⊤ m eλ Ymb 1 , i b i xi β * x⊤ i % * xi xi β xim Jβ ⊤ xm x⊤ m : S γβ T γβ S γβ − T γβ S γβ $ # + 9D =+9> 2 # 455A 9D =29> 66 655 9 'DN 9 * BL S * - Y 4A 4C * =5 9 # B5 = * > = = ! " > # * 4?O ?BO 6 @5 ◦ * L < = = :* = $ > > ! > # = 9 S * B > 63 ' : > 3 Y = : 64 65 > * 6 > ' k 63 44 6F 45 66 6B 4@ 4B 4B ! " + &." 6 C6 6 BB 6 @6 6 36 6 @? 6 3B 6 C3 6 F@ 6 3B 9D =+9> 6 C@ 6 B3 6 F? 6 @6 6 F4 6 F4 6 3C 6 ?@ 6 @A 6 33 6 BA 6 ?C 6 @A 6 F@ 6 F@ 6 @6 6 ?A 6 @? 6L " b &." P , '" 5 6B 5 6C 5 6? 5 5? 5 6B 5 5F 5 5? 5 5? 5 5? + 2 # &." 6 B4 6 C@ 6 C? 6 34 6 F6 6 B? 6 3@ 6 @A 6 @A k 63 6? 44 46 65 63 F 6A 44 9D =29> 6 C4 6 AA 6 F6 6 F5 6 ?C 6 @C 6 @F 6 ?A 6 AF 2 # # '" P P 9D 6 CA 6 AA 6 FA 6 FF 6 ?C 6 @3 6 F5 6 ?B 6 ?@ 9D '" 5 46 5 4? 5 6F 5 4C 5 65 5 6B 5 4@ 5 6C 5 63 =k P > # =ni > 6 b , , , = , # # # "J =B 6B>> − . − . , # a= 4> , # b I 5 AKF J * a 6 BA =4 B> 6 * 6 ?A * * * * - b +9 29 χ '" * J : b * J * * . , =X * p b* b b # . > 1 : , ( b : 8 b 6B + 9D =+9> &." '" 4F A4 6F @B 6B @3 5 F3 ?5 4@ @3 AC @F 5C 6 6F BA B6 @ @? B 3? 6 4A 4? ?C 65 55 ? 4C 5 F3 A 5F 3 4? 3 4F 5 F? ? 5@ 4 @C 4 ?@ 5 BA F 36 C C3 B FA 5 C@ B6 3F 4 35 4 ?? 5 FF 4@ 4F F 6F F AF 5 F4 ' ! " 4L " , P a # 2 # 9D =29> &." '" 34 6? 4? FC 4@ C4 6 4? A5 @5 3 5@ 3 3F 6 FC 3@ 43 A 5A F 3B 6 4B F@ 66 6C @@ 65 ?6 6 ?5 6B @F 4 36 4 AC 5 @4 6B 33 C ?C 3 ?C 5 3? A AC C F@ C ?C 5 F4 B? C3 6 @A C 56 6 4C A 3A 6 4F 5 ?6 6 B3 + 9D 2 # > ' 6 P 9D =&." P '" P % 8 * " 9- "' < ' < J = 1 > 2 & D D E F H , J. A. 6?A@ ' 433K43F , P. K. 455? & , Z., B # LQQ L K * * * , D Q Q # Q , P. 4565 + % ( 5@@664 , A., & F , , B. 6??F J , B., & K ! , B., M $ / 0 / , L R , S. J. 455C 0 / , J. R., & T 44BK4CB 1 . ". / 454 % , B. 4565 "; ( )**+ ,+-* ) , R., & J J I Q , I., & K , J. 455A 9 ' "' A?K6C4 , J L L H% D ?BK66C , M. 6??C - 6C * * K ! , W. S. 455C J , A. M., & I( , A. R. 455B ' @3K@A K ' , H., & C. , K. 455C ) 0 C@4KCFA P ! ' P *, J. N. 6?A6 ' / J # 2 # / & " / " ) ' * 43CK4@B , L. R. 6?@6 - T * , L. R., T * * T & * T * T , ' * " ' , G. D. 6??C 9 L 2 , L. P., & B , , , , R. P. 455B "; D 0 % / / 2 $" 6?6K45A L , H., L ! *, B. G., & W 2 P * 6?BK45? 1 * , R. A. J., W A56KA5C () FB4KFB3 !, I. P., & P , G. 6?AB - *, J. N. 6?AB S 3& ! , M. C. K. 6?AC , 4-567*8 " G. ., J. K., & R *, J. = . " 9 / 9 " L ( ' ( 63 * * @CK@A , > / : ' 0 9