
p7]c           @   s   d  Z  d d l Z d d l Z d e f d     YZ e j j e j j e	   Z
 d e f d     YZ d e f d     YZ d	 e f d
     YZ d e f d     YZ d S(   s-   
Test Results for discrete models from Stata
iNt	   Namespacec           B   s   e  Z RS(    (   t   __name__t
   __module__(    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR       s   t   Anesc           B   s#   e  Z d    Z d   Z e   Z RS(   c         C   s   d |  _  d S(   sI   r
        Results are from Stata 11 (checked vs R nnet package).
        i  N(   t   nobs(   t   self(    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   __init__   s    c       $   C   sy	  t    }  d |  _ d d d d d d d d	 d
 d d d d d d d d d d d d d d d d d d d d d d  d! d" d# d$ d% g$ } t j | 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 g$ } t j | d"d( d) |  _ t j t j j	 t
 dN   |  _ t j t j j	 t
 dO   |  _ d  |  _ dP |  _ dQ |  _ dR |  _ dS |  _ dT |  _ dU |  _ d#|  _ dW |  _ d& |  _ 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} g$ } 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 g$ } t j | d%d( d) |  _ d d g d d g d d g d d g d d g d d g g d d g d d g d d g d d g d d g d d g g d d g d d g d d g d d g d d g d d g g d d g d d g d d g d d g d d g d d g g d d g d d g d d g d d g d d g d d g g d d g d d g d d g d d g d d g d d g g g } 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 g  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 g  j |  _! t j d0d1d2d3d4g d5d6d7d8d9g d: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 g  j |  _" t j dSdTdUdVdWg dXdYdZ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 g  j |  _# t j dvdwdxdydzd{g d|d}d~dddg ddddddg ddddddg ddddddg ddddddg ddddddg g  j |  _$ t j ddddddg ddddddg ddddddg ddddddg ddddddg ddddddg ddddddg g  j |  _% t j ddddddg ddddddg ddddddg ddddddg ddddddg ddddddg ddddddg g  j |  _& t j ddddddg ddddddg d dddddg dddd	d
dg ddddddg ddddddg ddddddg g  j |  _' t j t j j	 t
 d dd |  _( |  S(&  Ni  gT-g13}?gA}˜.gz(?g'DHu?g:׿g?\gW?gOrgqTh,?g4?gg6 gzbJY?g~XOjgHK}g }z?g*qSg:Bpg@ku?gmXǁg`l?g)<?gDItgn8᷿g?|O'?gvqUg??gAI?g{L=g`¿g' @g(lQgh{O@n?gPS{?g	$6(i   it   ordert   FgbOu?g]?g@EWϹz?gTQֲ?g^)&?g'?'?g"I?g_P釵?gu5?gWaյ?g=Ж?g( Tl?gNb4?g$NK?g,Z 4?gvgP*?gB.5?guɦ1?ga
ڭk?gm?g.HO=?gd?g_,8Ӛ?gԙjӀ?g*q%?g/U?g.GU,?gõ?gÚL
?g)J?g;?g9[?gmk|?g̰2R?g;y?g/?s   yhat_mnlogit.csvs   phat_mnlogit.csvgޫזgcYg#@g% {S*g3?i   i$   i   gg@g2kĨ@gH2տg}/?:p	@gD !gR?gjuF?gƫgR.R!g}@gu@3%g"=D @g p0@g!DgCVHgcRa@g`[g:FS+ga#g?gAq [	g? g_v#@g<
g @g e	@grzggVr>
&@g[gtj@gpR0@guI giv
gX=7,@gC4g٩ME@g[B7@gR)&g;8uq]?g0a6u%X?g-xF!?g=?gKQ?gCv?go1O?gj3?gՉ(>o?g2MĚK?gi25?gʠmDj?gG͇-?gXp3?g6FOQ?gR?gE54?g{0Y?gP(?g LhD;gҷZ[?gq!H?gT?gY<gcwA0?gn]aɛ9g?gghH?g;?gz2ʑ<gG=u K?g~+?F6g39؂?ge-o:?geD>g9gE'g:-~?gj<?gbL?g_(Qgxxρ刿gtΛgnr?gg<)[?g8t̹g@R?ȍ?gL/ſgUgU?g5=nR?g@m`vgF@~gtQol?g,I?g
Xl8q?g&jo?gqQg%)g1˿g[6w?gE6F?g`4[K?ghge6#~?g̞LпggB?g3}ƅgԤo-?gu ]gG9MagVJ*ƿgEBwwgfbk?gj{?gQAəgSEu?ghm?gFB?go?gGc?gs|8"g\gjſgn~hg'
b?g" 89?g@7o͠ge:
whgy"*ĩ?gj3?gfI%d?gv ?g(5
In!gKrJgΒ̿gUCeحgg?gyZ~*@gOJgz?g=@,Z?g'l??g|C]?g<*I?g=],gCVz$gj1LEǁ?gt;S	gSf?gIfĝlOxg\a|?g$}  g YIag[ugCI,vg"I\pg윿giZmyPgRa86{?g;P?g!t[]gGH&wga(dAGge>dvg`žTA?g;Χ%Pge_ U?grhp??gQpG\?gXOa?g]$Yg905:?g(tFg09Zs?gacr%d?gR'Jɍb~g!w`X?g+D?gct_?gXv?go?g;,({?gBV=z)G?g5hxۀ?g0`?g sGo?g1l'%~?g4gH?gP^H?gfo`?gYi?g8=KEw?gꢡE?gd<|?gG^?gP_?g̿"j?g6bU9?g|Ѫq?g&+n0R?g vh?gD
J&v?g;MJC?gcz?g5^?gun-Um?g?Sw|?gَܬG?gXZ?g&Ob?gܩ(k?gr'?gjN3F?g,  ~?g2sa?g:ȸ?g{NϤg:$>cf?gW#A
gVn{g^Ն?g]y@gu'8bgdƂgeL|^gMDogȆmGgZ~g!SgE=>y?g~y,?gcg(!LgHSgg㮛C?g	#lg]}Jj?g|~r:H?g.t~?gΠJ,r?gKAwgH{?gU[BguH?g<UYIz?g=hFW}gT!(?g({7@?g^GW?g{ΞOw?g2Teq?gi+?g\3J?g)[?go+uc?g*ȉs?gSRFM?g-O?gy?gBGe?g36p?g|Z8?gr.!L?gEoj?g c?gLe?g%裓}?gFxA?g>"|w?gdh42Z?gbn?gؖD?gHv:H?gLÒ?g>c?gJ+<r?gߢ?g];M?g!&}?g7f?ga`jf?gL3?g"OsA?gÊy?ge\?g9)@~v?g	bg`U}h?gC}٣Lg4E6kgg?пg?dQl?g{lvgv]
_gpi66?gadig1jȬʿg<i?lzgo;tp?gvKgo[ц?grbf?ge2pĿg`	agng*r-Xgyڻrg
m:fE?gOgk]ŝC?g<v+<ngf7?g	96䚹(gz2x J?g;$.?gt0*Wa?gl$	4g-\Ng~cP>gȔ)8gQ]d?gyYcig9.VA?go܋3?gt?gOY;e?gs=Ks?gRphn?g޻E?gSWdF?gR+IT?g"qM_?gds?g-ARNo?g&?g{{oHH?g@+G?g[V^`?gj9C?g+Kni?gu?g碌+?PE?g>{?g*a]?g,<?gKt&/`?g&.t?g6P9?gwJ^q?gJ. TS?gr;T(?ge3~Wh?gz?g"@B?ga<Ȅz?gaW_?ga?gi(_`7l?gP?g#(+3F?gcDG ?g5&~sb?g}ã?g\"ʟ7j?gGP8?g3?D?gb+)
}?gk)4a?g0ݲ/?g(,+?g
g'?g9Lg#~gɢ	gC)/?g2}ƅؿgoSg20g"ӠgMgM!_[g&$gѿg!!ygΈo?g<?gd g"Ӫg`6!̿g+d0?ggjRDeGӆ?g~3PѿgiGhgZ7f?gl\bيW?gT2?gt6?gLH?g_aj?g(?gg/g^c/?gaG#?gsiq?gu$Ѫg
ox?g6f
 Tg8Ż?gGsϤ?ge0jK@g/ܛ?gi!?g+xyu?gOA6<?g20o?g֌?g@Wۚ?g&?g*2u?gr?gu/?g5L?gȮwrt?g7зU?gpo>:ܣ{?g?gfԔ?gY?g׬-?g$0?g2xB?gr-?g=?gdkv?g< Pj?g %%Hc?gz?gBJ?gNΕ?gd
?gD\P?g3QC]?gw^=w?gKOa,?gUB'?gÍ?g-z`-bK?g4Rʻ?gh-y?gKj?g9J?giؖ?s   mnlogit_resid.csvt	   delimitert   ,(   i   i(   i   ii  (   i   i(   i   i()   R    R   t   npt   reshapet   paramst   bset   loadtxtt   ost   patht   joint   cur_dirt   yhatt   phatt   Nonet
   cov_paramst   llft   llnullt   llrt
   llr_pvaluet	   prsquaredt   df_modelt   df_residt   Jt   Kt   aict   bict   zt   pvaluest   asarrayt   conf_intt   arrayt   Tt   margeff_dydx_overallt   margeff_dydx_overall_set   margeff_dydx_meant   margeff_dydx_mean_set   margeff_dydx_dummy_overallt   margeff_dydx_dummy_overall_set   margeff_eydx_dummy_overallt   margeff_eydx_dummy_overall_set   resid(   t   objR   R   R#   R$   R&   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   mnlogit_basezero   s   			!!																																																							!													!							(   R   R   R   R3   (    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR      s   		t
   DiscreteL1c           B   sY   e  Z d    Z d   Z e   Z d   Z e   Z d   Z e   Z d   Z e   Z RS(   c         C   s   d S(   s{   
        Special results for L1 models
        Uses the Spector data and a script to generate the baseline results
        N(    (   R   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR      s    c          C   s   t    }  t j } d d d d g |  _ d d g | | g d d g | | g g |  _ d | d	 | g |  _ d
 |  _ d |  _ d |  _ d | d | g | | | | g d | d | g | | | | g g |  _	 |  S(   s  
        Results generated with:
            data = sm.datasets.spector.load(as_pandas=False)
            data.exog = sm.add_constant(data.exog, prepend=True)
            alpha = 3 * np.array([0, 1, 1, 1])
            res2 = sm.Logit(data.endog, data.exog).fit_regularized(
                method="l1", alpha=alpha, disp=0, trim_mode='size',
                size_trim_tol=1e-5, acc=1e-10, maxiter=1000)
        g<^.ig        gn ?goRM"g!-J?gQ:gh'Nn?g}#@gEVǼ?i   gIE@g,F@gŅfK@g+/Sҿgu9?(
   R    R   t   nanR   R&   R   t
   nnz_paramsR!   R"   R   (   R2   R5   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   logit   s     
								c          C   sC   t    }  d d d d g d d d d g d	 d
 d d g g |  _ |  S(   si  
        Results generated with
            params = np.zeros((3, 4))
            alphas = np.array(
                    [[0.1, 0.1, 0.1, 0.1],
                        [0.4, 0.4, 0.5, 0.5], [0.5, 0.5, 1, 1]])
            model = sm.Logit(data.endog, data.exog)
            for i in range(3):
                alpha = alphas[i, :]
                res2 = model.fit_regularized(method="l1", alpha=alpha,
                                             disp=0, acc=1e-10,
                                             maxiter=1000, trim_mode='off')
                params[i, :] = res2.params
            print(params)
        gz$gbRC*@g+؈X?g8u @g&Ng+ ?g<ⴖ?gv!:?g?BigWY_M3?g       gm?(   R    R   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   sweep  s
    	c          C   s   t    }  t j } d d d d g |  _ d d g d d g d	 d
 g t j t j g g |  _ d d d t j g |  _ d |  _ d |  _ d |  _ d d d | g d d d | g d d d | g | | | | g g |  _	 |  S(   s  
        Results generated with
            data = sm.datasets.spector.load(as_pandas=False)
            data.exog = sm.add_constant(data.exog, prepend=True)
            alpha = np.array([0.1, 0.2, 0.3, 10])
            res2 = sm.Probit(data.endog, data.exog).fit_regularized(
                method="l1", alpha=alpha, disp=0, trim_mode='auto',
                auto_trim_tol=0.02, acc=1e-10, maxiter=1000)
        g7|gы ?g"'|J?g        gdڭ"g=V;qgݞ?g57@g*%縿g(?gn0Ky @gX?gqm?i   gX+3C@gZ~fE@g/H~.@g]E̲gc?zgy?g igSGUv?(
   R    R   R5   R   R&   R   R6   R!   R"   R   (   R2   R5   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   probit8  s     
								c          C   s  t    }  d d d d d d g d d d	 d
 d d g d d d d d d g d d d d d d g d d d d d d g d d d d d  d! g g |  _ d" d# g t j t j g d$ d% g t j t j g t j t j g d& d' g g d( d) g d* d+ g d, d- g d. d/ g d0 d1 g d2 d3 g g d4 d5 g d6 d7 g d8 d9 g t j t j g d: d; g d< d= g g d> d? g d@ dA g dB dC g dD dE g dF dG g dH dI g g dJ dK g dL dM g dN dO g dP dQ g dR dS g dT dU g g dV dW g dX dY g dZ d[ g d\ d] g d^ d_ g d` da g g g |  _ db dc dd de df dg g t j dh di dj dk dl g dm dn do dp dq dr g t j ds t j dt du dv g t j dw dx dy dz d{ g d| d} d~ d d d g g |  _ d |  _ d |  _ d |  _ |  S(   s  
        Results generated with
            anes_data = sm.datasets.anes96.load(as_pandas=False)
            anes_exog = anes_data.exog
            anes_exog = sm.add_constant(anes_exog, prepend=False)
            mlogit_mod = sm.MNLogit(anes_data.endog, anes_exog)

            alpha = 10 * np.ones((mlogit_mod.J - 1, mlogit_mod.K))
            alpha[-1, :] = 0
            mlogit_l1_res = mlogit_mod.fit_regularized(
            method='l1', alpha=alpha, trim_mode='auto', auto_trim_tol=0.02,
            acc=1e-10)
        gq$iP?gRQg,zgQ'kg.u߽g/`g        gB6"P?g-nؾ?gV?gq?g>?g˅ʿWgq[g:4g
ksg&c𬌿gϗSwg	c?g,2?g.?gb$l?g|6\Ʀ?ggP?g h?16?g5۪^]?gSLƸ?gƒ\6?gߧgѸ˭ g:gxgcYYpy"gIgjGi_?g&8pgrvgY.?gأ?g2 g ]Q?g$+¿gcpb?g8롿gHm}sg>33g!W?g~Wrs?g)5A?gZIg"B?g2ƿgeM?g cÿg[?g`.gIZ>?g\|g+Ư?g1.M:g<}Fпg	g҃ig?g@	?g B^?grOgI2?gוg)S^?g̀?g^/?gjS+gVL<gE<0쿿gAD?gr?gA'?gZegqٮg&އgڷ?g_eH?g ݬ?gtY,gsFgWi ſg)Qa2v>gep?gw0?gh;gw"!?gǖNeg/?g҉Z[?gu?g^,X5&g٧4gF剔$?g*B?gf1-?guh'V?gݣ@D?g:S?gꡦ+?gp?gʀ?g!]	0n?gf>R
˽?g`Δy?g{Wo?gNY?g
I?gE*k}?gsw~?gH٢7?g?gjOA?gXs~?gCλ?gJ3?g*AD?g%b?gX(K?gj>?g,Eч?g&{[?gJg?gk'F?g(P4?i    g	p@gNNIͨ@(	   R    R   R   R5   R&   R   R6   R!   R"   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   mnlogitW  s    																														(   R   R   R   R7   R8   R9   R:   (    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR4      s   								Zt   Spectorc           B   s;   e  Z d  Z d   Z d   Z e   Z d   Z e   Z RS(   s#   
    Results are from Stata 11
    c         C   s   d |  _  d  S(   Ni    (   R   (   R   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR     s    c       !   C   s  t    }  d |  _ d d d d g |  _ d d d d	 g d d
 d d g d d d d g d	 d d d g g |  _ 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+ 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 dE dF dG dH dI dJ dK dL dM dN dO dP dQ dR dS g  |  _ 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 g  |  _ t j	 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 g   |  _
 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 g   |  _ d |  _ d |  _ d |  _ dP|  _ d |  _ d |  _ d |  _ d |  _ d |  _ d d d d g |  _ d d g d d g d d g d d g g |  _ d d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _ d d d g |  _  d d d g |  _! d d d g |  _" d d d g |  _# d d d g |  _$ d d d g |  _% d d d g |  _& d d d g |  _' d d d g |  _( d d d g |  _) d d d g |  _* d d dg |  _+ dddg |  _, dddg |  _- dd	d
g |  _. dddg |  _/ dddg |  _0 dddg |  _1 dddg |  _2 dddg |  _3 ddd g |  _4 ddd g |  _5 d d d g |  _6 d d d g |  _7 ddd g |  _8 d!d"d#g |  _9 d$d%d&g |  _: d'd(d)g |  _; d d d*g |  _< d d d+g |  _= d d d,g |  _> d d d-g |  _? d.d/d0g |  _@ d1d2d3g |  _A d d d4g |  _B d d d5g |  _C d d d6g |  _D d d d7g |  _E d8d9d:g |  _F d;d<d g |  _G d=d>d?g |  _H d d@d g |  _I dAdBdCg |  _J d;d<d+g |  _K dDdEdFg |  _L dGdHdIg |  _M dddJg |  _N d!d"dKg |  _O d$d%dLg |  _P d'd(dMg |  _Q t j	 dNd g d dOg g  |  _R |  S(Q  Ni    g,@gq7@\?g<Gn@gN*
*gS4?gSq@碿geX^?g5>Kg&?g{?g)ֿg1!?gngJeQ8@g:5?gh~Br?g6t?g>)@g\c|&ſgDпg^Řj޿g͎TĿgd?g\hSȿg%ſg]Ϳgs
ֿg#E?g lEĿgJX޿gVgYX @g9f4g*8ڔƿg?=xοg#pӤɿg]=0g 9Ep?g^
]пgg	?g%gfYJ3gIȵ?gvLݕ?go1=?g624COgg-?gFkV?gm)g?$D@gڴ{ɵͿgb.jֿg[tWgFISͿg1
Ƿ?g'qѿg|#6ͿgM3Կgǘ߿g X`?g0~̿g&g]	V9g(Z&?gyxρHg,̰Ͽge\*?տg>pѿg.R\gi 3!?gF[kֿg9Gs?gtg(GgJf?gH`Z?g]y?gc%jg\K?gQlu?g?Vg)% @gC,cC7g|(Bvg"ǿgg߆?gdءgLT#gmegݳrg\	?ggAGZQȿgs]Կg8*7?gu׿g)s螿gG;ntg=	lgo>;g?grѬlgo[?g\TbοgFXTDg<i?g'?gs(CU?gmmyӿg 7C?gx@e?g>XƆgor?gBB7?gJ  v?gV!?gn `?g?<?gGء?g#?g:   e?g  r?g>1?g ?g  Q?g  ?gS$?gu?g=|?g t?g.   ?g   <?g  (%?gel?g$  @?g&  f?gTD?g  F?g?g]U?gx?g=?g*  @:@?g   ?gf_l?gMg)  @Sg|gg?g`
g?*g   KgB  _ g  {!?gB  {g  O3g   g}  Rg4Fg   (gg%  ?	g8  `B?g+  iV?gg  @gIg!  |?gv?:T?g:Tg$  @?g  \g\C?gg  @gٽ?g\ g	~)g{4i   i   gP.@g=x>?g.;TX?g@@g:}3C@gX@g?g#8@gi3ggs?g4@gnڌUǿg!?g"?g\@g+j0ï6gxZ
g{'Gؙ?g[?g%I?g++5?gمK4?g܈/ ?gK?g?զV?g;8?gKр0?g_+_?gnbҴ*h?gzu?g=bBW?g>/ˍݚ?g	R/?g
\$?g-b6>?gˬ?gDҵ@?g
|끋?gf·?gM :>g:B>gz">g1E2>g@|>g-C6>g(w?gNS?gt?g?F>%?g=!+.?gS[E?g NbUO?gy<?gHnF)?g%4 ?gck?gmѸ&	?g0g瀇?gs	@?i    g)W?gNg\W?gt?g?g&?gx"?gԷ?gCÖ?g呙
 @gy3?gMo?gu7?g#?gXF4/?gc@gՓeV?g(a)@ghr1?g+n?gԺj?gݛ@gknA=\?gt @g5?gqmM_t?g"~?g!t@gʰ?g"n+T?g9]@gU?gB!J?g[ET@gS??g(?g@
@go
@gyщ?gBr?@gc~9?gj'@g@g?gt?g
3a?g|&?grǆ?g-]4?g݇Xp9?gO?g6VD?gt'7?go^?gLA޸?g <?gu*-?g=yE7?g(,?gSk?gB?g홅sH?gعX?gSe?g@|k7?gm?g0N?gؚ`?gx?g1h'u4?g$?g=?gF= ?gY1׿6?gX_?g&a}y7g?g?g!O;ٚ?gh'u4?g|?g`G?g~6K?g^}y7g?gMWE7?g.V?g)dݲٚ?g9`d,?g/{UT?gS"?gҾ{T?gP1?i   i   i   (S   R    R   R   R   R   t   resid_pearsont	   resid_devt   resid_generalizedR   R'   R   R   R   R   R   R   R   R   R   R!   R"   R#   R&   R$   t   margeff_nodummy_dydxt   margeff_nodummy_dydx_set   margeff_nodummy_dydxmeant   margeff_nodummy_dydxmean_set   margeff_nodummy_dydxmediant   margeff_nodummy_dydxmedian_set   margeff_nodummy_dydxzerot   margeff_nodummy_dydxzero_set   margeff_nodummy_dyext   margeff_nodummy_dyex_set   margeff_nodummy_dyexmeant   margeff_nodummy_dyexmean_set   margeff_nodummy_dyexmediant   margeff_nodummy_dyexmedian_set   margeff_nodummy_dyexzerot   margeff_nodummy_dyexzero_set   margeff_nodummy_eydxt   margeff_nodummy_eydx_set   margeff_nodummy_eydxmeant   margeff_nodummy_eydxmean_set   margeff_nodummy_eydxmediant   margeff_nodummy_eydxmedian_set   margeff_nodummy_eydxzerot   margeff_nodummy_eydxzero_set   margeff_nodummy_eyext   margeff_nodummy_eyex_set   margeff_nodummy_eyexmeant   margeff_nodummy_eyexmean_set   margeff_nodummy_eyexmediant   margeff_nodummy_eyexmedian_set   margeff_nodummy_eyexzerot   margeff_nodummy_eyexzero_set   margeff_nodummy_atexog1t   margeff_nodummy_atexog1_set   margeff_nodummy_atexog2t   margeff_nodummy_atexog2_set   margeff_dummy_dydxt   margeff_dummy_dydx_set   margeff_dummy_dydxmeant   margeff_dummy_dydxmean_set   margeff_dummy_count_dydx_mediant"   margeff_dummy_count_dydx_median_set   margeff_dummy_eydxt   margeff_dummy_eydx_set   margeff_dummy_eydxmeant   margeff_dummy_eydxmean_set   margeff_count_dydxt   margeff_count_dydx_set   margeff_count_dydxmeant   margeff_count_dydxmean_set   margeff_count_dummy_dydxoverallt"   margeff_count_dummy_dydxoverall_set   margeff_count_dummy_dydxmeant   margeff_count_dummy_dydxmean_set   margeff_dummy_atexog1t   margeff_dummy_atexog1_set   margeff_dummy_atexog2t   margeff_dummy_atexog2_set
   pred_table(   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR7     sh   																$c       !   C   s  t    }  d |  _ d d d d g |  _ d d d d	 g d d
 d d g d d d d g d	 d d d g g |  _ d d d d g |  _ d |  _ d |  _ d |  _ d |  _ d |  _	 d |  _
 d |  _ d |  _ d |  _ d d d d  g |  _ d! d" g d# d$ g d% d& g d' d( g g |  _ d) d* d+ d, g |  _ 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 g  |  _ t j 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 g   |  _ 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 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 d d d d d d d d d g  |  _ t j d d g d d g g  |  _ |  S(   Ni    gQ?gZ5{7|?gA?gUI-g)Qs?gn^gc?guKg/n|?g՚Z?d?g\Nv湿gƿ?gBUÈgdoN@g,I4?gҴy?g


?g=hV@g|0:)g{4i   i   gdz/@g/ƙ(?g]6|W?g	r>@@g1C@g_|@gA5@g?g
:%-@g}ηrgW<H?g{,}@gIqټgՊTԪ?g
|$?g!X4@g;ǀ(g=)gMuϕ?gw3?g?gK;Զk?gEQW?gS^-?g<GO?gQ§2?g?g݋yg?g xq?gҦ?g6
ڻ?gH8?gf%}?gdvqv?g7h?gQ?gL+H?g 
~?gu+:H>e?g({?g&?g"<?gm\!?g~Z?g_[C{?gO?g_?gjIG?g?g2DQ?gj?g\٪x?gC5%YG?gI_?gߣ gg  g+  t g??g   g? g_2gg  ?g"g  gVݿg|g_׿g5  `5 g/  g?gzB??gK  `.?gj   ڞg?g   WMgg   h?g  `h?g `gZ?g  cܿg?+?g  ?ge @S?g!  g]ȿgji>#տggQȿgL_?g@<οg&ȿgۙSӿgsX~*޿gDNc?gPۆQǿg_g!AHg-?gnx	g"_H=ʿgOEӿgпg&wgLK1R?g{ֿgxͫ:?g:0oIgCgqD?g5C(^%?g0?gcg%P6?g^?g"07g7/N\ @gyY|Eg!=g!oտgl<bg?gQhYg?gc3g9ʿgG6u?g0kgi:տgݯg`-?g`g`̖pgj	g0eg)rg4+ۇ?g[Og, ?goDIݿgH'g?gֈ`?g1?g
<g:vP?g
<I?gs9>Zg"4q?i   i   i   (   R    R   R   R   R   R   R   R   R   R   R   R   R!   R"   R#   R&   R$   R   R   R'   R   R=   R   R<   R>   Ry   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR9     s    																									$(   R   R   t   __doc__R   R7   R9   (    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR;     s   		 		^t   RandHIEc           B   s   e  Z d  Z d   Z d   Z e   Z d   Z e   Z d   Z e   Z d   Z e   Z d   Z e   Z d   Z	 e	   Z	 d   Z
 e
   Z
 d	   Z e   Z d
   Z e   Z d   Z e   Z RS(   s(   
    Results obtained from Stata 11
    c         C   s   d |  _  d  S(   NiN  (   R   (   R   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR   0  s    c          C   s  t    }  d |  _ d d d d d d d d	 d
 d g
 |  _ d  |  _ d d d d d d d d d d g
 |  _ t j t j	 j
 t d  d d } | d  d   d f |  _ | d  d   d f |  _ d |  _ d |  _ d |  _ |  j |  j d |  _ d |  _ d |  _ d |  _ d  |  _ d! |  _ d" d# d$ d% d& d' d( d) d* d+ g
 |  _ d, d- g d. d/ g d0 d1 g d2 d3 g d4 d5 g d6 d7 g d8 d9 g d: d; g d< d= g d> d? g g
 |  _ d@ dA dB dC dD d dE dF dG d g
 |  _ dH dI dJ dK dL dM dN dO dP g	 |  _ dQ dR dS dT dU dV dW dX dY g	 |  _ dH dZ d[ dK d\ dM d] d^ d_ g	 |  _ dQ d` dS dT da dV db dc dd g	 |  _ t j t j	 j
 t de  d d |  _ |  S(f   NiN  gj媿gIϿg !&?g$gc?gv`?g7eg}WE?g&a?g.oJi?g1)g?gpjo?gRl]?glZ?glZ?g̈B?g?g۲Z?g 0?gMs܆?s   yhat_poisson.csvR	   R
   i    i   gcGrzg 1rEi	   gؗ@g_=?g2ZGղ{@gei礀@g^O$T72giЮE7gCM3@gLSp5gN36@gñN@gzągJ"?@gM~_@g؎H^O@gʭg gO@n6%ѿguih̿g O;?gZY@?gIyRgNugzyy?gp?giZZϠ?g'?g|Sg>Âv?gn?gL(?gf7?g|N|?g?p'?gڸЄ?g~8c0gx.5'gLuG9q0.gŚ^r*g9mdUi)g9-p?g|z&;?gG۟|<g|'Ǩ#<ÿg_AMgk29׹?g+@Qg?gѲڸ?gݽjUygb:?g@ԬS?gN}?g־?g<Rѝ}u?g??Tfr?gb׊A ?glQnH[?g;̀?gHl?g
@We?g:og29׹?gj?gO&g "?gw?gD&?>?gdA ?gk?gt~l?g7X@?s   poisson_resid.csv(   R    R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R   R!   R"   R#   R&   R$   t   margeff_dummy_overallt   margeff_dummy_overall_set   margeff_nodummy_overallt   margeff_nodummy_overall_seR1   (   R2   t   predict(    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   poisson3  sv    																									c          C   s  t    }  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 g |  _ d d d d d d d  d! d" d# g
 |  _ d$ d% d& d' d( d) d* d+ d, d- g
 |  _ d. d. d. d. d. d/ d/ d/ d/ d/ g
 |  _ d0 |  _	 d1 |  _
 d2 |  _ d3 d4 g d5 d6 g d7 d8 g d9 d: g d; d< g d= d> g d? d@ g dA dB g dC dD g dE dF g dG dH g g |  _ dI |  _ dJ |  _ dK |  _ dL |  _ dM |  _ dN |  _ |  S(O   NiN  gkL9g!o#ѿg?gUv4Ȇg/5?g2?g1ngfa9?gwIk?g;?gJb?g>G.z?g+wbq?g,'x?g^L'?g]	jp?g)7l?gP9aN?gqaUW?gǐzj?g3Y?gN?gRᢙ?gWr?gHc#gi;T'gkLR$@g9gT$&g;~P!@gfF:@g>}gңV?g	1~r6@gNI4:@gWͳd;g5*c^9gX;gW>w8+:g[B%<g[c?g?gn8f0+?g.-I?gtd#Zge«?g*H97?g!/U0@iN  i	   g8y㱿g.~	ˎg+(ӿg".l̿gymۜ?gg)YNB?gծgK	wg}G£?gi\?gI]M?g4?gwR	\g8sg٢@pgpf?gRy=?g>۴c?g+ ?gCfx?glf?gfvG?gG5@gF=-g|@g!/.g      "@g        (   R    R   R   t   lnalpha_std_errt   lnalphaR   R#   R$   t   fittedvaluesR!   R   R   R&   R"   R   R   R   R   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   negativebinomial_nb2_bfgs  s`    																																		c          C   s  t    }  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 g |  _ d d d d d d d  d! d" d# d$ g |  _ d% d& g d' d( g d) d* g d+ d, g d- d. g d/ d0 g d1 d2 g d3 d4 g d5 d6 g d7 d8 g d9 d: g g |  _ d; |  _ d< |  _	 d= |  _
 d> |  _ d? |  _ d@ |  _ dA |  _ dB |  _ dC dD dE dF dG dH dI dJ dK dL d@ g |  _ dM dN g dO dP g dQ dR g dS dT g dU dV g dW dX g dY dZ g d[ d\ g d] d^ g d_ d` g da db g g |  _ dc dc dc dc dc dd dd dd dd dd g
 |  _ |  S(e   NiN  g7ˡgVwoҿgZw?g!q=fgm5fg?gGzQ?gZAgJ ֨gȇ)џ?gx(;?gO'e@gƧE?gK K?gfˊu?gtH	1,?gsRÓk?g'+|h?g<=^?g,I
Q?g<a?gㄞ?g8߅?g6RQ?g. ?gnO^(g{.g3(@gfaH%g9ؙ	 @gBs'=@gw?YOgL8ga @gύ[A@gџ
M@g>ehgB<gjտgLPyпg?g;q+s?gmgqbŚgm#t?gake?gzfw?gRc?g8ⰿgK?<g²	\gL^?g.%?s?gs .?gQ@gHrh@gW]2?gQʪ?gK52!gF=-g}%Ĝ@gK528#@gkߨ(@g        g      "@g     @g*T^8g2R-d5gJ>18g:gc9<g/rvgfF1F'?gӉϺ?g:9ݤ?g:0'gC,ig/g|#_jտg>myпgZ/ɷ?g|#u?g'kgIg^`n?gQ`?gV
l){?gae?g0氿gTI~?g{`_g'3K?gV{Wr?g&7?gv?gSf?gQI@g@g!(?gn&|9?(   R    R   R   R   R   R   R#   R&   R   R   R   R!   R"   R   R   R   R$   R   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   negativebinomial_nb1_bfgs  st    																																							c          C   s  t    }  d |  _ d d d d d d d d	 d
 d g
 |  _ d d d d d d d d d d g
 |  _ d d d d d d d d d d g
 |  _ d  d! d" d# d$ d% d& d' d( d) g
 |  _ d* |  _ d+ |  _ d, |  _ d- |  _	 d. |  _
 d/ |  _ d0 |  _ d1 |  _ d2 d2 d2 d2 d2 d3 d3 d3 d3 d3 g
 |  _ d4 d5 g d6 d7 g d8 d9 g d: d; g d< d= g d> d? g d@ dA g dB dC g dD dE g dF dG g g
 |  _ |  S(H   NiN  g$HggӇoѿgMX?g,#yng4?g
x?g9S-gB^&Ǔ?g-?gS6E?g%6hv?gb`?gfn?g;Ui?gL2;?gGPmU?gcU?gmd?ghS?g}a
?g	E:g:Q>&8g,QH]:gO8gwۇ\ >;gΔnPgκ[t?gjY?gdR{?gYR:FCgpͺ$g^L{)gb3K%@gE!(gBjEa8#@gmc<@ggGv?gέ֩@gO|=@gb	C@g+H@g      "@g     @g>hBg-Tg`yHX@i    g.\6x?gɏk?gpb"gh9ugӿgBK̿gR%)?g?ǧDgvyejTgjDԢ?gջˇ?gȷl$?g*ˤ?gpu+gd];~g,Ԍ5g7rY}\?gꢕ9e?gàe@?g@)?gf{?(   R    R   R   R   R$   R#   R!   R"   R   R   R   R   R   R   R   R&   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   negativebinomial_geometric_bfgs!  sJ    																							c          C   s  t    }  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 g |  _ d |  _ d |  _ d d g d d  g d! d" g d# d$ g d% d& g d' d( g d) d* g d+ d, g d- d. g d/ d0 g d1 d2 g g |  _	 d3 |  _
 d4 |  _ d5 |  _ |  S(6   NiN  gJ'g=gzѿgScn?gԨ+CgP gWy?g%2o?g;bKg`}Tg5]t?gQ*:?g9-
n?gKl?gNgYZy?gd?g5$S.r?gqem?gT
?gǕ@?[?gs?gE2?go{?g`?g{ڥw?i	   g=J-)@g!㨲g0giGyԿg[ZzͿg^r?gy?g4ԧgJ³g?v?gW4?gWNΒâ?g4<l?gE.۱gK2z>̌g	쳿gLJI(K?gGg7;)?gm6VA0?g.?2{?g<?g[Ihd?gDE-Y.@gMWEVgEuْ@(   R    R   R   R   R   R   R   R   R!   R&   R"   t   wald_pvaluet   wald_statistic(   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   generalizedpoisson_gp2P  s8    																					c          C   s   t    }  d |  _ d d d d d d g |  _ d |  _ d	 d
 d d d d g |  _ d d g d d g d d g d d g d d g d d g g |  _ d |  _ d |  _ |  S(   NiN  g6 w?gbXgL
g?g斱TpgC·?gOg?]?ga]A?gr{?gw^?gGe"W?gL0k?g;֙t?g6uy?gxog@Cgu-kgtE)!XUg֦t?g??gEg˚g8Kr?g<?gww@gffffn@(   R    R   R   R   R   R&   R!   R"   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   zero_inflated_poisson_logits  s     										c          C   s   t    }  d |  _ d d d d d d g |  _ d |  _ d	 d
 d d d d g |  _ d d g d d g d d g d d g d d g d d g g |  _ d |  _ d |  _ |  S(   NiN  gN.!z߯?g=glz	g)?g{Un3pgyȷ?ggngs?go#?g|ҔS?gw^?gGe"W?g=P?gZf5?gO[?gagg*fD"g)(igydbRgUt?gX	?gEg&:gԕr?gl?g@gy@(   R    R   R   R   R   R&   R!   R"   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   zero_inflated_poisson_probit  s     										c          C   s   t    }  d |  _ d d d d d d g |  _ d |  _ d	 d
 d d d d g |  _ d d g d d g d d g d d g d d g d d g g |  _ d |  _ d |  _ |  S(   NiN  ghP~z?gcAaPQg gΟg;Ѻ?ghzjg/oj?g
ףpkg2?gcpp?g2?g̴X^?gch+3W?g`?gיtPӶ?gpĤ	?gT8T
gΨ*gk%Qg1iBI7g:?gD;?g{Y}a蟿gVW9g !J?g:?g̴l@g̬o@(   R    R   R   R   R   R&   R!   R"   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   zero_inflated_poisson_offset  s     										c          C   sd   t    }  d |  _ d d d d d d d g |  _ d	 |  _ d
 d d d d d d g |  _ d |  _ |  S(   NiN  ggC@g=1g Ac]˿gV?gag6vꭁ?g?g3333Mgfӓ?gAJi{@ghW!'?g#EdXk?gxĬg?gJ4?gT[r?iT (   R    R   R   R   R   R!   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt!   zero_inflated_generalized_poisson  s    				c          C   sO   t    }  d d d d d g |  _ d |  _ d d d	 d
 d g |  _ d |  _ |  S(   Ng{gUI$?g}KЏ$g)5ʿgje/5?ggW?gQgN@a?g?ܵ?gZM?goG8-xч?gVE?g\\@(   R    R   R   R   R!   (   R2   (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   zero_inflated_negative_binomial  s    			(   R   R   Rz   R   R   R   R   R   R   R   R   R   R   R   (    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyR{   ,  s,   		K		A		\		-		!										(   Rz   R   t   numpyR   t   objectR    R   t   abspatht   dirnamet   __file__R   R   R4   R;   R{   (    (    (    sR   lib/python2.7/site-packages/statsmodels/discrete/tests/results/results_discrete.pyt   <module>   s    y