U
    Kvf-                    @   s|  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
 d
ZG dd dZG dd dZG dd dZG dd dZG dd dZG dd dZG dd deZG dd deZG dd dZG dd deZG dd  d ZG d!d" d"ZG d#d$ d$eZG d%d& d&eZG d'd( d(ZG d)d* d*ZG d+d, d,Z G d-d. d.Z!G d/d0 d0Z"G d1d2 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                   @   s   e Zd ZdZdd ZdS )LongleyzM
    Longley used for TestGlmGaussian

    Results are from Stata and R.
    c                 C   s~  t dddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgd	d	d	d	d	gd
d
d
d
d
gdddddgdddddgdddddgdddddgdddddgdddddgg| _d| _t dddddddg| _t dddddddg| _d | _d!| _d"| _d#| _	d$| _
d| _d%| _d&| _d'| _d(| _d)| _d*| _t d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:g| _d S );Ngޭpp@gnWg2$G@g}աygn[s@gy)ogdgTW\*go{',@g<6Nv|@g2u6hD1gCg:5<bqcgUjUg{.{^u@gR@iig}⹭.@g|S7Vg<=n) gg|X=Q*gT@g
ףPIJg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self r    P/tmp/pip-unpacked-wheel-2v6byqio/statsmodels/genmod/tests/results/results_glm.py__init__   s     !                 zLongley.__init__N__name__
__module____qualname____doc__r"   r    r    r    r!   r      s   r   c                   @   s   e Zd ZdZdd ZdS )GaussianLogzC
    Uses generated data.  These results are from R and Stata.
    c              j   C   sl  t dddddgdddddgdddddgdddddgd	d	d
