
    M/Ph*-                    z   d Z ddlZddlZddlZddlmZ ddlm	Z	 dZ
dZdZdZdZ G d	 d
          Z G d d          Z G d d          Z G d d          Z G d d          Z G d d          Z G d d          Z G d de          Z G d de          Z G d d          Z G d de          Z G d d           Z G d! d"          Z G d# d$e          Z G d% d&e          Z G d' d(          Z G d) d*          Z G d+ d,          Z  G d- d.          Z! G d/ d0          Z" G d1 d2          Z#dS )3z
Results for test_glm.py.

Hard-coded from R or Stata.  Note that some of the remaining discrepancy vs.
Stata may be because Stata uses ML by default unless you specifically ask for
IRLS.
    N)add_constant)glm_test_resids            c                       e Zd ZdZd ZdS )LongleyzM
    Longley used for TestGlmGaussian

    Results are from Stata and R.
    c                    t          j        g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg          | _        d| _        t          j        g d          | _        t          j        g d          | _        d| _        d| _        d| _        d| _	        d| _
        d| _        d| _        d| _        d| _        d| _        d| _        d| _        t          j        g d          | _        d S ) N)ޭpp@r   r   r   r   )nWr   r   r   r   )2$G@r   r   r   r   )}աyr   r   r   r   )畃n[s@r   r   r   r   )y)or   r   r   r   )dr   r   r   r   )璄TW\*r   r   r   r   )o{',@r   r   r   r   )<6Nv|@r   r   r   r   )2u6hD1r   r   r   r   )Cr   r   r   r   ):5<bqcr   r   r   r   )UjUr   r   r   r   ){.{^u@r   r   r   r   )R@ir   r   r   r   i)g}⹭.@g|S7Vg<=n) gg|X=Q*gT@g
ףPIJ)gQ$:U@gl%?gsFdA?g'<;Um?gC@t?gBw|@g㥛h,+Ag8uIgm@gz?'-@gAk)AgB@g8uIg[g;Z^)A   	   gZ@g?      ?)gyR@g T7r @g ϖ[@g@g "@gy	2@g !@g@{#@g ĦK@g ɱt@gM@gDt?@g̨@g@@g ф@gO \F@)nparrayresidsnull_devianceparamsbseaic_R	aic_Statadeviancescalellf	bic_Statadf_modeldf_residchi2	prsquaredprsquared_cox_snellfittedvaluesselfs    l/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/statsmodels/genmod/tests/results/results_glm.py__init__zLongley.__init__   s   h) ) )                              ?      B 'h  -  -  - . . 8 , , , - - (
*)'
&&* &	#'  H &  &  &  ! !    N__name__
__module____qualname____doc__r4    r5   r3   r
   r
      s2         
E! E! E! E! E!r5   r
   c                       e Zd ZdZd ZdS )GaussianLogzC
    Uses generated data.  These results are from R and Stata.
    c                    t          j        g g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg d g d!g d"g d#g d$g d%g d&g d'g d(g d)g d*g d+g d,g d-g d.g d/g d0g d1g d2g d3g d4g d5g d6g d7g d8g d9g d:g d;g d<g d=g d>g d?g d@g dAg dBg dCg dDg dEg dFg dGg dHg dIg dJg dKg dLg dMg dNg dOg dPg dQg dRg dSg dTg dUg dVg dWg dXg dYg dZg d[g d\g d]g d^g d_g d`g dag dbg dcg dd          | _        de| _        t          j        g df          | _        t          j        g dg          | _        dh| _        di| _        dj| _        dk| _	        dl| _
        dm| _        dn| _        do| _        dp| _        t          j        g dq          | _        d S )rN)H0!5?r?   gzJ}L?r?   r?   )\uxcJ?r@   g#070a?r@   r@   )(M[Q3rA   g)8IrA   rA   )˪V?rB   gMW_m?rB   rB   )mRvbrC   gWDwrC   rC   )7N9R?rD   ǵ/f?rD   rD   )ƼqվrE   gpP rE   rE   )&:,B%rF   gɇE9rF   rF   )b%HrG   gS|\rG   rG   )L R"rH   glw4rH   rH   )*Ŏơ~GrI   g YrI   rI   )RirJ   g-Ѵ/rJ   rJ   )痀N5?rK   gDe?F?rK   rK   )hM=nrL   X
7krL   rL   )ͮ{+B?rN   g9mKNR?rN   rN   )h-6?rO   g'ZTE?rO   rO   )bHz?rP   gĭ(?rP   rP   )SIrQ   gˤrlXrQ   rQ   )5?rR   gwH1@	D?rR   rR   )4W_KrS   gॊXrS   rS   )?m1:rT   gYFrT   rT   )on="?rU   g1?74e/?rU   rU   )V0KrV   ge@WrV   rV   )<>rW   g/>rW   rW   )TrX   g	"rX   rX   ))H4"F?rY   gV%Q?rY   rY   )+|:?rZ   g>ckdC?rZ   rZ   )熺Ÿ́p>r[   gov>r[   r[   )mC{dJr\   g_--Sr\   r\   )(Kr]   g\ܐhSr]   r]   )\d=P?r^   g]2%V?r^   r^   )Q~3a@r_   g)yFr_   r_   )ٙQ?r`   grMdW?r`   r`   )sYTO?ra   g"S?ra   ra   )Q^?rb   g1HVb?rb   rb   )wS%%_rc   gǷwbrc   rc   )!J?rd   gEO?rd   rd   )*!P?re   gwۓR?re   re   )}	^Q?rf   g%FYR?rf   rf   )Mp9rg   gh<;rg   rg   ){QP?rh   gѵP?rh   rh   )ս_Fri   g34Fri   ri   )]c;Jrj   gcM*Jrj   rj   )~bU]rk   g?߾\rk   rk   )xC?rl   笢zB?rl   rl   )-0?rn   gf-?rn   rn   )lxz,CLro   gX$Hro   ro   )_S$6rp   gSb3rp   rp   )rM   rM   g<PrM   rM   )%F?rq   g*ЬA?rq   rq   )Tb8rr   gOk2rr   rr   )z rs   giZrs   rs   )3:{mZrt   gmrSrt   rt   )%LbZS1?ru   g}^(?ru   ru   )珠xs]]?rv   gyABT?rv   rv   )gej"rw   g4;rw   rw   )}!3rx   g]FV#`)rx   rx   )G[hNry   gyCi7Cry   ry   )_M?rz   g_ .B?rz   rz   )Ճ8řUr{   g}Ir{   r{   )ot{!r|   g[mDr|   r|   )M?r}   gf%}@?r}   r}   ){gUId/r~   gO r~   r~   )퐬_pG?r   gDXP8?r   r   )صݒ?r   g6I>r   r   )熥&fb_g?r   g>+NV?r   r   )ηO؏^>?r   gi)',?r   r   )>Ph2?r   g滥 ?r   r   )玽љj?r   gjۋ>r   r   )KEb7r   gZa/$r   r   )WmURr   gDs)">r   r   )G6uU?r   g! @?r   r   )<kWr   g!8)4Ar   r   )f:*0&r   govr   r   );,G#Pr   gDXP8r   r   )a4oX?r   g@?r   r   )	Br   gcKEr   r   )8ӅXA?r   g:m/'?r   r   )Rdar   g}3֣Er   r   )\+-eq>r   g)>r   r   )GR?r   g{эU5?r   r   )]0??r   g+ٱ!?r   r   )%?=Ar   gO"r   r   )]a&r   gŭr   r   )7h>^r   g1cH>r   r   )Vrz(+?r   g-C6
?r   r   )L!FLr   g D*r   r   ){r   gՃ8řr   r   )ƍA?r   gM?r   r   )3$Ir   gm[$r   r   )nq]r   g{/h7r   r   )xr+R?r   g;s	+?r   r   )?.q%Br   g
^r   r   )ĭH?r   g,gJF)|!?r   r   )焛*ø;?r   gFƴ?r   r   )f@O4r   g^W*^r   r   )@'SEr   g<ړr   r   )>U?r   '	*?r   r   ).9.Ir   g**,
r   r   )[v-e[?r   g{C/?r   r   gXL@)gnO?gT-xg¶8)g^n"?gu?]0K9>gl7%>g^F?ghD&g7K?g{m>g^F_@gb9{r   a   gVA)dg׾@gs[P@gzӈ@g~v@gƨV\
@gm'Ġ@g+HQ7@g~D@g;-i@gp(f@g٠@g7@>@gPn @gP:0| @g#$= @gjz̹?g7?g1+?gea?gqWM?g?g/Y?gN+?g;>?g=Ќ;i?gr+O?g^2O,?g7mˉ?g?g/,d?g̉?gtx1xB?g9?gj_PQ*?g䨥Y?gBj%?g<4?gP?gZJnEU?gE3!?gRR??gt/r2?gr}?g#?gܓ0?gd	?g\z?g;u!/1A?gSmlv?gge;?g0Ϗ?g)![1?g4Zx?g<1?g#H|?g;Vh?gXݰ>?gk?g&|?gҫXT?gwԜI?g|B?g|^%?gV?g
J?g*?g?gF?g[~̝"?g3sy3?g(nK?gYAj?g+?gsڷ?g>}1?g܎cw!?g%^?gqe?gOuC?gXmؐ5?g&hR/?g@W ?gR<?gwR?gGTu?ge*?g3?gwv?g5o?gcmt?g7وP?gwWS?g-/G?g?g>r2?g4[?g/:?vO?g4qS?gОӄ?g@S
?r   r    r!   r"   r#   r$   r%   r&   r'   r(   r)   r*   r+   r,   r-   r0   r1   s    r3   r4   zGaussianLog.__init__j   s	   h H 0? ? ?H 0/ / /H 0
0 0 0H 0/ / /H 00 0 0H 0/ / /H 00 0 0H 00 0 0H 0"0 0 0#H 0&0 0 0'H 0*0 0 0+H 0.0 0 0/H 02/ / /3H 060 0 07H 0:/ / /;H 0>/ / /?H 0B/ / /CH 0F0 0 0GH 0J/ / /KH 0N0 0 0OH 0R0 0 0SH 0V/ / /WH 0Z0 0 0[H 0^/ / /_H 0b0 0 0cH 0f/ / /gH 0j/ / /kH 0n/ / /oH 0r0 0 0sH 0v0 0 0wH 0z/ / /{H 0~0 0 0H 0B/ / /CH 0F/ / /GH 0J/ / /KH 0N0 0 0OH 0R/ / /SH 0V/ / /WH 0Z/ / /[H 0^0 0 0_H 0b/ / /cH 0f0 0 0gH 0j0 0 0kH 0n0 0 0oH 0r/ / /sH 0v/ / /wH 0z0 0 0{H 0~0 0 0H 0B0 0 0CH 0F/ / /GH 0J0 0 0KH 0N0 0 0OH 0R0 0 0SH 0V/ / /WH 0Z/ / /[H 0^0 0 0_H 0b0 0 0cH 0f0 0 0gH 0j/ / /kH 0n0 0 0oH 0r0 0 0sH 0v/ / /wH 0z0 0 0{H 0~/ / /H 0B/ / /CH 0F/ / /GH 0J/ / /KH 0N/ / /OH 0R/ / /SH 0V0 0 0WH 0Z0 0 0[H 0^/ / /_H 0b0 0 0cH 0f0 0 0gH 0j0 0 0kH 0n/ / /oH 0r0 0 0sH 0v/ / /wH 0z0 0 0{H 0~/ / /H 0B/ / /CH 0F/ / /GH 0J0 0 0KH 0N0 0 0OH 0R0 0 0SH 0V/ / /WH 0Z0 0 0[H 0^0 0 0_H 0b/ / /cH 0f0 0 0gH 0j0 0 0kH 0n/ / /oH 0r0 0 0sH 0v/ / /wH 0z/ / /{H 0~0 0 0H 0B0 0 0CH 0F/ / /GH 0J0 0 0KH 0N/ / /OH 0 H1 H1T 0h  >  >  > ? ?8LLLMM(
+,+
$+%	H "&6 "&6 "&6 "7 "7r5   Nr6   r;   r5   r3   r=   r=   f   s2         {7 {7 {7 {7 {7r5   r=   c                       e Zd ZdZd ZdS )GaussianInversezG
    This test uses generated data.  Results are from R and Stata.
    c                    t          j        g g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg d g d!g d"g d#g d$g d%g d&g d'g d(g d)g d*g d+g d,g d-g d.g d/g d0g d1g d2g d3g d4g d5g d6g d7g d8g d9g d:g d;g d<g d=g d>g d?g d@g dAg dBg dCg dDg dEg dFg dGg dHg dIg dJg dKg dLg dMg dNg dOg dPg dQg dRg dSg dTg dUg dVg dWg dXg dYg dZg d[g d\g d]g d^g d_g d`g dag dbg dcg dd          | _        de| _        t          j        g df          | _        t          j        g dg          | _        dh| _        di| _        dj| _        dk| _	        dl| _
        dm| _        dn| _        do| _        dp| _        t          j        g dq          | _        d S )rN)"R@r   gRt@?r   r   )_Q+r   r   r   r   )i :w&r   gŷ$?r   r   )_?r   g?'[r   r   )\{QFr   rm   r   r   )InsS?r   g[RjOr   r   )T,Xr   gHaR?r   r   )yT5?r   gP G-r   r   )Xr4r   g^!+?r   r   )缅d8?r   gtI.r   r   )祌??r   gm)2r   r   )ڔu>Wr   gC I?r   r   )+ٱ?r   gavtr   r   )]6:?r   g..*r   r   )p7r   gk$	%?r   r   )ſF}3r   g ?r   r   )Nr   gQݷ8?r   r   )^gS?r   gf(ϼ<r   r   )<*6F?r   g`Z8Z/r   r   )M~0Ur   g;s	;?r   r   )Pl-?r   gZ1r   r   )qڝhIQ?r   g:O_3r   r   )(Yi0?r   g>$Er   r   )[+:?r   g?yr   r   ).ԫv[r   g3:?r   r   )vqvRH?r   gQq)%r   r   )1QaJ+?r   g34r   r   )X=Jr   g²lh$?r   r   )uC?r   gF%r   r   )n~hP?r   gAȘ&r   r   )6-2?r   g/bDGr   r   )9y@r   g=?r   r   ),-~1?r   g
