
p7]c           @   s  d  Z  d d l Z d d l Z d d l m Z d d l m Z d d l	 m
 Z
 d d l Z d d l m Z m Z d d l Z e j d d	 k  Z d
 Z d	 Z d Z d Z d Z d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ d e f d     YZ  d e f d     YZ! d  e! f d!     YZ" d" e f d#     YZ# d$ e f d%     YZ$ d& e$ f d'     YZ% d( e$ f d)     YZ& d* e f d+     YZ' d, e f d-     YZ( d. e f d/     YZ) d0 e f d1     YZ* d2 e f d3     YZ+ d4 e f d5     YZ, d S(6   s   
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.
iN(   t   asbytes(   t	   NP_LT_113i   (   t   glm_test_resids(   t   add_constantt   categoricali    i   i   i   t   Longleyc           B   s   e  Z d  Z d   Z RS(   sM   
    Longley used for TestGlmGaussian

    Results are from Stata and R.
    c         C   s)  t  j 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  j d d d d d d d g  |  _ t  j d d d d d d d g  |  _ d  |  _ d! |  _ d" |  _ d# |  _	 d$ |  _
 d |  _ d% |  _ d& |  _ d' |  _ d( |  _ t  j d) d* d+ d, d- d. d/ d0 d1 d2 d3 d4 d5 d6 d7 d8 g  |  _ d  S(9   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^)Ai   i	   gZ@gyR@g T7r @g ϖ[@g@g "@gy	2@g !@g@{#@g ĦK@g ɱt@gM@gDt?@g̨@g@@g ф@gO \F@(   t   npt   arrayt   residst   null_deviancet   paramst   bset   aic_Rt	   aic_Statat   deviancet   scalet   llft	   bic_Statat   df_modelt   df_residt   chi2t   fittedvalues(   t   self(    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyt   __init__"   sx    																		(   t   __name__t
   __module__t   __doc__R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR      s   t   GaussianLogc           B   s   e  Z d  Z d   Z RS(   sC   
    Uses generated data.  These results are from R and Stata.
    c      i   C   s  t  j 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/ d0 d/ d/ g d1 d1 d2 d1 d1 g d3 d3 d4 d3 d3 g d5 d5 d6 d5 d5 g d7 d7 d8 d7 d7 g d9 d9 d: d9 d9 g d; d; d< d; d; g d= d= d> d= d= g d? d? d@ d? d? g dA dA dB dA dA g dC dC dD dC dC g dE dE dF dE dE g dG dG dH dG dG g dI dI dJ dI dI g dK dK dL dK dK g dM dM dN dM dM g dO dO dP dO dO g dQ dQ dR dQ dQ g dS dS dT dS dS g dU dU dV dU dU g dW dW dX dW dW g dY dY dZ dY dY g d[ d[ d\ d[ d[ g d] d] d^ d] d] g d_ d_ d` d_ d_ g d d da d d g db db dc db db g dd dd de dd dd g df df dg df df g dh dh di dh dh g dj dj dk dj dj g dl dl dm dl dl g dn dn do dn dn g dp dp dq dp dp g dr dr ds dr dr g dt dt du dt dt g dv dv dw dv dv g dx dx dy dx dx g dz dz d{ dz dz 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 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 gd  |  _ d |  _ t  j d d d g  |  _ t  j d d d g  |  _ d |  _ d |  _ d |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d |  _ t  j d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d ddddddddd	d
ddddddddddddd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?gzB?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'	*?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{i   ia   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   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   k   s   																																																																																																																																																																																																																																															(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   g   s   t   GaussianInversec           B   s   e  Z d  Z d   Z RS(   sG   
    This test uses generated data.  Results are from R and Stata.
    c      i   C   s  t  j 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/ d0 d/ d/ g d1 d1 d2 d1 d1 g d3 d3 d4 d3 d3 g d5 d5 d6 d5 d5 g d7 d7 d8 d7 d7 g d9 d9 d: d9 d9 g d; d; d< d; d; g d= d= d> d= d= g d? d? d@ d? d? g dA dA dB dA dA g dC dC dD dC dC g dE dE dF dE dE g dG dG dH dG dG g dI dI dJ dI dI g dK dK dL dK dK g dM dM dN dM dM g dO dO dP dO dO g dQ dQ dR dQ dQ g dS dS dT dS dS g dU dU dV dU dU g dW dW dX dW dW g dY dY dZ dY dY g d[ d[ d\ d[ d[ g d] d] d^ d] d] g d_ d_ d` d_ d_ g da da db da da g dc dc dd dc dc g de de df de de g dg dg dh dg dg g di di dj di di g dk dk dl dk dk g dm dm dn dm dm g do do dp do do g dq dq dr dq dq g ds ds dt ds ds g du du dv du du g dw dw dx dw dw g dy dy dz dy dy 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 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 gd  |  _ d |  _ t  j d d d g  |  _ t  j d d d g  |  _ d |  _ d |  _ d |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d |  _ t  j d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d ddddddddd	d
ddddddddddddd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+g'	*?gi :w&gŷ$?g_?g?'[g\{QFgzB?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{i   ia   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   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   m  s   																																																																																																																																																																																																															(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   i  s   t   Star98c           B   s   e  Z d  Z d   Z RS(   s0   
    Star98 class used with TestGlmBinomial
    c      0  C   s-  df|  _  dg|  _ d+ |  _ d, |  _ d- |  _ d. |  _ d/ |  _ d0 |  _ d1 |  _ d2 |  _	 d3 |  _
 d4 |  _ d5 |  _ d6 |  _ t j |  _ t j 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g/ |  _ d  S(h  NgP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  i   g$U9@g[r3@g)ԔF@gnU9mgm9mg      ?gKק@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?(   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?(   R
   R   R	   t   df_nullR   R   R   R   R   R   R   t	   llf_StataR   t   pearson_chi2R   t   star98_residsR   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   ^  s         	     													(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   Z  s   t   Lbwc           B   s   e  Z d  Z d   Z RS(   s\   
    The LBW data can be found here

    https://www.stata-press.com/data/r9/rmain.html
    c         C   s  t  j j t  j j t  j j t   d  } t s9 t r t | d   } t	 j
 |  } Wd  QXt	 j d    } | | d  | d <n t j |  j d t  } t | d d d t } | j |  _ t	 j | d | d	 | d
 | d | d | d | d | d f  } t | d t |  _ d |  _ d |  _ d# |  _ d$ |  _ d% |  _ d& |  _ d' |  _ d( |  _ d) |  _ d* |  _  d+ |  _! d, |  _" d- |  _# d. |  _$ t% j& |  _' t	 j( d/ d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d: d; d< d= d= d> d? d@ dA dB dC dD dE dF dG dH dI dJ dK dK dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dY dZ d[ d\ d] d^ d_ d` da db dc dd de df dg 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 du d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d 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(   Ns   stata_lbw_glm.csvt   rbc         S   s   |  j  t d   S(   Ns   "(   t   stripR    (   t   x(    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyt   <lambda>  t    t   racet   indext   colt   dropt   aget   lwtt
   race_blackt
   race_othert   smoket   ptlt   htt   uit   prependg;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@i   gU.Yg$ <>9@g0)Um@g5@"0i   i   i   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_:?(	   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?(*   t   ost   patht   joint   dirnamet   abspatht   __file__R   t   PY2t   openR   t
   recfromcsvt	   vectorizet   pdt   read_csvt
   to_recordst   FalseR   t   Truet   lowt   endogt   column_stackR   t   exogR
   R   R   R   R   R   R   R   R	   R   R   R   R   R    R   t
   lbw_residsR   R   R   (   R   t   filenamet   datafilet   datat   vfunct   design(    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s    $	  	  													(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR"     s   t   Scotvotec           B   s   e  Z d  Z d   Z RS(   s2   
    Scotvot class is used with TestGlmGamma.
    c      !   C   s   d< |  _  d= |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _	 d |  _
 d |  _ d |  _ t j |  _ t j 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; g   |  _ d  S(>   Ng~Ӏ
?g
~`?g`gC)\BR?g[WgZy@g{Q\ľgR1g: >ggioA?g6Oo>gO$F?g[aј>gO}/?gN>gcj?gYw'?i   gە_?i   i   gjQ3Nf@g6oq%@gsTg*ܹqdgϬwH\m?gJ0?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@(   g~Ӏ
?g
~`?g`gC)\BR?g[WgZy@g{Q\ľgR1(   g: >ggioA?g6Oo>gO$F?g[aј>gO}/?gN>gcj?(   R
   R   R	   R   R   R   R   R   R   R   R   R   R    R   t   scotvote_residsR   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s2     	 												(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRN     s   t   Cancerc           B   s   e  Z d  Z d   Z RS(   s`   
    The Cancer data can be found here

    https://www.stata-press.com/data/r10/rmain.html
    c         C   s   t  j j t  j j t  j j t   d  } t j t | d   } | j	 |  _
 t j | j | j f  } t | d d d t } t j | d d d } t | d t |  _ d  S(   Ns   stata_cancer_glm.csvR#   R*   i   R+   t   axisR4   (   R5   R6   R7   R8   R9   R:   R   R=   R<   t	   studytimeRE   RF   R,   t   drugR   RC   t   deleteR   RB   RG   (   R   RI   RK   RM   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   B  s    $	(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRP   <  s   t	   CancerLogc           B   s   e  Z d  Z d   Z RS(   s+   
    CancerLog is used TestGlmGammaLog
    c      5   C   s  t  t |   j   t j 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/ d0 d1 d2 g d3 d4 d5 d6 d7 g d8 d9 d: d; d< g d= d> d? d> d@ g dA dB dC dD dE g dF dG dH dI dJ g dK dL dM dN dO g dP dQ dR dS dT g dU dV dW dX dY g dZ d[ d\ d] d^ g d_ d` da db dc g dd de df dg dh g di dj dk dl dm g dn do dp dq dr g ds dt du dv dw g dx dy dz 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 d d d g g0  |  _ d |  _ t j d d d d g  |  _ t j d d d d g  |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ t j d d dddddddddddd	dd
dddddddd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(#  Ng27{Hg{9yQg_FҞCg`gr"g [gu? ?gx~:gfgE.8g,4k_gHߠ=g6XDgOggyܭgҾ^gB&9;g&RgFg\g"N'g6FPegZ_$	Cg\,JgCg蹅DͿgklexϿgg;#5Ͽg4gM3ȿgBʿg~W[gbȠʿg8Q?gv2=ֿg=yX5ڿgoi5g '*ڿg#ӡgXT$[g ﭿgLgﭿg^a=$߿g;?g2g?g@"?g|g?g(@ ?g
%U ˿gV5RͿg;%5giTNͿgE4~5gh>ӿg ?rֿg7dDg_q^Rkֿgfg3P>タg;r30gR0g)g`^6giR
?g}?gOz12@g,J	z?g @gy)?gE?gص@g˂?g?gؼm1?gT?gʡA@g9}?g9@gWm?gS	hy?g>xEG@gn<՜r?g)yu@gh}?gW\wN?g_LIO@g-4?gL~@g9?q`?ga/?gK+\@gjkD0?gL
F%+@g!Va@?giH?goM`@g{?gm~)@g[Nֿg҆hٿgc('U=gVٿg+Lk	g=#ٿg
Ϳ޿g,+MJADgG6Ws޿g`JgBk}߿g&p{gwT\WgS66GggXg<"ٿged_hd޿g(rVgB:}Q޿gנ/IgB0X{gA$gPkw`kgLϓg =%g]wgHgjgZ_$E)gpUjgS[ g	E-n?gM*?gt	1@gv?gaL{?g5,vοgWnпgߡ(AWg3kпg.g?gaR||Bv?gjWVDO@gDPg?gY(@gdmlvgu䚧ߢg=$P,g{ӖߢgqtgyS4?g_FE?gp:ǲP@gBus?gKu/@gu?gWJ;?gEP\@g1Z웠2?gp @g z*L?g7*?gDR*$@gǞP?gF?g5:L?gH|o@g4y;?gE>+@g2"rg_g	c~gh?RDUg_Ahj3gJsugjg?ߊ{g%T=Og?o*Ra0gZZ Ŀgqx(ƿgSrPg+]&ƿg'󏾉
g!rhmۿgag<,Ԛ}g1kVg׻,gk?gP'?ghC@gPYJ?g+~w?gǮY˿g=(Ϳg$jgeDXCͿgm},g345g2 ֐g&" %g]U"ڿg7-:?g#?g1殩[@gP+?gn@g~nȿg51ʿgH}omghI3̿ʿgq g/ CE?g7t?g+ebg@gnr?q?gzM@gֲx?gy]?g`"p@gDP?gW])@g-qU?g~J
?gH}8q@g}}m?g)*.@gcz?g^kz?gM*k@g.hk(?g?g}i?glV}
r@gPI?gTH9#@g˜;@g_gE.?a?gXln?g:׃@gm:,?gq5?gjO?g)?gx[>t@g_YǉK@g ,0@g}Z?gUeg4Eci   i,   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=@(   t   superRU   R   R   R   R   R	   R
   R   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   Q  s    																																																																																																									(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRU   M  s   t   CancerIdentityc           B   s   e  Z d  Z d   Z RS(   s:   
    CancerIdentity is used with TestGlmGammaIdentity
    c      5   C   s  t  t |   j   t j 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/ d0 d1 d2 g d3 d4 d5 d6 d7 g d8 d9 d: d; d< g d= d> d? d> d@ g dA dB dC dD dE g dF dG dH dI dJ g dK dL dM dN dO g dP dQ dR dS dT g dU dV dW dX dY g dZ d[ d\ d] d^ g d_ d` da db dc g dd de df dg dh g di dj dk dl dm g dn do dp dq dr g ds dt du dv dw g dx dy dz 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 d d d g g0  |  _ t j d d d d g  |  _ t j d d d d g  |  _ d |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ t j d d dddddddddddd	dd
dddddddd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(#  Ng27{Hg{9yQg_FҞCg`gr"g [gu? ?gx~:gfgE.8g,4k_gHߠ=g6XDgOggyܭgҾ^gB&9;g&RgFg\g"N'g6FPegZ_$	Cg\,JgCg蹅DͿgklexϿgg;#5Ͽg4gM3ȿgBʿg~W[gbȠʿg8Q?gv2=ֿg=yX5ڿgoi5g '*ڿg#ӡgXT$[g ﭿgLgﭿg^a=$߿g;?g2g?g@"?g|g?g(@ ?g
%U ˿gV5RͿg;%5giTNͿgE4~5gh>ӿg ?rֿg7dDg_q^Rkֿgfg3P>タg;r30gR0g)g`^6giR
?g}?gOz12@g,J	z?g @gy)?gE?gص@g˂?g?gؼm1?gT?gʡA@g9}?g9@gWm?gS	hy?g>xEG@gn<՜r?g)yu@gh}?gW\wN?g_LIO@g-4?gL~@g9?q`?ga/?gK+\@gjkD0?gL
F%+@g!Va@?giH?goM`@g{?gm~)@g[Nֿg҆hٿgc('U=gVٿg+Lk	g=#ٿg
Ϳ޿g,+MJADgG6Ws޿g`JgBk}߿g&p{gwT\WgS66GggXg<"ٿged_hd޿g(rVgB:}Q޿gנ/IgB0X{gA$gPkw`kgLϓg =%g]wgHgjgZ_$E)gpUjgS[ g	E-n?gM*?gt	1@gv?gaL{?g5,vοgWnпgߡ(AWg3kпg.g?gaR||Bv?gjWVDO@gDPg?gY(@gdmlvgu䚧ߢg=$P,g{ӖߢgqtgyS4?g_FE?gp:ǲP@gBus?gKu/@gu?gWJ;?gEP\@g1Z웠2?gp @g z*L?g7*?gDR*$@gǞP?gF?g5:L?gH|o@g4y;?gE>+@g2"rg_g	c~gh?RDUg_Ahj3gJsugjg?ߊ{g%T=Og?o*Ra0gZZ Ŀgqx(ƿgSrPg+]&ƿg'󏾉
g!rhmۿgag<,Ԛ}g1kVg׻,gk?gP'?ghC@gPYJ?g+~w?gǮY˿g=(Ϳg$jgeDXCͿgm},g345g2 ֐g&" %g]U"ڿg7-:?g#?g1殩[@gP+?gn@g~nȿg51ʿgH}omghI3̿ʿgq g/ CE?g7t?g+ebg@gnr?q?gzM@gֲx?gy]?g`"p@gDP?gW])@g-qU?g~J
?gH}8q@g}}m?g)*.@gcz?g^kz?gM*k@g.hk(?g?g}i?glV}
r@gPI?gTH9#@g.g*nP@g240@g׋{C@g}(?gΑy@g.1@gl|b @gFAiFt@gX<@gW0.@gB?g˜;@geeg͝gci   i,   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ʄ?;@(   RV   RW   R   R   R   R   R
   R   R   R   R   R   R	   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s    																																																																																																									(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRW     s   t   Cpunishc           B   s   e  Z d  Z d   Z RS(   s   
    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 j d d d d d d  d! d" d# d$ d% d& d' d( d) d* d+ g  |  _ d  S(.   Ngٜ.C1?gc?gҼgMg>"= ?gh@g7l[ٸ2gQ,4g2?g#ۖTS?g%x?g]D ?g:Kڞj?g[!"@g_@gNTa@i   gu2@i
   i   gLu贲vS@gQ@gz!{#g`@gie?i   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^?(   gٜ.C1?gc?gҼgMg>"= ?gh@g7l[ٸ2gQ,4(   g2?g#ۖTS?g%x?g]D ?g:Kڞj?g[!"@g_@(   R
   R   R	   R   R   R   R   R   R   R   R   R   R   R    R   t   cpunish_residsR   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   j  s,     	 														(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRX   d  s   t   Cpunish_offsetc           B   s   e  Z d  Z d   Z RS(   sS   
    Same model as Cpunish but with offset of 100.  Many things do not change.
    c         C   s)   t  t |   j   d |  _ d |  _ d  S(   Ng9m4&gٜ.C1?gc?gҼgMg>"= ?gh@g7l[ٸ2gJ+@g1?gjMS?g)Wx?g      ?gBfj?gA`"@(   g9m4&gٜ.C1?gc?gҼgMg>"= ?gh@g7l[ٸ2(   gJ+@g1?gjMS?g)Wx?g      ?gBfj?gA`"@(   RV   RZ   R   R
   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s      	 (   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRZ     s   t   InvGaussc           B   s   e  Z d  Z d   Z RS(   s   
    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   s;  t  j j t  j j t  j j t   d  } t | d  & } t j | d d d t	 d } Wd  QX| d  d  d f |  _
 | d  d  d d   f |  _ t |  j d	 t |  _ d|  _ d|  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ t j |  _ t j 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dddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddddd ddddddddd	d
ddddddddddddd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g |  _ d  S(  Ns   inv_gaussian.csvR#   t	   delimitert   ,t   dtypei   i  i    R4   gb0?gsF9=пg|H?g\s)?g+?g->䱏?gWJR@i  g̢?@i  i   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?(   gb0?gsF9=пg|H?(   g\s)?g+?g->䱏?(    R5   R6   R7   R8   R9   R:   R<   R   t
   genfromtxtt   floatRE   RG   R   RB   R
   R   R	   R   R   R   R   R   R   R   R   R   R    R   t   invgauss_residsR   R   R   (   R   RI   t   fdRK   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s   $	%													(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR[     s   t   Medpar1c           B   s   e  Z d  Z d   Z RS(   s_   
    The medpar1 data can be found here.

    https://www.stata-press.com/data/hh2/medpar1
    c         C   s   t  j j t  j j t  j j t   d  } t j |  j   } | j	 |  _
 t j | j | j f  } t | d d d t } t j | d d d } t | d t |  _ d  S(   Ns   stata_medpar1_glm.csvR*   i    R+   i   RQ   R4   (   R5   R6   R7   R8   R9   R:   R?   R@   RA   t   losRE   R   RF   t   admitypet   codesR   RC   RT   R   RB   RG   (   R   RI   RK   RM   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s    $	(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRc     s   t   InvGaussLogc           B   s   e  Z d  Z d   Z RS(   s5   
    InvGaussLog is used with TestGlmInvgaussLog
    c      ]  C   s,  t  t |   j   t j j t j j t j j t   d  } t	 j
 | d d d d  j |  _ d |  _ t j d d d g  |  _ t j d	 d
 d g  |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ t j d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d  d! d" d" d" d" d" d  d  d" d" d# d  d! d" d  d d$ d" d  d" d$ d" d  d d d" d% d" d" d d" d d d" d" d" d d" d d# d" d d" d& d% d" d" d# d% d" d" d" d" d" d" d$ d" d" d# d d! d" d" d d" d d" d d$ d  d! d" d! d" d" d d! d" d" d" d" d" d d" d  d  d" d" d d" d d" d! d" d" d  d" d% d" d" d% d# d" d" d  d d% d" d! d" d" d" d% d" d" d" d" d! d" d# d" d! d! d" d d# d# d d# d" d" d" d" d" d" d" d  d d" d d" d" d" d$ d" d" d" d" d d d" d$ d" d d" d" d! d% d" d d" d$ d$ d  d# d" d" d  d" d" d# d" d  d d" d$ d  d" d" d! d  d" d" d" d" d$ d" d d# d  d" d% d$ d" d! d" d d  d d" d d# d" d% d" d" d d  d d d" d" d$ d% d" d" d" d" d  d  d" d" d d" d! d  d" d! d! d" d" d d d d d! d" d# d" d d" d" d" d! d# d" d d d d" d" d  d d% d d# d" d! d d" d" d" d# d d$ d" d" d" d d  d" d" d" d" d" d" d" d# d  d d" d d# d" d" d d" d  d" d$ d" d  d# d# d" d  d d" d d  d" d# d! d" d  d# d d  d$ d" d d d! d" d" d" d" d" d# d" d! d" d" d  d" d" d" d" d  d  d" d  d" d" d" d" d  d! d" d  d" d  d" d" d" d" d! d  d d$ d" d" d% d" d! d" d" d" d" d# d# d# d" d" d" d" d d% d! d" d# d" d" d" d" d" d# d" d" d" d" d# d" d# d d" d d" d" d" d" d" d" d" d" d d  d" d" d" d" d" d" d" d  d d  d" d d" d# d" d" d" d# d% d" d" d# d" d" d" d" d" d" d" d" d" d" d" d" d" d d" d" d" d" d d  d# d" d d d# d" d" d! d$ d" d" d# d" d d# d" d" d d" d! d  d  d% d d" d! d" d d# d d d" d" d  d# d" d" d" d" d" d" d" d d! d# d d! d d# d" d" d" d" d% d" d" d" d" d" d# d d  d d% d" d# d" d d" d" d# d" d# d" d" d" d" d d# d d  d$ d" d! d" d d" d  d" d" d" d d$ d" d" d" d# d d" d  d d  d" d" d  d" d" d" d$ d  d" d" d d d  d" d" d! d% d d" d" d" d" d" d  d" d" d# d d" d d" d  d$ d" d  d" d" d" d  d" d" d" d" d" d" d! d" d% d" d  d d  d" d d" d# d! d" d" d" d" d  d" d# d d# d" d d" d" d" d" d d" d d" d" d" d" d" d" d" d" d# d& d" d" d" d" d d" d d" d" d# d" d" d" d" d" d" d% d% d" d" d" d# d" d  d! d  d! d! d d" d d" d d" d" d" d" d" d% d" d# d" d d" d# d" d" d" d" d" d  d" d! d  d! d" d# d" d d  d" d" d! d# d d" d" d" d! d" d" d  d d! d% d" d d% d" d" d% d% d" d! d! d" d" d# d# d d" d# d d" d# d  d" d d" d  d" d! d# d" d" d! d  d" d" d" d" d" d" d  d% d  d$ d# d" d d" d  d" d$ d" d" d" d" d% d" d  d d! d d# d" d" d  d# d d! d" d" d" d" d" d# d" d" d" d" d" d" d  d d" d" d  d" d% d" d d" d d" d  d" d" d" d" d  d$ d" d" d$ d d" d" d$ d! d d! d" d" d" d" d" d" d! d d  d" d" d" d" d  d# d" d  d" d$ d$ d" d d" d d" d! d$ d" d" d$ d d" d" d% d" d# d" d" d% d# d d d" d# d! d" d! d" d% d! d d d" d" d d d d" d" d# d# d% d" d" d$ d# d" d$ d d d! d" d d" d" d" d" d" d" d  d" d" d" d" d" d" d" d" d d" d" d d" d" d d d" d  d! d" d d# d! d" d$ d! d  d" d" d" d" d% d  d" d d" d! d" d" d d! d" d" d" d" d" d! d! d% d d" d% d% d" d# d  d% d d% d d d" d" d! d d" d! d" d d! d$ d" d" d" d d" d" d" d" d d  d! d" d# d d d" d! d" d" d% d" d  d" d" d" d d& d" d" d d  d# d" d$ d" d# d" d" d# d$ d" d d! d" d  d" d" d! d" d" d" d# d d" d d" d d d# d" d" d" d  d" d" d  d" d" d" d" d" d" d" d" d" d" d d d# d" d! d d" d d d% d# d" d  d" d" d" d" d! d  d! d" d% d" d! d! d d" d d! d" d  d d d  d" d# d  d! d! d d" d# d" d" d d  d% d# d" d  d" d" d  d" d! d  d d" d" d  d d" d" d" d  d  d" d% d" d" d d" d" d" d! d" d d" d" d" d# d! d  d" d d# d  d! d" d  d  d$ d" d" d" d" d" d  d" d d" d" d" d! d" d% d# d" d" d  d" d" d" d" d" d" d" d& d" d# d  d d" d" d# d" d d$ d" d d$ d d" d" d" d! d" d d# d$ d d" d& d d! d  d  d% d! d! d% d! d  d! d  d# d! d" d" d" d" d# d" d" d" d! d" 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('   Ns   medparlogresids.csvt   sepR]   t   headergrlvt@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@i   iY  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@(   RV   Rg   R   R5   R6   R7   R8   R9   R:   R?   R@   t   Nonet   valuesR   R	   R   R   R
   R   R   R   R   R   R   R   R   R   R   R   (   R   RI   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s   $	!										(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRg     s   t   InvGaussIdentityc           B   s   e  Z d  Z d   Z RS(   sN   
    Accuracy is different for R vs Stata ML vs Stata IRLS, we are close.
    c      ]  C   s6,  t  t |   j   t j d d d g  |  _ t j d d d g  |  _ t j j	 t j j
 t j j t   d  } t j | d d	 d
 d  j |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ d |  _ t j d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d  d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d  d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d! d" d# d$ d% d% d% d% d% d# d# d% d% d& d# d$ d% d# d! d' d% d# d% d' d% d# d! d" d% d( d% d% d" d% d! d! d% d% d% d" d% d" d& d% d! d% d) d( d% d% d& d( d% d% d% d% d% d% d' d% d% d& d! d$ d% d% d" d% d! d% d" d' d# d$ d% d$ d% d% d" d$ d% d% d% d% d% d" d% d# d# d% d% d! d% d! d% d$ d% d% d# d% d( d% d% d( d& d% d% d# d! d( d% d$ d% d% d% d( d% d% d% d% d$ d% d& d% d$ d$ d% d" d& d& d" d& d% d% d% d% d% d% d% d# d! d% d! d% d% d% d' d% d% d% d% d" d" d% d' d% d" d% d% d$ d( d% d" d% d' d' d# d& d% d% d# d% d% d& d% d# d" d% d' d# d% d% d$ d# d% d% d% d% d' d% d! d& d# d% d( d' d% d$ d% d! d# d" d% d! d& d% d( d% d% d! d# d" d" d% d% d' d( d% d% d% d% d# d# d% d% d" d% d$ d# d% d$ d$ d% d% d! d" d! d" d$ d% d& d% d" d% d% d% d$ d& d% d" d" d! d% d% d# d! d( d! d& d% d$ d" d% d% d% d& d" d' d% d% d% d" d# d% d% d% d% d% d% d% d& d# d" d% d" d& d% d% d" d% d# d% d' d% d# d& d& d% d# d! d% d" d# d% d& d$ d% d# d& d! d# d' d% d! d" d$ d% d% d% d% d% d& d% d$ d% d% d# d% d% d% d% d# d# d% d# d% d% d% d% d# d$ d% d# d% d# d% d% d% d% d$ d# d! d' d% d% d( d% d$ d% d% d% d% d& d& d& d% d% d% d% d" d( d$ d% d& d% d% d% d% d% d& d% d% d% d% d& d% d& d! d% d" d% d% d% d% d% d% d% d% d" d# d% d% d% d% d% d% d% d# d! d# d% d! d% d& d% d% d% d& d( d% d% d& d% d% d% d% d% d% d% d% d% d% d% d% d% d! d% d% d% d% d! d# d& d% d" d" d& d% d% d$ d' d% d% d& d% d! d& d% d% d" d% d$ d# d# d( d! d% d$ d% d" d& d" d" d% d% d# d& d% d% d% d% d% d% d% d! d$ d& d" d$ d! d& d% d% d% d% d( d% d% d% d% d% d& d" d# d! d( d% d& d% d! d% d% d& d% d& d% d% d% d% d" d& d" d# d' d% d$ d% d" d% d# d% d% d% d! d' d% d% d% d& d" d% d# d" d# d% d% d# d% d% d% d' d# d% d% d" d" d# d% d% d$ d( d! d% d% d% d% d% d# d% d% d& d" d% d! d% d# d' d% d# d% d% d% d# d% d% d% d% d% d% d$ d% d( d% d# d" d# d% d" d% d& d$ d% d% d% d% d# d% d& d" d& d% d! d% d% d% d% d! d% d" d% d% d% d% d% d% d% d% d& d) d% d% d% d% d" d% d! d% d% d& d% d% d% d% d% d% d( d( d% d% d% d& d% d# d$ d# d$ d$ d" d% d" d% d! d% d% d% d% d% d( d% d& d% d" d% d& d% d% d% d% d% d# d% d$ d# d$ d% d& d% d! d# d% d% d$ d& d! d% d% d% d$ d% d% d# d" d$ d( d% d! d( d% d% d( d( d% d$ d$ d% d% d& d& d! d% d& d! d% d& d# d% d! d% d# d% d$ d& d% d% d$ d# d% d% d% d% d% d% d# d( d# d' d& d% d! d% d# d% d' d% d% d% d% d( d% d# d! d$ d" d& d% d% d# d& d" d$ d% d% d% d% d% d& d% d% d% d% d% d% d# d! d% d% d# d% d( d% d" d% d" d% d# d% d% d% d% d# d' d% d% d' d! d% d% d' d$ d! d$ d% d% d% d% d% d% d$ d" d# d% d% d% d% d# d& d% d# d% d' d' d% d! d% d! d% d$ d' d% d% d' d" d% d% d( d% d& d% d% d( d& d" d" d% d& d$ d% d$ d% d( d$ d" d! d% d% d! d! d! d% d% d& d& d( d% d% d' d& d% d' d! d! d$ d% d! d% d% d% d% d% d% d# d% d% d% d% d% d% d% d% d" d% d% d" d% d% d" d" d% d# d$ d% d" d& d$ d% d' d$ d# d% d% d% d% d( d# d% d" d% d$ d% d% d" d$ d% d% d% d% d% d$ d$ d( d" d% d( d( d% d& d# d( d" d( d! d" d% d% d$ d" d% d$ d% d" d$ d' d% d% d% d" d% d% d% d% d" d# d$ d% d& d! d" d% d$ d% d% d( d% d# d% d% d% d! d) d% d% d! d# d& d% d' d% d& d% d% d& d' d% d" d$ d% d# d% d% d$ d% d% d% d& d" d% d" d% d" d! d& d% d% d% d# d% d% d# d% d% d% d% d% d% d% d% d% d% d! d! d& d% d$ d! d% d" d" d( d& d% d# d% d% d% d% d$ d# d$ d% d( d% d$ d$ d" d% d" d$ d% d# d! d! d# d% d& d# d$ d$ d! d% d& d% d% d" d# d( d& d% d# d% d% d# d% d$ d# d" d% d% d# d! d% d% d% d# d# d% d( d% d% d! d% d% d% d$ d% d" d% d% d% d& d$ d# d% d! d& d# d$ d% d# d# d' d% d% d% d% d% d# d% d" d% d% d% d$ d% d( d& d% d% d# d% d% d% d% d% d% d% d) d% d& d# d" d% d% d& d% d! d' d% d! d' d" d% d% d% d$ d% d" d& d' d" d% d) d! d$ d# d# d( d$ d$ d( d$ d# d$ d# d& d$ d% d% d% d% d& d% d% d% d$ d% 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{+?s   igaussident_resids.csvRh   R]   Ri   grlvt@i[  g+bs@iY  i   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 #@(    RV   Rl   R   R   R   R
   R   R5   R6   R7   R8   R9   R:   R?   R@   Rj   Rk   R   R	   R   R   R   R   R   R   R   t	   llf_stataR   R   R   R    R   (   R   RI   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s   $	!													(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRl     s   t	   Committeec           B   s   e  Z d    Z RS(   c         C   s}  t  j 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/ d0 d1 d2 g d3 d4 d5 d6 d7 g d8 d9 d: d; d< g d= d> d? d@ dA g dB dC dD dE dF g dG dH dI dJ dK g dL dM dN dO dP g dQ dR dS dT dU g dV dW dX dY dZ g d[ d\ d] d^ d_ g d` da db dc dd g g  |  _ de |  _ t  j df dg dh di dj dk dl g  |  _ t  j dm dn do dp dq dr ds g  |  _ dt |  _ du |  _ dv |  _ dw |  _	 dx |  _
 dy |  _ dz |  _ d{ |  _ d| |  _ d} |  _ t  j d~ d d d d d d d d d d d d d d d d d d d g  |  _ 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@i   i   g=)@gD7.>@g9{5@gv@g\Y@g׀JM~k@gFt@g9lS#Ѵ[@g}xTQ\@gG{k@g[0lL@g&KeN@g\@g"ǵ%R@gM'N@gHT"5@g|3aE@gJK@gnf5?gfх%2@(   R   R   R   R	   R
   R   R   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s    																																																			(   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRn     s   t   Wfsc           B   s   e  Z d  Z d   Z RS(   sN   
    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/ 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 gF |  _ d  S(g   Ngd3n'@gPG?gݼ?gيG?g6G?g\BZ?gvǼ?g^7i?g_8c/?gUgj}ӿgXg#8?gL"%?gq76?gc058?gD^?gd{Z?g^nT?g163?g;d?gNǭD?gf8TBf?g##Q@g|	91@gKQ@g      ?g_bEFBogÓfi
   i;   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   t
   wfs_residsR   R	   R
   R   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s<    										(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRo     s   t   CpunishTweediePower15c           B   s   e  Z d  Z d   Z RS(   s  
    # 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/ d0 d1 d2 d3 g } d4 d5 d6 d7 d8 d9 d: d; d< d= d> d? d@ dA dB dC dD g } | |  _  | |  _ | |  _ | |  _ dE dF g |  _ dG dH g |  _ dI |  _ dJ |  _ dK |  _ dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ d[ 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@i   i   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?(
   t   resid_responset   resid_deviancet   resid_pearsont   resid_workingR
   R   R   R   R   R   (   R   t
   resid_respt	   resid_devt
   resid_peart
   resid_work(    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   A  sT    																								(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRq   4  s   t   CpunishTweediePower2c           B   s   e  Z d  Z d   Z RS(   s  
    # 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/ d0 d1 d2 d3 g } d d d d d d d d d	 d
 d d d d d d d g } | |  _  | |  _ | |  _ | |  _ d4 d5 g |  _ d6 d7 g |  _ d8 |  _ d9 |  _ d: |  _ d; d< d= d> d? d@ dA dB dC dD dE dF dG dH dI dJ dK g |  _	 d  S(L   Ng_<@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/@i   i   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?(
   Rr   Rs   Rt   Ru   R
   R   R   R   R   R   (   R   Rv   Rw   Rx   Ry   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     sX    																							(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyRz     s   t   CpunishTweedieLog1c           B   s   e  Z d  Z d   Z RS(   s  
    # 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/ d0 d1 d2 d3 g } d4 d5 d6 d7 d8 d9 d: d; d< d= d> d? d@ dA dB dC dD g } | |  _  | |  _ | |  _ | |  _ dE dF g |  _ dG dH g |  _ dI |  _ dJ |  _ dK |  _ dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ d[ 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@i   i   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   (   R   Rv   Rw   Rx   Ry   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR     s\    																					(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR{     s   t   FairTweedieLog15c           B   s   e  Z d  Z d   Z RS(   sH  
    # 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/ d0 d1 d2 d3 d4 d5 d6 g } d7 d8 d9 d: d; d< d= d> d? d@ dA dB dC dD dE dF dG dH g } | |  _  | |  _ | |  _ | |  _ dI dJ dK g |  _ dL dM dN g |  _ dO |  _ dP |  _ dQ |  _ dR dS dT dU dV dW dX dX dY dZ d[ d\ d] d\ d^ d_ d` dT g |  _	 d  S(a   NgPNЉ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#~?igI#!ؿg4"N^(?g`ڸgT^gV?g{x"z?gx\L?gGztA@i   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?(
   Rr   Rs   Ru   Rt   R
   R   R   R   R   R   (   R   Rv   Rw   Rx   Ry   (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR   2  sT    																											(   R   R   R   R   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyR|     s   (-   R   t   numpyR   t   pandasR?   t   statsmodels.compat.pythonR    t   statsmodels.compat.numpyR   R'   R   R5   t   statsmodels.apiR   R   t   syst   version_infoR;   t	   DECIMAL_4t	   DECIMAL_3t	   DECIMAL_2t	   DECIMAL_1t	   DECIMAL_0t   objectR   R   R   R   R"   RN   RP   RU   RW   RX   RZ   R[   Rc   Rg   Rl   Rn   Ro   Rq   Rz   R{   R|   (    (    (    sK   lib/python2.7/site-packages/statsmodels/genmod/tests/results/results_glm.pyt   <module>   sX   K `_#!    )    E/KLN