d	d	gdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgddd ddgd!d!d"d!d!gd#d#d$d#d#gd%d%d&d%d%gd'd'd(d'd'gd)d)d*d)d)gd+d+d,d+d+gd-d-d.d-d-gd/d/d0d/d/gd1d1d2d1d1gd3d3d4d3d3gd5d5d6d5d5gd7d7d8d7d7gd9d9d:d9d9gd;d;d<d;d;gd=d=d>d=d=gd?d?d@d?d?gdAdAdBdAdAgdCdCdDdCdCgdEdEdFdEdEgdGdGdHdGdGgdIdIdJdIdIgdKdKdLdKdKgdMdMdNdMdMgdOdOdPdOdOgdQdQdRdQdQgdSdSdTdSdSgdUdUdVdUdUgdWdWdXdWdWgdYdYdZdYdYgd[d[d\d[d[gd]d]d^d]d]gd_d_d`d_d_gdddaddgdbdbdcdbdbgdddddeddddgdfdfdgdfdfgdhdhdidhdhgdjdjdkdjdjgdldldmdldlgdndndodndngdpdpdqdpdpgdrdrdsdrdrgdtdtdudtdtgdvdvdwdvdvgdxdxdydxdxgdzdzd{dzdzgd|d|d}d|d|gd~d~dd~d~gdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddggd| _d| _t dddg| _t dddg| _d| _d| _d| _d| _	d| _
d| _d| _d| _d| _t ddddddddddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;gd| _d S (<  NgH0!5?gzJ}L?g\uxcJ?g#070a?g(M[Q3g)8Ig˪V?gMW_m?gmRvbgWDwg7N9R?ǵ/f?gƼqվgpP g&:,B%gɇE9gb%HgS|\gL R"glw4g*Ŏơ~Gg YgRig-Ѵ/gN5?gDe?F?ghM=ngX
7kgͮ{+B?g9mKNR?gh-6?g'ZTE?gbHz?gĭ(?gSIgˤrlXg5?gwH1@	D?g4W_KgॊXg?m1:gYFgon="?g1?74e/?gV0Kge@Wg<>g/>gTg	"g)H4"F?gV%Q?g+|:?g>ckdC?gŸ́p>gov>gmC{dJg_--Sg(Kg\ܐhSg\d=P?g]2%V?gQ~3a@g)yFgٙQ?grMdW?gsYTO?g"S?gQ^?g1HVb?gwS%%_gǷwbg!J?gEO?g*!P?gwۓR?g}	^Q?g%FYR?gMp9gh<;g{QP?gѵP?gս_Fg34Fg]c;JgcM*Jg~bU]g?߾\gxC?笢zB?g-0?gf-?glxz,CLgX$Hg_S$6gSb3g<Pg%F?g*ЬA?gTb8gOk2gz giZg3:{mZgmrSg%LbZS1?g}^(?gxs]]?gyABT?ggej"g4;g}!3g]FV#`)gG[hNgyCi7Cg_M?g_ .B?gՃ8řUg}Igot{!g[mDgM?gf%}@?g{gUId/gO g퐬_pG?gDXP8?gصݒ?g6I>g&fb_g?g>+NV?gηO؏^>?gi)',?g>Ph2?g滥 ?gљj?gjۋ>gKEb7gZa/$gWmURgDs)">gG6uU?g! @?g<kWg!8)4Agf:*0&govg;,G#PgDXP8ga4oX?g@?g	BgcKEg8ӅXA?g:m/'?gRdag}3֣Eg\+-eq>g)>gGR?g{эU5?g]0??g+ٱ!?g%?=AgO"g]a&gŭg7h>^g1cH>gVrz(+?g-C6
?gL!FLg D*g{gՃ8řgƍA?gM?g3$Igm[$gnq]g{/h7gxr+R?g;s	+?g?.q%Bg
^gĭH?g,gJF)|!?g*ø;?gFƴ?gf@O4g^W*^g@'SEg<ړg>U?'	*?g.9.Ig**,
g[v-e[?g{C/?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Ag׾@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   r   r   r    r    r!   r"   j   s                                                                                                          K                                                                    zGaussianLog.__init__Nr#   r    r    r    r!   r(   f   s   r(   c                   @   s   e Zd ZdZdd ZdS )GaussianInversezG
    This test uses generated data.  Results are from R and Stata.
    c              j   C   sn  t dddddgdddddgdddddgdddddgd	d	d
d	d	gdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgddd ddgd!d!d"d!d!gd#d#d$d#d#gd%d%d&d%d%gd'd'd(d'd'gd)d)d*d)d)gd+d+d,d+d+gd-d-d.d-d-gd/d/d0d/d/gd1d1d2d1d1gd3d3d4d3d3gd5d5d6d5d5gd7d7d8d7d7gd9d9d:d9d9gd;d;d<d;d;gd=d=d>d=d=gd?d?d@d?d?gdAdAdBdAdAgdCdCdDdCdCgdEdEdFdEdEgdGdGdHdGdGgdIdIdJdIdIgdKdKdLdKdKgdMdMdNdMdMgdOdOdPdOdOgdQdQdRdQdQgdSdSdTdSdSgdUdUdVdUdUgdWdWdXdWdWgdYdYdZdYdYgd[d[d\d[d[gd]d]d^d]d]gd_d_d`d_d_gdadadbdadagdcdcdddcdcgdededfdedegdgdgdhdgdggdididjdidigdkdkdldkdkgdmdmdndmdmgdododpdodogdqdqdrdqdqgdsdsdtdsdsgdududvdudugdwdwdxdwdwgdydydzdydygd{d{d|d{d{gd}d}d~d}d}gdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddggd| _d| _t dddg| _t dddg| _d| _d| _d| _d| _	d| _
d| _d| _d| _d| _t dddddddddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<gd| _d S (=  Ng"R@gRt@?g_Q+r*   gi :w&gŷ$?g_?g?'[g\{QFr)   gInsS?g[RjOgT,XgHaR?gyT5?gP G-gXr4g^!+?gd8?gtI.g??gm)2gڔu>WgC I?g+ٱ?gavtg]6:?g..*gp7gk$	%?gſF}3g ?gNgQݷ8?g^gS?gf(ϼ<g<*6F?g`Z8Z/gM~0Ug;s	;?gPl-?gZ1gqڝhIQ?g:O_3g(Yi0?g>$Eg[+:?g?yg.ԫv[g3:?gvqvRH?gQq)%g1QaJ+?g34gX=Jg²lh$?guC?gF%gn~hP?gAȘ&g6-2?g/bDGg9y@g=?g,-~1?g
FgɵWՁb*?gb:g>Җ)SgnBA#?gD &_Fg`=I?g  B?gx+Ng"/l(@?g!
g_L5?gʹ gF<IgU?g(>g&Jw¾gr$g-C6>g,?gy;>1;gu5hFg&(4
?gbHzIgb:?gpCP?gt><Kga2U0*S?g(M[QgaHIg_eG|?g}6gZQf>g/ز!S?gMbg;TYgr?gM?g/\cgy`[>gTM>gZ/-Sg
^
?gV8Q?gPgfLg1gkNu>g<ܴIgtR#>gV Y?gÔ/g[d8^g+,?gs'[?g<ړg!C}4g$z>gZ_&P?gAg1{vA?gs>gDėgۢ>gF+R%Igv!>gR[SF?gg% "g@cG_>g
OmI?gYVPhgՋ)0?g~վg3kB?g,g,ZgH?gtR#gGW:Kg
h`?>g4F&;aggg?g	hǲ)MgM >gQ?g	\?gkB](g<>gŭ?gBY[%gXCQ?giUMugSV#M?g]FV#`gK8,g`>gL¶ZgOne.>gν,,?g};hžgXf=2g>g<-?p'@g34>gVF?(gJ֗[>gyxρIg!>@>gOЁ@I?gl(1,߾gkZs2?gqdqAƾgrnOQ?g{fdg+4?gʔA{ƾgpDkLg,z>g@I(gTqs*>g3G5?g!+lěNžg7 >?gΗ;gߏ61g~>gI	fU?g8{gk0h|Ug?̔>gB^j7?gtݛ?,gCgP3bn>g'ZTEg>>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γuYCDAgp:;?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,   r   r    r    r!   r"   l  s                                                                                                          J                                                                                zGaussianInverse.__init__Nr#   r    r    r    r!   r-   h  s   r-   c                   @   s   e Zd ZdZdd ZdS )Star98z0
    Star98 class used with TestGlmBinomial
    c              1  C   sV  d| _ d| _d| _d| _d| _d| _d| _d| _d	| _d
| _	d| _
d| _d| _d| _d| _d| _tj| _tddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>g/| _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   Zstar98_residsr   r   r   r   r   r    r    r!   r"   ]  s                                                                                                                                                                                                                                                     zStar98.__init__Nr#   r    r    r    r!   r.   Y  s   r.   c                   @   s   e Zd ZdZdd ZdS )Lbwz\
    The LBW data can be found here

    https://www.stata-press.com/data/r9/rmain.html
    c                 C   sb  t jt jt jtd}t|}tj|j	ddt
d}tj||gdd}|j| _|ddd	d
ddddg }t|dd| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _tj| _ t!"dddd d!d"d#d$d%d&d'd(d)d*d+d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdcdwdxdydydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddNddddddddddddg| _#d S )Nzstata_lbw_glm.csvraceFprefix
drop_firstdtyper   )ZaxisageZlwtZ
race_blackZ
race_otherZsmokeZptlZ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{?gri?gRS?go/)7?gn?g?g9?g|}?g#L?gWL:m?gybx?g'k5K?ga-?gal?gt)*V?gtܩ?g]?gQw36?gyeJ?g˒2?gK?gQ-d&?gLd?gRگ?g[?g&v?gM/]E?g_7?g˽Wڲ?g`5Q?g(6Il-?g?g!0A3?g=՝?g?gmB?g$;?gq?gT?g k\(?g*ѫdt?g`8$9	?g6 ?gE7?gi^?gxJ?g?giG7'?gg=M{?g?>Z?g] T?g\tL?g凌f?gme?g;>?gМq"?gTi0	?g7?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#?g'na?g=$cRİ?g-b1?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?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_dummiesr4   floatconcatlowendogr   exogr   r   r   r   r   r   r   r   r   r   r   r   r0   r2   r   Z
lbw_residsr   r   r   r   r   filenamedataZdummiesZdesignr    r    r!   r"     s   
                                                                                                                                                         zLbw.__init__Nr#   r    r    r    r!   r3     s   r3   c                   @   s   e Zd ZdZdd ZdS )Scotvotez2
    Scotvot class is used with TestGlmGamma.
    c              "   C   s   d| _ d| _d| _d| _d| _d| _d| _d| _d	| _d
| _	d| _
d| _d| _d| _d| _tj| _tddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/g | _d S )0N)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   r0   r   r   r   r   r   r   r   r   r2   r   r   r   Zscotvote_residsr   r   r   r   r   r    r    r!   r"     sd                            zScotvote.__init__Nr#   r    r    r    r!   rQ     s   rQ   c                   @   s   e Zd ZdZdd ZdS )Cancerz`
    The Cancer data can be found here

    https://www.stata-press.com/data/r10/rmain.html
    c                 C   sx   t jt jt jtd}tt|d}|j	| _