Fr   r   )ɵWՁb*?r   gb:r   r   )>Җ)Sr   gnBA#?r   r   )D &_Fr   g`=I?r   r   )  B?r   gx+Nr   r   )"/l(@?r   g!
r   r   )_L5?r   gʹ r   r   )F<Ir   gU?r   r   )(>r   g&Jw¾r   r   )r$r   g-C6>r   r   )羒,?r   gy;>1;r   r   )u5hFr   g&(4
?r   r   )bHzIr   gb:?r   r   )pCP?r   gt><Kr   r   )a2U0*S?r   g(M[Qr   r   )aHIr   g_eG|?r   r   )}6r   gZQf>r   r   )/ز!S?r   gMbr   r   );TYr   gr?r   r   )玷M?r   g/\cr   r   )y`[>r   gTM>r   r   )Z/-Sr   g
^
?r   r   )V8Q?r   gPr   r   )fLg1r   gkNu>r   r   )<ܴIr   gtR#>r   r   )纜V Y?r   gÔ/r   r   )煳[d8^r   g+,?r   r   )s'[?r   g<ړr   r   )!C}4r   g$z>r   r   )Z_&P?r   gAr   r   )1{vA?r   gs>r   r   )Dėr   gۢ>r   r   )F+R%Ir   gv!>r   r   )R[SF?r   gr   r   )% "r   g@cG_>r   r   )
OmI?r   gYVPhr   r   )Ջ)0?r   g~վr   r   )3kB?r   g,r   r   ),ZgH?r   gtR#r   r   )GW:Kr   g
h`?>r   r   )4F&;ar   ggg?r   r   )	hǲ)Mr   gM >r   r   )痰Q?r   g	\?r   r   )kB](r   g<>r   r   )ŭ?r   gBY[%r   r   )XCQ?r   giUMur   r   )SV#M?r   g]FV#`r   r   )K8,r   g`>r   r   )L¶Zr   gOne.>r   r   )ν,,?r   g};hžr   r   )Xf=2r   g>r   r   )<-?p'@r   g34>r   r   )VF?(r   gJ֗[>r   r   )yxρIr   g!>@>r   r   )OЁ@I?r   gl(1,߾r   r   )kZs2?r  gqdqAƾr  r  )rnOQ?r  g{fdr  r  )+4?r  gʔA{ƾr  r  )pDkLr  g,z>r  r  )@I(r  gTqs*>r  r  )3G5?r  g!+lěNžr  r  )7 >?r  gΗ;r  r  )ߏ61r  g~>r  r  )I	fU?r	  g8{r	  r	  )k0h|Ur
  g?̔>r
  r
  )B^j7?r  gtݛ?,r  r  )Cr  gP3bn>r  r  )'ZTEr  g>>r  r  gHj_?<@)gڶAP?g7ame?gSgP?)g9=?g2,]?g>q>J>g68}gp`Ȁ&gi?gVŨ>g68}@geg9{r   r   gγuYCDA)dgp:;?gT?gAMm?gV?gH\ךs1?g&Kq?g&x?g&.?g>.%?gkAfa?g?ggpB?gQ ?gm1Sf?gݯ?gy=?g 	O?gq::?gWn?g%c?g?gS3?g&&|?gy$I?ghd$?g|6%!?gJH'?g97?geP?gwcřq?gQ?gx*?gC?giI?gԮϽ?g]?g<9?gڋea?ggo?g^d?gc=?g(7H?gr??g&MB?gy:\?gNO?gl X?g~,?gh[?g<8?gr?gԵ?g"?g=~J?gQ?gYV1#?g:ѼZ?gzC?g.,?gaq?gp?g}ʵ?g*B^	?g^?gZ?gy?g8+?g-?g~qS?gZ@?gD?g*_'?gap?gJ\n?gN~?gd
?gNU\ھ?gl8?g=?g倸?gCp?gήR?g#U?gQ[؄ͺ?g?I?gvvɹ?gT~L?gSuӸ?gl]?gMZP?gظz?g"?gpgϣ?g~q<?g@7׵?g/{~u?g%O?gtف)?gJENW^?g=?r   r1   s    r3   r4   zGaussianInverse.__init__l  s	   h H 1@ @ @H 10 0 0H 1
0 0 0H 1/ / /H 10 0 0H 1/ / /H 10 0 0H 1/ / /H 1"0 0 0#H 1&/ / /'H 1*/ / /+H 1.0 0 0/H 12/ / /3H 16/ / /7H 1:0 0 0;H 1>0 0 0?H 1B0 0 0CH 1F/ / /GH 1J/ / /KH 1N0 0 0OH 1R/ / /SH 1V/ / /WH 1Z/ / /[H 1^/ / /_H 1b0 0 0cH 1f/ / /gH 1j/ / /kH 1n0 0 0oH 1r/ / /sH 1v/ / /wH 1z/ / /{H 1~0 0 0H 1B/ / /CH 1F/ / /GH 1J0 0 0KH 1N0 0 0OH 1R/ / /SH 1V/ / /WH 1Z/ / /[H 1^0 0 0_H 1b/ / /cH 1f0 0 0gH 1j/ / /kH 1n0 0 0oH 1r0 0 0sH 1v/ / /wH 1z/ / /{H 1~0 0 0H 1B0 0 0CH 1F/ / /GH 1J0 0 0KH 1N/ / /OH 1R0 0 0SH 1V0 0 0WH 1Z/ / /[H 1^0 0 0_H 1b0 0 0cH 1f/ / /gH 1j0 0 0kH 1n/ / /oH 1r0 0 0sH 1v/ / /wH 1z/ / /{H 1~0 0 0H 1B0 0 0CH 1F/ / /GH 1J0 0 0KH 1N/ / /OH 1R/ / /SH 1V/ / /WH 1Z/ / /[H 1^0 0 0_H 1b0 0 0cH 1f0 0 0gH 1j/ / /kH 1n0 0 0oH 1r/ / /sH 1v/ / /wH 1z/ / /{H 1~0 0 0H 1B0 0 0CH 1F/ / /GH 1J0 0 0KH 1N0 0 0OH 1R0 0 0SH 1V0 0 0WH 1Z/ / /[H 1^/ / /_H 1b/ / /cH 1f/ / /gH 1j0 0 0kH 1n0 0 0oH 1r/ / /sH 1v/ / /wH 1z0 0 0{H 1~/ / /H 1B0 0 0CH 1F/ / /GH 1J0 0 0KH 1N0 0 0OH 1 H2 H2R -hEEEFF8LLLMM(
+.+
%+%	H &L &L &L M Mr5   Nr6   r;   r5   r3   r   r   h  s7         jM jM jM jM jMr5   r   c                       e Zd ZdZd ZdS )Star98z0
    Star98 class used with TestGlmBinomial
    c                 >   d| _         d| _        d| _        d| _        d| _        d| _        d| _        d| _        d	| _        d
| _	        d| _
        d| _        d| _        d| _        d| _        d| _        t           j        | _        t'          j        g d          | _        d S )N)gP7gJY'S?gHMp,gBK)gURI?g:y?g<?g;<gR<aտg	 ſg++#t?gH_SmgR,Ԍg̉gpgL#' pgz?g02?g㴩w?gtg-?g8%mbg0ȫ@)g*|jp<?g0t?C?gp@Y]H?g$"o<?g[j?gOM09A?g?g~9F?g"i^?g(*<?gAS-$T?g7J-?gܯ[4_?gm:z ??gBO?gi?geT\~?g,ԤXX?g*V>goe6?gJU?g+@i.  g\ݯ@i     g$U9@g[r3@g)ԔF@gnU9mgm9mr   gKק@g0'?(/  g[}?gm{?gj?gޝ$Y?g$?g,:a?gO?gPٸ?gba?g_vj.7?gz?gen:ٸ?gtr}?g\آE?g$:`?gh#?gJM?g05)?gF?g?gfl
\E?g̈́?g^?gdr?g[?gVjb?g=ș?g]o?gb[0?g?gR5?gTWOU?g-?gf&C?gIQ?g7D&*e?g_D?gc?gTk]g?g=?g?g?AtQ?gl%?g J]?g?gX.?gtI?g?g$o?g:Vg4?gS7?g;?gֱ#?g_F?gk?gpB?gS#%?gdĪ?g6{?g=WY?gH{t?g\K:?gЄ`4?g0`?g'wLyK?gܦ?gq]P?g<?g'?gV1?gϰ,C?gp)^?g!I?g.yݦ?g$v?g^+;?gE9?gyp ?gU?g
k߳`?gpn?gwE?gpsF?g<W"?gMIR?gŃ?gDě?g.@6$Ik?go?g֤?g0?gL?gD0P?gRi?grO?gM)?g??gW-~?g"?gI@T3S?g9=?gZq?gs?gtH'?gX<w)?gٶz?gRA?g1?g@߾?gY>?gR<EID?g?gѿ?gVWp5?gbѶ?g&3X?g_w?ge;:	 ?g
>,?g$Wh?g=]?g76?g8H?gMD?gkQُ?g}Up?gpzӀ?g
;4?gN/D?g)Y?g:a_?gY	?g滄?g	1(U?gբO?g-5;?gj
*L?g V,L?g^oc?gd"<?gwv?g!|@?g=RfZM?goA?gw;3?gL%?gfgn?g{skʼ?gSix?gHt?gêI	D?gwԘ?g,c)>?gV;/!?gK?ge?gfT2c?g@Q_y?g8p\cc?g´!?giW[q?gY.?gVC?g#?gU=?gG<?gD Ւ?g;tYФ?g|El?gNWk?gh$?g
Ϯ?g]?g&.o?gZN?gi>X?gL@]?gZ<?gޣs?gs$~b2?g>?gE^?gꬻ?ga+?gH+9?gsJh?g?gy?g0uB?gOz?g^`?g0B/?gJZ?g|;?gps?g	n;?g܃?gdqʍ?gf =G?gP6?g. ?gԇ?ggBD?gL9?g&	?gЙ0?gw7?gL?goVH?g W?g{h?g	[???gÀFe;?g1Iu?guh`"?g@}?gZAuW?g	v?g?XH?g"dN?gƁ?gd,?gjƛ_?gU7i?g]
U?g7= ?g(),ȁ?gQ@p	?g72?gb?g>@1?gJ-M?g9e?g=|6?gt6?gL0k?g&*L!?geWG?gI\Di?gٺim?g,V?g{i?gIB?gGp1?gƣ?g?guo=O?g@+r?gBrt?g&L >N?gP??gp	~#ӄ?gM4?gݎx8?g+_a?g>O1?gXxPW?g|ѢgHO?g]C?gw ?gK&?g3?g%7Y?gXh?gC^|?gr?g۪=l?g;?gnH?g7]+?g7Tq?g>}_cE?g2R?g|z?g?g9?g՝F?g`&5r?g~J?gh$ ?g<&1?g5?g	v?gҴv?g d?gB~pg?g~?gܸ?g^2^5?gf,&?gd  W?g>?g+ ?g}p^q?gA?gBv?g}??g`9	'?g(<?gUl]?g<U?gHȀC?g{9?)r#   r$   r"   df_nullr'   r,   r+   r%   r&   r*   r)   	llf_Statar(   pearson_chi2r.   r/   r   star98_residsr!   r   r    r0   r1   s    r3   r4   zStar98.__init__]  s    ) +%"
**$+
'#) %3H =&2 =&2 =&2 =3 =3r5   Nr6   r;   r5   r3   r  r  Y  s2         [3 [3 [3 [3 [3r5   r  c                       e Zd ZdZd ZdS )Lbwz\
    The LBW data can be found here

    https://www.stata-press.com/data/r9/rmain.html
    c                    t           j                            t           j                            t           j                            t
                              d          }t          j        |          }t          j        |j	        ddt                    }t          j        ||gd          }|j        | _        |g d         }t          |d          | _        d	| _        d
| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        t<          j        | _         tC          j"        g d          | _#        d S )Nzstata_lbw_glm.csvraceFprefix
drop_firstdtyper   )axis)agelwt
race_black
race_othersmokeptlhtuiprepend)	g;V/gBdgһ_3?gk'?go
?gDV?ggE/Q?gԚE?g.*?)	gB?g0^|?g~)O?g?g!?g6ʓ(?g?d!?g9>%af?gPE?gA$Unk@geޓ?gU.i@gU.Yg$ <>9@g0)Um@g5@"0         g
8f@)gZ%?g`?g\v?gA@ ?g4E?ghRa?g"#S?gC?g_(3R-?gC5?gԗJ':?g_	 
?g]8J!?gr{?ri?r,  gRS?go/)7?gn?g?g9?g|}?g#L?gWL:m?gybx?g'k5K?ga-?gal?gt)*V?tܩ?r-  g]?gQw36?gyeJ?g˒2?gK?gQ-d&?gLd?gRگ?g[?g&v?gM/]E?g_7?g˽Wڲ?`5Q?r.  g(6Il-?g?g!0A3?g=՝?g?gmB?$;?gq?gT?g k\(?g*ѫdt?g`8$9	?g6 ?E7?r0  gi^?gxJ?g?giG7'?gg=M{?g?>Z?g] T?g\tL?g凌f?gme?g;>?gМq"?gTi0	?7?g]''Ƹ?gx틶?gcT?g'a^)?g'_[?g`p?g굨?g~M'?g^-?g;p w?gMM|?gnQ5?g-#qf&?gH ?gJֿ?g)eRp?glh?gkj?gSYV#?r1  g'na?g=$cRİ?-b1?r2  g6$2ô?gfN?gb\?g?gQ+%?g<?g_6u8?g,3M;?g͓)0?gΕ{`K8?g^ 
ú?gD?g- n?go?g:"`?gSX9?gΙ?g;@t?g's7?g??gk?g WX?gm-و?gwic?go]R1?g|ܹ?g'u?gji?g<Jg?g`~?gDW?gS4?gqج?gS9?ggB?gL?gA>?gks?gM]?gpjA?gϾ.?gEUM?gwR?g,uW?gB9?gv+jl?g\&?gb+&>?gߩQ?g* E?g|5d?g-?g`ٳ?g.?gg?g 
fL?gִ ΢?g8G?g+F?gjjӂ?g{j/?g! V?g&s$~?gd?gH֧=.?g3?g7/+u?g
X?g2NZ?gs?ghEJ?gCNC?gB9?gW24?g<s?goWhֈ?gWb?r/  g<?gZ:Q??gHGA	W?gU?g+U?g,?g@iT?g$?g :o)?g0Xb?g|3[	?g.^w_:?)$ospathjoindirnameabspath__file__pdread_csvget_dummiesr  floatconcatlowendogr   exogr#   r$   r%   r&   r'   r(   r)   r-   r"   r*   r,   r+   r  r  r   
lbw_residsr!   r   r    r0   r2   filenamedatadummiesdesigns        r3   r4   zLbw.__init__  sU   7<<0I0I J J 35 5 {8$$.6e',. . .y$q111X
 * * * + 777	"& #
*)
%%	-+,%0H &&> &&> &&> &? &?r5   Nr6   r;   r5   r3   r  r    s2         
I? I? I? I? I?r5   r  c                       e Zd ZdZd ZdS )Scotvotez2
    Scotvot class is used with TestGlmGamma.
    c                 0   d| _         d| _        d| _        d| _        d| _        d| _        d| _        d| _        d	| _        d
