
p7]c           @   s   d  d l  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 S(   iNc         C   s^   t  j |   }  t t  j t |     } t  j t  j |  j | |  d d d d d d S(   sq   
    A convenience function to make the SAS covariance matrix
    compatible with statsmodels.rlm covariance
    it   axisi   i    (   t   npt   asarrayt   intt   sqrtt   lent   rollt   reshape(   t   arrt   side(    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyt   _shift_intercept   s    t   Huberc           B   s   e  Z d  Z d   Z d d d d d d d d	 d
 d d d d d	 d d g Z e e  Z d d d d d d d d d d d d d d d d g Z e e  Z d d d d d  d! d" d# d d" d$ d% d& d# d% d' g Z	 e e	  Z
 RS((   s   
    c         C   s  t  j d d d d g  |  _ t  j d d d d g  |  _ d	 |  _ t  j 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" g  |  _ d# |  _ d$ |  _ t  j d% d& d' d( g d& d) d* d+ g d' d* d, d- g d( d+ d- d. g g  |  _	 t  j d/ d0 d1 d2 d3 d4 d5 d5 d6 d7 d8 d9 d: d; d< d= d> d? d@ dA dB g  |  _
 t  j dC dD dE dF g  |  _ d  S(G   Ng;?g 5?gw|\]g
^eDg45Pv?g۵	k?g=~?g\#@g=@g      ?g%%?gy(?g@?gf@gޞ gov@gX4@g\JeVg
g̩Kg</Sfgg8		g@K?g	E?g;v^0?ggwFg!@g$j@?g7.퓨gެ?gnO?gza?g_Ky!i   i   gkS\?gYnlgh,_@FgV?g.n^?g 3\ܔ>gq}^gc?gHܾeĿgPE*@g浯yC@g8,C@gH(Qi@@gr5@gʥT6j3@g_Ay4@gzɜ44@g:A8A1@gx]	+@g_#\(@g9Q"#'@gy*@g'Nn)@g(#@g\iRw@g>!@gntB@g};y!@g0*@gͥ<7@g3Q@gm%P7k@g'Ԩ4gp(   R   t   arrayt   paramst   bset   scalet   weightst   residt   df_modelt   df_residt   bcov_unscaledt   fittedvaluest   tvalues(   t   self(    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyt   __init__   sD    	-							g}8gW@gD?g3VCܿg\r),gm{?gAh:;?gKqUg3fsglxz,Cܿg =~?gh㈵>g<,gm{?g=UT@gM~T?gPkwοg^2gjMS?g}9?gCfڞg3.lguοgD4ؙ?gm2[gd;Og"?gׁsQ@g,g~5gvq-gb(goŏ1g?gPyg7T7`g)Ǻ?g#Gkg(geN?(   t   __name__t
   __module__t   __doc__R   t   huber_h1R
   t   h1t   huber_h2t   h2t   huber_h3t   h3(    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR      s   	(t   Hampelc           B   s   e  Z d  Z d   Z d d d d d d d d	 d
 d d d d d	 d d g Z e e  Z d d d d d d d d d d d d d d d d g Z e e  Z d d d d d  d! d" d# d$ d" d% d& d' d# d& d( g Z	 e e	  Z
 RS()   s   
    c         C   s  t  j d d d d g  |  _ t  j d d d d g  |  _ d	 |  _ t  j 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  g  |  _ d! |  _ d" |  _ t  j d# d$ d% d& g d$ d' d( d) g d% d( d* d+ g d& d) d+ d, g g  |  _	 t  j d- d. d/ d0 d1 d2 d3 d3 d4 d5 d6 d7 d8 d9 d: d; d< d= d> d? d@ g  |  _
 t  j dA dB dC dD g  |  _ d  S(E   Ngr?g]D?gif¿gF<Dgm A?gGZD?g#48 ?gVK'@gMĴ@g      ?g/?g_xF\@g3 gMst@g*{;%@gfq{87gb9[gRg.zrg۞1g9Q?g56C#@gv0@g)gK(xпg*@gu]>L?gE1g&K0׿gk,ܿg};"}?gGli   i   gkS\?gYnlgh,_@FgV?g.n^?g 3\ܔ>gq}^gc?gHܾeĿgPE*@gz2wC@gQ:ÚC@goQ@@gR!66@gs3@g,gkT4@gY*W5@gԸy32@gX)@gr2/7'@gHs%@g8@)@gEc(@gj@gQT-xv@g)"&U"@g1_	 @g^c"@g0[0+@gt[6@g@'@g.RC
@g.v`gᇵ9(   R   R   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR   P   sF    	0							gSa@gG`?gO@age6g|?5^?ga?g	m9⪢g X4}gEg6XO?gh㈵>gbX9ge1?g-`@gN?g9]׿g}gK?ggs?gwn\g	m9wgV-׿g:M?gggoʡgKqUٗ?gx`@g:?g(~k	grmogS㥫?g'IL?g=
YgXjpg~jtg67'?gQ,zg;Ong؞Y?(   R   R   R   R   t	   hampel_h1R
   R   t	   hampel_h2R   t	   hampel_h3R!   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR"   M   s   	&t   BiSquarec           B   s   e  Z d    Z d d d d d d d d d	 d d
 d d d d d g Z e e  Z d d d d d d d d d d d d d d d d g Z e e  Z d d d d d d d d d d d d  d d d  d! g Z e e  Z	 RS("   c         C   s  t  j d d d d g  |  _ t  j d d d d g  |  _ d	 |  _ t  j 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/ d0 d1 d2 d3 g  |  _ d4 |  _ d5 |  _ t  j d6 d7 d8 d9 g d7 d: d; d< g d8 d; d= d> g d9 d< d> d? g g  |  _	 t  j d@ dA dB dC dD dE dF dF dG dH dI dJ dK dL dM dN dO dP dQ dR dS g  |  _
 t  j dT dU dV dW g  |  _ d  S(X   NgAw?gOk`?g -gb1$EgyMl?g߁*A?gW?g@w)#@g8|ObMA@gIrLX?gfQ?gߗK?g}?g@G?gz"?g(z?gM?g3c?gj*?g9˪?g$W51?gהT;?g>{E?gy]?glN?g叒?g:?g6w?g" ?g_àDb?g@gxAg1T|@g[@g2@oAgjUgQKg;FW.ٿggMg!b{KοgPx?gѯ)g?gBg'Rtg@gQ筲]?gQ>׿gS$?gsEu&?g
8?g:D+$i   i   gkS\?gYnlgh,_@FgV?g.n^?g 3\ܔ>gq}^gc?gHܾeĿgPE*@g!0AӾC@g7*4C@gU:@@g`C5@gX3@gWj4@g^d4@g|o޴0@gQ,uz,@gzt*@gr,a)@gQ9Dn,@ge\*@g %4;@g
CI@g @g.S@g= @g^*@g"o9@gQY;+!@g,@g]Yl?g5bϾ(   R   R   R   R   R   R   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR      sL    										gŏ1wV@g?g0ڿg;fg9m4?gsƇ?gխޗgvꭁrgeaڿgQi>"?gh㈵>g1w-!g9(a?gB=P@g]j~^?gfҿgwJg P\?gL֝gK8lg˞6?gMbPgL?g~rH@g&z<?gV`Ŀg: g3k) ?g:fg+~7egWۼ?gt^cZg۠[;?(
   R   R   R   t   bisquare_h1R
   R   t   bisquare_h2R   t   bisquare_h3R!   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR&      s   	)t   Andrewsc           B   s   e  Z d    Z d   Z d d d d d d d d d d d	 d
 d d d
 d g Z e e  Z d d d d d d d d d d d d d d d d g Z e e  Z d d d d d d d d d d d d d d d d g Z	 e e	  Z
 RS(    c         C   s!  d d d d 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  g |  _ d! d" d# d$ d% d& d' d( d) d* d+ d, d- d. d/ d0 d1 d2 d3 d4 d5 g |  _ d6 d7 d8 d9 d: d; d< d= d> d? d@ dA dB dC dD dE dF dG dH dI dJ g |  _ d  S(K   Ngݵ|г?g??g>W[g/$%Eg^)?gݓ?g_v?g6<R"@gv=@g      @g      1@gn?P@g3g/^p@gA
zŻ@g|vJgnMd-Vg|IgO&ٿg4ƱFg+YCMϿgX?g.W\b?gu(g7g͔@gpW?gcӷlw׿gNyC?gjq;?g^/ ?gsP$gϡ4f?gCNg#?g6''S@g!6gQr*g&lg5Kԛƿg2+Ajg|ugf[_®?gA ?gG
F1gd)!gr6?giDk?gLĿg>eWG>?g.E?g&Sk?gfdSgJ}g?gtd>?gNd>?gV?)?g񺒥=?g "j?gxi?gWCe?g,\?gyb%>?gMnW?g?g?g:VA?g[hnR?gp_}?giR?gF?gۅ Fr?g.?i    (   R   R   R   R   R   R   t   sresidsR   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR      s,    			c         C   s   d	 d
 d d g S(   Ng镲?gV-?gڊe?gX9v?g<,Ԛֿgb=y?gbX9PNgS㥛7(   g镲?gV-?(   gڊe?gX9v?(   g<,Ԛֿgb=y?(   gbX9PNgS㥛7(    (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyt   conf_int   s    	gHU@gY!?g|ٿgeS.g
ףp=
?gs,&6gQ,rgrar?g.>g?ghUMP@gl/ں?g?Կg*Pgދ/?g]o%gH}mgNP?g_LEgCup?g]OH@goI?gϽϿgV	g Z+ڌ?g;Dgvӂgg¹?gC8
Qg.u׃I?(   R   R   R   R,   t
   andrews_h1R
   R   t
   andrews_h2R   t
   andrews_h3R!   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR*      s   		t
   HuberHuberc           B   s   e  Z d    Z d   Z RS(   c         C   s  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 g |  _ t |  j  |  _ 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/ 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 g |  _ dI dI dI dJ dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dK g |  _ dW |  _ dX |  _ dT |  _	 dU |  _
 dV |  _ d  S(Y   NgSt$\@gc]K?g"gnض(gX_"?g׻?gxĬwgKTo?gE/ӑ>gu=?g	hY@g[v?g͒ 5lտg`g<HO?g+Hgqg6U?gnagOyt#,*?gHV@g,d?ghW!gQg *?gdg"fgX zR&5?g3krgѓ2?gՕ%E@g x	g%HZ@gB+@g+-zg^чPgJgXs6 gΰ:lgIga[?gd'Р?g2(?):?g牷3? g9gr@gd<k3%?gXϙgHgiUrgV<?gg* g/A?g$^ɞg\?g~ q?gPGgܷggk{tӿg0gU$?g׭ۗ?gH]"?g{cԧg'-	ǿg1,f?gbAO4?g-'ԿgR]g ED"yVg)>?g~ڢi   gh!1?gD㢽n?g rh?g(\?gDJgfjDg:M?gv/?gH.?gQIf%@gz6>W[
@i   i   (   g rh?g(\?gDJgfjD(   g:M?gv/?gH.?gQIf%@(   R   R
   R   R!   R   R+   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR      s<    -				c         C   s   d	 d
 d d g S(   Ng|гY?g6<R?g:M?gFx$?g48EGڿgCl?gOgY 4(   g|гY?g6<R?(   g:M?gFx$?(   g48EGڿgCl?(   gOgY 4(    (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR,   %  s    (   R   R   R   R,   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR0      s   	't   HampelHuberc           B   s   e  Z d    Z d   Z RS(   c         C   s  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 g |  _ t |  j  |  _ 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/ 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 g |  _ dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dI dJ g |  _ dV |  _ dW |  _ dS |  _	 dT |  _
 dU |  _ d  S(X   Ng7[ ob@g3.?g?RDUgOjMgiq0'h?g9#J{g~tgy~gjh?g&|>g֩=#?g~ja@gj=&R?g4Opؿg82gE
)??gbEagVF?xg.s/?g5)^hgV
\?g(~k`@g鷯?g_5!gp:ǀgGoܚ?gKgTpgnʆ5?g_Q{gJ_?g?l|	@gր- gA!?@g_)@gjbgK_g>g|;<gP2׋g[Eg?g6@gqnF:@g*˚gJƿgL@ggF?g:0g:?ٿg`TE߿gT57b
?g68gBOx?g~|gl?gv0y;v?gΝ=agpmz2g8dg]h|ڿg	J4ggl?gݲR'1?gWQ;1?g[VRK޿g66go?g88	c?gWۿg,ľg+ήeÿgE:s?g#GJi   g{A?gBfj?gc]F?g	c¿gmV}"Dg`TR'?g]Fx?gY8m?gaTR'I(@g47B
@i   i   (   gBfj?gc]F?g	c¿gmV}"D(   g`TR'?g]Fx?gY8m?gaTR'I(@(   R   R
   R   R!   R   R+   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR   +  s>    3				c         C   s   d	 d
 d d g S(   Ng oŏ?gǺ?gX9v?gz):?gQ|ݿg0*?gԚPgcZB>x0(   g oŏ?gǺ?(   gX9v?gz):?(   gQ|ݿg0*?(   gԚPgcZB>x0(    (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR,   S  s    (   R   R   R   R,   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR1   *  s   	(t   BisquareHuberc           B   s   e  Z d    Z d   Z RS(   c         C   s  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 g |  _ t |  j  |  _ 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/ 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 g |  _ dI dJ dK dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ d[ d\ d] g |  _ di |  _ dj |  _ df |  _	 dg |  _
 dh |  _ d  S(k   Ng؁sF>`@g0~?ge`TR'geag=#?g*gt^czg&OYM׿?g+<->gg%?gw/q[@g>Ȳ`?g+MJAпg(~kg'?gGŧ gY4rg"J?g{jgƟlX?g˜.V@gCVzNzg6?g-	PSgfHū?g~@gj>"Dbgc!:?gZӼ}g?g'CwG@g{- ghc@g54>@gYtg&IKg gK)I)g=ڹgz?gө?gLaHC?gc( g0'8)go@gP?ggް{`gu.͢gΤ?gH: gb?gQgNy?g(m~?grugjrxgVdŇvg| :ӿgGgo\l?gf׉?g8?gNg9l0ǿgߠ7?g'dk?g3b'yԿgrg=+J	g(N=|?ghg?g'V?g.w$fg?g?geZG8?gxp?g[$>?gh{n?g;?gޢM?gDd?g=3?gڍ̐?goX?gYl?gɆ?gEʳ?gE!:?gFi?gzl?g[t?gea?gN@a?gJ+g_LrDgaTR'?g&S?g"lxz,?g镲&@g镲q
@i   i   (   gea?gN@a?gJ+g_LrD(   gaTR'?g&S?g"lxz,?g镲&@(   R   R
   R   R!   R   R+   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR   Y  s@    				c         C   s   d	 d
 d d g S(   Ngfc]F?gX9v?g7A`?g:pΈ?gK=Uۿg{Gz?gzOg46<2(   gfc]F?gX9v?(   g7A`?g:pΈ?(   gK=Uۿg{Gz?(   gzOg46<2(    (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR,     s    (   R   R   R   R,   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR2   X  s   	(t   AndrewsHuberc           B   s   e  Z d    Z d   Z RS(   c         C   s  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 g |  _ t |  j  |  _ 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/ 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 g |  _ dI dJ dK dL dM dN dO dP dQ dR dS dT dU dV dW dX dY dZ d[ d\ d] g |  _ di |  _ dj |  _ df |  _	 dg |  _
 dh |  _ d  S(k   Ng|a2=`@g`6?g%g~@gsd?ga)gt^czgȮԿ?gOڟ+>g؁sF?gS㥛p[@gJ\Ǹ?g$~пgdzg?,Օ?g^2᧿gvꭁrg_5?g&KjgaO;5Y?gTƿV@gXni5$g
)??gn/igOGɫs?g)sgHPsbg:>Z1?gǵb|g![?gnd5R@g:ӧ g/E@g7B@gmgPHg?t g~4f'gE᲻g4?g /?g,U?go" g2+1gCF@gԠe?g|1>gƿgrgH~?g˝dnh gRom?gվ3Bg5;B~?g\?g#v0ge~0Qg_Ov^sg3ӿgG񏖁g|g?gAMv?g %ž?gK3gz0ǿg?gd0MS?g|vӁԿg =8gx27g.ԦY|?ghߐbg?gйa/?go@mV?g-?giL84?gf}??g?gM,?gGD4?giU?g^i7?g~f?gY-֎?g M?g儦x?gLa?grx	?gK;?gU\Ӿ?gCi?g
s_?g^?g^)?go_gpΈpDgaTR'?gTt$?g"lxz,?g$(~&@gHPs
@i   i   (   g^?g^)?go_gpΈpD(   gaTR'?gTt$?g"lxz,?g$(~&@(   R   R
   R   R!   R   R+   R   R   R   R   R   R   (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR     sD    				c         C   s   d	 d
 d d g S(   Ng$C?g#~j?gzG?gͪ?gK=UۿgPkw?gjMSOg<,Ԋ2(   g$C?g#~j?(   gzG?gͪ?(   gK=UۿgPkw?(   gjMSOg<,Ԋ2(    (   R   (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR,     s    (   R   R   R   R,   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyR3     s   	*(   t   numpyR   R
   t   objectR   R"   R&   R*   R0   R1   R2   R3   (    (    (    sK   lib/python2.7/site-packages/statsmodels/robust/tests/results/results_rlm.pyt   <module>   s   	
=;<9-..