tjtj|jdddd}t|j|ft}t|dd| _d S )	Nzstata_cancer_glm.csvrbcategory)r8   T)r7   Fr;   )r@   rA   rB   rC   rD   rE   r   Z
recfromcsvopenZ	studytimerL   rF   rH   ZSeriesZdrugcolumn_stackr9   astyperI   r   rM   rN   r    r    r!   r"   8  s    zCancer.__init__Nr#   r    r    r    r!   rU   2  s   rU   c                       s    e Zd ZdZ fddZ  ZS )	CancerLogz+
    CancerLog is used TestGlmGammaLog
    c              6      s  t t|   tdddddgdddd	d
gdddddgdddddgdddddgdddddgdd d!d"d#gd$d%d&d'd(gd)d*d+d,d-gd.d/d0d1d2gd3d4d5d6d7gd8d9d:d;d<gd=d>d?d>d@gdAdBdCdDdEgdFdGdHdIdJgdKdLdMdNdOgdPdQdRdSdTgdUdVdWdXdYgdZd[d\d]d^gd_d`dadbdcgdddedfdgdhgdidjdkdldmgdndodpdqdrgdsdtdudvdwgdxdydzd{d|gd}d~dddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgg0| _d| _tddddg| _tddddg| _d| _	d| _
d| _d| _d| _d| _d| _d| _d| _tdd dddddddddddd	dd
dddddddddddddddddddddddddddd dd!dd"g0| _d S (#  N27{H{9yQ_FҞC痮`r" [u? ?x~:fE.8,4k_Hߠ=6XD糕OgyܭҾ^B&9;繍&RF\"N'6FPeZ_$	C\,J琢C蹅DͿklexϿ;#5Ͽ4M3ȿBʿ~W[bȠʿ8Q?v2=ֿ=yX5ڿoi5 '*ڿ#ӡXT$[ ﭿLﭿ^a=$߿;?2g?@"?|g?(@ ?
%U ˿V5RͿ;%5iTNͿE4~5獛h>ӿ ?rֿ7dD_q^Rkֿf3P>タ;r30R0)g`^6iR
?}?Oz12@,J	z? @y)?E?ص@˂??ؼm1?T?ʡA@9}?9@Wm?S	hy?>xEG@n<՜r?)yu@h}?W\wN?_LIO@-4?L~@9?q`?a/?K+\@jkD0?L
F%+@!Va@?iH?oM`@{?m~)@[Nֿ҆hٿc('U=Vٿ+Lk	=#ٿ
Ϳ޿,+MJADG6Ws޿`JBk}߿&p{wT\WS66GgX<"ٿed_hd޿(rVB:}Q޿נ/IB0X{A$Pkw`kLϓ =%]wgH糁jZ_$E)pUjS[ 	E-n?M*?t	1@v?aL{?5,vοWnпߡ(AW3kп硅.?aR||Bv?jWVDO@DPg?Y(@dmlvu䚧ߢ环=$P,{ӖߢqtyS4?_FE?p:ǲP@Bus?Ku/@箕u?WJ;?EP\@1Z웠2?p @ z*L?7*?DR*$@ǞP?F?5:L?H|o@4y;?E>+@2"r_	c~h?RDU_Ahj3Jsuj?ߊ{%T=O?o*Ra0ZZ Ŀqx(ƿSrP+]&ƿ'󏾉
!rhmۿa<,Ԛ}1kV׻,k?P'?hC@PYJ?+~w?ǮY˿=(Ϳ$jeDXCͿm},3452 ֐&" %]U"ڿ7-:?#?1殩[@P+?n@~nȿ51ʿH}omhI3̿ʿq / CE?7t?+ebg@nr?q?zM@ֲx?y]?`"p@DP?W])@-qU?~J
?H}8q@}}m?)*.@cz?^kz?M*k@.hk(??}i?lV}
r@PI?T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@gj5#@gwʯ@gGY\oݭ@gh=nM$@g- @gӚIH@gM7@g`	@g8'@g"^4&@gEg!@gb2n@gMIs;%@g8\7: @gk"@gÛ[%$@g婻p+@gE$.@g(I04@g`fL(@g+v
'@g·3@gFe/@g:4~6@g]Z2@gHnj9@g5e:@gMDQ4@g^*f@@gh86@g?@gUм5=7@g<bA@g.M8@g=~NŶ2@gR=@)superr[   r"   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   	__class__r    r!   r"   G  s                                                   a                                       zCancerLog.__init__r$   r%   r&   r'   r"   __classcell__r    r    rK  r!   r[   C  s   r[   c                       s    e Zd ZdZ fddZ  ZS )CancerIdentityz:
    CancerIdentity is used with TestGlmGammaIdentity
    c              6      s  t t|   tdddddgdddd	d