| _	        d| _
        d| _        d| _        d| _        d| _        t          j        | _        t%          j        g d          | _        d S )N)g~Ӏ
?g
~`?g`gC)\BR?g[WgZy@g{Q\ľgR1)g: >ggioA?g6Oo>gO$F?g[aј>gO}/?gN>gcj?gYw'?   gە_?      gjQ3Nf@g6oq%@gsTg*ܹqdgϬwH\m?gJ0?g~jt?g2Y?) gH˻L@gyJ@g? HI@g{%>*M@g̝Q@g2gqL@g\fP@g @.P@gKL@gxeQO@gzJ@gȹN@gPM 3P@g$ܼO@gq6YN@g+qR@g*zHL@gEWu'R@g٩=P@gH^w!J@g;!	P@gQ@gfF@g1؉rwK@gqP@gZ6,J@g!БL@gγM@gr;MIP@gVƯ>N@gmhtR@gg&^gQ@)r#   r$   r"   r  r'   r,   r+   r%   r&   r*   r)   r(   r  r.   r/   r   scotvote_residsr!   r   r    r0   r1   s    r3   r4   zScotvote.__init__  s    JD &&%
!" 
'#* %5H &C &C &C D Dr5   Nr6   r;   r5   r3   rH  rH    s2         D D D D Dr5   rH  c                       e Zd ZdZd ZdS )Cancerz`
    The Cancer data can be found here

    https://www.stata-press.com/data/r10/rmain.html
    c                    t           j                            t           j                            t           j                            t
                              d          }t          j        |          }|j        | _	        t          j
        t          j        |j        d          d          }t          j        |j        |f                              t"                    }t%          |d          | _        d S )Nzstata_cancer_glm.csvcategory)r  T)r  Fr'  )r3  r4  r5  r6  r7  r8  r9  r:  	studytimer?  r;  Seriesdrugr   column_stackr  astyper<  r   r@  rB  s        r3   r4   zCancer.__init__8  s    7<<0I0I J J 68 8{8$$^
