U
    ÃmœdK  ã                   @   s„   d dl Zdd„ ZG dd„ dƒZG dd„ dƒZG dd	„ d	ƒZG d
d„ dƒZG dd„ dƒZG dd„ dƒZG dd„ dƒZ	G dd„ dƒZ
dS )é    Nc                 C   s@   t  | ¡} tt  t| ƒ¡ƒ}t jt j|  ||¡ddd�ddd�S )zq
    A convenience function to make the SAS covariance matrix
    compatible with statsmodels.rlm covariance
    éÿÿÿÿé   )Zaxisr   )ÚnpZasarrayÚintÚsqrtÚlenZrollZreshape)ZarrZside© r   úe/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/statsmodels/robust/tests/results/results_rlm.pyÚ_shift_intercept   s    
r
   c                   @   sœ   e Zd ZdZd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d*S )+ÚHuberú
    c                 C   s(  t  ddddg¡| _t  ddddg¡| _d	| _t  d
d
ddd
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
dg¡| _t  ddddddddddddddddddd d!d"g¡| _d#| _d$| _t  d%d&d'd(gd&d)d*d+gd'd*d,d-gd(d+d-d.gg¡| _	t  d/d0d1d2d3d4d5d5d6d7d8d9d:d;d<d=d>d?d@dAdBg¡| _
t  dCdDdEdFg¡| _d S )GNg”Òù;Šê?g 5Ÿ­¢í?gw²‚|\]À¿g
”^eƒDÀg•45¿Pv¼?gÛµÎ	kÓ?gÚÞ=~À?gË¹œ\�#@gª=ô”†@ç      ð?g%´ü%é?g˜Äöæy(à?gò@…Ÿ�×?g¥¡÷f@g‚‰ÃÂÞž ÀgÃÍo¹vµ@g¢X4¾«@g›\Je£Vú¿gü
Êý“Àgž—Ì©K³û¿g</™S—fç¿gg8Ð	Â	Àgñ@KôÅÞ?g	Eã…ú?g;ví^0àö?g�‹ç•ÀgwÒëäFé¿g¿®²¸!@g$µÔjÝ@ì?g7.ˆí“¨ê¿gš¡ÈÅÞ¬¯?gnO�ØîÞÐ?g˜z¢a°ø?g_›KyÆÕ!Àé   é   ç¦ºÇ÷kS\?ç‹YÞßÆnl¿çhÏ,_‰@F¿çVßæ­þ›?ç‰–ã.n^Š?ç 3\Ü”µ°>çq}’^¹¯¿çö÷€ðc?çøHÛÜ¾eÄ¿çPEËËç*@gæµ¯��yC@g˜8,ìí‰C@gHÒ(Qi@@g×ér•5@gÊ¥T6j¥3@g_Ay£’4@gzÉœº4»4@g:A8A1@gxø¥]Ð	+@g_—#\¯¼(@g9Q"ôù#'@gäâÁyå†*@g'½Nn‘)@g (ù¦#ï@g\i¥Räw@gã‚Ø>‰ª!@g½ntB¦À@g…};‰y!@g­0¿Ëóé*@g¯Í¥<ãê7@g3QÆð¿Ö@gm%ÏP7k@g'Ô¨4Àï¿gãp¤»À©r   ÚarrayÚparamsÚbseÚscaleÚweightsÚresidÚdf_modelÚdf_residZbcov_unscaledZfittedvaluesZtvalues©Úselfr   r   r	   Ú__init__   sÀ     ÿ                  ý                û ÿÿÿÿù	               ú ÿzHuber.__init__g¿}8gøW@g©ÐDØðÈ?g3Šå–VCÜ¿g\rÜ),ò¿gmÅþ²{òÈ?gAš±h:;‰?g…¶œKqU™¿g3ùf›Ós¿glxz¥,CÜ¿gç Á=~·?çñhãˆµøä>g±¿ìž<,ò¿gmÅþ²{ò�?g=›UŸ§T@gMø¥~ÞT´?g±PkšwœÎ¿gò^µ2á—î¿gjMóŽS´?g}è‚ú–9�?gC�ƒfÚž¿g•ñï3.l¿gušžÎ¿gD4ºƒØ™º?gì†m‹2[¿gd;ßO�—î¿g§"Æ‚Œ?gõÛ×�sŠQ@g,g~5§¿gvqà-€¿gñ�˜õb(é¿go�Å�1§¿gÚá¯Éõ�?gðPèy¢¿g7TŒó7¡`¿g)ËÇº¸½?g�Ï#„Gk¿g¸¯çŒ(é¿gº ¾eN‡?N)Ú__name__Ú
__module__Ú__qualname__Ú__doc__r%   Zhuber_h1r
   Úh1Zhuber_h2Úh2Zhuber_h3Úh3r   r   r   r	   r      sp   (             ý             ý             ýr   c                   @   sœ   e Zd ZdZd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d+S ),ÚHampelr   c                 C   s(  t  ddddg¡| _t  ddddg¡| _d	| _t  d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
d
dg¡| _t  ddddddddddddddddddddd g¡| _d!| _d"| _t  d#d$d%d&gd$d'd(d)gd%d(d*d+gd&d)d+d,gg¡| _	t  d-d.d/d0d1d2d3d3d4d5d6d7d8d9d:d;d<d=d>d?d@g¡| _
t  dAdBdCdDg¡| _d S )ENgÅr�Çó¶ç?g]âDËì™ó?gif¹®� Â¿gÅF¬Ä<DÀgëm úAÁ?gGZÇDŒ×?gûó#48 Ä?gV÷áƒKÉ'@g¯MŠÄ´@r   g¬½Ã/Íé?g_xF×\€@g¥3¬© ÀgMàs³�t@g¸*�{;%@gf�q{87ü¿gb9[£’èÀg½ÏÈR¹Àg.zŸ‘¥rõ¿g�ÆÍÛž1Àgµ9­ãQñ?gœ56C#@gávì�0Æ@g�ÏÇ)øú¿gK£°(xÐ¿g‘¨Æ*‘@gu]•—>Lì?gE1©Êö¿g&øK0×¿gÈÃk,›Ü¿g}Þ;"}ö?gÑÑGl¹£Àr   r   r   r   r   r   r   r   r   r   r   r   gz˜‹2úwC@gàQ:ÃšŠC@göƒ‘ÉoQ@@gRµ!±66@g·‡sÃ3@g,gkTý4@g£÷Y*W5@gÔ¸yÛ32@géÞXŠÃÕ)@gú˜r2/7'@gHâÄÛs%@gòù8¿@)@g…EÁc€(@g¸ò«œj·@gQT-xv@g)"&U™Á"@g1Á_‚	¸ @g^cÙÄà"@g0„¿[0+@gtô[î¨6@g@Þ'‚ü@g.R¬C£
@g.v¹�`Íí¿g¶á‡µ¾9Àr   r#   r   r   r	   r%   P   sÀ     ÿ                   þ               ú ÿÿÿÿù	              ù ÿzHampel.__init__gÙÎ÷Sã©a@gÅ™Gþ`Ò?gO@aÃÓä¿gÛe6Èú¿g‘í|?5^Ò?g¦aøˆ˜’?g	m9—âª¢¿g£¯ ÍX4}¿g¢E¶óýÔä¿g�6ŽX‹OÁ?r&   gbX9´Èú¿g™Óe1±ù˜?g-²�ïç`@g’–ÊÛNÇ?g—9]›×¿güŠ}¢ù¿g²�ï§ÆKÇ?g¡g³êsµ•?gwÖn»Ð\§¿g¶Ö	m9w¿gV-²�×¿g£’:MÄ?gàœ¥½Ág¿goƒÀÊ¡ù¿g·œKqUÙ—?g¶óýÔx`@gÄÎ:¯±«?g(~Œ¹k	™¿grm¨çoø¿gÙÎ÷Sã¥«?g‘'I×L¾™?g›=Ð
Y­¿gXÎüjp¿gú~j¼t“˜¿gòÍ67¦'È?g¢—Q,·´z¿g;ßO�—nø¿gØžY ¦–?N)r'   r(   r)   r*   r%   Z	hampel_h1r
   r+   Z	hampel_h2r,   Z	hampel_h3r-   r   r   r   r	   r.   M   sp   &             ý             ý             ýr.   c                   @   s˜   e Zd Zd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
d$S )%ÚBiSquarec                 C   s(  t  ddddg¡| _t  ddddg¡| _d	| _t  d
ddddddddddddddddddddg¡| _t  dd d!d"d#d$d%d&d'd(d)d*d+d,d-d.d/d0d1d2d3g¡| _d4| _d5| _t  d6d7d8d9gd7d:d;d<gd8d;d=d>gd9d<d>d?gg¡| _	t  d@dAdBdCdDdEdFdFdGdHdIdJdKdLdMdNdOdPdQdRdSg¡| _
t  dTdUdVdWg¡| _d S )XNgýAw®í?g¯Ok`ÌÒä?g -ô¹Á¼¿g�Îb1ƒ$EÀg°yMöl©»?gß�*AßÒ?gWƒÀ?g©@w“)#@g8|ObMA@gIrLX’ì?gíÿófšQì?gÛß—Ké?gñ�×í}Õ?g@ÎÐGî?gú™zÝ"Ðì?g(zúëî?gŠ³àMÐéï?gó3—Ècî?gãj*¬Ð÷ï?g9À•Ëªï?gÕ$W51òï?g¾×”T;ë?g>{EªÝî?g–‘Ýy ]í?g‡lÿNƒ—ï?gåå�’ìï?g:—€÷æï?g½à6Çw‘ï?g"”ï ì¯î?gÝ_ÚÃ Db?g¯ïýìË@gÁx£âAÓÀg�1®úT|@g–[ñº@g2ˆÕ@oAü¿g˜…¸jUÀgéŽQà•Kö¿g¤;F�W.Ù¿gÅgºðæMû¿g!b{K�Î¿g©žPxž²è?g¿Ñ¯)gÚÓ?gBå¹Àg'ãRt¨ßö¿g¾µ–‰@gQ¯ç­²]ë?gŸ€Q>�×¿gSí$ÇÂÚ?gsEïuÎ&ì?gÝâ
½8 ø?g:DÔ+ß$Àr   r   r   r   r   r   r   r   r   r   r   r   g!0AÓ¾C@gŒ7*4ÍC@gèU°:¸@@g`é¹C5@gƒXôÄ3@g³W­j4@gï^¹d4@g|¦oÞ´0@g‰íQ,uz,@gözÖt*@gr�²Æ,a)@gQ9�Dn,@ge\ŠõÛ*@g÷ %ò4;@g
CªI”@gôˆò�¼ @g.±�ÓS@g©¡“= @g¤£^èøë*@g"ê�o9@gQYâ;+!@g°Ð,…§@g]Yýl?´ì¿g¢5ÔbÏ¾Àr   r#   r   r   r	   r%   ‰   sÀ     ÿ                û               ú ÿÿÿÿù	              ù ÿzBiSquare.__init__gÅ�1w•V@g�¨ŒÇ?g©0¶ä Ú¿g;äf¸ñ¿g9ÖÅm4€Ç?g»ÕsÒûÆ‡?g‘Õ­ž“Þ—¿çv‰ê­�­r¿gþe÷äa¡Ú¿gQƒi>"¶?r&   g�1w-!ñ¿g9(a¦í�?gB•š=ÐôP@g]j„~¦^·?gfÚþ••Ò¿gwJëÿè¿g �¨PÝ\Œ?g’LàÖ�¿gK°8œùÕl¿gËž6çà¹?gü©ñÒMbP¿g¡L£ÉÅˆ?g—ÿ�~ûrH@g&z‹‡÷<?gªÔì�V`Ä¿g�Ñ:ªš á¿g3k) í�?gœ§:äf¸¡¿g+~©Ÿ7e¿gÔîW¾Û¼?gt^c—¨ÞZ¿g­Û ö[;�?N)r'   r(   r)   r%   Zbisquare_h1r
   r+   Zbisquare_h2r,   Zbisquare_h3r-   r   r   r   r	   r/   ˆ   sn   )             ý             ý             ýr/   c                   @   s    e Zd Zdd„ Zd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d#S )$ÚAndrewsc                 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/d0d1d2d3d4d5g| _d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHdIdJg| _d S )KNgÝµ„|Ð³í?g¼?Æä?g«>W[±¿¼¿g/Ý$�%EÀgšž^)»?gö—Ý“‡…Ò?gŸ«­Ø_v¿?g6<½R¶"@gvàœ¥=@g      @g      1@g¼n?PÖ@g3ìßÀgÙ/þ‹^p@g©A
zÅ»@g|¿ž�vJü¿gnMd-ØVÀgáÿÓ|±Iö¿gƒÿOóÅ&Ù¿g4üáÆ±Fû¿gÑ+Y¦CMÏ¿g‹ïöæX‚è?g.W\°bÓ?góëuõ(Àg»ˆ¥7’òö¿g¶šÍ”‡@g¯pŽ›Wë?gcÓ·lw×¿gNyCÅðÔÚ?g›€jqÂ;ì?g^/ãã ›ø?g×s§Pä$ÀgÏ¡4f‘ñ?gôãCÑNò¿gÕŸ#òø?g6''Á¸S@g!Õ6ãÐè¿gQšr*¾öð¿g&làŒã¿g5KÔ›éÆ¿g2+Ajíç¿g·ù|u»¿gf[_Â®Õ?gçï¸A¡ Á?gG¨È
F1õ¿gd«)›!ä¿gƒró6ÃÀî?gïiDkû×?gèýL«‹•Ä¿g>eWG>‰Ç?g–.EýÃØ?g&S‡k•å?g—üfšdSÀg€J}gˆì?g—t¡d¼>ì?g¹NÅd>é?g’V?)èŒÕ?góñº’¥=î?g "jÂì?gxÿi´æî?gà–WCeéï?g,²Ðà\î?gybè%>÷ï?gŠMnWªï?gý÷Ûü‘òï?g†‹¦¼àë?g:VAþÕî?g[hÃn“Rí?gpÚÓ_}•ï?gÅiRìï?gå›úFæï?gÛ… FrŽï?g.ó‰¼í©î?r   )r   r   r   r!   r"   r    Úsresidsr   r#   r   r   r	   r%   Å   sŽ                   ú	               ú	                ûzAndrews.__init__c                 C   s   ddddgS )N)g±áé•²ç?g“V-ò?)gºÚŠýe÷´?gX9´Èvó?)gíž<,ÔšÖ¿g·bÙ=yÀ?)gbX9´PNÀgøSã¥›ô7Àr   r#   r   r   r	   Úconf_intã   s    ÿzAndrews.conf_intgç§èHâU@g§Y Ý!ÅÆ?g’®™|³ÍÙ¿geS®ð.—ð¿g
×£p=
‡?gsº,&6—¿g˜Q,·´r¿graŠrµ?gçó·ë.Û>g�½ÞýñŽ?gÇhUM P@gµl­/Úº?g¼?ÆÔ¿g*Æù›Pˆè¿g•óÅÞ‹/Š?gö]üo%›¿g’ËH¿}m¿gNœÜïP¸?g”ö_˜LE¿gCá³up°‡?g]›�OH@goIØÕä¡?gòÏ½‡Ï¿gõÖÀV	á¿gÆÝ Z+ÚŒ?g÷;úDž¿g¦òv„Ó‚g¿gä×±ÁÂ¹?gCÅ8
Q¿g.u�×ƒI�?N)r'   r(   r)   r%   r3   Z
andrews_h1r
   r+   Z
andrews_h2r,   Z
andrews_h3r-   r   r   r   r	   r1   Ä   sp                ý             ý             ýr1   c                   @   s   e Zd Zdd„ Zdd„ ZdS )Ú
HuberHuberc                 C   sH  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/d0d1d2d3g| _d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHg| _dIdIdIdJdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdKg| _dL| _dM| _dN| _	dO| _
dP| _d S )QNg�St$—Ÿ\@gñcÌ]KÈÍ?g«ì»"øßà¿gÁnØ¶(³õ¿g½ÄX¦_"Ž?g²×»?ž¿gx¹ˆïÄ¬w¿g¯±KTo¼?gE/Ó‘„Æá>gŸu�–=”?g•Ô	hÒY@gÃï¦[vˆÃ?gÍ’ 5µlÕ¿gÔ`†�ˆó¿g·³¯<HO‘?gœ¼è+H£¿gôÃáÑÆq¿g“6U÷ÈæÀ?g¹nÀç‡a¿gOyt#,*’?gŽðHðV@g,�Ïd°?gŒ¡œhW!½¿gàóÃáQñ¿gÀ éÓ*ú“?g„d¨¿gõó¦"Æf¿gX zR&5Ä?g±3…Îkìr¿gÑ“2©¡�?g„Õ•%óE@g x	«ÎÀg%H�Z‰@gB+ø¹§@g’+¡-zÞû¿g—ñ^Ñ‡PÀgJgXs6 Àg”Î°æ:lð¿g�¿îI¿Àgìú»a[æ?gÝd'Ð þ?g2(?):ý?gç‰·3ë? ÀgŒœ9¯§ùâ¿gÙãr@gd<k3%Ôì?g¿XùÏ™ð¿gH‘££–’·¿gáiUò�‚r¿gV÷�<Œ÷?ggÿ«*ÿ Àg/Ì°Aì?gÕÅ$^Éžå¿g�\«ûÏô?g¢¥~ qîý?gP½G“ëà¿g™ûÜ·ë¿gö®À®ã¿gîÊk{tðÓ¿g·ß§�0é¿göU—ä$Ë?g×­â°Û—â?gH—€]"¾á?g{c£Ô§ºã¿g'æÃ-Ë	Ç¿g1°,fæ?gbAµO4€Ñ?g€½¦-Ì'Ô¿gRê]¼žœ¿g E’D"yV¿g¤þ)>—Ü?g~¼«¬Ú¢Àr   g¹h!1Üç?gDã¢½n¯à?)gÅ °rh‘é?gÃõ(\�Âð?gÐDØðôJÁ¿gÏfÕçj‹DÀ)g“:M„¿?gÉv¾Ÿ/Õ?g§èH.ÿÁ?gƒQI�€f%@gz6«>W[
@r   r   ©r+   r
   r,   r-   r    r2   r   r   r   r   r!   r"   r#   r   r   r	   r%   þ   sú                 ý             ý             ý                û               ú                   þzHuberHuber.__init__c                 C   s   ddddgS )N)g…|Ð³Yõá?g6<½R–ð?)g£’:M„Ù?gÜFx$û?)g×4ï8EGÚ¿g–C‹lçûÁ?)g�äòÒOÀgY·Ñ 4Àr   r#   r   r   r	   r3   %  s     ÿzHuberHuber.conf_intN©r'   r(   r)   r%   r3   r   r   r   r	   r4   ý   s   'r4   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚHampelHuberc                 C   sH  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/d0d1d2d3g| _d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHg| _dIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdIdJg| _dK| _dL| _dM| _	dN| _