gdddddgdddddgdddddgdddddgdd d!d"d#gd$d%d&d'd(gd)d*d+d,d-gd.d/d0d1d2gd3d4d5d6d7gd8d9d:d;d<gd=d>d?d>d@gdAdBdCdDdEgdFdGdHdIdJgdKdLdMdNdOgdPdQdRdSdTgdUdVdWdXdYgdZd[d\d]d^gd_d`dadbdcgdddedfdgdhgdidjdkdldmgdndodpdqdrgdsdtdudvdwgdxdydzd{d|gd}d~dddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgg0| _tddddg| _tddddg| _d| _d| _	d| _
d| _d| _d| _d| _d| _d| _d| _tdd dddddddddddd	dd
dddddddddddddddddddddddddddd dd!dd"g0| _d S (#  Nr\   r]   r^   r_   r`   ra   rb   rc   rd   re   rf   rg   rh   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   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r  r  r  r  r  r  r	  r
  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r   r!  r"  r#  r$  r%  r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r0  r1  r2  r3  r4  r5  r6  r7  r8  r9  r:  r;  r<  r=  r>  r?  r@  rA  rB  rC  rD  rE  rF  rG  g.g*nP@g240@g׋{C@g}(?gΑy@g.1@gl|b @gFAiFt@gX<@gW0.@gB?rH  geeg͝gcr   rI  g1=2}I@g<@g@X?@gz/"@gT	&@gK!@g\X"@g=~@goH@gv/N)@gq9;(@gm=;"@g]@gc$('@g-J\ @g,s%)"@gIqe%@g`,@g]t.@gV+ 3@g9])@g:(J(@gCdn2@g |c/@gUo34@gr1@g鋭9@g(عH%,:@gV6@g{<@gNUF8@gO<R<@gu8@gd=@g59@gag5@gUʄ?;@)rJ  rO  r"   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   rK  r    r!   r"     s                                                   b                                         zCancerIdentity.__init__rM  r    r    rK  r!   rO    s   rO  c                   @   s   e Zd ZdZd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                 C   s   d| _ d| _d| _d| _d| _d| _d| _d| _d	| _d
| _	d| _
d| _d| _d| _tj| _tdddddddddddddddddg| _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   r0   r   r   r   r   r   r   r   r   r   r2   r   Zcpunish_residsr   r   r   r   r   r    r    r!   r"   `  sD                 zCpunish.__init__Nr#   r    r    r    r!   rP  Z  s   rP  c                       s    e Zd ZdZ fddZ  ZS )Cpunish_offsetzS
    Same model as Cpunish but with offset of 100.  Many things do not change.
    c                    s   t t|   d| _d| _d S )N)g9m4&rQ  rR  rS  rT  rU  rV  )gJ+@g1?gjMS?g)Wx?g      ?gBfj?gA`"@)rJ  rY  r"   r   r   r   rK  r    r!   r"     s    zCpunish_offset.__init__rM  r    r    rK  r!   rY  {  s   rY  c                   @   s   e Zd ZdZd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                C   sBM  t jt jt jtd}t|ddd}tj|dt	ddd  }W 5 Q R X |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| _d| _tj| _tdddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐdddddddddddddddddddddddddddddddd	d 	d	d	d	d	d	d	d	d	d		d
	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d 	d!	d"	d#	d$	d%	d&	d'	d(	d)	d*	d+	d,	d-	d.	d/	d0	d1	d2	d3	d4	d5	d6	d7	d8	d9	d:	d;	d<	d=	d>	d?	d@	dA	dB	dC	dD	dE	dF	dG	dH	dI	dJ	dK	dL	dM	dN	dO	dP	dQ	dR	dS	dT	dU	dV	dW	dX	dY	dZ	d[	d\	d]	d^	d_	d`	da	db	dc	dd	de	df	dg	dh	di	dj	dk	dl	dm	dn	do	dp	dq	dr	ds	dt	du	dv	dw	dx	dy	dz	d{	d|	d}	d~	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	dÐ	dĐ	dŐ	dƐ	dǐ	dȐ	dɐ	dʐ	dː	d̐	d͐	dΐ	dϐ	dА	dѐ	dҐ	dӐ	dԐ	dՐ	d֐	dא	dؐ	dِ	dڐ	dې	dܐ	dݐ	dސ	dߐ	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d	d