.49J!G!G!G,02 2 2$(G!455<<UCC 777			r5   Nr6   r;   r5   r3   rO  rO  2  -         
8 8 8 8 8r5   rO  c                   "     e Zd ZdZ fdZ xZS )	CancerLogz+
    CancerLog is used TestGlmGammaLog
    c                 (   t                                                       t          j        g g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg d g d!g d"g d#g d$g d%g d&g d'g d(g d)g d*g d+g d,g d-g d.g d/g d0          | _        d1| _        t          j        g d2          | _        t          j        g d3          | _        d4| _        d5| _	        d6| _
        d7| _        d8| _        d9| _        d:| _        d;| _        d<| _        t          j        g d=          | _        d S )>Ng27{Hg{9yQg_FҞCg`gr"g [gu? ?gx~:gfgE.8g,4k_gHߠ=g6XDgOggyܭgҾ^gB&9;g&RgFg\g"N'g6FPegZ_$	Cg\,JgCg蹅DͿgklexϿgg;#5Ͽg4gM3ȿgBʿg~W[gbȠʿg8Q?gv2=ֿg=yX5ڿgoi5g '*ڿg#ӡgXT$[g ﭿgLgﭿg^a=$߿g;?g2g?g@"?g|g?g(@ ?g
%U ˿gV5RͿg;%5giTNͿgE4~5gh>ӿg ?rֿg7dDg_q^Rkֿgfg3P>タ;r30gR0rh  g)g`^6giR
?g}?gOz12@g,J	z?g @gy)?gE?gص@g˂?g?gؼm1?gT?gʡA@g9}?g9@gWm?gS	hy?g>xEG@gn<՜r?g)yu@gh}?gW\wN?g_LIO@g-4?gL~@g9?q`?ga/?gK+\@gjkD0?gL
F%+@g!Va@?giH?goM`@g{?gm~)@g[Nֿg҆hٿgc('U=gVٿg+Lk	g=#ٿg
Ϳ޿g,+MJADgG6Ws޿g`JgBk}߿g&p{gwT\WgS66GggXg<"ٿged_hd޿g(rVgB:}Q޿gנ/IgB0X{gA$gPkw`kgLϓg =%g]wgHgjgZ_$E)gpUjgS[ g	E-n?gM*?gt	1@gv?gaL{?g5,vοgWnпgߡ(AWg3kпg.g?gaR||Bv?gjWVDO@gDPg?gY(@gdmlvgu䚧ߢg=$P,g{Ӗߢgqt濩gyS4?g_FE?gp:ǲP@gBus?gKu/@gu?gWJ;?gEP\@g1Z웠2?gp @g z*L?7*?gDR*$@r}  gǞP?gF?g5:L?gH|o@g4y;?gE>+@g2"rg_g	c~gh?RDUg_Ahj3gJsugjg?ߊ{g%T=Og?o*Ra0gZZ Ŀgqx(ƿgSrPg+]&ƿg'󏾉
g!rhmۿgag<,Ԛ}g1kVg׻,gk?gP'?ghC@gPYJ?g+~w?gǮY˿g=(Ϳg$jgeDXCͿgm},g3452 ֐g&" %r  g]U"ڿg7-:?g#?g1殩[@gP+?gn@g~nȿg51ʿgH}omghI3̿ʿgq g/ CE?g7t?g+ebg@gnr?q?gzM@gֲx?gy]?g`"p@gDP?gW])@g-qU?g~J
?gH}8q@g}}m?g)*.@gcz?^kz?gM*k@r  g.hk(?g?g}i?glV}
r@gPI?gTH9#@˜;@)g_gE.?a?gXln?g:׃@)gm:,?gq5?gjO?g)?gx[>t@g_YǉK@g ,0@g}Z?gUeg4Ecr   ,   g'(<cB@)0gj5#@gwʯ@gGY\oݭ@h=nM$@- @gӚIH@gM7@`	@r  r  r  羯8'@g"^4&@gEg!@r  gb2n@gMIs;%@r  g8\7: @r  k"@gÛ[%$@婻p+@gE$.@g(I04@g`fL(@g+v
'@煭·3@r  r  gFe/@r  g:4~6@g]Z2@Hnj9@g5e:@MDQ4@g^*f@@gh86@g?@r  gUм5=7@<bA@g.M8@r  g=~NŶ2@r  gR=@)superr4   r   r    r!   r"   r#   r$   r%   r&   r'   r(   r)   r*   r+   r,   r-   r0   r2   	__class__s    r3   r4   zCancerLog.__init__G  s   h ` 00 0 0` 00 0 0` 0
0 0 0` 00 0 0` 00 0 0` 00 0 0` 00 0 0` 00 0 0` 0"0 0 0#` 0&/ / /'` 0*0 0 0+` 0.0 0 0/` 020 0 03` 06/ / /7` 0:/ / /;` 0>/ / /?` 0B/ / /C` 0F/ / /G` 0J/ / /K` 0N/ / /O` 0R0 0 0S` 0V0 0 0W` 0Z0 0 0[` 0^0 0 0_` 0b0 0 0c` 0f0 0 0g` 0j/ / /k` 0n0 0 0o` 0r/ / /s` 0v0 0 0w` 0z/ / /{` 0~/ / /` 0B/ / /C` 0F/ / /G` 0J0 0 0K` 0N0 0 0O` 0R0 0 0S` 0V0 0 0W` 0Z/ / /[` 0^0 0 0_` 0b0 0 0c` 0f/ / /g` 0j0 0 0k` 0n/ / /o` 0r/ / /s` 0v/ / /w` 0z/ / /{` 0~/ / /` 0 `1 `1B /h  =  =  = > >8NNNOO'
**(
 &+%	H &C &C &C D Dr5   r7   r8   r9   r:   r4   __classcell__r  s   @r3   rY  rY  C  sT         BD BD BD BD BD BD BD BD BDr5   rY  c                   "     e Zd ZdZ fdZ xZS )CancerIdentityz:
    CancerIdentity is used with TestGlmGammaIdentity
    c                 (   t                                                       t          j        g g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg dg d g d!g d"g d#g d$g d%g d&g d'g d(g d)g d*g d+g d,g d-g d.g d/g d0          | _        t          j        g d1          | _        t          j        g d2          | _        d3| _        d4| _        d5| _	        d6| _
        d7| _        d8| _        d9| _        d:| _        d;| _        d<| _        t          j        g d=          | _        d S )>Nr[  r\  r]  r^  r_  r`  ra  rb  rc  rd  re  rf  rg  ri  rj  rk  rl  rm  rn  ro  rp  rq  rr  rs  rt  ru  rv  rw  rx  ry  rz  r{  r|  r~  r  r  r  r  r  r  r  r  r  r  r  r  r  r  )g.g*nP@g240@g׋{C@)g}(?gΑy@g.1@gl|b @gFAiFt@gX<@gW0.@gB?r  geeg͝gcr   r  g1=2}I@)0g<@g@X?@gz/"@T	&@K!@g\X"@g=~@oH@r  r  r  v/N)@gq9;(@gm=;"@r  g]@gc$('@r  g-J\ @r  ,s%)"@gIqe%@`,@g]t.@gV+ 3@g9])@g:(J(@Cdn2@r  r  g |c/@r  gUo34@gr1@鋭9@g(عH%,:@V6@g{<@gNUF8@gO<R<@r  gu8@筮d=@g59@r  gag5@r  gUʄ?;@)r  r4   r   r    r!   r#   r$   r%   r&   r'   r(   r"   r)   r*   r+   r,   r-   r0   r  s    r3   r4   zCancerIdentity.__init__  s(   h ` 00 0 0` 00 0 0` 0
0 0 0` 00 0 0` 00 0 0` 00 0 0` 00 0 0` 00 0 0` 0"0 0 0#` 0&/ / /'` 0*0 0 0+` 0.0 0 0/` 020 0 03` 06/ / /7` 0:/ / /;` 0>/ / /?` 0B/ / /C` 0F/ / /G` 0J/ / /K` 0N/ / /O` 0R0 0 0S` 0V0 0 0W` 0Z0 0 0[` 0^0 0 0_` 0b0 0 0c` 0f0 0 0g` 0j/ / /k` 0n0 0 0o` 0r/ / /s` 0v0 0 0w` 0z/ / /{` 0~/ / /` 0B/ / /C` 0F/ / /G` 0J0 0 0K` 0N0 0 0O` 0R0 0 0S` 0V0 0 0W` 0Z/ / /[` 0^0 0 0_` 0b0 0 0c` 0f/ / /g` 0j0 0 0k` 0n/ / /o` 0r/ / /s` 0v/ / /w` 0z/ / /{` 0~/ / /` 0 `1 `1D h  >  >  > ? ?8 ? ? ? @ @ (
 +*(
.
 &+%	H &C &C &C D Dr5   r  r  s   @r3   r  r    sT         GD GD GD GD GD GD GD GD GDr5   r  c                       e Zd ZdZd ZdS )Cpunishz
    The following are from the R script in models.datasets.cpunish
    Slightly different than published results, but should be correct
    Probably due to rounding in cleaning?
    c                 "   d| _         d| _        d| _        d| _        d| _        d| _        d| _        d| _        d	| _        d
| _	        d| _
        d| _        d| _        d| _        t          j        | _        t#          j        g d          | _        d S )N)ٜ.C1?c?繮ҼgM>"= ?h@7l[ٸ2gQ,4)g2?g#ۖTS?g%x?g]D ?g:Kڞj?g[!"@g_@gNTa@   gu2@
   r   gLu贲vS@gQ@gz!{#g`@gie?r   g8@)gqA@g.vVfd @gͼJ?gm"}@gG @gN3^
?gY?g[*MV@gh?g6("@gI} 7@gM޶
@g?&J7@gI[pN?gb<V?g>B(J?g^?)r#   r$   r"   r  r'   r,   r+   r%   r&   r*   r-   r)   r(   r  r   cpunish_residsr!   r   r    r0   r1   s    r3   r4   zCpunish.__init__`  s    86 -*%
*+%	%
'%4H &: &: &: ; ;r5   Nr6   r;   r5   r3   r  r  Z  s-         
; ; ; ; ;r5   r  c                   "     e Zd ZdZ fdZ xZS )Cpunish_offsetzS
    Same model as Cpunish but with offset of 100.  Many things do not change.
    c                 d    t                                                       d| _        d| _        d S )N)g9m4&r  r  r  r  r  r  )gJ+@g1?gjMS?g)Wx?g      ?gBfj?gA`"@)r  r4   r#   r$   r  s    r3   r4   zCpunish_offset.__init__  s/    -r5   r  r  s   @r3   r  r  {  sB         	- 	- 	- 	- 	- 	- 	- 	- 	-r5   r  c                       e Zd ZdZd ZdS )InvGaussz
    Usef

    Data was generated by Hardin and Hilbe using Stata.
    Note only the first 5000 observations are used because
    the models code currently uses np.eye.
    c                    t           j                            t           j                            t           j                            t
                              d          }t          |d          5 }t          j        |dt                    dd          }d d d            n# 1 swxY w Y   |d ddf         | _
        |d ddd f         | _        t          | j        d	
          | _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        t4          j        | _        t          j        g d          | _        d S )Nzinv_gaussian.csvzutf-8)encoding,)	delimiterr  r   i  r   Fr'  )gb0?gsF9=пg|H?)g\s)?g+?g->䱏?gWJR@i  g̢?@i  r   gt7Kkó@gp!#?g\xg?̔fKgY?g~c@(  gPeb?g^Hp?gӀ
L?g?gp'č?gt?gkN?g6?gN,*f?g%1,?gO]m?gW2t@?gUjrG?gf@j
?gG?g)Y&?gQTOIJ?g4?g{?gʉΤ?gCR?gF$e*?g*0Kت?gw!?g=܂C?ga={?gd㠨?gw#ro?g
+*?gw?gCn'?gzx	-oW?gqFʮ
?gy0?g7Z?g|}k/?gM'Fx?gpkz?g-?gJ?gFS$x?g+n7_?gQB?g 1?gKYE?g%6?gۨF	?gx.?g??g+/u?gD-?g?;fl?g)Vj?gYe?gil@?ge?gO/!x?g~?ggae?gά\?g}
|?gQzk??g8ڵ?g)/?g\Y?ghV:Hm?gZ!0?gAR?gn?g8:?g|HF?g7J?g g?gL=0?gtYb?g70?g_#?gO j h?gD]?g
?g%UN?gx^?gsڎ?gߟ]n?g_Yi?g-Nx?gQ芌?gj$t]?g?gԼU?gtK?gܝ ?gqPJ?gW?gأG8?g4Ix?g
F?gڤhF?g=(?gʽ"?g!8F?gOiDY?g=((E+?g_Vq?gRr_?g};??g.UT?g,
h?g(e`?g,͔sW?g r$?g5_`??gt?gu?g7w]?g0-a?goۥX?g.KC?glN7?ggyY?gy?gGl?gE~?-?gm&O?gZP11?gSy>^?gT&NQ?gsʊ?g>1?g~+ϊ?gI?gH쇖D?gǥt?ggK)?gbu?gz+r=?g1?gv%?gx%Wk?gom?gQ?gyh?gz̷?ggQ?gf 8?gC?g$}?g<\?gc?g_?gE!?gz0?gC1?g56?gcG?gv*?gH0?g#&s$^?glL?g+zo?giBR?gi?g;?gu{	?g F0?gm,BU?g]3T ?g&?gT#9b?g0(?g[K?gBhs<?g>I|?gxSx?g|/?g9CD?gL^?gޜB7S?gj?gd2u?gr}i?gЄ^j<?gJȅ5+?gN^H?g(sԕ?g3@;?g!/?gd1?gim{?g;جD?g1r?gRX?gZR?gb%I?goh/?g},TiQ?g52\?gq3F?g5<?g^mZc?gmSy?gů?g4i?gL\dl!?g	1?g+?g?g{m?g@U}?ggn ?g7?g-BJ?gSqd?g]!b9(?gx6o?g%[?gtm`J?gsP?g_Ɓ?gKHT?gư5?g:?g*?g/+I?gHy{J?g窳9?gk2
Fc?gb/Q?gOq?gK<?g)E^?gph?glÀ?g~]<rZ?g?gc >/?g7ZH?g"XH?g>
Sk?ga*?g$b@Q?gߖ??gA|?gP?gO|sf?g!7&?g'`?gv.?gA"?g?4mT?g`Jst?gLY?g7ê͡?gIW?gdv?g/};?go~&?gcU)?gk?gd?ger=?g6A&:?g
?gf u?g*`k?gHACA!?gUL?goY`?gŽ?g?*b?gX?g[%F?g;
B??g`q?g?g?gH9?gY?gbs?g5F~9?gKNE?g)RS?g?g{xk?go"?g1ñ?g"?gѹ0?g"0?g!x?g1 ?g4H^?gN?g\(Z?gsl?geJN?gf?[?g'O?g/p?g@0/?g1Up?g%o?g
:1Va?gD>?gq<?g`v.q?gN4?gr-?g7?gsY?gg35C?gZ?g@
C)?gh5?g`?gz+?g68`?g "	?gA?g9#J{?gJM?gO[?g`
?gK?gK!KI?g?g{G?g$cf?gkXa?g'?g?g-t?gʐ ^?go?g"r^@?girR?g<?gU=?gƅD?g'+?gjK?gv^?gË?g!~?gfEC?gk$\?g}#/?g ]?g^C>?gASD?ge?g?gO5|7?gu?gmS]<h?gݢy?g1?g5Ɣ?gu|q?gѮc8?gu?J?gQ?gO-O?g>C?geDqx?gm-?g݊$?g?Uf?g[8p1?gSu8?ga0M*?g85%?g
G?gZxP#?g_6w?gsdq?gxp?gM$?gaH?g9 Wnr?gq{{?gl?gf?gme[Mh?gF]0d?gje%Z?ga?g{`lo?gE{aכ=?gA?g%ylӑ?gqrxs?g
\?g,n?gVk?gu?g˘_?g=<SA?gPE?gC~f?gj}?g4H?gg?g8ѣ?g*e`?g,cSͩ?gɺ?g;v^0p?gPdr?gP*mA?g	
?g3J]Sz?g9?g3vxU-?gR ?gRu!?gjUdJ?g7?gh?gx<?g1?g4Su?glZ)r?g"b?g&?gh05׃?glG?g.]w'?g >p?g,L"?gH>?g4)=?g
?gipK?gt?gN8;?gBDbi-?g.q?g=6c?gܳP?g}ZfQ?g}v?gkt2?g,G?gb?gYvg?g	v|D?g5	#?gtk?gH08?gM?g]cY?gL#?g|ǁ?gR+~?g2?g\0#?gMr.n?g./?gr0q?g=^L+?gݠG?gABK:?g!Sb,?g~@{?gy?gvO?gz-?g@sD?gMb?gly?g=\#?gX6VSd?g?g̦+?g'AmO?g| ?gLE2?g5qh0d?gmP-?gYɠ?g?*('?g?g		*?g4/?g8πz#?g+?g{mO"?g-`2?gT?gV]{n?gm2pa?gّ?g)Hj?g;z?g}?gfQ?g(ݯe?gv?gQ2?gQh?gil?g	ڜs?g¬iAE?gvr>?gM?gG?gn?g1VX?g+3Q?g]?g޹SZ?gLu?g=#?gqPi?g[?gbҥW3?gNm?gf?gKog?gR`?gwp"?g54?g
|T,?gZ\?g{
?gC?gvQR?g[?g3?g[j?g=7?g<st?goK?gM?g!?goB(?g7S^?gM~:?g&
M?gtqej?g9F@?gv?gBGE?g׎Vs?gF~l?g4U?gBF.?gi&k?gR	\?gJTA?gu?gБ6?g-r?gAPɔjz?gL[+^?g~Z?g?gs,	?gzJ?gPI?g?g?f?gܯ
?gOS?g)w,S?g"?g8l?g5?gX֨$?geKR?gl<?gvw?gxY+_?gYS۱?gi)3?g9?g
(?gUzx?g+?gU>r?gxO"?gau?gYH?gAY?g0\c?gN?gċ^?g'S?gݏ>z?g,T'%?gV>?g< }Q/?ggy?gǺU?g=2?g
j?gy@#?gQM((?g?gf?g	^{?g<҅C_?g+h%?gzhH)?g]1'?g;f?gͲV?gϴX?gwc ?g[8s|?gըA?gu.bm?gS)i9?g2{?g'ф?g?g,e?g?Օ֪?guB!?gMט|?g?gL~@?gJ ̣?gd\?g~T~Ol?g}r$?gIU?gDg?g8?g]BO2!?gR?gu#fk?g{nc?gY~t?g^gu?gK+?g~Y??goD:?g@ݻ?gZ??g?gɫ6?gW?g2IÝi?gD@?g柍s?g2-n?g}O?g20#o?g b}?g0FL?g#u:?gLNJ?g6?gm?g?gh[(?gߌ?g!$_B?g4J4J?g+_?gs4?g2hH?geU!qѡ?ge?g㗏O?g{?g,8Xh?g+t5en?g}F?g3'7	?g:Y?gl?gbd8?g\<v?giDN?g?g|]?g;y`?g`0/?gNG?g>M?gw?gz:??gu%V?gkD?g33?g2?gm̩lq?g?g?!v?gDW}	p?g=X|?gS9h?gu?gY>fF?g*?gO)?gfm?g䷮G?gpV+S?gGHy?gjf?g6ˡ?g?gVX:?g$<e?gW?g!??g>Z?gE?g`YH?gmM?g)A<?g}?g@xEl?g}v?g )}?gTO?g}?gV?gl?gSUg?gx?gTqsz?gB;+DI?gL`B?g0?gE?gGV?gg?gXL?gF_R?g;?g%?g$\?gIw<P?gyt?gT?gK3?g?gK?gq.0h?g<"2n?g|o?gX)?g>%?g4Gw?gzJG?g")?gP6
?g^|veR?g]	?gt\X?g%#78?g&.d?g_q ^?g򦍰v?gcK?gc?gW,?g&#Ll?g	?gWU?g^pݟ?gW^?g5?ggi<K/?gKE?g6~T?g=?g%?gɿ|?g@?gb^x?gR<W)w?gwh?gm?g[A?g:m2O?g=Ⱥ*?gRr?gY:??gꪂI?g&Y?g\?g!?g68?gh?g|oF?gh?g/#?gIVUZw?gV?ghRr?g]m|?gNJ?g{$ʃy?g2s}@|?g(߿[?g~6P?g?g!HV?g$k%?g	&X?g}6Y@?g%<?gdi?g>ܼ?gISNz?gϦ?g 28?gד
t?gw?g:`?gH!=r?g$c"?g#i`?g8F?g0J?g5B`a?g&[e|?g-i#?g%}?g֩M7?g?gjF?gĺ+R?g;?g]4?giOMrU3?gå[];?g?gRʹ?g͚Po"?gfF`ۧ?g䉖?gaOQ?g[?g0g2?g2?gVA?g(??g2?gG7?gU?g#3%F?gBԖ6?g+?g@?gUxE?g;?gJ;?g'?gŜY?g6]I?gOuY+?gy(Γ#?g)fo?g:j3?g	<?g dY?g4Op2?g[?gS@?g`#a!?gO#x?g-³k?gb?g[MӈU?gZR?gwYX?g<aPC?gk#?g	jv?gP=?g> ?gV9?gg
J?gS?gVJVx?g-_?g?g6Y?gX>?g\zn?g ih!?g	?gu?g DU#f?g9z>x5?gˎs?gS3?gXq޹?gT#?g8o?gE?gqRi?gID?g'z*?g#K?g _U?gHg?gJLH?gL?gM?g{O?g.?gW1dh?gW?g)3?gJ}8 ?gp\J?g9~,?g؉?gM-?gYε;?gW˳?gs?g% |?gN	?glڕ\?gȄox|?gR^
s?gcbG?gSm$?gwo1?gJ[$?g&?gϢV?g@w?gP.ީ??g?Ne(9?gQ_?gCO?gH"i?g:Հ?g&?g|Iԅ?gAwQ^?goLhB?g~]d?g`3#0T?gMD"?g?g)Pr$?g;6Zim?g*?g1g?g`A@?g?gcF?,?g[v?g{,;?g!@)?gdqʭ?gKm6?g>RL?gt$c?g7!?gq- U?gL'	?g5?gR2?g~F?g)s?gCO߶?gzW5?gî^?gP{?g_ :|?gjc$?gSpat?gّ?gX?gàx?g>.˙?g޹A?g:Q?g$?ghB?go?gUԉ?g='+?g`?gYx"i?gR?g^*?g_t?g(P?gD=:?g7?gY8A?g0J?gN`?gkE?gtZX?g) qw?g%"xS?gd9v?gэU?g:e?g?gi?gc#?g$X?g?g|$?g?gKq?g{?gӼ?gLgY?gS$1?g3xNI?g`E?gF	v?gK?gb)?grH?g?g"Rpk?gӟ:?g@?g4ȣ?g)N?g}us?gd]?g(Lk*?gIlU?gSK?gʹ=?gฌ?gi:w9?glP?g/#w?gg(\?gqĴ?gӖq4?gvHx&?g0?gIB,?gkr7?gm`}?g>Z?gh%J?gZ3\?g1dX?gSp?g0@?g%y?g `
IE?g޺{?g6E?gePmp"?g<?gLW<?g양83s?g>Ak?g'?g?gSG}H?g;o?gkM!?g;?g@$"?g*is?g殩(?gPC+b6?g%)q?gU?gJ?g7ؽ?g&7?gif@?gz:\?g?g+>p[?grWD'?gǢ#?g8_?gL?gEj5?gS[Ɩ?gđZ?geW2?g7Pg?ga?geb?g<	"^?gW0W?g,ǥ?g
c?ge+?g6?gg`o?g`?gudҹ?g3Z?ʬ?gigK?ghH?ga?gJ@?ghBa?gA I?g,q?g4?gP1;ͺ?g͡5b?g<l?gF'?g:R=E?g]?g0?gu(} ֓?g3?gI8?g5v?gkO:?glG?gX9?gCY?gB ?g(q?g2z?gXwH?gDG?g|f>?g*?g?g ,z?gL۴?gtŵo(?gݛ?gÔ+?gI:?g(r?gFs?g??g  ?g_--٧?gt?gI?gv?gph$?guX?g ?g  ?g0.n?gؓeNv?gh?g4?g`?g@?g:*M?gnٞ?g)%j?gA[?g&6?gˍ?gw?g+f)?g1O?g]aӑ?gʒFK?gDS\?g}Fj?gvzg?g'?g8?g/?g kB?gv?g{+T/?gK_?gI7?gXMv?g s?g,{Z?g?gU9/2?g&\X}?gRp?gI?go$R?g`u?gg.?g,Jh?gR2=G?gƐsͪ?g|?g0/?g1ƌ?gGI?g;3?gb?g0Gh ?g&?g-F?gN?guA؄?gT;`q?g@?g,h5h?gf}?g%M?gČvR?g}"?g?]@?gO+d?g#.+?gu3n?gm<W?g^?g]}?gm72?g* iG	?g\f?gRr?g>[?g&d޴?g\?׮?g-F$B?gAu?gؘN?g^L;?gynV?g~<?g@*?gIsX?ghpcI?gdX?gԼa ?gJ?gξ?gl8 ?g<E7h?g>>A7?g:sQ%?g>c?gQ^?gm7?g_6l!?gD)?gE?g|9껬?g%T=Or?g0@z\?gH?g{Fl&9?gܘ?gw)u?gӰ?gJr/?gA?gQ?gbCy?gTm?g7^?g]I@?gX&?g}F-i?g,*Xg?g~?gZR	?g$M?gR䵷?g`}U~	?gl#?gm?gyK?g3?gBDbi}?gB~.?g!?g~yuL?gt?gFD?g(Y?g:e?g	g?g"?g*I?guV?gEx3?g:T?g㱗T
-?g$cu?gUI?gpH?g$ I"?gOB[?gjDF?gֹ7?g0%?gFt?g81p?gVPj̘?gzH?gWv?gi(?g/*?gC`?gm>?g9"?g牄B?gk?g؄ ?g ?g~qi?gPx?gLP8?gO[~?gt =?gu?gsI?gOr?gT@Ҏ?gM]?gBqKԀ?g?gjwY?g?ES$?g\?g`
?gXTn`z?gT qx?g_*k?gtO(6?g&J?gϘa&?g=}%
?gy?g," ?g ?gkf?g#*?g("?gvZФ?gz9?g3<Qm?gCq詥?gCWJ?g\?g`3?g@իȨ?g1$?gkim?gkٌ8?g<	?g(Ӊ?gt]?gb7Cb?gx
u?g3b?gN?gb?<?g<I?gZ6g?g_?gX?g	MW?gJ~Z"?g-?g/?g1x?gH?gcp?gbE.J?g1N?ǧDy?g	<*?g }?g%?g[?g2JDn?g8Ϟb?g`?g 7n8?gn7ԭ?ghb 4,?g(?g.a?gn?gTL?g?iyL?gUEZ?gQˋ?gl%m	?gf;?gxe<?g!E?gU?gE;ٓ?gdEh?gE^?gkC?gיkfg?gQ?gz,ޫ?gq?g[WPA?ga?g0wHs?g}1[E?g@Y*?g'B>?gG9?gAB?gwR6Ep?g?:?gf3]K?gCO#`?g|?g_L?g`O?g>
@s?gc՜u?g@ov?gi?gh]ݐj?guD?g'm?gR:_?gB{1?g3.S?g܁?g,^6L<?g ?gl.?g}t?g)l?gkY?g#V?gV?g	^?gaK?g?gE>T?gm<?g+
?g^#͈?gm̼?gGkZN?g>.?gF,?g;?{?gР*?gq?gM5?g*'oD?g:?g?gq`: ?gユ>f?gZ*-?g*r#?gR7z?g㑲?gWė&h?gG"?gF ?gphd?gM:=?gB=?ga-y?g>4^?g<SZg?gah&?gK?gRla?g%6?gݔL ?gUq?g	V?gC?g?/`T?g_?gc0?g5v&<?g {?g*j.<?g׳?gtIj?gqT:?gJ?g
yNlS?gP?g>Z?g?g@;?g<?gu?gl9]7!?g
L|?g̷?g2s?g.,?gpk4,?g]U`?grq?g:?gַ#O?g!la?g>d-?g&ET?gq&"??g1_?gV\?gUhy?g)g9?g7?g?g@!H?gT?g
7ʏr?gW>?gVlĺ?gZ(-?gmf?gTV(?g~$~^?g"{|?ga&?g{?gyj8U?g5>?gI?g2?gaԜ?gܼ8?gá?g·FH?giF?gcE{?gԨU.3?gEt=?g<p?gܶ.?g6(?g=[j16?gi&WU?gC;N?gƘ?gɋ܁?g1O?gvӂ7?gǁW?g޿;?gn}Mb?gNm9?gZ.UH?gp;?g49JH?gl{%9?g"%?g5t?ga?g>b?g-І?grJ?gxEڃ?g#?gD.?ga}_?gʲP ?g[z?g+^?g۱?gHp?g}-;?ghÍc}?ge?ga?gߍ?g#UF?g J?gt>?gUe!?gO~f?g?gm1?g)G}i:?gLَ?gcȘXl?gv?gQ>?g'Q\?gJ|?g1WiA?gO?g3P?g/	k B?g@M'?g9(?g7&<?g~e7?gOH?g$})?gEu&2?g|5G?gJyG?g߀{?g(u`?gv?gc#?gܕh?gyvu#?gfS?goXb?gny?gb$?gDn$?g%ggG?g{9\?g?gm?gd?gZ:-?gI>]z?gBE>L?g{@?g-?g0q@?gՑဏ?g߻0?gCt+{?gzER?g(m}I?g@edw.?guc_?gz?g̀?gʊ]?g
!?g;lj4?gv$$2?g W?g*&%~R?g,?gL(+#?gtS%?g::?g!?gr?gҰ?gVQ?gOuY+?g)R夅?gAа?g'KJ?g[Ia?g r?gȉ?g;R[U?g}"E?g?gS!;S?gůn<?g9A?g֗,R?gEY?gZW.?gdh?gpi,H?g⒡2;?guU?g4?gYݨ?gm܇S?gFѝ?g<9?g!eYjQ?g̡=Ef?g(T?gҽV[?gZ<f<?gg?gdmH'?g?g$wR?gj$t?gos?gL
?gfZj֜?g*Y,?gh4?gm6Vb?gߓ^v?g@?g3_?g.9IN?gZ44?g&kdL?gI?gp6ʗ?g~}o?gHFR?gܨ?g!bLu?g)fM?gVWyo?gv|?gnu4?g%?g)2f?gaHVv?gSM,?ghm~?gwb?g7?g&?gsd?g	+?g|?g(|:G?g &?gp~?gfQ?g.*?g5?g޽7ͯ?g?klT?g](?gx3o?gȯb?g4`^?g0`,&?gVM|?g}\z?gĻo1?gQ`?g	?g~3Rk?g5D?g;:?gӿ@?gF?gr^Zt?g2^+?g&Ğn#?gM/?gF?&j?gQ!?g;c%(?g?gʼ/?gڱ	?gO*?gk>~]?gS?g8?gU?gev?gcea?giߘ?gKC?gT*4R?gv?g/Nv?g?g*-
X>?gP=?gq?g-^L?gyvT0?g#;?g	?gz<!?gjB?g+D,?gYua?g?gC?gQ?g%HN?g\^
?g̰A,?gf\?g P;?gqBX?gÏ Z?g`_ ?g	ͥ?gɱm?gOjy/F?gzc?g3_Z?gN%?g~lPA?g?(a?gd/%?gx4?g'8i?gVe!?g\PE?g*?gG66J?g4;6?gPH?g)\?g+?g}?g!͒ߦ?gK @?g2b&?gJC?g~eR-?gU*,?gH?gUᲬ?g%w?gU]?gS?go^?g yf?gwjT?gjR?gȮD?g4	?g@?g?gka?gs|Y?g%R?g0?g.?g`G5O"`?g@uH{?go\ux?gӬ?gA?gCe?g?gra?g=?gi?g!iJ?g3K?gf?g݃?g7oV+?gE2?gzX?g5E?ge!?gPm?g<7+h{?g"%`?gO?g2Y?g!Թ`?ghY?gE?gv-sj?g1l$?g=8gb?ghp?g(m}G?glqS$۰?gLoC?g k8ı?g0?gC5?ga%?gS?gOH?gi?gLH^?gx{?gXӫH?g㏚]t?gȾV?gk\?g nz
?gJ?g֍q?g+?g˴d?gBD!9?g:?g?g^?gf^ sK?gVf?gPQb?gVH5?g_H$ѩ?g*$r?g1d+?g2)?g6 ?gai3?gl,?g|YY?gd@?gޭɦ}?gAb?gRwe?g5x¿}?g@5?go?gPG?gǪ?ghNE!?g"I?gcKUt?g\k?gtd+?g*?gBS ?g뙐?gq?g
ҜB?g]-J6?g?gV*?g|E?g_\?gG?gL0އj?gZ+59?gId?g@#?gߣ?gQ뢍?gMŊ;?gan?gҁCZ?gR s?gpqX?gM	z?g/
?gL?gU !̌?g"=?gE?gs
}?ge?gP?g.R(?g,,v?g.Q?gpX%?gUX>?g>*z~?gWqp
?go?gʳ?2?g )?g(2?gu)?g?g=`?g<SQ?gE"Z?g1Ti?gam?gYF?g90b?gw?g Yf[?g﨔K?g!?g1??g[71fCU?gð	?g"r]?gy=?gK}?gWZ?g?ms?gFa?g<?g"u1?g\?X!@!?gN?g[ֳF?gCDD?gr?gK?g$Z4?g@8
w?g?gr5?g=q2?gh#? ?g,^n(#?g_Q:&?g-%;??gdl?gaX?gv?g8#`?g U+?g`?gN?g,B*[?g?gR?g'?g*pp?g^?gjR!?gƂ%?g_pr?g^6t1?gK]'?g"e ?gc?g_?g<S?g?gb[?gLX?gUa?gA ?g#T?gx>1?g2|])?gIW?gnR\?g\)W?g݂|;O?gsxM?glK߭?gUv+?gn?lM?g?MV?g̅@e?g0O,b?g	ސF?g4|^?g3]+3?g"?gg?g1Ĺz?gOnS?g;d1?gqw`M?gxwl?gd73?gKesI?gnG?gښc/?ghR+?g/E?gd?g%T?g?geb?gg!?gq)z?gH@"?gPPH:?g/LQ?gjH_?gr,?get?gu?g?gO%%?g"x2?g
?gc>?gO?gũ$?gxn?g:yO%?gഁ?gB-?g-)<?g|8?g|h?g,.ک?gy?gMt?gf?g}~R?gE-(gC?g<h?gE?gAs?gMS?gcWb?g)=?g'T?gLX?gy?gNN	g;?g.0.?g-:J.?gc?g>)co?gY};?gGT,	?g|#?gf/W?g#?gk) `?g.};D?gFD?gیA?gfT}?gCMW?g')?g$?gU"`?g"?gL?gxbi?gb?gX?gk?gKi?gp+?g2Ќ?gx.A?g#޶ҝ?g8~?g)?gbDs?g۪^] ?gęxjB?g0h?g???g&i?g\gK?gM(?gsQ?g薩	?gzEC?g?g>+;?g!4:?g+´?g)?g] q?gj"}?gTMË?gtө?g*0?g7 ?glȖ?gqOr?g?g?g*1-S?g}s?g/M?g@0s?gC?g?g~vfX?gХn=?g!q˄?g$iP?g㪪~?g)Lo^?gL.%?gة[r?gs6?gnL?g$C?gW0?gkJO+?gI|0?gh_W?g}?g^R?gSBM#?gh)o?g,Uf?g-S(??g2hH?gcVu?gdC?gձl?g/ã?g*g?g΂y]?g_$?g(q?g 9{;?gbOco?g[wc?g_~?g˚ܨ?g.r.I)?g%%C?g
 J[?gQ?gj(?gYc?g~M9s?g6g?gv6?g8
9d?gHSk?g50$?gOpVC?g:hjv?gGD2v?gs(?g"j?g?gq9?g9
P?g'oRbr?gU4X7?g#7j?g4U?gר=?gʸ5+V?gwZD?gg?g: Y?gw[d?gӃ'AA?gU>6,?gEzz:?gv?gwM?gb8K?g2b?g:?g74q?g&?g,h5?gh:d?gpr`c?g	D?go3?gc[#?g5V?g?Sm??gJpe?gN
-P?gvE[[?g:mS]<S?gg"Sٹ?gTz&?go?g^?g\3(?gc&Kj?g+E[?gaU?gi,?gIF?gŅY=?gŜy?gɅ0߮?gX	?gyxS?gߙ;$?gNĎ?g	N?gFa?g^gE?g]`Ѿ?g#2?gˈ#N?gAX?guH?gUQ`D?g`el?g^P|?g>ʒ?gFKj?gW!>?g:LI?g mu9?gJ?gy:?g?g`ng?g=?g:Y?gㅶZ?gID?g7A?g!P?gU?g+w?gwY?gmVcd?gkR,g?gҼF?g[J>(?g-~3O?gX?gAWl?g8`4?gV@?gZdx?g@N?gA>N^?g;,#`?g?gj?gZW?gyj?cp?gr֎?g9?gXH+&?gxZ?gGrz?gyɜ?g2^or*?gZt?gcO-l?g u}?g~.d8?gr?gC._?g2ʎ?g@*?g@?gm?g7WI;?g,-¼?g\<"?g2?gu?gA?g̩	G?gD54?g"'(c?gFn?g?gG?^P?g:8?ge?gjlmQ?g\xz}?g=[xj?g3V?g%?gǟ?g6.?g}"O?gP<?gCC?gXb?gB?g+?gjvr?g\km?gW&?g]DW?gX?gQi?g?잟V?gb+˄?g\?y/)?gc&?g(J?g>2?g5ui?g	\0?gWe?g"?g%?g'	?g\d?gx:?g![к?gֺ.?g]HMܦ?gNZ?gOm?gKh`?gduFkw?g?gͳ?g@<L?g?gG	]?gLy߾?g۴B&S?gTZ?gGlV]?gxH?gcTA?g?gOS?g5?gɤLv?g`Է?g7|rk?gX?gV;?gUΈ?g k?gs7?g%}?gȽ?g:?g1HUn?g\{?g"Ƚ?goO?g2R9X?gbeU cH?gjr?g]S?gVժKv?g,?gHJ@m?gsS?g3 t?g:U?gV?g ?gLI?g{V?glZ?g4;^?g+3?gP
?g^?g{*@~?gXe?g>|p?gs?g/L#?goiow?g"h?gO·vò?gJ?g??g~{	?gz_?g6gsp?g;D?gJJN?gƙ?g'P?gm?gm)?gm-?g3.=6?gܕC?g#[l?gIFRBB?g<p?go@m?gڜ?gJx̗?gC?g~'X?g`@iC?gw޾?gOs'?gH?P?g^?g>9
?gzbf?g?gGER?g#5i*?g<4$4?gMM?g,y?gu%G?gtJb?gB]gf?g?g#A$?g)G>?gc[˱?g`?gVz?g
SN?g.Ӕ?g`,h?gwFh?gm5?g[ P?g]vr?g9?gHt?gJ?g6^?g@ #?gvTN?gu?g}Gd?g$y?gn;8d?glEq?gY??g샏5?gU8?gjПOGt?gP?g;/?g#s?g[5e]?gL.?g<;K?g
h]a	?gIvX?g
H.Ip?ghMn=?g0P?gw7?g-p0'*?gG?gm-8(?gg'AR?gܙҡ?g f?gJ{Y?g~Ӊ>V?g)I?gÌ ٱR?g$/)݃?gw?gT /{d?gToe?grG}m?gxQ6?g/*0?gF?gf?g&V??g78K(?g拤z?gzcl?g>?g\?g}6`F?gt?g<`G?gޘ?gm&?gY?gc ,K?g;gj|?gX}O?gDz@q?g?gtf?gxW?gV*?g3P|)?g >e3?gʇj?gH?gіƎ?g3F(|?gTC?g(?g7bH?gx]?gtWLy?g璁?g<-S?g@?gTI?g_#u?gs?gP?gyo?g!I?goQ?gQj?g;hZ?gTK?g.?gjwby?grK?gݽbh?gnK.?gGZ?gaxX+F?giC/x?g?g	?g)I?gk eݢ ?g18?g= ?g,eņQ?g[?gxD"?gҌ?grD)?g]I$?g8N
g?g[?gIv?gÖ`ig?g)/:?gsQ?gv-c;?g^]T?gܭUP ?gƿK?ghc?gT!?g=*~|W?grh?g%"WE?g.f_%?gN?gCm?gy3P?gbV?gCJyΒ?g Q?gb?g֯?gqL#9?gnQ+?gLLv?gLA^u?g@?g3t?gZ즭D?gI?gkbj?gUOAٟ?gB0?g"<g{?g?gh)E?g{k2E?g` B;?gIW)|2?gB½S?g43
?g^@"?g,?gȴW?gUj@+0?gA26?g5r?gmX?g9\Y+?gXV?goZ"Z?gBZ?gꡍq?gB{8?gث<??gU\?g
#bo?gڋe?gZGv?g.b+v4?g7)p?gR?g'?gU
&?gܤ\f?gR
q?gc/?ga<?gSC?gUw?g<?gt0rx?gL!AȤ?gwX?gH`Z?ghLl?gjW?g)??g+тB?gw#^d#\?gc<1?ghby?g);!?g@+Qr?g"f?g3~?gQ#W?gч_E?gAa?gn?9	?gP?g@r 7
?g	N3>?g]O?g?g92?g5Z?gWc"?gK~Z1?glC\?g!Et?gUS?gA15?g%&?gb*?g<h4?g%%|?g?g8*ܗ?g-Yd?ghTz,?gO9A?g?g4)?gTr?g1mS?g7
Y?g|
?g9 ]&?gKfie`?g|؇?gV?g1?gSFr;?g?g;?gQ?gi?gy?#?g5F<?g<?gUK?gQ"3T?g-Zҧ?gxk
?gP?g{?gɣD?gp?g8E˪;?g,'?gΨu?gm.?gB쯉?g^zG?gEN?g$?gM/RҚ?gaO?gN6m.?gU9?g"k?gul&?gF y?gAQ?gnSs?g? g1?gN?gs?gY*V*?g0+?gIxEϱ?g4-[?gމ?gk]?goa?gJn*q?gd?gV94c?g*?gxEf?gбC?gkShF?g54?gy2?g5ug?gԑ\?gI2?gn6z?gi *?g
?gPÖY?gfO!2?g$O?g)uM?gk ?gSޟގ?gه6?g6"?gީF?g|N'6?gC%*+?g$0?gT#r?ggp)!7g?g}?rƛ?g3IM?g7?g}?gQ?g>rR5͝?gQ`	a?g4?gyIA?g7+g?g0?g4 ?gH-?g¡xx?gT(rh?gr&S+?g0R ?gX)o?gA?g
}d?g\^j?g|p?gL<AI?g3~?ge_>?gsʫ?gN?g^<s?g2mb5?gA?ggC?g9?gt?g,~ذ?gBq?gD?g?=G?g".?g*?gux?g
1:?gB??g"t?g"b?gm?gW98?g&x?gM;(?g2?g4?gd?g\e(?g@E8?gch?gY	?g`X`j?g@,?g3?g^y?g3ٲ\?ge.?g?hp?g&M5?g(I0?gwcO?g^1e?gf2Dq?gz~B?g"k.?gY?g/kbb?gk>)G?gne׶?g".?gݕ<"?g&,w?ghҿ$?gdH?glc?g\@?gbc?gNS?gd1wN?gm[J?gs8?gj/n^?g\@
?gCV3?gP|?gf"?gyB
?gbr?g4?g81
Aq?gH?gf?gZØV?g@T?g;+]?gr?g*lf?g^ח?gaX?gC ?g+8?ghfk:?gj/	k`?gid?gۀ2:?gP?g; ?gɋ=?gqx9ӡ?gtzA?g TF?gyc?gM?g]b?gQ,z?g]㦄?g"r5?g|'$?gHPW?gb'p?g,bEZ?g?g_D1u?gfPe?gUdjͮ?gC&?g>?gU?g p?gY?g­T?g6߱m?g-T?gӺӾ5'?g6DL?gd[]?g;'?gS<.E$?gk>8?gi,2E?gYv;?g!a?g9ij?g)v4?g6 ?gN4?gԝ?gv\?gx4?g(Y5H?g;KVfƒ?gCޓ?gPS?g}
?g"HJ?g\Nt#M?g#[w?g9#1?g˧S?g xߥ?g~A?gOR}?g-w?ga"ןB?gvb?gg/?gDsZ?g@K?g[?g@!L?g^sX+?grM?gr	?go1F?gf12?gjQُx?g?^u?gG7Wd?g|?g?c?gMd?g;%?gvfR?gF*v*?g,8Sb?g,?gU?gz"}?gloN?gX?g.PZI?gn'C??g{þ?g8! z?g	?g.t;?g\WV?gs^ku?gFIs?gd@?gRq?glw0f<?g{@?gDRJ?g,?gkwt?g,K}?gYAu
?g|I?g80Qd?gN-Ds?gܜm?g]E\?giL8$?g*3gl?g6=b?g,(K?gn4p?g^ꌴ?g6"?gHQ?g1?g(?gz>n?ga\?gCn?g	W?gJbA8?gzґu?gB|?gư?gMU?giN?gAҗ?g&
Oq?g<?gα?gʞ(?gtx6?g⑙#?g>/?g5?D?g?g5?T?g,Wy?g8ɟ?g얪?gor	I?g$"e?gtڸ.?gnskN?gi?gfPe_?gT\i?gEsl?gp6^?g?g3?g Іp?g
}Q?gר=?g9\x?gd?gD,?g?g͐?g;16?g?g]?gR?gg|?g?gN|8?gu1?g߶?giM?g?gh>LO?g<?g{Jm?gAݴ?g3z?gɯ?gEX+m?gyϻG?gr?g!q?gM?g%#?g;@؇?gz5cz?gf-y?g@Q??gVQ?g`q=?g[?g<?g[Ԃr6?ga?g'1s8?g=o?g1?g2+E?gL}>Y?gN,lX?g	d<k3e?g%{7b?gη˟[?gGP{?gr~?g-T_?gޕ4?gC_rJ?g)U?gzͣɢ?gW-?g1h?gxǭr?g&~?g:l?gZ5C?g w[?g`7?gw?g>Q0L?gȤ(?g	HR?gb?g̎\PB?gI#?gD6稾?g~?g?g
V?gd?g`d?g˰3?gu?g,z?g_\ ;?g.I3Pp?g !?g,?gH[?gr!?g3~}?ghf=4?g,J|?g)HR?g5R?g?gaԜ?gp?g֢˼?gZ?gŃw
?g^`^u?g#"?gB?gt3?gHC:g?geL7?gWVֲ?gJұ?g]?gkMNH?g4;y&?g:?g5;N?gA?gZB=?gq/q?g{Q\?g~v?gN?g}Ct?gsC?gh\?gRcGg?gNC=?gnaf?g6#?g\n?gkd9?g3DVm?gv?gyumh?gYKZF?gNdg?g+n?gU	?gc;J'?gGru?g#?gKtT=R?g?gױ:?g^+;?g'C?g5M
Y?gG?g^X?gB?gg?g!`?gyrY?g蕙?gg..]?gr?gJ"?gA'?gbc?g?gt@.?gSK7?g&?gǔrk?gg36?gj	JT?goJf?gש?g	|z??gT5"?g?g+P?gF?gs
?g__D:?g3 Z/?gxhM?gZ?g {y;?gH?g<y&n?g=?gx4?gsaI^B?g ?g*?gN?g&9?g9*4y?g[E'?gtK?gɟk?gd]u?gehf:?g5r^\?g]
Y'x?g&?gL0y?gǲ?g(?gs"x?g~ߞkH?gm5~?gzR?g<kkQ?gf`}?g6U[?grMh?g䏚]?g_-˞?gU<?gZ?gfd?gmj\O?g+?gN?g~M?g0ie9?gu?gQ<?g 0GM?g?g|V%ә?g}d?gT`u,?goѰJR?gwu;?g'A?gߴǨ~?g<w?gox`<?gfNf?g|A?g5I?gf?gXWj1?gԢӍ?g|S&o_?g(-0T?g<?g?gvH9 ?g[ë?gЪ\*?g)R?g҄q`?g?g6?g-f[86?gF?gz5@i?gf!?g,\)?gԱU?gh~?g?P?gn?gW?g<.h)?gL| ?<?g^cv?gMM?gT7?gFQ?gĄ-;?gҾ-?g?g
?g\?g#G?g(Y삣?g(Z$?g[x?gX
wr?gS@3F?g#QG?g 2Yc"?g3g?g?gGb8?gǲa\?gn.?gz?g.?gn\?g뀹P؉?g8L?gD?gpǂ?g#3Z+?gY~Ǎ?gI<e;?g>i?g/u,?g*}?g8 c?gv?gI?g}bSD?gz?g=ʫK?gg?g!?g}-?gӜ?gG R?gR?g;cM?gh_?gy#M?g܌<m?g`m?gĻ<?gZz	R?gJ{Y?g<4)%?gSQ-:?g ?g7M?gQU?g021l?ga	Z7?gWY]?gהu?gSs?gc?gTD]֊?gMiˡ?gbSBd?gz!U?g	by?g)e)e?gb?gx1g?ghUN?gS 2?gq?g1
?geāg-?g-2I?g&H)3?g	?gb2z?gch'?gi?gwz#?gP-)-?gǲGE?g	?gpn_?g*/?g?ge:?g](ç?g뉮?gL.Σ?g9s?gͦ?gGx)?ghr?gR?gZ?g\?gF0'˽>?gA:}?gR`Ch}?g+?gwsl?ggD?gX0_b?gv?g|?geH?g?g!;WO?g2G?gn?gKzr?gPWo?g =?gl57?g	?g	|y?gic?gUp?g[= ?g??gU+?g+8?gy}?g{a:xa?g
&?g
T?gjptZ?g1WN?gt?gPc?g$z?gZ?gnH';?g%eA~t?gS?g	r?g~Sm?gJ?g4$1&?g]pF?g<Z](?gM?gT<q?gܶRl?gǒ#?gU?g#?g[?g/ a]?gg@?gHn)?g)?g '5?gV?g<?gZJ7[?g 1?gY?Kz?gʍ\?g|y?g't?g=X?g+)b?g^&?g=wHRN?gQ Ia??gռ?gtv?gxR'?gx!T-?gjB?g_`t?gY8 t?gq_?g˽P$?gu?gZ(?g(iB*?gTɾ3?g:Ԓ?g ?g=Lu?gLg?gl0D-?gsfN5?gq ?g[fN?gb^j?gFv?g(B?geM,?gOa|YP?g<,?g{l?gXF=?g=P~g?gGe?g_6C?gO)D?gί?gGQ-C2?gW_?gXF=?gS?F?gZoW?gM
sݫ?gP6B?gX?g?gyH\?gDY?gN=?gHB2?gk?gQȡU?g9a?gʙ;?g"6wM?gדȭz?g-?gT.4?g(?g+d?gtsj?gF)?gDH#5?g|6?g[[?g8hL?gɖl?g\?g6?gɵ?gm;-?gD̺׍?gg?g<ZŰ=?g b'?gQ?g7s?g2?g}$?g=?gq;n?gu? ??g2cn?g?g
?g?u(:?gl?gt~?gguI?gWwn?g|?gӥ<?gF:wB?gKudY?g0r?g*	"z?gT'?gs%?gmGB?g\p?gSeR"?gj)?gصb$ ?gzbI?ggB?g
B?g~@	??gIj?gQ?gzO?g.
/?g0I?gy?g?gM.!z?g-?g,x?g<?g_+?gHJ4,?gG?g<y?gD?g+x~fL?g)lL?gjus?g2
?gG?g,3o?gG?g>:?g#,?gYo?gz6?gHx3?g Jj?gHʫ?gO?gn*.?gC0ª?g*,hs?g_w?gxA?gJLnA?gШ,?gD?g+kPf?g=7?ghԶ@{?ghT	L?gW+,?g=?g,"HM:?gxSP?gj%2ȩ?gAB?gJbA?gGӰJ?gDtc{%?gjt{?gz?gddZ?g?ѧI?gr?gQx\?gi}$?glMLr?gJ@J?gZI?g>M)1?g' ?g/k?g Bذ?go)kaC?gQ?gT?g?gXuV?g>Η8?g$\?gI;=X?g7.Oa?giv?gge?g
[`Mm?g,[=?g=rB?g?goMO?gd-?geAp?gZ1ȉ?gˮŸ,?g6Fm?gEc?g.-U/?g+F?g1u4?gP75*?ghѧ?gT|ϊ ?gmT?g{OĦ?g\mo"?gA2?gN7?gpŤ?gLpd?g;t?g{?gP!}?g}xhq?gt:?gGORG?gvUN?g^?giAVJ?g:?g?gQiwX?g@qJH]?g]Q?gњ77?gD=?g`|ʸR?g
<;?g:y/?gH?g)bi?gS4??g4c4?gC0?gXm
+?g_mX2?g)M?g:?gGUg?g)p?g/?g {?gMK?g
i?g?g|"?g~t?gx\TX?gg=M?g?gdsv9?gJ?gו";gr?g]i
?g^u?gN!0~?g,K?gb?gD]?g1?g?ge;?g3c,?g.e?g@-?g[M?gZ׮\?g@o?gڟ?g]\?g?gXhs?gq?gP3?grGIYE?gm?gQl"?go5>?gSA^r?gفh?g}PcM?g@/ܹ0?g
?g-?gs=p?gb%`?gڥ?g?g\ً?gSf?gwb
?ggQ-?gt,+?g(7?gх?g֦p?g{A ?gF?g6?g|.	?g?gjM5%Z?gmN?g=U$?g|[4?g$O'?gӧѻ?gҎS|?g}?gis?g==?g^D.?gSU\$?g('w?g̀#o?gJqfr?gT[$l?g?<FZ?g90?gWԻ?g~8ī?g鲘?g,7B?g
8|c?gV ?g<nM?gS?g<QX?gv?g"2J?go?g؁=?gtm`?gH?gj'l?gP5#?gLy0?g }r?g^c0?gv4G.?g3 sL?g]n0a?g==?gSwD	#?gN?grD9?gpb?gAgt,?g8"-?g񈆜)?glS#	?g{?gۖu<V?g?g.P9`@?g(٠e?g2?g.9j?g\!?ggy]?g?gZ?gI?g[\_?g_l?g^ŘjX?g1?g3E%?gm^?gi*?gvJ\H?gԾ ?gE?g,vI?ginK?g?gʈ,.?gA}|?g?g
h?g2m?g;R?gq&\~?g7?gu?gs?g3WJf?gbv?g`?g}?g
T?gCc;B:?g?g-b?g2Sؾ?g?gQ9y?gʙyiS?g+7?g[Պj?guFM?g5?g3?gC@	]?g'm?ggzG?gl7[?g&$|M?gR'ތ?g{?g!`
?g2`a?g$c!?g5?gu1?gc?gD	?g8z>?gI?gO?g"?gŰ?`?g.[?g{6?gC?gBJP?g|?gԤ32?g'?gc?g|Cxj$?g]?g_vd?gJ0?gGF?g`G$g?g5E?gt/?g?g!x?g-0T?g5V?g~PK?gP87}?g|L{)?g/Eyj?gjx?gL"?g;~?gNa,?grZ*N9O?gwW.Y?gKԠ??g?gg+H?gaZ=l?g?gB ?gW?gaP?g;?glF9
?g]*,c?g\2?g ?gt?g?g%ۧ?g1!?giQ?g*Lc$?gmW?g>H?gп4wt?g,X?gF?gfZP?g3?g?g?g1/I?g'5b?gk?gښ!<?gWc֩?g?gKV=?gf
s?g[??u?gdsDY?gY?g8J=z]?gx"*?g)4<?g:|??gXR>?g©?grC .?gu5?ghĒ?g]%?g.z?gN|?g([
?gh?g"O=?g Wq?g/JP?g̓$?g,宩?gVV?g`^g?gEfr?g!)?g6?gZ[?gߎh?gș?g~[?gx:"[M?gvm?gv?g&?g!S?g21:?g{?gx?g0效c?g4?gW]?g$L?gZ^?g,iԴI?g1j+?gj ?g8<dG?gz{%?gRg?g,֔?gF&v?g698?g4?g~(Z?gܼ?gzB?g2A?ge[?gSk?g$r?g?gmYg\?g2/p?gfT,?g	#?g'L?g}.?gcƳz?gHQo?gfe?gꪂI?gJ&?g+?gܦ6?gR?g|B~?gtd?g,G2?g:п47?gƩZ?g e?gB>?gO"6?gw$?giVte?gE?gbŜ?g/ϔ?gn2?gEe?gg?gF'?g^:?gQ8?gJ?gN?g@?g~?g,o?gP?g+f?ge9^?g+qW?goX?g,?gyܪ?g'?gLB?g@D!Uv?gA?gȐ#-?g^ۙ;?gN-?gIK?gIR?g?gs#0?gcu/ 	?g~rC[?gI?gdE?g?g+#V?gagLy?g	q?gW)T?gs?gyQM/?g\
z3?gvAb?g1Ӵڮ?gЏLH?ghH)?gKos2?g-	PS?g7yj?gVe?g-fW?geuj%?gs:D2?gd/gG?gPTA?gx!?gw?gF7go?g?]?gYQ9?gYG[?gE)?gVuQ?gMp?g$@?gǸ$j?g?G?gne?gD!{X?g+ņ?gdF.?gM?gd?guf`?gSCG?g= ?g1MC&?grH]?g$?g#a;?gH_B?g-j +2?gid?gܕh?g:/ృ?g]?gS?gBH+?g1W	?gbib?g{b8?giF?gK;:K?g?g耣?gÔ/?g'?g]>?gSI+?g5&|?g֟P;?gJ5?ghpL?g"0?g]8)?g-Ykc?g;ٓ?g	v?gȟ	>?gl亙?gi?g2C&a?g̨C6?g0s<,?gܨA?g+?gOX?g?gBg?g{BQ?gK?grg~?g¸|?g|_?g"Isx?g(ç?gA1?g*CM\?g/k?gߙo?ge152]B?g0]o85?g`72
K?gKD?gY*l?gϋ??g?gmo?g'?gďֽ?gR?g7Gnm?gi?gB)?gq ?g/9*?gbEp?gE!?g-?gт?gX[1?g!~M?g57F?g]7K?g-x?ghJ4?g(	?g"?g4۽?gfA(?g
V>?g[]?gEpV?g1`E?g1	?gZ?g%3?gG^R^?g3z:?ge$?gT0?gjj9?gw/	?gf^dܻ?gB@ߒ0?g8fx?ghd?gu_?ghi]YX?gzoc?g0j?g7@7K?g/{k?g~įXe?gj?gGxZo?gW?gh4K?g$yK?gE?g;Y?gFH&?gQ?gQX8?go+6?gk?gFƎ?gr	?g^Đv?g}wLY]9?g\wa?g
\?gov4g?g!6b?g21?gНԺ5?g[gts?gчW'?gwwP!?g?gr|K?g:yO%?gQ?gǆ?g	?ggHOp5?g袾?g&FY?gݲ?gm]?gv_?gf~?gU?g W?g+D?gs{?g̽IO?g5c?g* ?gk?g|4?gt=W?g(?g?g5G?gz4?gyM?gӴ>?gJa?g|?g?gg?g!far?g?gѬI?gPW%+?g;?gL?gcm!b?glv?gє|pg$?gk?gklDj?gI|K?g	Dc?gK*"u?g1?g(?g^	XW?gUJ?gUǽ?g3Oϗ?gҼF??gx?g{Ѽ?g{# ?gָ?g ]?gvhXj?gχD?g]B-s?gF{zt<)?g4s?g1V)?gN?gAt?g{z?gf(n?gźeE?gg@?gHOQ]?gsT?gI?gz݅Z@?gG%86?g$C?g䂯?gLuد?gikG?g&̋?g *?g::?gMن?g]Ε?gm"?g`r?g~?gԒ?g4DR`?gS,?g4IG?gZ1T?gw
G?gJf/z?g?g}?gƸ7^?gu ʛ4?g8B?g/B%lG?g7%I?gy?g$?g2c<?g9HD?g|A?g4 &?g,y=5?gI9k?glEr?gϼg^?gq&?g|Mjüb?gL?gr%+T?gi0Oy?gqi?g4?gքb?g/naE?g(?g	?gs& ?g,|wL?grsq?g$.?gz\?gq?gĕ?gdhL?g
?gFE?ge*̔?g[_?gĎ?g$ƾS?gg7Yl?gsۄ{%?go"
l%?gHLg?gaD.z?g5W?gOb?gBN?gؼh?g?g+?g65-?getK?g}e?gL
?g{{|?gB	,?g&y$3?g*В?gAx0?gW(?grD?ge?gr?g풞?g@'G?g,Z͝?gj0?g?:?gUb"(?g?g?ׄմ?g̵G?g`f?g^zGb?g,h1S?gP'@?gjD?g?g#w;T?gi.d?gqV"I?gWZV?gj{?gF?gxx?gړu?gt+?g*o?g7ܘN?gH^Bh?gre?g^1u?gƳ6P?g-7!c?g OZ?g^:?gI@?gκ?gK(3?g,[o?gN|}?g@l?gMr?g$?gɪ'?gN
[X?gȧ5>?g!B?g'y)?gm̈^m?g΀?g%zr+?g@?g+|}n?g[{?gCY?gyX?g0?g\}?gV&?g?g?g|F?g$-2h#?gF5?giM?g;_"?gպ? ?guw[?g#,(?g
U?gL~?gv?gsj,?g-e?g~aC?g_?g&7?gR9^?gz0S?gl?gڸ?g\?g[?g/?gt,Hw?g&Zd?g*%D?gS?g\
l?g??X?g*"?gP&^3s?gAhCu?g>D)?g02?g?/	?g5tU?gIc?g3?gH?g@X]?g+P?g~Zk?g$e,?g9#p?gT'V\?gr?gW卛?gTu4
?g ?o?g?gMTc??gks
?gO}N,?gp^?gݙd?f?gI0?gvi~;&?g?gi?g`Bl?g4G?g1"M?g
?gfh~'?gs9?gU?gHk?g?g'ķ?g1xw?g#g?g{Q?g7?gH?gJ5?g	\xV?g-Բ?g#?gJ?ge[M?g

N?gv9?g6?g)5e?gM-?gA_a^ӥ?go~?gu?g^q?g
q?gMԶ?gg?gZ4?g}?g}-Ŝ?ggC}O?g,NL?g+ެ!?gcdJ?geW\??gh?g[27B?g^?g?g
y?gxlc?gk?gzNt?gƿY?gVZ?g8fv?g*[?g<2V?gu?g6??g`hC?gSz?gM?gY?gIS[?gU?g'GB?g~}?g*ɣ~?g3%?g[F?gZZ?gam?g`GC[x?g'pB?gc?gI(?glC?gI??g/yL?g*6Pb?g~?gZ2C
?gu?U?gMku?gR?g@'?ggP?g|?r&?g!?g»}?g܃)?g-T?gn,j?g|r?gf@?g"`?g͎i?g
~?gx(e?g׽^.?g;*?g503?gGt?g-V2?gM׮;?g)D~?g{u;?gh <?gC"+}?ggR?g8?gH?gTK.?gpL?gp޵?g.}?g:T[o?g?g YYD?gd}'?gd?gƦ?gI?g#\Tͤt?g	H)L?gfOጊ-?g!4?gyOD?gCZ}S?g&v?g&D+W?gB~O~?gG*F?gSl?gFj V?gظz?g ?gIc;?gE?gh?g<`J4?g!!?g<@H?gx)TO?g.7?g17k2ۜ?gEAt1
?g)w?gy]?g<pD?gV?gL\dl?gg'?g:v^0?g5fg(?gV}e?g>l{9?ga6}-?gJ7;?g84#?g4?gD#?g.;@?gˏzg?g#!`ު?g1}_?gv"?gHpϏ?goF?g9?gWx4?gUS6?g~|?gR]?gtvb?gBЀ?gzuyw?gs~c?gHè?gFˁ*?gd?gJ#?g%i*?ggۉ?g2Ib?grS=x?g}qJ?g3ެ?g¢#?gP?gߔ(?gK!?gK$?g#{?gF{zt<i?g&?c*?g[v-e?gc@^?g!:?gi񵩂?g}?gmH?g?gp~?gb?g'?gW?g(ă?g!¡?gK5?gZA?gt?gȁӒM?g5ZP?guf;2?gW?gUqL?guUE]?g}|Rt?g2"jZ?g0?gקV>O?g{r$?g#!?g1uxJ?g{YJk?g^?gۧۀ?gNk?gO?g4?g%gQ?g^m9?g7?g8?g+Xx6b?g݋7K,?go_?gY$.?g7#3?g9>Y?gŅg?g=?gL'[?g?g?z?g2kM?gmsAD?gH?g/rm$?g?g+G?g> TX5?g?Q]J.?g$?ga_D6+?gJ(>?gOz|?g&ke"?g#֠7N?gE"?gyye?gx**?gpH^G?gф7?g~R)?g۶Z=?gzl?g1?g?g5v?gw?g-?g<?gޠ(?g&N?g:^?giՀv?g?gf_?g1b]䷋?g9F?g@L?gYԪ$?g=?g𽅖?gq1?gmfL8?gA?g*"?g:{>?gc?gD?g9bNL"?gC{6?gʳ*$?g<'^?go=S?gFXM?g?g]]_&?gv
g?g\O?gT?gM6X;j?g4,gkT"?g.D?gTX?gfT?g?g?gpGQ?gH똈?g&Y\?g?gw{8?g/ܹ0ҋ?g<a0?gv@V?gku+?g<C0?gkS?g8%(	?gjx?gp4?g1piʈ?gʨ2?g0ѥ?g!|o?g{W?g
x?gŎc?gw?gQg?gac?gq^?gg?gQ]
Y?g+Ļ?gm?g!?g[?g$P?gmq*u?g^?gmD?g]YxG?g랆?ggC/?gս_?g6?gO?g+J?ghqL?g r?gR	.L?gm?g\tL%?g8`|?g\?g5`ʿ?gDjV0?gģ\S?g^p?gb+?g@`Pl&?gAF^R?gW3N?gbX?g@O$2?gPPH:?g[B?gQwв?gϪ"?g{?g.:u?g \E?gت,#?gP?gy?g~%}?g4+(?gJm5-?g}Zc ?gJ4?gq?g#?gSY:	?gtv~l?g$0u?gc?g[Vi}?gw?g'Jd?gt_?gCKGΌ?gѬ?ge$?gl8T!?gn>?gTj ?g~q?gC?gDIWx?ghW?gils?gj[?gw< ;?g`ҭ@?gD?g.쒂x?gڬN?g
`x?gV7?g{y?gd?g=}%Z?gͳw?g9ފ?g{?g2>j?gCEq~?g9ѮBʏ?g?g#)/?gy?gÒzh?gL<?g@!?gדJ?g3?gS,+?gm;?gv!i$?g,T'?g:G?g!o]>?gLC?g:P1?ga#{
?gٿL ?gҀ?gCi5E?g-1v?g*Yg?g,~SX?gZQsfy?gf8?g#;<?g 7?g@ӧ?g-s?gkO-?gO*?g8?g'?g[?gȨ:?g-Fۄ?gN~?g0?g0<aO\?gROF?gJsy6?ggefs?)r3  r4  r5  r6  r7  r8  openr   
genfromtxtr<  r?  r@  r   r#   r$   r"   r  r'   r,   r+   r%   r&   r*   r)   r(   r  r   invgauss_residsr!   r    r0   )r2   rC  fdrD  s       r3   r4   zInvGauss.__init__  s   7<<0I0I J J 24 4(W--- 	E=s%@@@DD	E 	E 	E 	E 	E 	E 	E 	E 	E 	E 	E 	E 	E 	E 	E%4%(^
$O	 E:::	
 97.)%
*+'
'%5H h&L h&L h&L hM hMs   1%B""B&)B&Nr6   r;   r5   r3   r  r    s7         :DM DM DM DM DMr5   r  c                       e Zd ZdZd ZdS )Medpar1z_
    The medpar1 data can be found here.

    https://www.stata-press.com/data/hh2/medpar1
    c                    t           j                            t           j                            t           j                            t
                              d          }t          j        |                                          }|j	        | _
        t          j        |j        ddt                    }t          j        |j        |f                              t                    }t%          |d          | _        d S )Nzstata_medpar1_glm.csvr  Tr  Fr'  )r3  r4  r5  r6  r7  r8  r9  r:  
to_recordslosr?  r;  admityper<  r   rU  codesrV  r   r@  rB  s        r3   r4   zMedpar1.__init__  s    7<<0I0I J J 79 9{8$$//11X
.v$',. . .$*g!677>>uEE 777			r5   Nr6   r;   r5   r3   r  r    rW  r5   r  c                   "     e Zd ZdZ fdZ xZS )InvGaussLogz5
    InvGaussLog is used with TestGlmInvgaussLog
    c                 v   t                                                       t          j                            t          j                            t          j                            t                              d          }t          j	        |dd           j
        | _        d| _        t          j        g d          | _        t          j        g d          | _        d| _        d| _        d	| _        d
| _        d| _        d| _        d| _        d| _        d| _        t          j        g d          | _        d S )Nzmedparlogresids.csvr  sepheaderrlvt@)g!j?gSȿg'0/?)g<4x?g|܎?g(I21+?gru@gIJy@g_9Ys@gCG7?g\gWa&gFx@r   Y  (\  f!@r  r  r  7h@r  r  e'y@r  r  r  r  r  r  r  Z@r  r  甝ev@r  r  r  r  r  r  r  X@r  r  r  r  r  r  r  r  r  4@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  :@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  }m	@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  F@:<&@9l@K7j>@79@r  r  r  r  r  r  r  r  Z!)r:@r  r  r  r  r  MɕQ/@r  r  r  r  r  r  r  r  r  Eaʚ	@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  d@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  )r  r4   r3  r4  r5  r6  r7  r8  r9  r:  valuesr!   r"   r   r    r#   r$   r%   r&   r'   r(   r)   r*   r-   r+   r,   r0   r2   rC  r  s     r3   r4   zInvGaussLog.__init__  s   7<<0I0I J J 57 7k(DAAAH.hEEEFF8BBBCC'
**(
%+%	H `& `& `& ` `r5   r  r  s   @r3   r  r    sK         r r r r r r r r rr5   r  c                   "     e Zd ZdZ fdZ xZS )InvGaussIdentityzN
    Accuracy is different for R vs Stata ML vs Stata IRLS, we are close.
    c                    t                                                       t          j        g d          | _        t          j        g d          | _        t          j                            t          j        	                    t          j        
                    t                              d          }t          j        |dd           j        | _        d| _        d| _        d| _        d	| _        d
| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        d| _        t          j        g d          | _        d S )N)g3C>?gq!ɋg۪A@)gTy}?g/pjҳ?gm{+?zigaussident_resids.csvr  r  r  i[  g+bs@r  r   g"X@g|7'z@g&g og{@g镲gୗu8?gt\8w@(\  4Sd@r  r  r  %+@r  r  ?@r  r  r  r  r  r  r  x@r  r  D@r  r  r  r  r  r  r  [{@r  r  r  r  r  r  r  r  r  S:ܳ@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  '5@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  hG
@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  SX@@By`@}@ϻ0(@r  r  r  r  r  r  r  r  :{@r  r  r  r  r  J_@r  r  r  r  r  r  r  r  r  far@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  G #@r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  )r  r4   r   r    r#   r$   r3  r4  r5  r6  r7  r8  r9  r:  r  r!   r"   r  r'   r,   r+   r%   r&   r*   	llf_statar-   r)   r(   r  r0   r  s     r3   r4   zInvGaussIdentity.__init__  s%   hEEEFF8BBBCC7<<0I0I J J 8: :k(DAAAH.*'
*++%	(
-H `& `& `& ` `r5   r  r  s   @r3   r  r    sK         u u u u u u u u ur5   r  c                       e Zd Zd ZdS )	Committeec                    t          j        g dg dg dg dg dg dg dg dg d	g d
g dg dg dg dg dg dg dg dg dg dg          | _        d| _        t          j        g d          | _        t          j        g d          | _        d| _        d| _        d| _        d| _	        d| _
        d| _        d| _        d| _        d | _        d!| _        t          j        g d"          | _        d S )#N)gdق(gU&go1Tgg@,{)goͿgOXпg1*kg{YHgֈ`)g%VFc?gt:D?g@;~@g]P2@gσv6@)g ms>ҿgjVԿg(0g8Kr
g۷a)g<l?gȥV0?gJ+@g9ahVv?g.Q5@)g~p>uٿg#fyr޿gp=
gg[Ɏ<V)gſgevqvRƿgg)yugJ)ggK~ۿgQngQ^ggYG)g%?gj(?gC@g$?gM֨GB@)gG2;??gw$[?gV-Ru@gYV^V?gMbX@@)g8^ud?gHڍ>?g㥛Ġt@g{%9`W?gctv2(8@)g?gӡ.a?gL7A@gߢ?gZ/r;@)g<,?gME?g@g[/?gGZ*oG;@)g	?J?g1zZ?g1Zq@gr-Z @gS@)g<8g6gbX9]gݘģgUDA)g=	ObgP4gͪ+qg\>g9#*)gn`zտg2|ؿgvg~֤RgD -)gD+g?gK1?g`vO@g#?gpn(@)gI5o?gg?gZ-?gT#9?gP黶?)gzچVIgl@rgpgXXp?Agz]zK,gŠ;@)gL5Pug;?g I@gvEֿg{?gVԿg)=;)gĻy?g@9?g\	c?g]ʿ?gUx/W?gZӖ?gݯ?gvdMk@gM؊%@g(Jv@gn?gvdMUYgN`MUYg R@g2yf@r      )g=)@gD7.>@g9{5@gv@g\Y@g׀JM~k@gFt@g9lS#Ѵ[@g}xTQ\@gG{k@g[0lL@g&KeN@g\@g"ǵ%R@gM'N@gHT"5@g|3aE@gJK@gnf5?gfх%2@)r   r    r!   r"   r#   r$   r%   r&   r'   r(   r)   r  r*   r-   r+   r,   r0   r1   s    r3   r4   zCommittee.__init__  s   h0 0 00 0 0/ / /0 0 0/ / /0 0 00 0 00 0 0/ / // / // / // / // / // / /0 0 00 0 00 0 0/ / // / /0 0 0O( 1 (2 (2R /h  @  @  @ A A 8      (
*)(
&++%	H &F &F &F G Gr5   N)r7   r8   r9   r4   r;   r5   r3   r  r    s-        AG AG AG AG AGr5   r  c                       e Zd ZdZd ZdS )WfszN
    Wfs used for TestGlmPoissonOffset

    Results are from Stata and R.
    c                     t           j        | _        d| _        g d| _        g d| _        d| _        d| _        d| _        d| _	        d| _
        d	| _        d