dO| _d S )PNg7À[ ob@g•ñï3.Ó?g?RD†U¼å¿gOjMóû¿giqÆ0'h“?gìÀ9#J{£¿g~têÊgy~¿gjh°Â?g¢&ú|”ç>gÖ©ò=#š?gÛù~j¼¤a@gÏj�=&RÈ?g4ôOp±¢Ø¿g82�üÁÀú¿gE»
)?©–?gábE¦a¨¿gVñFæ‘?x¿g.sž±/Å?g5)Ý^Òh¿g›V
�\â˜?g(~Œ¹kÑ`@gé·¯ç¬?gÄ_“5ê!š¿gp”¼:Ç€ù¿g¶Go¸�Üš?gâKº ®¿gT©Ù­Àp¿gnü‰Ê†5É?g_ÎQÚ{¿gJÐ_è£—?gÒÌ?l|	@gâÖ€- ÀgAö!?€¼@g_ï)ù@g’«jßbéû¿gÏKë�_öÀg>�ˆÞÀg|;<¼õ¿gPÚ2×‹Àg×�[Egò?g6µŠŠÚ@gqnFß:­@g*ÀËšÅÜø¿g…›Jñ²Æ¿gLôïš¯@gÿgFŸŒŠì?g:‹0Éáö¿g�Š§:¾?Ù¿g`T…EÕß¿gT57b
‡ö?gé6œŠ8ÀgBúð£Oxî?g‡~–|¶ã¿gl…õî«œõ?gv0y;vü?gÎ�=ðaá¿gpmãzÔ2í¿g8dÎö’ýæ¿g]’h¡|Ú¿gŽ	J4õèí¿gg‚ƒëÁlÖ?gìÝ²R'1è?gW˜Q;â1é?g[VR“KÞ¿g6å6¹¨«¿goàåÄæ?g88	Œ’cÑ?g¯ñWÐÍáÛ¿gÌò,Ä¾¿g”+Î®eÃ¿gˆE¶:sÛ?g#˜´óGJÀr   gö·{îAþë?)gBÏfÕçjç?gˆc]ÜFô?g	ŠcîÂ¿gmV}®¶"DÀ)gª`TR' Á?g]ÜFxØ?gY†8ÖÅmÄ?gaTR' I(@g4¢´7øB
@r   r   r5   r#   r   r   r	   r%   +  sú                 ý             ý             ý               ú               ú                   þzHampelHuber.__init__c                 C   s   ddddgS )N)g o�Å�Ý?gÇº¸�ð?)gX9´Èvà?gzÇ):’Ëÿ?)gÎQÚ|Ý¿g¾0™*Å?)gÔšæ§PÀgcîZB>x0Àr   r#   r   r   r	   r3   S  s     ÿzHampelHuber.conf_intNr6   r   r   r   r	   r7   *  s   (r7   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚBisquareHuberc                 C   sH  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/d0d1d2d3g| _d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHg| _dIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]g| _d^| _d_| _d`| _	da| _
db| _d S )cNgØ�sF”>`@gˆ0~÷æÐ?geª`TR'ã¿gþe÷äa¡ø¿gªò=#‘?g—¨ÞØ*¡¿çt^c—¨Þz¿g¯&OYM×¿?góÜÝ+<-ä>g–˜g%­ø–?gw¾Ÿ/q[@g“È>È²`º?g+MJA·—Ð¿g(~Œ¹k‰ô¿g' ‰°áé•?gG¬Å§ ¨¿gY4��Žr¿gí»"øßJÄ?g{Úá¯Éj¿gÚÆŸ¨lX“?gËœ.‹‰óV@gCV·zNz·¿g�—6ÊúÅ?g-	PSËÖð¿g¡fHÅ«œ?gå~‡¢@Ÿ°¿gj>"¦Db¿gÖýc!:Ê?g†ZÓ¼ã}¿gèº�ê�?g‡'�CwG@g{-Ê Àgh¨Ëcë@g�54>Á@gœ’ˆýYtû¿g&ïIKØÀg¦â”ó¤ ÀgKÅ)çI)ð¿gœ�=Ú¹ÀgÙz•Ðæ?g¢ßÓ©æ¢þ?gL¡˜aHCý?g¢cÀå( Àg0§'Ó8)ã¿gå®Õo@g²žP„ë?gô‘ç¹êð¿gÞ°îÑ{`¿¿g—uœ.Í¢¿gÎ¤¦�…‰÷?gü¿H:Ý Àg£‚©·bí?g°€¶“�Qä¿gÖÓNšôyô?g¹(Äm~õý?gr�uòœà¿gjrx°Áê¿gVdÅ‡vã¿g| :ƒ�Ó¿gøÄGé¿goô\ŸlœË?g˜Êf×‰â?gßÁ‘8µá?gœÏÅN”Žã¿g9çßâl0Ç¿gß 7Á©âå?gù'dk¦Ð?g3ªbã'yÔ¿gàrÛ�ü¢¿g=�+J	Á†¿g(¤N£=|Ü?g‰½ÏÑhÀg‰ÿ°¯¹–í?g³'VÑÕî?g.w$fgë?g›¥ã¸ó—æ?geZG8ï?g�xëpþí?g™‡[$>îî?gh{þønºï?gôáªÉ;î?gÞ¢³MÝï?gD–»…dï?g=å3ï?gÚ�Ì�ëî?goX¡è„çï?g©Yîl’¦î?gá�É†ŽÍï?g€Eƒ¶Ê³ï?gE!:ùþï?g¿¦Fièÿï?gz–ÑÜËlï?g[¸t¬ß?)g‹ýe÷äaé?g¥N@aÃð?gJ+‡Á¿g_˜LŒrDÀ)çaTR' ‰À?g&S£’Ö?ç"lxz¥,Ã?g±áé•²Ì&@gâé•²q
@r   r   r5   r#   r   r   r	   r%   Y  sú                 ý             ý             ý                û                û                ûzBisquareHuber.__init__c                 C   s   ddddgS )N)gfˆc]ÜFá?gX9´Èv¾ð?)g7‰A`åÐÖ?gû:pÎˆÒû?)çKÈ=›UÛ¿g{®GázÄ?)g¼–�zžOÀg46<�2Àr   r#   r   r   r	   r3   �  s     ÿzBisquareHuber.conf_intNr6   r   r   r   r	   r8   X  s   (r8   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚAndrewsHuberc                 C   sH  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/d0d1d2d3g| _d4d5d6d7d8d9d:d;d<d=d>d?d@dAdBdCdDdEdFdGdHg| _dIdJdKdLdMdNdOdPdQdRdSdTdUdVdWdXdYdZd[d\d]g| _d^| _d_| _d`| _	da| _
db| _d S )cNgÚ|a2=`@gå`6†åÐ?g„�çÞÃ%ã¿gå~‡¢@Ÿø¿gsdå—Á‘?ga†�ˆ)¡¿r9   gÖÈ®´ŒÔ¿?gOÀÚŸ+ä>gØ�sF”ö–?gøSã¥›p[@g„J\Ç¸âº?g$—ÿ�~ûÐ¿gñ×d�zˆô¿g–?ß,Õ•?gþñ^µ2á§¿r0   gö´Ã_“5Ä?g–&¥ ÛKj¿gaO;ü5Y“?gTÆ¿ÏôV@g£Xni5$¶¿g»
)?©öÃ?gƒn/iŒÖð¿gOGÉ«sœ?gÜ)¬ÿs°¿gHPüs×b¿g:>Zœ1ÌÉ?gŠÍÇµ¡b|¿g![–¯Ëð�?gnd5•ùR@gì:ÕÓ§¾ Àg¡/E¹ð@gÞï7BÁ@gšÝÀmû¿g§PHøÿÀg?t³ Àg~è4f'ð¿gE¹Çá²»ÀgÖÈÙ4éæ?g /Ïó±þ?g†÷·,¨Uý?gßçš÷o" Àg2+1ðã¿gCFŽ”@gÔ ¹e¿‚ë?g|1ù>óð¿gàÕ®³Æ¿¿grÂó¤¿gHÈ~‹÷?gË�dnhÙ ÀgRo²°mí?gòÕ¾3Bä¿g5†;BÆ~ô?g\áóý?g#´v¿0—à¿ge~ì0Q½ê¿g_Ov^sã¿gýÓØ3‹Ó¿gGñ�–�é¿g|Ñøg¸Ë?gÉÚAMv‘â?g ñ%Å¾á?gKÒ3…ã¿gÍzô0ûÇ¿gõú�áå?gd0ýMS¤Ð?g|¬vÓ�Ô¿g« =ßß8£¿g™x�“27ˆ¿g….Ô¦Y|Ü?ghâÇß�bÀg¯³ßâúŠí?g¾Ð¹a/Òî?go@»«mVë?gš-Ž±Œæ?giL±8Ö4ï?gf}•?‰öí?gÌäåúˆéî?géM¯,¹ï?gœGD4î?gi¡ÛUÜï?gú^i7ÿï?gù~¢fÓï?gçY-ÖŽçî?g�ð Mçï?gå„¦xž î?gæLùa›Ìï?gr x¯	²ï?gKÁµ;íþï?gU\Ó¾äÿï?gžÞCÄøiï?g
sã_ÉÊß?)gšž^é?gž^)ËÇð?gÒo_ÎÁ¿gpÎˆÒÞpDÀ)r:   gTt$—ÿ�Ö?r;   g$(~Œ¹Ë&@gHPüs
@r   r   r5   r#   r   r   r	   r%   ‡  sú                 ý             ý             ý               ú               ú                ûzAndrewsHuber.__init__c                 C   s   ddddgS )N)gÝ$�•Cá?g#Ûù~j¼ð?)gáz®GáÖ?g„žÍªÏÕû?)r<   gò°PkšwÄ?)gjMóŽSœOÀgíž<,ÔŠ2Àr   r#   r   r   r	   r3   ±  s     ÿzAndrewsHuber.conf_intNr6   r   r   r   r	   r=   †  s   *r=   )Únumpyr   r
   r   r.   r/   r1   r4   r7   r8   r=   r   r   r   r	   Ú<module>   s   
=;<9-..