d 
d
d
d
d
d
d
d
d
d	
d

d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d 
d!
d"
d#
d$
d%
d&
d'
d(
d)
d*
d+
d,
d-
d.
d/
d0
d1
d2
d3
d4
d5
d6
d7
d8
d9
d:
d;
d<
d=
d>
d?
d@
dA
dB
dC
dD
dE
dF
dG
dH
dI
dJ
dK
dL
dM
dN
dO
dP
dQ
dR
dS
dT
dU
dV
dW
dX
dY
dZ
d[
d\
d]
d^
d_
d`
da
db
dc
dd
de
df
dg
dh
di
dj
dk
dl
dm
dn
do
dp
dq
dr
ds
dt
du
dv
dw
dx
dy
dz
d{
d|
d}
d~
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
dÐ
dĐ
dŐ
dƐ
dǐ
dȐ
dɐ
dʐ
dː
d̐
d͐
dΐ
dϐ
dА
dѐ
dҐ
dӐ
dԐ
dՐ
d֐
dא
dؐ
dِ
dڐ
dې
dܐ
dݐ
dސ
dߐ
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
dd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddÐdĐdŐdƐdǐdȐdɐdʐdːd̐d͐dΐdϐdАdѐdҐdӐdԐdՐd֐dאdؐdِdڐdېdܐdݐdސdߐddddddddddddddddddddddddddddddddd ddddddddd	d
dddddddddddddddddddddd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfdgdhdidjdkdldmdndodpdqdrdsdtdudvdwdxdydzd{d|d}d~ddddddddddddddddddddddddddddddddddg| _d S (  Nzinv_gaussian.csvrzutf-8)encoding,)	delimiterr8   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?)r@   rA   rB   rC   rD   rE   rX   r   Z
genfromtxtrI   rL   rM   r   r   r   r   r0   r   r   r   r   r   r   r   r   r2   r   Zinvgauss_residsr   r   r   )r   rO   fdrP   r    r    r!   r"     sN'   "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       zInvGauss.__init__Nr#   r    r    r    r!   rZ    s   rZ  c                   @   s   e Zd ZdZdd ZdS )Medpar1z_
    The medpar1 data can be found here.

    https://www.stata-press.com/data/hh2/medpar1
    c                 C   sp   t jt jt jtd}t| }|j	| _
tj|jddtd}t|j|ft}t|dd| _d S )Nzstata_medpar1_glm.csvr4   Tr5   Fr;   )r@   rA   rB   rC   rD   rE   rF   rG   Z
to_recordsZlosrL   rH   ZadmityperI   r   rY   codesrZ   r   rM   rN   r    r    r!   r"     s    zMedpar1.__init__Nr#   r    r    r    r!   r`    s   r`  c                       s    e Zd ZdZ fddZ  ZS )InvGaussLogz5
    InvGaussLog is used with TestGlmInvgaussLog
    c              ^     sj  t t|   tjtjtjtd}t	j
|dd dj| _d| _tdddg| _tdd	d
g| _d| _d| _d| _d| _d| _d| _d| _d| _d| _tddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd d!d!d!d!d!ddd!d!d"dd d!ddd#d!dd!d#d!dddd!d$d!d!dd!ddd!d!d!dd!dd"d!dd!d%d$d!d!d"d$d!d!d!d!d!d!d#d!d!d"dd d!d!dd!dd!dd#dd d!d d!d!dd d!d!d!d!d!dd!ddd!d!dd!dd!d d!d!dd!d$d!d!d$d"d!d!ddd$d!d d!d!d!d$d!d!d!d!d d!d"d!d d d!dd"d"dd"d!d!d!d!d!d!d!ddd!dd!d!d!d#d!d!d!d!ddd!d#d!dd!d!d d$d!dd!d#d#dd"d!d!dd!d!d"d!ddd!d#dd!d!d dd!d!d!d!d#d!dd"dd!d$d#d!d d!dddd!dd"d!d$d!d!ddddd!d!d#d$d!d!d!d!ddd!d!dd!d dd!d d d!d!ddddd d!d"d!dd!d!d!d d"d!dddd!d!ddd$dd"d!d dd!d!d!d"dd#d!d!d!ddd!d!d!d!d!d!d!d"ddd!dd"d!d!dd!dd!d#d!dd"d"d!ddd!ddd!d"d d!dd"ddd#d!ddd d!d!d!d!d!d"d!d d!d!dd!d!d!d!ddd!dd!d!d!d!dd d!dd!dd!d!d!d!d ddd#d!d!d$d!d d!d!d!d!d"d"d"d!d!d!d!dd$d d!d"d!d!d!d!d!d"d!d!d!d!d"d!d"dd!dd!d!d!d!d!d!d!d!ddd!d!d!d!d!d!d!dddd!dd!d"d!d!d!d"d$d!d!d"d!d!d!d!d!d!d!d!d!d!d!d!d!dd!d!d!d!ddd"d!ddd"d!d!d d#d!d!d"d!dd"d!d!dd!d ddd$dd!d d!dd"ddd!d!dd"d!d!d!d!d!d!d!dd d"dd dd"d!d!d!d!d$d!d!d!d!d!d"dddd$d!d"d!dd!d!d"d!d"d!d!d!d!dd"ddd#d!d d!dd!dd!d!d!dd#d!d!d!d"dd!dddd!d!dd!d!d!d#dd!d!dddd!d!d d$dd!d!d!d!d!dd!d!d"dd!dd!dd#d!dd!d!d!dd!d!d!d!d!d!d d!d$d!dddd!dd!d"d d!d!d!d!dd!d"dd"d!dd!d!d!d!dd!dd!d!d!d!d!d!d!d!d"d%d!d!d!d!dd!dd!d!d"d!d!d!d!d!d!d$d$d!d!d!d"d!dd dd d dd!dd!dd!d!d!d!d!d$d!d"d!dd!d"d!d!d!d!d!dd!d dd d!d"d!ddd!d!d d"dd!d!d!d d!d!ddd d$d!dd$d!d!d$d$d!d d d!d!d"d"dd!d"dd!d"dd!dd!dd!d d"d!d!d dd!d!d!d!d!d!dd$dd#d"d!dd!dd!d#d!d!d!d!d$d!ddd dd"d!d!dd"dd d!d!d!d!d!d"d!d!d!d!d!d!ddd!d!dd!d$d!dd!dd!dd!d!d!d!dd#d!d!d#dd!d!d#d dd d!d!d!d!d!d!d ddd!d!d!d!dd"d!dd!d#d#d!dd!dd!d d#d!d!d#dd!d!d$d!d"d!d!d$d"ddd!d"d d!d d!d$d ddd!d!dddd!d!d"d"d$d!d!d#d"d!d#ddd d!dd!d!d!d!d!d!dd!d!d!d!d!d!d!d!dd!d!dd!d!ddd!dd d!dd"d d!d#d dd!d!d!d!d$dd!dd!d d!d!dd d!d!d!d!d!d d d$dd!d$d$d!d"dd$dd$ddd!d!d dd!d d!dd d#d!d!d!dd!d!d!d!ddd d!d"ddd!d d!d!d$d!dd!d!d!dd%d!d!ddd"d!d#d!d"d!d!d"d#d!dd d!dd!d!d d!d!d!d"dd!dd!ddd"d!d!d!dd!d!dd!d!d!d!d!d!d!d!d!d!ddd"d!d dd!ddd$d"d!dd!d!d!d!d dd d!d$d!d d dd!dd d!ddddd!d"dd d dd!d"d!d!ddd$d"d!dd!d!dd!d ddd!d!ddd!d!d!ddd!d$d!d!dd!d!d!d d!dd!d!d!d"d dd!dd"dd d!ddd#d!d!d!d!d!dd!dd!d!d!d d!d$d"d!d!dd!d!d!d!d!d!d!d%d!d"ddd!d!d"d!dd#d!dd#dd!d!d!d d!dd"d#dd!d%dd ddd$d d d$d dd dd"d d!d!d!d!d"d!d!d!d d!d!d!d"dd"d dd"d!d!d!d#d!d!d!d!d!d!d!d"d$dddg\| _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  gf!@g7h@ge'y@gZ@gev@gX@g4@g:@g}m	@gF@g:<&@g9l@gK7j>@g79@gZ!)r:@gMɕQ/@gEaʚ	@gd@)rJ  rb  r"   r@   rA   rB   rC   rD   rE   rF   rG   valuesr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   rO   rK  r    r!   r"     s                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    zInvGaussLog.__init__rM  r    r    rK  r!   rb    s   rb  c                       s    e Zd ZdZ fddZ  ZS )InvGaussIdentityzN
    Accuracy is different for R vs Stata ML vs Stata IRLS, we are close.
    c              ^     s|  t t|   tdddg| _tdddg| _tj	tj
tjtd}tj|dd d	j| _d
| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _d| _tdddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd d!d"d#d$d$d$d$d$d"d"d$d$d%d"d#d$d"d d&d$d"d$d&d$d"d d!d$d'd$d$d!d$d d d$d$d$d!d$d!d%d$d d$d(d'd$d$d%d'd$d$d$d$d$d$d&d$d$d%d d#d$d$d!d$d d$d!d&d"d#d$d#d$d$d!d#d$d$d$d$d$d!d$d"d"d$d$d d$d d$d#d$d$d"d$d'd$d$d'd%d$d$d"d d'd$d#d$d$d$d'd$d$d$d$d#d$d%d$d#d#d$d!d%d%d!d%d$d$d$d$d$d$d$d"d d$d d$d$d$d&d$d$d$d$d!d!d$d&d$d!d$d$d#d'd$d!d$d&d&d"d%d$d$d"d$d$d%d$d"d!d$d&d"d$d$d#d"d$d$d$d$d&d$d d%d"d$d'd&d$d#d$d d"d!d$d d%d$d'd$d$d d"d!d!d$d$d&d'd$d$d$d$d"d"d$d$d!d$d#d"d$d#d#d$d$d d!d d!d#d$d%d$d!d$d$d$d#d%d$d!d!d d$d$d"d d'd d%d$d#d!d$d$d$d%d!d&d$d$d$d!d"d$d$d$d$d$d$d$d%d"d!d$d!d%d$d$d!d$d"d$d&d$d"d%d%d$d"d d$d!d"d$d%d#d$d"d%d d"d&d$d d!d#d$d$d$d$d$d%d$d#d$d$d"d$d$d$d$d"d"d$d"d$d$d$d$d"d#d$d"d$d"d$d$d$d$d#d"d d&d$d$d'd$d#d$d$d$d$d%d%d%d$d$d$d$d!d'd#d$d%d$d$d$d$d$d%d$d$d$d$d%d$d%d d$d!d$d$d$d$d$d$d$d$d!d"d$d$d$d$d$d$d$d"d d"d$d d$d%d$d$d$d%d'd$d$d%d$d$d$d$d$d$d$d$d$d$d$d$d$d d$d$d$d$d d"d%d$d!d!d%d$d$d#d&d$d$d%d$d d%d$d$d!d$d#d"d"d'd d$d#d$d!d%d!d!d$d$d"d%d$d$d$d$d$d$d$d d#d%d!d#d d%d$d$d$d$d'd$d$d$d$d$d%d!d"d d'd$d%d$d d$d$d%d$d%d$d$d$d$d!d%d!d"d&d$d#d$d!d$d"d$d$d$d d&d$d$d$d%d!d$d"d!d"d$d$d"d$d$d$d&d"d$d$d!d!d"d$d$d#d'd d$d$d$d$d$d"d$d$d%d!d$d d$d"d&d$d"d$d$d$d"d$d$d$d$d$d$d#d$d'd$d"d!d"d$d!d$d%d#d$d$d$d$d"d$d%d!d%d$d d$d$d$d$d d$d!d$d$d$d$d$d$d$d$d%d(d$d$d$d$d!d$d d$d$d%d$d$d$d$d$d$d'd'd$d$d$d%d$d"d#d"d#d#d!d$d!d$d d$d$d$d$d$d'd$d%d$d!d$d%d$d$d$d$d$d"d$d#d"d#d$d%d$d d"d$d$d#d%d d$d$d$d#d$d$d"d!d#d'd$d d'd$d$d'd'd$d#d#d$d$d%d%d d$d%d d$d%d"d$d d$d"d$d#d%d$d$d#d"d$d$d$d$d$d$d"d'd"d&d%d$d d$d"d$d&d$d$d$d$d'd$d"d d#d!d%d$d$d"d%d!d#d$d$d$d$d$d%d$d$d$d$d$d$d"d d$d$d"d$d'd$d!d$d!d$d"d$d$d$d$d"d&d$d$d&d d$d$d&d#d d#d$d$d$d$d$d$d#d!d"d$d$d$d$d"d%d$d"d$d&d&d$d d$d d$d#d&d$d$d&d!d$d$d'd$d%d$d$d'd%d!d!d$d%d#d$d#d$d'd#d!d d$d$d d d d$d$d%d%d'd$d$d&d%d$d&d d d#d$d d$d$d$d$d$d$d"d$d$d$d$d$d$d$d$d!d$d$d!d$d$d!d!d$d"d#d$d!d%d#d$d&d#d"d$d$d$d$d'd"d$d!d$d#d$d$d!d#d$d$d$d$d$d#d#d'd!d$d'd'd$d%d"d'd!d'd d!d$d$d#d!d$d#d$d!d#d&d$d$d$d!d$d$d$d$d!d"d#d$d%d d!d$d#d$d$d'd$d"d$d$d$d d(d$d$d d"d%d$d&d$d%d$d$d%d&d$d!d#d$d"d$d$d#d$d$d$d%d!d$d!d$d!d d%d$d$d$d"d$d$d"d$d$d$d$d$d$d$d$d$d$d d d%d$d#d d$d!d!d'd%d$d"d$d$d$d$d#d"d#d$d'd$d#d#d!d$d!d#d$d"d d d"d$d%d"d#d#d d$d%d$d$d!d"d'd%d$d"d$d$d"d$d#d"d!d$d$d"d d$d$d$d"d"d$d'd$d$d d$d$d$d#d$d!d$d$d$d%d#d"d$d d%d"d#d$d"d"d&d$d$d$d$d$d"d$d!d$d$d$d#d$d'd%d$d$d"d$d$d$d$d$d$d$d(d$d%d"d!d$d$d%d$d d&d$d d&d!d$d$d$d#d$d!d%d&d!d$d(d d#d"d"d'd#d#d'd#d"d#d"d%d#d$d$d$d$d%d$d$d$d#d$d$d$d%d"d%d#d!d%d$d$d$d&d$d$d$d$d$d$d$d%d'd"d!d"g\| _d S ))Ng3C>?gq!ɋg۪A@gTy}?g/pjҳ?gm{+?zigaussident_resids.csvr]  rc  rf  i[  g+bs@rg  r   g"X@g|7'z@g&g og{@g镲gୗu8?gt\8w@g4Sd@g%+@g?@gx@gD@g[{@gS:ܳ@g'5@ghG
@gSX@g@gBy`@g}@gϻ0(@g:{@gJ_@gfar@gG #@)rJ  rj  r"   r   r   r   r   r@   rA   rB   rC   rD   rE   rF   rG   rh  r   r   r0   r   r   r   r   r   r   Z	llf_statar   r   r   r2   r   ri  rK  r    r!   r"     s                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    zInvGaussIdentity.__init__rM  r    r    rK  r!   rj    s   rj  c                   @   s   e Zd Zdd ZdS )	Committeec                 C   s  t dddddgdddd	d
gdddddgdddddgdddddgdddddgdd d!d"d#gd$d%d&d'd(gd)d*d+d,d-gd.d/d0d1d2gd3d4d5d6d7gd8d9d:d;d<gd=d>d?d@dAgdBdCdDdEdFgdGdHdIdJdKgdLdMdNdOdPgdQdRdSdTdUgdVdWdXdYdZgd[d\d]d^d_gd`dadbdcddgg| _de| _t dfdgdhdidjdkdlg| _t dmdndodpdqdrdsg| _dt| _du| _dv| _dw| _	dx| _
dy| _dz| _d{| _d|| _d}| _t d~dddddddddddddddddddg| _d S )Ngdق(gU&go1Tgg@,{goͿgOXпg1*kg{YHgֈ`g%VFc?gt:D?g@;~@g]P2@gσv6@g ms>ҿgjVԿg(0g8Kr
g۷ag<l?gȥV0?gJ+@g9ahVv?g.Q5@g~p>uٿg#fyr޿gp=
gg[Ɏ<VgſgevqvRƿgg)yugJggK~ۿgQngQ^ggYGg%?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Ag=	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   r1   r   r   r   r   r   r   r    r    r!   r"     s
                       )                        zCommittee.__init__N)r$   r%   r&   r"   r    r    r    r!   rk    s   rk  c                   @   s   e Zd ZdZdd ZdS )WfszN
    Wfs used for TestGlmPoissonOffset

    Results are from Stata and R.
    c              F   C   s  t j| _d| _dddddddd	d
ddg| _dddddddddddg| _d| _d| _d| _d| _	d| _
d| _d| _d| _d | _d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]d^d_d`dadbdcdddedfgF| _d S )gNgd3n'@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ÓfrX  ;   gv~6F@g"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   Z
wfs_residsr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r    r    r!   r"     s                                                                             zWfs.__init__Nr#   r    r    r    r!   rm    s   rm  c                   @   s   e Zd ZdZd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                 C   s  ddddddddd	d
dddddddg}ddddddddddddddd d!d"g}d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3g}d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDg}|| _ || _|| _|| _dEdFg| _dGdHg| _dI| _dJ| _dK| _dLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\g| _	d S )]NgZ<@gn5I?g}A@g2+gvT*Qgi?gǉ?gQtg8gZX'ggXgg3gs/ggXAW^g:p%"ٿg9tog4@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տgy7gH~@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   r   r   Z
resid_respZ	resid_devZ
resid_pearZ
resid_workr    r    r!   r"   7  s                                     	           

        zCpunishTweediePower15.__init__Nr#   r    r    r    r!   ro  *  s   ro  c                   @   s   e Zd ZdZd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                 C   s  ddddddddd	d
dddddddg}ddddddddddddddd d!d"g}d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3g}ddddddddd	d
dddddddg}|| _ || _|| _|| _d4d5g| _d6d7g| _d8| _d9| _d:| _d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJdKg| _	d S )LNg_<@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1gV4@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   rp  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?rq  rv  r    r    r!   r"     s                                   
           

        zCpunishTweediePower2.__init__Nr#   r    r    r    r!   rw  u  s   rw  c                   @   s   e Zd ZdZd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                 C   s  ddddddddd	d