| _        d| _        d| _        g d| _        d S )Ngd3n'@)gPG?gݼ?gيG?g6G?g\BZ?gvǼ?g^7i?g_8c/?gUgj}ӿgX)g#8?gL"%?gq76?gc058?gD^?gd{Z?g^nT?g163?g;d?gNǭD?gf8TBf?g##Q@g|	91@gKQ@r   g_bEFBogÓfr  ;   gv~6F@)Fg"v@g>"D3@gL7A`@@go%;@@gCR'@gbg
;@gYni5$A@gqB@g}˜.P@gQkw[@gK4Y@geI)A@gvۅ8@gR@gk+#J@gSqC@g6TA@g䃞ͪY@gMbZ@gk	D@goʋh@g$~t@gCi@g-ęE@g0_^D@gr))X@gBO@gŏ1w->@g}iƢQ@g鷯e@gHY@g3·E@g|?5^dv@g//@g-f@gW>:@gqhDO@gDla@g{<J@g~4*@g6>W[\@gV-j@gvӂV@gqS2@g˄@gh |@gqa@gQi@gj+[@g'hX@guL@g=~o3@gSt$O`@gsAb@gZkQ@gb *@gҍ@gfjz@g;pΈ`@gcZB>"@gUNr@g(+f@gp>?4G@g
|@gZӼt@g\mt@gU@gJ+Ö@gGz0|@g-P@)r   
wfs_residsr!   r"   r#   r$   r%   r&   r'   r(   r)   r*   r+   r,   r-   r0   r1   s    r3   r4   zWfs.__init__  s    %0*  H H H "
*'
%+ &	= = =r5   Nr6   r;   r5   r3   r  r    s-         
&= &= &= &= &=r5   r  c                       e Zd ZdZd ZdS )CpunishTweediePower15a  
    # From R
    setwd('c:/workspace')
    data <- read.csv('cpunish.csv', sep=",")

    library(statmod)
    library(tweedie)

    summary(glm(EXECUTIONS ~ INCOME + SOUTH - 1,
            family=tweedie(var.power=1.5, link.power=1),
            data=data))
    c                     g d}g d}g d}g d}|| _         || _        || _        || _        ddg| _        ddg| _        d	| _        d
| _        d| _        g d| _	        d S )N)gZ<@gn5I?g}A@g2+gvT*Qgi?gǉ?gQtg8gZX'ggXgg3gs/ggXAW^g:p%"ٿg9to)g4@g/Pa?g|[S( @g̣IgUg.W㻶<?g#by.?g@1g&R^g=]g@rĊ3g61ܿg(w`g(
Rgo߿g^<rտgy7)gH~@gI,?g9.Q@gczgodvg*ۿ?gC!g?gM%YgUg6Rg9]!]߿g<a@ٿgY^g~Lg-'ܿgPR$ӿg#B
߿)gK<@g$"I?gYlNA@g
0gJSg$y?g ?g1
tgkg@7"*g|Vgg<3gUk2gD~g^gv]pٿgԑ;nog'Ƌl;?gv@g9jY>g ;@g[,B@r      )gZ$0 @gȤl @gȜ?gǕ@gw*Q@g_{K?gL:?g(}: @gEh	 @gX'@gST,3?g-u?g./@g]4xF@gz+J/?gF]d?gRn?
resid_responseresid_devianceresid_pearsonresid_workingr#   r$   r'   r+   r,   r0   r2   
resid_resp	resid_dev
resid_pear
resid_works        r3   r4   zCpunishTweediePower15.__init__7  s    8 8 8
4 4 4	     
4 4 4
 )'''#\2 ,/ ( 
	 	 	r5   Nr6   r;   r5   r3   r
  r
  *  s-         ; ; ; ; ;r5   r
  c                       e Zd ZdZd ZdS )CpunishTweediePower2a  
    # From R
    setwd('c:/workspace')
    data <- read.csv('cpunish.csv', sep=",")

    library(statmod)
    library(tweedie)

    summary(glm(EXECUTIONS ~ INCOME + SOUTH - 1,
            family=tweedie(var.power=2, link.power=1),
            data=data))
    c                     g d}g d}g d}g d}|| _         || _        || _        || _        ddg| _        ddg| _        d| _        d	| _        d
| _        g d| _	        d S )N)g_<@gNþ]?gr;@ggqEg!!y?g1~?g1mQgqxW,gZ,gu<<œgZgg6<g')ug=@g:gxmٿgѪ徂)gO[4C @g4Tu?g?g`)\ggW悇?g9/?g7փz'gsg!g8fgp]ڿgyₔggg\MݿgPԿg#p1)gV4@gZt?g@lA]`@gΙOBc޿gg^goh9?g61B?g^u"(gSFgOt6gd
;ۿg#%S׿g4|g9gܩٿgь1gҿg<*ۿg1<Ej?gZp к@gio>ge@gňnq/@r   r  )g71 @gT# @g6?gW@gqE@g!?g]?gVOڶ( @g9+e @gZ,@g:I?gx|V-?gg6<@gS@gRC?g^[lu?gWr_A?r  r  s        r3   r4   zCpunishTweediePower2.__init__  s    4 4 4
4 4 4	4 4 4
4 4 4
 )'''+-=>(*:; ) 
	 	 	r5   Nr6   r;   r5   r3   r  r  u  s-         < < < < <r5   r  c                       e Zd ZdZd ZdS )CpunishTweedieLog1a  
    # From R
    setwd('c:/workspace')
    data <- read.csv('cpunish.csv', sep=",")

    library(statmod)
    library(tweedie)

    summary(glm(EXECUTIONS ~ INCOME + SOUTH - 1,
            family=tweedie(var.power=1, link.power=0),
            data=data))
    c                     g d}g d}g d}g d}|| _         || _        || _        || _        ddg| _        ddg| _        d	| _        d