dddddddg}ddddddddddddddd d!d"g}d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3g}d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDg}|| _ || _|| _|| _dEdFg| _dGdHg| _dI| _dJ| _dK| _dLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\g| _	d S )]Ng4$<@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߿gegWٔ@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   rp  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?
rr  rs  ru  rt  r   r   r   r   r   r   rv  r    r    r!   r"     s                                   
         	

        zCpunishTweedieLog1.__init__Nr#   r    r    r    r!   rx    s   rx  c                   @   s   e Zd ZdZd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                 C   s  ddddddddd	d
ddddddddg}dddddddddddddd d!d"d#d$g}d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3d4d5d6g}d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHg}|| _ || _|| _|| _dIdJdKg| _dLdMdNg| _dO| _dP| _dQ| _dRdSdTdUdVdWdXdXdYdZd[d\d]d\d^d_d`dTg| _	d S )aNgPNЉ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`_;?g|f>RO?gM%3?gbKh?g#?ggy?g1?g|0?gkM?g>"7v?gnZ?gpr?ry  rv  r    r    r!   r"   (  s                                                                 zFairTweedieLog15.__init__Nr#   r    r    r    r!   rz    s   rz  )$r'   r@   Znumpyr   ZpandasrF   Zstatsmodels.apir   Z statsmodels.genmod.tests.resultsr   Z	DECIMAL_4Z	DECIMAL_3Z	DECIMAL_2Z	DECIMAL_1Z	DECIMAL_0r   r(   r-   r.   r3   rQ   rU   r[   rO  rP  rY  rZ  r`  rb  rj  rk  rm  ro  rw  rx  rz  r    r    r    r!   <module>   sn   N   rbR% 
 !        -     ~      E/KLN