| _        d| _        g d| _	        d S )N)g4$<@gӿgK@gx_z
gQjg$>?g‚!!?gu	UghegcgWezg!*a}g\SgЃGgqA53gy3[gGn'uZgDy$)g-|v8@ggޠ.Z@g]xqgkgzFM#?g;?gX} g$[ug'gE;g'^qgF	g"z^g^W*_g(g+=.)gBą#@geɹgt@ghc:gjPgWU?g$p?gW9Fg	gOgxaB^)gdygѼg)!gzg
V4߿ge)gWٔ@gC砿gʞ@gpqܿgDsl7g²WF?gAƎhc?gQgJRdzK1gZ}g^ݿg0`$ۿg	S
gR}gn0u	ܿgEiؿg._4ݿg;_>gKWh?go>g9	VB9?g̀/|W@r   r  )gō @gX"@gZur ?gy/=@gQj@gxxp?g,f~?g:ꄪ @gq31 @gWez@go?gh.e?gЃG@g њ@g&-t?ga:-?g<Qm?
r  r  r  r  r#   r$   r'   r+   r,   r0   r  s        r3   r4   zCpunishTweedieLog1.__init__  s    4 4 4
3 3 3	  
4 4 4
 )'''+-=>(*;< ) 
	 	 	r5   Nr6   r;   r5   r3   r  r    s-         > > > > >r5   r  c                       e Zd ZdZd ZdS )FairTweedieLog15aH  
    # From R
    setwd('c:/workspace')
    data <- read.csv('fair.csv', sep=",")

    library(statmod)
    library(tweedie)

    model <- glm(affairs ~ rate_marriage + age + yrs_married -1, data=data,
             family=tweedie(var.power=1.5, link.power = 0))
    r <- resid(model, type='response')
    paste(as.character(r[1:17]), collapse=",")
    r <- resid(model, type='deviance')
    paste(as.character(r[1:17]), collapse=",")
    r <- resid(model, type='pearson')
    paste(as.character(r[1:17]), collapse=",")
    r <- resid(model, type='working')
    paste(as.character(r[1:17]), collapse=",")
    paste(as.character(model$coefficients[1:17]), collapse=",")
    s <- summary(model)
    paste(as.character(sqrt(diag(s$cov.scaled))), collapse=",")
    s$deviance
    paste(as.character(model$fitted.values[1:17]), collapse=",")
    c                     g d}g d}g d}g d}|| _         || _        || _        || _        g d| _        g d| _        d| _        d| _        d	| _        g d
| _	        d S )N)gPNЉg%@gEce=?gur)?g?i]5@g!E@gaj?g$vJ?gr@gwxkg?gu4?g|֯o?gaпg 3?gx8@gZ?gG~?g)g@qg)5@gEH ?gDhmR?g;o@g
@gR7??go@gL<@gstEc?g^Ti?g-P?g^ξ=׿g^m?gsa<}@gࢄ3?gs+?g )g	Tg+A%@gLk?ghdm?g5HZ"@gb@O@g6%@g@gYPv@g*鲤?g? 'n?g3OQ@gԿgY@g;@gnU&?g}[(=?g )g)h9g/[@g ܅?gWO?gl"3'@g 2@g*G@g<ț!@g33u@gz~e?gF^?g%䌃	@g(U;ֿgu@guf@gs\<?gd#~?)gI#!ؿg4"N^(?g`ڸ)gT^gV?g{x"z?gx\L?gGztA@r   i  )g6&Wa?gmC6?穘?g`_;?g|f>RO?gM%3?bKh?r#  g#?ggy?g1?|0?gkM?r$  g>"7v?gnZ?gpr?r"  r  r  s        r3   r4   zFairTweedieLog15.__init__(  s    F F F
D D D	D D D
4 4 4
 )'''I I II I I
 ! 
	1 	1 	1r5   Nr6   r;   r5   r3   r  r    s-         0<1 <1 <1 <1 <1r5   r  )$r:   r3  numpyr   pandasr9  statsmodels.apir    statsmodels.genmod.tests.resultsr   	DECIMAL_4	DECIMAL_3	DECIMAL_2	DECIMAL_1	DECIMAL_0r
   r=   r   r  r  rH  rO  rY  r  r  r  r  r  r  r  r  r  r
  r  r  r  r;   r5   r3   <module>r.     s    
			         ( ( ( ( ( ( < < < < < < 					K! K! K! K! K! K! K! K!\7 7 7 7 7 7 7 7DnM nM nM nM nM nM nM nMb_3 _3 _3 _3 _3 _3 _3 _3DO? O? O? O? O? O? O? O?d"D "D "D "D "D "D "D "DJ8 8 8 8 8 8 8 8"FD FD FD FD FD FD FD FDRKD KD KD KD KDV KD KD KD\; ; ; ; ; ; ; ;B- - - - -W - - - bM bM bM bM bM bM bM bMJ!8 8 8 8 8 8 8 8"v v v v v' v v vry y y y yw y y yxBG BG BG BG BG BG BG BGJ,= ,= ,= ,= ,= ,= ,= ,=^H H H H H H H HVI I I I I I I IXK K K K K K K K\U1 U1 U1 U1 U1 U1 U1 U1 U1 U1r5   