Ë
    âQ(hýT ã                   ó  — d Z ddlZddlZddlZddlmZ ddlmZ ddlmZm	Z	 ddl
mc mZ ddlmZ ddlmZ ddlZddlmZ dd	lmZmZmZmZmZmZmZ dd
lmZ ddlmZmZ ddl m!Z! ddl"m#Z#m$Z$  G d„ d«      Z%dBd„Z&dBd„Z'e! G d„ d«      «       Z(e! G d„ d«      «       Z) G d„ d«      Z* G d„ d«      Z+e! G d„ d«      «       Z,e! G d„ d«      «       Z- G d„ d«      Z. G d „ d!«      Z/e! G d"„ d#«      «       Z0e! G d$„ d%«      «       Z1d&„ Z2d'„ Z3 G d(„ d)«      Z4d*„ Z5d+„ Z6d,„ Z7d-„ Z8e! G d.„ d/«      «       Z9 G d0„ d1«      Z: G d2„ d3«      Z; G d4„ d5«      Z< G d6„ d7«      Z= G d8„ d9«      Z> G d:„ d;«      Z? G d<„ d=«      Z@e! G d>„ d?«      «       ZA G d@„ dA«      ZBy)Cz?
Tests for the stats.mstats module (support for masked arrays)
é    N)Únan)ÚmaskedÚnomask)Ústatsé   )Úcheck_named_results)Úraises)Úassert_equalÚassert_almost_equalÚassert_array_almost_equalÚassert_array_almost_equal_nulpÚassert_Úassert_allcloseÚassert_array_equal)Úsuppress_warnings)Ú_mstats_basicÚ	_stats_py)Úskip_xp_invalid_arg)ÚSmallSampleWarningÚtoo_small_1d_not_omitc                   ó   — e Zd Zd„ Zy)ÚTestMquantilesc                 óÎ   — t        j                  g d¢g d¢g d¢g d¢g d¢g d¢g d¢g d¢g d	¢g d
¢g d¢g«      }g d¢g d¢g d¢g}t        j                  |dd¬«      }t	        ||«       y )N)ç      @ç      @ç      ð?)g     €G@ç      .@ç       @)g     €H@ç      B@ç      @)r   ç     €C@ç      @)ç      E@ç      D@ç     8�À)ç     €D@r&   r%   )r   r%   r%   )r!   r%   r%   )g     €E@r%   r%   )r$   r%   r%   )r   r%   r%   )g3333333@g333333-@g333333÷?)r$   g     ÀB@ç      @)gffffffE@gfffffD@gffffff@r   )r   é2   )ÚaxisÚlimit)ÚnpÚarrayÚmstatsÚ
mquantilesr   )ÚselfÚdataÚdesiredÚquantss       úa/var/www/skyplay_api_hub/venv/lib/python3.12/site-packages/scipy/stats/tests/test_mstats_basic.pyÚtest_mquantiles_limit_keywordz,TestMquantiles.test_mquantiles_limit_keyword   sg   € ä�x‰xšÚ'Ú'Ú'Ú*Ú*Ú+Ú,Ú,Ú,Ú,ð
.ó 
/ˆò &Ú$Ú&ð(ˆô ×"Ñ" 4¨a°wÔ?ˆÜ˜F GÕ,ó    N)Ú__name__Ú
__module__Ú__qualname__r4   © r5   r3   r   r      s   „ ó-r5   r   c                 ó|   — t        j                  | ||¬«      }t        |||¬«       t        |j                  |«       y ©N)r)   Údtype©Úrtol)r-   Úgmeanr   r
   r<   ©Ú
array_liker1   r)   r<   r>   Úxs         r3   Úcheck_equal_gmeanrC   /   s.   € ä�‰�Z d°%Ô8€AÜ�A�w TÕ*Ü�—‘˜%Õ r5   c                 ó|   — t        j                  | ||¬«      }t        |||¬«       t        |j                  |«       y r;   )r   Úhmeanr   r
   r<   r@   s         r3   Úcheck_equal_hmeanrF   6   s.   € Ü�‰�J T°Ô7€AÜ�A�w TÕ*Ü�—‘˜%Õ r5   c                   ó†   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zej                  j                   eed«       d¬«      d	„ «       Zd
„ Zy)ÚTestGeoMeanc                 óT   — g d¢}t        j                  dd«      }t        ||d¬«       y )N©r   é   é   é   é   ç      Ð?ç›+¡†›„=r=   )r+   ÚpowerrC   ©r/   Úar1   s      r3   Útest_1dzTestGeoMean.test_1d>   s#   € ÚˆÜ—(‘(˜7 EÓ*ˆÜ˜!˜W¨5Ö1r5   c                 óÌ   — t        j                  g d¢«      }d}t        ||«       t        j                  g d¢g d¢¬«      }t        j                  dd«      }t        ||d¬	«       y )
N©
é
   é   é   é(   r(   é<   éF   éP   éZ   éd   gs©ÑÅ¤F@rJ   ©r   r   r   r   ©Úmaské   çUUUUUUÕ?rP   r=   )Úmar,   rC   r+   rQ   rR   s      r3   Ú
test_1d_mazTestGeoMean.test_1d_maC   sL   € ä�H‰HÒ>Ó?ˆØˆÜ˜!˜WÔ%ä�H‰H’\ªÔ5ˆÜ—(‘(˜5 %Ó(ˆÜ˜!˜W¨5Ö1r5   c                 ój   — t         j                  j                  g d¢g d¢¬«      }d}t        ||«       y )NrV   ©
r   r   r   r   r   r   r   r   r   r   ra   gdåÝq_¼D@©r+   re   r,   rC   rR   s      r3   Útest_1d_ma_valuezTestGeoMean.test_1d_ma_valueM   s.   € ä�E‰E�K‰KÒAÚ;ð ó =ˆàˆÜ˜!˜WÕ%r5   c                 ób   — t         j                  j                  g d¢«      }d}t        ||«       y )N)
rW   rX   rY   rZ   r(   r[   r\   r]   r^   r   r   ri   rR   s      r3   Útest_1d_ma0zTestGeoMean.test_1d_ma0T   s$   € ä�E‰E�K‰KÒ?Ó@ˆØˆÜ˜!˜WÕ%r5   c                 óÔ   — t         j                  j                  g d¢«      }t         j                  }t        j                  d¬«      5  t        ||«       d d d «       y # 1 sw Y   y xY w)N)
rW   rX   rY   rZ   r(   r[   r\   r]   r^   éÿÿÿÿÚignore)Úinvalid)r+   re   r,   r   ÚerrstaterC   rR   s      r3   Útest_1d_ma_infzTestGeoMean.test_1d_ma_infZ   sI   € ä�E‰E�K‰KÒ@ÓAˆÜ—&‘&ˆÜ�[‰[ Ô*ñ 	*Ü˜a Ô)÷	*÷ 	*ñ 	*ús   ÁAÁA'Úfloat96úcannot find float96 so skipping©Úreasonc                 óÚ   — t        j                  g d¢g d¢¬«      }t        j                  dd«      j	                  t        j
                  «      }t        ||t        j
                  d¬«       y )NrJ   r`   ra   rc   rd   rP   )r<   r>   )re   r,   r+   rQ   Úastypers   rC   ©r/   rS   Ú
desired_dts      r3   Útest_1d_float96zTestGeoMean.test_1d_float96a   sE   € ô �H‰H’\ªÔ5ˆÜ—X‘X˜e UÓ+×2Ñ2´2·:±:Ó>ˆ
Ü˜!˜Z¬r¯z©zÀÖFr5   c                 óÔ  — t        j                  g d¢g d¢g d¢gg d¢g d¢g d¢g¬«      }t        j                  g d¢«      }t        ||dd¬«       t        j                  t        j                  d	d
«      t        j                  dd«      t        j                  dd«      g«      }t        ||dd¬«       g d¢g d¢g d¢g}d}t        t        j                   j                  |«      |«       y )NrJ   ©r   r   r   r   ©r   r   r   r   ©r   r   r   r   ra   r   rP   ©r)   r>   rN   rO   rc   ç      à?rM   rn   ©rW   rX   rY   rZ   ©r(   r[   r\   r]   ©r^   r_   én   éx   g/,$»qJ@)re   r,   r+   rC   rQ   rR   s      r3   Ú
test_2d_mazTestGeoMean.test_2d_mah   s·   € Ü�H‰H’l¢L²,Ð?Ú'ª²|ÐDôFˆä—(‘(š<Ó(ˆÜ˜!˜W¨1°5Õ9ä—(‘(œBŸH™H W¨eÓ4ÜŸH™H S¨%Ó0ÜŸH™H S¨%Ó0ð2ó 3ˆô 	˜!˜W¨2°EÕ:ò Ò/Ò1DÐEˆØˆÜœ"Ÿ%™%Ÿ+™+ a›.¨'Õ2r5   N)r6   r7   r8   rT   rf   rj   rl   rr   ÚpytestÚmarkÚskipifÚhasattrr+   r{   r‡   r9   r5   r3   rH   rH   <   s[   „ ò2ò
2ò&ò&ò*ð ‡[�[×Ñ™G B¨	Ó2Ð2Ø@ð ó BñGóBðGó
3r5   rH   c                   ón   — e Zd Zd„ Zej
                  j                   eed«       d¬«      d„ «       Z	d„ Z
y)ÚTestHarMeanc                 ó  — t        j                  g d¢g d¢¬«      }d}t        ||d¬«       t        j                   j                  g d¢«      }d}t        ||«       t        j                   j                  g d¢g d	¢¬«      }d
}t        ||«       y )NrJ   r`   ra   ç/ºè¢‹.ú?rP   r=   rV   g¶=b¹#A@rh   g¹´óÜOÐ?@©re   r,   rF   r+   rR   s      r3   rT   zTestHarMean.test_1d{   sq   € Ü�H‰H’\ªÔ5ˆØ+ˆÜ˜!˜W¨5Õ1ä�E‰E�K‰KÒAÓBˆØˆÜ˜!˜WÔ%ä�E‰E�K‰KÒAÚ;ð ó =ˆàˆÜ˜!˜WÕ%r5   rs   rt   ru   c                 ó¼   — t        j                  g d¢g d¢¬«      }t        j                  dt        j                  ¬«      }t        ||t        j                  ¬«       y )NrJ   r`   ra   r�   ©r<   )re   r,   r+   Úasarrayrs   rF   ry   s      r3   r{   zTestHarMean.test_1d_float96‰   s9   € ô �H‰H’\ªÔ5ˆÜ—Z‘ZÐ 9ÄÇÁÔLˆ
Ü˜!˜Z¬r¯z©zÖ:r5   c                 ó4  — t        j                  g d¢g d¢g d¢gg d¢g d¢g d¢g¬«      }t        j                  g d¢«      }t        ||dd¬«       g d	¢}t        ||d
d¬«       g d¢g d¢g d¢g}d}t        t        j                   j                  |«      |«       y )NrJ   r}   r~   r   ra   r   rP   r€   )gº…ëQ¸þ?g433333@gš™™™™™ù?rn   r‚   rƒ   r„   gmôW¶UC@r�   rR   s      r3   Útest_2dzTestHarMean.test_2d�   s|   € Ü�H‰H’l¢L²,Ð?Ú'ª²|ÐDôFˆä—(‘(š<Ó(ˆÜ˜!˜W¨1°5Õ9âLˆÜ˜!˜W¨2°EÕ:âÒ/Ò1DÐEˆØˆÜœ"Ÿ%™%Ÿ+™+ a›.¨'Õ2r5   N)r6   r7   r8   rT   rˆ   r‰   rŠ   r‹   r+   r{   r•   r9   r5   r3   r�   r�   y   sE   „ ò&ð ‡[�[×Ñ™G B¨	Ó2Ð2Ø@ð ó Bñ;óBð;ó
3r5   r�   c                   ó   — e Zd Zd„ Zy)ÚTestRankingc                 óž  — t        j                  g d¢«      }t        t        j                  |«      g d¢«       t
        |ddg<   t        t        j                  |«      g d¢«       t        t        j                  |d¬«      g d¢«       t        j                  g d	¢«      }t        t        j                  |«      g d
¢«       t        j                  g d¢g d¢g«      }t        t        j                  |«      g d¢g d¢g«       t        t        j                  |d¬«      g d¢g d¢g«       t        t        j                  |d¬«      g d¢g d¢g«       y )N)
r   r   r   r   rK   rL   rM   é   r™   rc   )
r   rL   rL   rL   r™   rc   é   ç      !@r›   rW   rL   rM   )
r   r'   r'   r   r   rM   r™   ç      @rœ   é   T)Úuse_missing)
r   r'   r'   ç      @rŸ   rM   r™   rœ   rœ   r�   )
r   r   r™   r   rK   rM   rL   r™   r   rc   )
r   rL   r›   rL   r™   rš   rc   r›   rL   rW   )r   r   r   r   rK   )rL   rM   r™   r™   rc   )r   rL   rL   rL   r™   )rc   rš   r›   r›   rW   r   ©r)   )r   rK   ç      @r¡   r™   r   ©r   r   r   r   r   ©rK   rK   rK   rK   rK   )re   r,   r   r-   Úrankdatar   ©r/   rB   s     r3   Útest_rankingzTestRanking.test_rankingŸ   sð   € Ü�H‰HÒ+Ó,ˆÜœFŸO™O¨AÓ.Ú6ô	8äˆˆ1ˆQˆ%‰ÜœFŸO™O¨AÓ.Ú9ô	;äœFŸO™O¨A¸4Ô@Ú=ô	?ä�H‰HÒ+Ó,ˆÜœFŸO™O¨AÓ.Ú6ô	8ä�H‰H’k¢<Ð0Ó1ˆÜœFŸO™O¨AÓ.Ú(Ò*:Ð;ô	=äœFŸO™O¨A°AÔ6Ú(ª/Ð:ô	<äœFŸO™O¨A°1Ô5Ú(ª,Ð7õ	9r5   N)r6   r7   r8   r¦   r9   r5   r3   r—   r—   ž   s   „ ó9r5   r—   c                   ó   — e Zd Zd„ Zd„ Zd„ Zd„ Zej                  j                   e
j                  «       dk(  d¬«      d„ «       Zej                  j                   e
j                  «       dk(  d¬«      ej                  j                  d	„ «       «       Zd
„ Zej                  j!                  dd«      ej                  j!                  dd«      d„ «       «       Zd„ Zej                  j&                  d„ «       Zed„ «       Zy)ÚTestCorrc           
      óÚ  — t        j                  d«      }t        j                  «       5  t        j                  dt
        «       t        t        j                  ||«      d   d«       t        t        j                  ||d d d…   «      d   d«       t        j                  |d¬«      }t        j                  ||«      }t        |d   t        u «       t        |d	   t        u «       d d d «       t        j                  g d
¢«      }t        j                  g d¢«      }t        j                  ||«      \  }}t        |t        j                  d«      dz  «       t        |d«       g d¢}t        j                  g d¢|¬«      }t        j                  g d¢|¬«      }	t        j                  ||	«      \  }}t        |t        j                  d«      dz  «       t        |d«       y # 1 sw Y   ŒüxY w)NrW   Úerrorr   r   rn   ç      ð¿Tra   r   )r«   ç        r   )r   r   rL   rL   rK   rd   )FFFT)r«   r¬   r   g     ÀX@)r   r   rL   rn   )re   ÚarangeÚwarningsÚcatch_warningsÚsimplefilterÚRuntimeWarningr   r-   Úpearsonrr,   r   r   r+   Úsqrt)
r/   rB   ÚprÚx1Úy1ÚrÚprb   Úx2Úy2s
             r3   Útest_pearsonrzTestCorr.test_pearsonrµ   sm  € ä�I‰I�b‹MˆÜ×$Ñ$Ó&ñ 	%ô ×!Ñ! '¬>Ô:ä¤§¡°°1Ó 5°aÑ 8¸#Ô>Ü¤§¡°°1±T°r°T±7Ó ;¸AÑ >ÀÔEä—‘˜ Ô&ˆAÜ—‘  AÓ&ˆBÜ�B�q‘EœV�OÔ$Ü�B�q‘EœV�OÔ$÷	%ô �X‰XÒ&Ó'ˆÜ�X‰X’iÓ ˆÜ�‰˜r 2Ó&‰ˆˆ1Ü˜AœrŸw™w q›z¨!™|Ô,Ü˜A˜uÔ%ò +ˆÜ�X‰XÒ,°4Ô8ˆÜ�X‰X’m¨$Ô/ˆÜ�‰˜r 2Ó&‰ˆˆ1Ü˜AœrŸw™w q›z¨!™|Ô,Ü˜A˜uÕ%÷5	%ð 	%ús   ªB<G!Ç!G*c                 óˆ  — t         j                  j                  g d¢g d¢¬«      }t         j                  j                  g d¢g d¢¬«      }t        j                  g d¢«      }t        j                  g d¢«      }t	        j
                  ||«      \  }}t        j
                  ||«      \  }}t        ||«       t        ||«       y )N©r   rK   rL   rM   r™   rc   )r   r   r   r   r   r   ra   )é	   r�   rš   rc   r™   r¾   )r   r   r   r   r   r   )r   rM   r™   rc   )r¾   rc   r™   r¾   )r+   re   Úmasked_arrayr,   r-   r²   r   r
   )	r/   ÚmxÚmyrB   ÚyÚmrÚmpr·   r¸   s	            r3   Útest_pearsonr_misaligned_maskz&TestCorr.test_pearsonr_misaligned_maskÔ   s�   € Ü�U‰U×ÑÒ 2Ò9KÐÓLˆÜ�U‰U×ÑÒ 2Ò9KÐÓLˆÜ�H‰H’\Ó"ˆÜ�H‰H’\Ó"ˆÜ—‘  RÓ(‰ˆˆBÜ�~‰~˜a Ó#‰ˆˆ1Ü�R˜ÔÜ�R˜Õr5   c                 ó   — g d¢g d¢}}t        t        j                  ||«      d   d«       ddddt        j                  gd	d
d
dt        j                  g}}t        j                  |«      t        j                  |«      }}t        t        j                  ||«      d   d«       g d¢}g d¢}t        t        j                  ||«      d   d«       dddddddddddddddt        j                  g}dddd d!d"d#d$d%d"dd&dd'dt        j                  g}t        j                  |«      t        j                  |«      }}t        t        j                  ||«      d   d«       t        t        d(«      «      }t        t        d(«      «      }|d)   |d   c|d<   |d)<   |d*   |d+   c|d+<   |d*<   |d,   |d-   c|d-<   |d,<   t        t        j                  ||«      d   d.«       t        j                  ||«      }d/}t        ||d0¬1«       y )2N)ç333333@ç      @ç®Gáz®	@çHáz®G@)çffffffú?ç…ëQ¸@rÌ   çÍÌÌÌÌÌ@r   gIœQ=ä¿rÇ   rÈ   rÉ   rÊ   rË   rÌ   rÍ   ©r   ç33333³G@r#   çš™™™™™%@çÍÌÌÌÌN@ç333333û?ç      P@çÍÌÌÌÌŒO@r   çffffffö?çš™™™™™@ç333333Ó?ç333333@r×   çÍÌÌÌÌÌ@©çš™™™™™6@çš™™™™™ @ç333333F@çÍÌÌÌÌÌ'@çš™™™™™8@ç333333ã?çÍÌÌÌÌÌ@çÍÌÌÌÌÌD@r¬   rà   rÙ   çffffff@r   ç333333ó?rÕ   gàˆ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ä   éÐ  r¾   i²  rW   iå  i³  gV-²�ïï?©ÚcorrelationÚpvalueT©re   )
r   r-   Ú	spearmanrr+   r   re   Úfix_invalidÚlistÚranger   ©r/   rB   rÂ   ÚresÚ
attributess        r3   Útest_spearmanrzTestCorr.test_spearmanrÞ   só  € â'Ò)>ˆAˆÜœF×,Ñ,¨Q¨qÓ1°!Ñ4°jÔAØ˜˜T $¤r§v¡vÐ.°°T¸$¸tÄBÇFÁFÐ/KˆAˆÜ—.‘. Ó#¤R§^¡^°AÓ%6ˆAˆÜœF×,Ñ,¨Q¨qÓ1°!Ñ4°jÔAò0ˆò0ˆäœF×,Ñ,¨Q¨qÓ1°!Ñ4°iÔ@Ø�$˜˜d D¨#¨t°TØ�#�s˜C  c¨3´·±ð8ˆà�3˜˜d D¨#¨s°DØ�#�s˜C  c¨3´·±ð8ˆä—.‘. Ó#¤R§^¡^°AÓ%6ˆAˆÜœF×,Ñ,¨Q¨qÓ1°!Ñ4°iÔ@ô
 ”�t“ÓˆÜ”�t“ÓˆØ�q‘T˜1˜Q™4ˆ
ˆˆ!‰ˆa�‰dØ˜#™  "¡ˆˆˆ"‰ˆq�‰vØ˜D™' 1 S¡6ˆˆˆ#‰��$‘ô 	œF×,Ñ,¨Q¨qÓ1°!Ñ4°eÔ<ô ×Ñ˜q !Ó$ˆØ.ˆ
Ü˜C °Ö5r5   c                 óR  — g d¢}g d¢}d}t        j                  ||«      \  }}t        ||«       t        |d«       t        j                  ||d¬«      \  }}t        ||«       t        |d«       t        j                  ||d¬«      \  }}t        ||«       t        |d	«       d
}t        j                  dd|«      }d|z  t        j
                  j                  |«      z   }t        j                  ||«      \  }}t        j                  ||d¬«      \  }	}
t        |
|dz  «       t        j                  ||d¬«      \  }}t        |d|dz  z
  «       ||	cxk(  r|k(  sJ ‚ J ‚t        j                  t        d¬«      5  t        j                  ||d¬«       d d d «       y # 1 sw Y   y xY w)NrÎ   rÚ   g+þ;
æ?g œ·‚r?Úgreater©Úalternativeg œ·‚b?ÚlessgšæcH}íï?r_   r   r™   çš™™™™™¹?rK   r   zalternative must be 'less'...©Úmatchz	ekki-ekki)
r-   rê   r   r+   ÚlinspaceÚrandomÚrandrˆ   r	   Ú
ValueError)r/   rB   rÂ   Úr_expr·   r¸   ÚnÚstat1Úp1Ústat2Úp2Ústat3Úp3s                r3   Útest_spearmanr_alternativez#TestCorr.test_spearmanr_alternative  s„  € ò0ˆò0ˆð #ˆä×Ñ  1Ó%‰ˆˆ1Ü˜˜5Ô!Ü˜Ð,Ô-ä×Ñ  1°)Ô<‰ˆˆ1Ü˜˜5Ô!Ü˜Ð,Ô-ä×Ñ  1°&Ô9‰ˆˆ1Ü˜˜5Ô!Ü˜˜?Ô+ð ˆÜ�K‰K˜˜1˜aÓ ˆØ�‰E”B—I‘I—N‘N 1Ó%Ñ%ˆä×$Ñ$ Q¨Ó*‰	ˆˆrä×$Ñ$ Q¨°yÔA‰	ˆˆrÜ˜˜B ™FÔ#ä×$Ñ$ Q¨°vÔ>‰	ˆˆrÜ˜˜A  Q¡™JÔ'à˜Ô& Ò&Ð&Ñ&Ð&Ð&ä�]‰]œ:Ð-LÔMñ 	<Ü×Ñ˜Q ¨{Õ;÷	<÷ 	<ñ 	<ús   Å;FÆF&Úppc64lezfails/crashes on ppc64leru   c                 ór	  — t        j                  t        j                  g d¢«      «      }t        j                  t        j                  g d¢«      «      }ddg}t        t        j                  t        j                  ||«      «      |«       t        j                  t        j                  d«      «      }t        j                  t        j                  d«      «      }ddg}t        t        j                  t        j                  ||«      «      |«       t        t        t
        j                  ||d¬«       |d	   }|d
   |d	<   ||d
<   ddg}t        t        j                  t        j                  ||«      «      |«       |d   }|d   |d<   ||d<   ddg}t        t        j                  t        j                  ||«      «      |«       t        j                  t        j                  d«      «      }t        j                  t        j                  d«      d d d…   «      }ddg}t        t        j                  t        j                  ||«      «      |«       |d	   }|d
   |d	<   ||d
<   ddg}t        t        j                  t        j                  ||«      «      |«       |d   }|d   |d<   ||d<   ddg}t        t        j                  t        j                  ||«      «      |«       t        j                  ddddt        j                  g«      }t        j                  ddddt        j                  g«      }t        j                  ddddt        j                  g«      }t        t        j                  t        j                  ||«      «      dd g«       t        t        j                  t        j                  ||d!¬«      «      dd"g«       t        t        j                  t        j                  ||«      «      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%t        j                  g«      }t        j                  d%d)d)d)dd*d'd%d+d,d-d)d)d)d)d)d)d%dd.t        j                  d%g«      }t        j                  ||«      }t        t        j                  |«      d/d0g«       d1}t        ||d2¬3«       y )4N)r¾   rK   r™   rc   )rM   rš   r¾   é   r¬   r   rW   g´žx·O~¢>Úbanana©Úmethodr   rK   g?é“>é“î?gaÆV¥ã×>r™   rc   g}Ò'}Ò'í?go‡›…&5ÿ>rn   r«   g?é“>é“î¿g}Ò'}Ò'í¿rÇ   rÈ   rÉ   rÊ   rË   g     €:@g¸…ëQ¸Àg×£p=
×@rÌ   rÍ   g‡¼Š1UUÕ?ç      è?Ú
asymptoticgälgNÍß?g~öÈ‹ñ†á¿gÂQM¤�ÔÑ?r   rX   r[   rZ   r]   é!   éC   é   é   é-   g‘ŒæXJÄ¿gv™PwTkÚ?ræ   Tré   )re   r,   r+   r   r“   r-   Ú
kendalltaur­   Úassert_raisesrý   rë   r   r   )r/   rB   rÂ   ÚexpectedÚbÚzÚresultrð   s           r3   Útest_kendalltauzTestCorr.test_kendalltau1  s
  € ô �H‰H”R—X‘XšlÓ+Ó,ˆÜ�H‰H”R—X‘XšmÓ,Ó-ˆð ˜�:ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJô �H‰H”R—Y‘Y˜r“]Ó#ˆÜ�H‰H”R—Y‘Y˜r“]Ó#ˆð Ð+Ð,ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJô 	”j¤&×"3Ñ"3°Q¸À(ÕKð ˆa‰DˆØ�‰tˆˆ!‰Øˆˆ!‰ð 'Ð(:Ð;ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJð ˆa‰DˆØ�‰tˆˆ!‰Øˆˆ!‰ð 'Ð(:Ð;ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJô �H‰H”R—Y‘Y˜r“]Ó#ˆÜ�H‰H”R—Y‘Y˜r“]¡4 R 4Ñ(Ó)ˆð Ð,Ð-ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJð ˆa‰DˆØ�‰tˆˆ!‰Øˆˆ!‰ð (Ð);Ð<ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJð ˆa‰DˆØ�‰tˆˆ!‰Øˆˆ!‰ð (Ð);Ð<ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?ÀÔJô �N‰N˜D $¨¨d´B·F±FÐ;Ó<ˆÜ�N‰N˜D $¨¨t´R·V±VÐ<Ó=ˆÜ�N‰N˜D $¨¨d´B·F±FÐ;Ó<ˆÜœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?Ø'¨Ð.ô	0äœBŸJ™J¤v×'8Ñ'8¸¸AÀlÔ'SÓTØ'¨Ð3ô	5äœBŸJ™J¤v×'8Ñ'8¸¸AÓ'>Ó?Ø'¨Ð3ô	5ô �N‰N˜A˜q ! Q¨¨B°°2°q¸"Ø  A r¨1¨b°!°Q¸¸1¸aÄÇÁðIó Jˆä�N‰N˜A˜r 2 r¨2¨r°2°q¸"¸bØ  B¨¨B°°B¸¸2¸rÄ2Ç6Á6È1ðNó Oˆä×"Ñ" 1 aÓ(ˆÜœBŸJ™J vÓ.°¸YÐ0GÔHð /ˆ
Ü˜F J°4Ö8r5   c                 ó@  — t        j                  dt        ¬«      }t        j                  |d«      }t        j                  dt        ¬«      }t        j
                  |dd  |d d f«      }t        t        j                  t        j                  ||«      d   «      «       y )Nrå   r’   iË  éè  r   )
r+   r­   Úfloatre   Úmasked_greaterÚconcatenater   Úisfiniter-   r  )r/   rB   rÂ   s      r3   Útest_kendalltau_largezTestCorr.test_kendalltau_large‰  sx   € ô �I‰I�d¤%Ô(ˆÜ×Ñ˜a Ó&ˆÜ�I‰I�d¤%Ô(ˆÜ�N‰N˜A˜d˜e˜H a¨¨ hÐ/Ó0ˆÜ”—‘œF×-Ñ-¨a°Ó3°AÑ6Ó7Õ8r5   c                 óZ  — t         t         dddddddddddgg d¢dddd	d
ddddt         ddt         gt         d	dddt         d	ddddddgg}t        j                  |«      j                  }t	        j
                  |«      }t        |d   dd«       t        |d   j                  d«      g d¢«       y )NrM   rK   é   é   r™   r   rL   ©rM   rL   r™   rL   rK   rš   rL   r   r   rK   rL   r™   rL   rc   é   r¾   r	  é   zglobal p-value (indep)gü©ñÒMb€?zseasonal p-value)g
×£p=
Ç?gö(\�Âõà?çš™™™™™É?g{®Gáz¤?)r   re   rë   ÚTr-   Úkendalltau_seasonalr   Úround)r/   rB   Úoutputs      r3   Útest_kendalltau_seasonalz!TestCorr.test_kendalltau_seasonal•  s¼   € ä”3˜˜1˜b " a¨¨A¨q°!°Q¸Ð:Ú4Ø��A�q˜"˜a  A q¬#¨q°!´SÐ9Ü�1�b˜!˜R¤ a¨¨A¨q°!°Q¸Ð:ð<ˆô �N‰N˜1Ó×ÑˆÜ×+Ñ+¨AÓ.ˆÜ˜FÐ#;Ñ<¸eÀQÔGÜ˜FÐ#5Ñ6×<Ñ<¸QÓ?Ú1õ	3r5   r  )Úexactr  rõ   ©ú	two-sidedró   rö   c                 ó¶  — t         j                  j                  d«       d}t         j                  j                  |«      }t         j                  j                  |«      }t         j                  j                  |«      dkD  }t	        j
                  ||¬«      }t	        j
                  ||¬«      }t        j                  ||||¬«      }	|j                  «       }
|j                  «       }t        j                  |
|||¬«      }t         j                  ||<   t         j                  ||<   t        j                  |||d|¬«      }t        |	|«       t        ||«       y )Nr   r(   r�   ra   )r  rõ   Úomit)r  Ú
nan_policyrõ   )r+   rû   Úseedrü   re   r,   r-   r  Ú
compressedr   r   r   )r/   r  rõ   rÿ   rB   rÂ   rb   Úx_maskedÚy_maskedÚ
res_maskedÚx_compressedÚy_compressedÚres_compressedÚres_nans                 r3   Útest_kendalltau_mstats_vs_statsz(TestCorr.test_kendalltau_mstats_vs_stats¡  s  € ô 	�	‰	�‰�qÔØˆÜ�I‰I�N‰N˜1ÓˆÜ�I‰I�N‰N˜1ÓˆÜ�y‰y�~‰~˜aÓ  3Ñ&ˆä—8‘8˜A DÔ)ˆÜ—8‘8˜A DÔ)ˆÜ×&Ñ&Ø�h v¸;ôHˆ
ð  ×*Ñ*Ó,ˆØ×*Ñ*Ó,ˆÜ×)Ñ)Ø˜,¨vÀ;ôPˆô —&‘&ˆˆ$‰Ü—&‘&ˆˆ$‰Ü×"Ñ"Øˆq˜¨FÀôMˆô 	˜
 NÔ3Ü˜ Õ0r5   c           	      óš   — dddddddddœ}|j                  «       D ]-  \  }}t        j                  |d   |d	   «      }t        ||«       Œ/ y )
Ngé�£»mä?g‹Ž:¹4Wã?g,
ù:ÅH r¬   gäKªØëç?gù°¼(Â6Ú?))r_   iY	  )ée   i„	  )éª   r   )é«   r   )rA  r   )é¬   r   )éÈ   iE&  )éÉ   i¸%  r   r   ©Úitemsr   Ú_kendall_p_exactr   ©r/   ÚexpectationsÚncr  rï   s        r3   Útest_kendall_p_exact_mediumz$TestCorr.test_kendall_p_exact_mediumÀ  se   € ð &<Ø%;Ø"8Ø"%Ø"8Ø"%Ø%;Ø%;ñ=ˆð )×.Ñ.Ó0ò 	/‰LˆB�Ü×0Ñ0°°A±¸¸1¹Ó>ˆCÜ  XÕ.ñ	/r5   c                 ó–   — dddddddœ}|j                  «       D ]-  \  }}t        j                  |d   |d	   «      }t        ||«       Œ/ y )
NgEë‡ß?g§U`Ux<å?g½åveÅÛ?gÓÂZ¸šá?gôÐ=¹ñê?gG÷ ’ÉÖ?))i�  i5˜  )i‘  i\š  )i   idd )i!  i™h )i@  i º	 )iA  i ž	 r   r   rE  rH  s        r3   Útest_kendall_p_exact_largez#TestCorr.test_kendall_p_exact_largeÏ  s`   € ð '=Ø&<Ø'=Ø'=Ø(>Ø(>ñ@ˆð )×.Ñ.Ó0ò 	/‰LˆB�Ü×0Ñ0°°A±¸¸1¹Ó>ˆCÜ  XÕ.ñ	/r5   c                 óL  — g 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‘d‘t         j                  ‘}t        t        j                  ||«      d   d d!«       t        j                  ||«      }d"}t        ||d#¬$«       y )%N)#r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   rn   gš™™™™™-@gš™™™™™+@gÍÌÌÌÌÌ(@g333333$@gffffff@gffffff@g333333@gffffff@g333333@r¡   gffffff
@gš™™™™™	@r    gffffff@r'   g333333@gffffff@gÍÌÌÌÌÌ @rÒ   g      ø?gÍÌÌÌÌÌô?rä   gš™™™™™ñ?çš™™™™™é?gffffffæ?rà   r�   r(  r÷   r   gÂ/õó¦"×?r™   ræ   Tré   )r+   r   r   r-   Úpointbiserialrr   rî   s        r3   Útest_pointbiserialzTestCorr.test_pointbiserialÞ  s‰  € ò
5ˆðBˆTð B�4ð B˜ð B˜tð B Sð B¨#ð B¨sð B°Cð B¸ð B¸cð BÀ3ð BÈð BØðBØðBØðBØðBØ!$ðBØ&)ðBØ+.ðBØ03ðBØ58ðBØ:=ðBØ?BðBØDGðBØILðBàðBàðBàðBàðBà!$ðBà&)ðBà+.ðBà03ðBà58ðBä:<¿&¹&ðBˆô 	œF×1Ñ1°!°QÓ7¸Ñ:¸GÀQÔGô ×#Ñ# A qÓ)ˆØ.ˆ
Ü˜C °Ö5r5   N)r6   r7   r8   r»   rÅ   rñ   r  rˆ   r‰   rŠ   ÚplatformÚmachiner  Úslowr!  r-  Úparametrizer=  rK  ÚxslowrM  r   rQ  r9   r5   r3   r¨   r¨   ´   s'  „ ò&ò>ò$6òL+<ðZ ‡[�[×ÑÐ(˜×(Ñ(Ó*¨iÑ7Ø9ð ó ;ñT9ó;ðT9ðl ‡[�[×ÑÐ(˜×(Ñ(Ó*¨iÑ7Ø9ð ó ;à‡[�[×Ññ9ó ó;ð9ò
3ð ‡[�[×Ñ˜XÐ'>Ó?Ø‡[�[×Ñ˜]Ð,LÓMñ1ó Nó @ð1ò:/ð ‡[�[×Ññ/ó ð/ð ñ6ó	 ñ6r5   r¨   c                   óB   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	d„ Z
d	„ Zy
)ÚTestTrimmingc                 ó^  — t        j                  d«      }t        t        j                  |«      g d¢«       t        j                  d«      }t        t        j                  |d«      g d¢«       t        j                  d«      }t        t        j                  |dd¬«      g d¢«       t        j                  d«      }t        t        j                  |dd	¬
«      g d¢«       t        j                  d«      }t
        x|ddg<   |d<   t        t        j                  |d«      g d¢«       t        j                  d«      j                  dd«      }dgdz  dgdz  z   dgdz  z   }t        j                  |dd	d ¬«      }t        |j                  j                  «       |«       t        j                  |dd	d¬«      }t        |j                  j                  «       |«       t        j                  |dd	d¬«      }t        |j                  j                  j                  «       |«       t        j                  d«      j                  dd«      }t
        |d<   dgdz  dgdz  z   dgdz  z   }t        j                  |dd	d ¬«      }t        |j                  j                  «       |«       t        j                  |dd	d¬«      }t        |j                  j                  «       |«       t        j                  |j                  dd	d¬«      }t        |j                  j                  j                  «       |«       y )NrW   ©
r   r   rK   rL   rM   r™   rc   rš   r�   r¾   )rK   r�   )
NNrK   rL   rM   r™   rc   rš   r�   N©FF©ÚlimitsÚ	inclusive)
NNNrL   rM   r™   rc   rš   NN)r÷   r(  T)r]  Úrelative)
Nr   rK   rL   rM   r™   rc   rš   NNé   r   rn   r™   )NNrK   rL   rM   Nrc   rš   r�   NNNr_   r   r\   rX   )r_  r)   r…   r	  )
re   r­   r
   r-   Útrimr   ÚreshapeÚ_maskÚravelr)  )r/   rS   rB   r  Útrimxs        r3   Ú	test_trimzTestTrimming.test_trimó  s`  € Ü�I‰I�b‹MˆÜ”V—[‘[ “^Ò%:Ô;Ü�I‰I�b‹MˆÜ”V—[‘[  5Ó)Ò+IÔJÜ�I‰I�b‹MˆÜ”V—[‘[ ¨%¸-ÔHÚ9ô	;ä�I‰I�b‹MˆÜ”V—[‘[ ¨)¸TÔBÚ3ô	5ô �I‰I�b‹MˆÜ!Ð!ˆˆ1ˆRˆ&‰	�A�a‘DÜ”V—[‘[  EÓ*ÚKô	Mô �I‰I�c‹N×"Ñ" 2 rÓ*ˆØ�3�r‘6˜Q˜C ™F‘? a S¨¡VÑ+ˆÜ—‘˜A˜y°4¸dÔCˆÜ�U—[‘[×&Ñ&Ó(¨(Ô3Ü—‘˜A˜y°4¸aÔ@ˆÜ�U—[‘[×&Ñ&Ó(¨(Ô3Ü—‘˜A˜y°4¸bÔAˆÜ�U—[‘[—]‘]×(Ñ(Ó*¨HÔ5ô �I‰I�c‹N×"Ñ" 2 rÓ*ˆÜˆˆ!‰Ø�3�r‘6˜Q˜C ™F‘? a S¨¡VÑ+ˆÜ—‘˜A˜y°4¸dÔCˆÜ�U—[‘[×&Ñ&Ó(¨(Ô3Ü—‘˜A˜y°4¸aÔ@ˆÜ�U—[‘[×&Ñ&Ó(¨(Ô3Ü—‘˜AŸC™C °TÀÔCˆÜ�U—W‘W—]‘]×(Ñ(Ó*¨HÕ5r5   c                 óš  — t        j                  d«      }t        t        j                  |«      j                  «       d«       t        t        j                  |d¬«      j                  «       d«       t        |dd t        j                  |«      }t        |j                  «       d«       t        |j                  d	gd
z  dgdz  z   d	gdz  z   dgdz  z   d	gd
z  z   «       t        |_        d|_
        t        t        j                  |«      j                  «       d«       t        t        j                  |«      j                  «       d«       y )Nr_   r[   r·   )Útailr]   r(   r\   é0   r   r#  r   é"   rX   é   )rW   rW   )re   r­   r
   r-   ÚtrimbothÚcountÚtrimtailr   rc  r   Úshape)r/   rB   re  s      r3   Útest_trim_oldzTestTrimming.test_trim_old  s÷   € Ü�I‰I�c‹NˆÜ”V—_‘_ QÓ'×-Ñ-Ó/°Ô4Ü”V—_‘_ Q¨CÔ0×6Ñ6Ó8¸"Ô=Üˆˆ"ˆRˆÜ—‘ Ó"ˆÜ�U—[‘[“] BÔ'Ü�U—[‘[ 1 # b¡&¨A¨3¨r©6¡/°Q°C¸±FÑ":¸a¸SÀ¹VÑ"CÀqÀcÈ"ÁfÑ"LÔMÜˆŒØˆŒÜ”V—_‘_ QÓ'×-Ñ-Ó/°Ô4Ü”V—_‘_ QÓ'×-Ñ-Ó/°Õ4r5   c                 óì   — t        j                  d«      }t        j                  |dd¬«      }t        j                  g d¢g d¢¬«      }t        ||«       t        |j                  |j                  «       y )NrW   )g333333Ã?gìQ¸…ëÁ?r[  r\  rZ  )
r   r   r   r   r   r   r   r   r   r   ra   )re   r­   r-   Útrimrr,   r
   rb   )r/   rB   r  r  s       r3   Ú
test_trimrzTestTrimming.test_trimr%  sR   € Ü�I‰I�b‹MˆÜ—‘˜a¨ÀÔOˆÜ—8‘8Ò:Ú!?ôAˆä�V˜XÔ&Ü�V—[‘[ (§-¡-Õ0r5   c                 óø   — t        j                  g d¢«      }t        t        j                  |d«      dd«       t        t        j                  |d«      dd«       t        t        j                  |d«      dd«       y )N©éM   éW   éX   ér   é—   éÒ   éÛ   éö   éý   i  i(  i+  i2  ix  i¬  i  iš  i  i3
  r÷   iW  r   )r÷   r÷   ©r(  r(  i  )re   r,   r   r-   Útrimmed_mean©r/   r0   s     r3   Útest_trimmedmeanzTestTrimming.test_trimmedmean-  sb   € Ü�x‰xò @ó AˆäœF×/Ñ/°°SÓ9¸3ÀÔBÜœF×/Ñ/°°YÓ?ÀÀaÔHÜœF×/Ñ/°°YÓ?ÀÀaÕHr5   c                 ó  — t         j                  j                  d«      }|j                  d¬«      }t        j                  |«      }t	        j
                  |g d¢¬«      }t        t        j                  |d«      |j                  «       «       y )Nl   l{ýf0j rX   ©Úsize©r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   ra   r÷   )
r+   rû   Údefault_rngÚsortre   r,   r   r-   Útrimmed_varÚvar©r/   ÚrngÚ	data_origr0   s       r3   Útest_trimmedvarzTestTrimming.test_trimmedvar4  ók   € ô �i‰i×#Ñ#Ð$7Ó8ˆØ—J‘J B�JÓ'ˆ	Ü�w‰w�yÓ!ˆÜ�x‰x˜ò $Bô Cˆäœ×*Ñ*¨9°cÓ:¸D¿H¹H»JÕGr5   c                 ó  — t         j                  j                  d«      }|j                  d¬«      }t        j                  |«      }t	        j
                  |g d¢¬«      }t        t        j                  |d«      |j                  «       «       y )Nl   œ/{Gzi— rX   r„  r†  ra   r÷   )
r+   rû   r‡  rˆ  re   r,   r   r-   Útrimmed_stdÚstdr‹  s       r3   Útest_trimmedstdzTestTrimming.test_trimmedstd>  r�  r5   c                 ó¶   — t        j                  g d¢«      }t        t        j                  |d«      dd«       t        t        j                  |d«      dd«       y )Nru  r  g>®ãL@r™   r(  )re   r,   r   r-   Útrimmed_stder�  s     r3   Útest_trimmed_stdezTestTrimming.test_trimmed_stdeH  sI   € Ü�x‰xò @ó AˆäœF×/Ñ/°°YÓ?ÀÈ1ÔMÜœF×/Ñ/°°SÓ9¸8ÀQÕGr5   c                 ót  — t        j                  g d¢«      }t        t        j                  |d«      j                  d¬«      dd«       t        t        j                  |dd«      j                  d¬«      dd«       t        |d<   t        j                  |«      }t        |j                  |j                  «       y )	Nru  r  r   ©Úddofgš™™™ÙÕ@r[  gffff¦7Ç@r™   )	re   r,   r   r-   Ú	winsorizerŠ  r   r
   rb   )r/   r0   Ú
winsorizeds      r3   Útest_winsorizationzTestTrimming.test_winsorizationN  s›   € Ü�x‰xò @ó AˆäœF×,Ñ,¨T°)Ó<×@Ñ@ÀaÐ@ÓHØ# Qô	(äÜ×Ñ˜T 9¨]Ó;×?Ñ?ÀQÐ?ÓGØ�Qô	ô ˆˆQ‰Ü×%Ñ% dÓ+ˆ
Ü�Z—_‘_ d§i¡iÕ0r5   c           
      ór  — t        j                  t        j                  t        j                  dddg«      }t	        t
        t        j                  |dd¬«       t        t        j                  |d«      t        j                  g d¢«      «       t        t        j                  |d	«      t        j                  t        j                  t        j                  t        j                  t        j                  t        j                  g«      «       t        t        j                  |dd
¬«      t        j                  t        j                  t        j                  dddg«      «       t        t        j                  |d	d
¬«      t        j                  t        j                  t        j                  dddg«      «       y )Nr   r   rK   )çš™™™™™©?rž  Úraise)r3  )çš™™™™™Ù?r   r£   )rO  rO  r2  )	re   r,   r+   r   r  rý   r-   rš  r
   r�  s     r3   Útest_winsorization_nanz#TestTrimming.test_winsorization_nanZ  s  € Ü�x‰xœŸ™¤§¡¨¨A¨qÐ1Ó2ˆÜ”j¤&×"2Ñ"2°D¸,Ø!(õ	*ô 	”V×%Ñ% d¨JÓ7Ü—X‘XšoÓ.ô	0ä”V×%Ñ% d¨JÓ7Ü—X‘XœrŸv™v¤r§v¡v¬r¯v©v´r·v±v¼r¿v¹vÐFÓGô	Iä”V×%Ñ% d¨JÀ6ÔJÜ—X‘XœrŸv™v¤r§v¡v¨q°!°QÐ7Ó8ô	:ä”V×%Ñ% d¨JÀ6ÔJÜ—X‘XœrŸv™v¤r§v¡v¨q°!°QÐ7Ó8õ	:r5   N)r6   r7   r8   rf  rp  rs  r‚  rŽ  r“  r–  rœ  r¡  r9   r5   r3   rX  rX  ð  s4   „ ò#6òJ5ò1òIòHòHòHò
1ó:r5   rX  c                   ó,  — e Zd Zg d¢Z ej
                  dddddej                  g«      Z ej                   ej                  g d¢g d¢g d	¢g d
¢g d¢g«       ej                  g d¢g d¢g d¢g d¢g d¢ge
¬«      ¬«      Zdddœd„Zd„ Zd„ Zd„ Zd„ Zd„ Zy)ÚTestMomentsrJ   g¤p=
×£ò?gæ?¤ß¾ä?gDúíëÀ9³?gª‚QI�€Ö?g}?5^ºIæ¿)gÂq2ÔÕÛª?gþ{•!3à?gr!faµë?gìØkî¦£×?gžØYÞ?)gSH*:/¡½?g#RÚ/Ô?g>yõ}zdÝ?gVß¯¸î<ç?g¼F¢íè?)ggË w¶©å?gžyZú§fí?g½�îœ=¹?g!@VyÈ§?g# 1öØ?)geY(€ÿ9®?g—£�D$Ó?gÁ;·ØÑÕ?gðJX<q0å?g#dˆŒœê?)gŠ÷�®¡­ä?g×çA¬Ä]î?g_.ÚÊèÓÑ?g¥�ÉÔ²?g2^ë9†“Ù?)TFFTF)TTTFT)FFFFF)TTTTT)FFTFFr’   ra   N©ro  r<   c                ó´   — t        j                  |«      }|�t        j                  ||«      }t        ||«       |€|j                  }|j                  |k(  sJ ‚y ©N)r+   r“   Úbroadcast_tor   r<   )r/   ÚactualÚexpectro  r<   s        r3   Ú_assert_equalzTestMoments._assert_equal�  sP   € Ü—‘˜FÓ#ˆØÐÜ—_‘_ V¨UÓ3ˆFÜ˜6 6Ô*Øˆ=Ø—L‘LˆEØ�|‰|˜uÒ$Ð$Ñ$r5   c                 ób  — t        j                  | j                  d«      }t        |dd«       t        j                  | j                  d«      }t        |d«       t        j                  | j                  d«      }t        |d«       t        j                  | j                  d«      }t        |d«       t        j                  | j                  g d	¢«      }t	        |g d
¢«       t        j                  | j                  d«      }t	        |d«       t        t        t         j                  | j                  d«       t        j                  | j                  g d¢«      }t	        |g d
¢«       t        j                  g «      }| j                  |t        j                  t        j                  ¬«       t        j                  t        j                  g t        j                  ¬«      «      }| j                  |t        j                  t        j                  ¬«       t        j                  t        j                  d«      d¬«      }| j                  |g dt        j                  ¬«       t        j                  g gd¬«      }| j                  |t        j                  dt        j                  ¬«       t        j                  g gddgd¬«      }| j                  |g d¬«       t        j                  d«      }t        j                  |d<   t        t        j                  |d«      t         j"                  «       y )Nr   r¬   rW   rK   ç      ô?rL   rM   ç     €@rJ   )r   r¬  r   r­  r   rä   )r   rK   rL   r"   r’   )r   r   r   r    )r   r¤  )r   )Úmomentr)   )rK   r   )ro  g      $@r¾   )r-   r®  Útestcaser   r   r  rý   rª  r+   r   Úfloat64r,   Úfloat32Úzerosr­   r
   re   r   )r/   rÂ   rB   s      r3   Útest_momentzTestMoments.test_momentŠ  s  € Ü�M‰M˜$Ÿ-™-¨Ó*ˆÜ˜A˜c "Ô%Ü�M‰M˜$Ÿ-™-¨Ó*ˆÜ˜A˜dÔ#Ü�M‰M˜$Ÿ-™-¨Ó*ˆÜ˜A˜cÔ"Ü�M‰M˜$Ÿ-™-¨Ó*ˆÜ˜A˜fÔ%ô �M‰M˜$Ÿ-™-ªÓ6ˆÜ˜Ò/Ô0ô �M‰M˜$Ÿ-™-¨Ó-ˆÜ˜˜3ÔÜ”j¤&§-¡-°·±ÀÔDÜ�M‰M˜$Ÿ-™-Ò)9Ó:ˆÜ˜Ò/Ô0ô �M‰M˜"ÓˆØ×Ñ˜1œbŸf™f¬B¯J©JÐÔ7Ü�M‰Mœ"Ÿ(™( 2¬R¯Z©ZÔ8Ó9ˆØ×Ñ˜1œbŸf™f¬B¯J©JÐÔ7Ü�M‰Mœ"Ÿ(™( 6Ó*°Ô3ˆØ×Ñ˜1˜b¨´B·J±JÐÔ?Ü�M‰M˜2˜$ QÔ'ˆØ×Ñ˜1œbŸf™f¨D¼¿
¹
ÐÔCÜ�M‰M˜2˜$¨¨1 v°AÔ6ˆØ×Ñ˜1˜b¨ÐÔ/ä�I‰I�c‹NˆÜ�v‰vˆˆ!‰Ü”V—]‘] 1 aÓ(¬"¯)©)Õ4r5   c                 ó\   — t        j                  | j                  «      }t        |dd«       y )Nçþí¿Å%ŸÜ?rW   )r-   Ú	variationr¯  r   )r/   rÂ   s     r3   Útest_variationzTestMoments.test_variation¯  s#   € Ü×Ñ˜TŸ]™]Ó+ˆÜ˜AÐ.°Õ3r5   c                 óx   — t        j                  g d¢«      }t        j                  |d¬«      }t	        |d«       y )N)r   rK   rL   rM   r™   r   r˜  gâ<<'�Ýà?)r+   r,   r-   r¶  r   )r/   rS   rÂ   s      r3   Útest_variation_ddofzTestMoments.test_variation_ddof³  s/   € ô �H‰H’_Ó%ˆÜ×Ñ˜Q QÔ'ˆÜ˜AÐ1Õ2r5   c           	      ó€  — t        j                  | j                  «      }t        |dd«       t        j                  | j                  d¬«      }t        |dd«       t        j                  | j                  «      }t        |dd«       t        j                  t        j                  g d¢«      t        j                  g d¢t        ¬	«      ¬
«      }t        t        j                  | j                  d«      |«       t        | j                  «      D ]'  \  }}t        t        j                  |«      ||   «       Œ) t        j                  t        j                  g d¢«      t        j                  g d¢t        ¬	«      ¬
«      }t        t        j                  | j                  dd¬«      |«       t        | j                  «      D ])  \  }}t        t        j                  |d¬«      ||   «       Œ+ t        t        j                  | j                  dd d …f   «      t        j                  | j                  dd d …f   «      «       y )Ng7l•*ÄÒ¿rW   r   ©Úbiasg¨õ2Û ùÛ¿r¬   )gÒ½•Šræ?r   g>HÃ›eÑ?r   gIåÝÇ®ª¿©FFFTFr’   ra   r   )g`ŽƒÙ¨ùú?r¬   gÂûé‡nÙ?r   g0¡·¿FrK   )r-   ÚskewÚtestmathworksr   r¯  re   r,   r+   Úboolr   Útestcase_2dÚ	enumerater   ©r/   rÂ   Ú
correct_2dÚiÚrowÚcorrect_2d_bias_correcteds         r3   Útest_skewnesszTestMoments.test_skewnessº  sµ  € Ü�K‰K˜×*Ñ*Ó+ˆÜ˜AÐ/°Ô3Ü�K‰K˜×*Ñ*°Ô2ˆÜ˜AÐ0°Ô4Ü�K‰K˜Ÿ™Ó&ˆÜ˜A˜c "Ô%ô —X‘XÜ�H‰Hò /ó 0ä—‘Ò<ÄDÔIô
ˆ
ô
 	œŸ™ D×$4Ñ$4°aÓ8¸*ÔEÜ × 0Ñ 0Ó1ò 	A‰FˆAˆsÜ¤§¡¨CÓ 0°*¸Q±-Õ@ð	Aô %'§H¡HÜ�H‰Hò ,ó -ä—‘Ò<ÄDÔIô%
Ð!ô
 	œŸ™ D×$4Ñ$4°a¸eÔDØ1ô	3ä × 0Ñ 0Ó1ò 	>‰FˆAˆsÜ¤§¡¨C°eÔ <Ø 9¸!Ñ <õ>ð	>ô
 	œŸ™ D×$4Ñ$4°Qº°TÑ$:Ó;ÜŸ
™
 4×#3Ñ#3°A²q°DÑ#9Ó:õ	<r5   c           	      ó�  — t        j                  | j                  ddd¬«      }t        |dd«       t        j                  | j                  dd¬«      }t        |dd«       t        j                  | j                  dd«      }t        |d«       t        j                  t        j                  g d¢«      t        j                  g d	¢t        ¬
«      ¬«      }t        t        j                  | j                  d«      |«       t        | j                  «      D ]'  \  }}t        t        j                  |«      ||   «       Œ) t        j                  t        j                  g d¢«      t        j                  g d	¢t        ¬
«      ¬«      }t        t        j                  | j                  dd¬«      |«       t        | j                  «      D ])  \  }}t        t        j                  |d¬«      ||   «       Œ+ t        t        j                  | j                  dd d …f   «      t        j                  | j                  dd d …f   «      d¬«       y )Nr   r   )Úfisherr¼  gOß»S@rW   gâx|ïN@g=
×£p=ú?)ç      ø¿ç      Àgã'„@=�÷¿r¬   gt‚
ÀQô¿r½  r’   ra   )rË  rÌ  g´?õ<þ¿r¬   g`‡œ@7Àà¿Fr»  rK   rM   )Únulp)r-   Úkurtosisr¿  r   r¯  re   r,   r+   rÀ  r   rÁ  rÂ  r   r   rÃ  s         r3   Útest_kurtosiszTestMoments.test_kurtosisÛ  sÐ  € ô �O‰O˜D×.Ñ.°¸!À!ÔDˆÜ˜A˜°Ô3ô
 �O‰O˜D×.Ñ.°q¸qÔAˆÜ˜AÐ0°"Ô5Ü�O‰O˜DŸM™M¨1¨aÓ0ˆÜ˜A˜tÔ$ô —X‘XœbŸh™hò (8ó 9ä#%§8¡8ò -4Ü;?ô$AôBˆ
ô 	"¤&§/¡/°$×2BÑ2BÀAÓ"FØ",ô	.ä × 0Ñ 0Ó1ò 	E‰FˆAˆsÜ¤§¡°Ó 4°jÀ±mÕDð	Eô %'§H¡HÜ�H‰HÒIÓJÜ—‘Ò<ÄDÔIô%KÐ!ô 	"¤&§/¡/°$×2BÑ2BÀAØ7<ô#>à";ô	=ô   × 0Ñ 0Ó1ò 	>‰FˆAˆsÜ¤§¡°¸%Ô @Ø 9¸!Ñ <õ>ð	>ô
 	'¤v§¡°t×7GÑ7GÈÊ1ÈÑ7MÓ'NÜ',§~¡~°d×6FÑ6FÀqÊ!ÀtÑ6LÓ'MØ,-ö	/r5   )r6   r7   r8   r¯  re   rë   r+   r   r¿  r,   rÀ  rÁ  rª  r³  r·  r¹  rÈ  rÏ  r9   r5   r3   r£  r£  i  s·   „ ò €HØ"�B—N‘N E¨6°6¸6À7Ø$&§F¡Fð$,ó -€Mà�"—(‘(Øˆ�‰ÒNÚNÚNÚNÚMð	Oó 	Pð
 ˆR�X‰XÒ8Ú6Ú:Ú5Ú9ð	;ð CGô	Hô
I€Kð 6:Àô %ò#5òJ4ò3ò<óB%/r5   r£  c                   ó   — e Zd Zd„ Zd„ Zy)ÚTestModec                 óÒ  — g d¢}t        j                  |d«      }t        j                  g d¢«      }t        j                  |d«      }t        j                  t        j                  |«      dkD  |«      }t        j                  |dkD  |«      }t        j                  |dk  |«      }t        j                  t        j                  |«      dk  |«      }t        t        j                  |d ¬«      d«       t        t        j                  |d¬«      d«       t        t        j                  |d ¬«      d	«       t        t        j                  |d ¬«      d«       t        t        j                  |d ¬«      d	«       t        t        j                  |d ¬«      d
«       t        t        j                  |d ¬«      d«       t        t        j                  |d¬«      g d¢gg d¢gf«       t        t        j                  |d¬«      g d¢gg d¢gf«       t        t        j                  |d¬«      dgdgdggdgdgdggf«       t        t        j                  |d¬«      dgdgdggdgdgdggf«       t        t        j                  |d¬«      ddggddggf«       t        t        j                  |d¬«      dgdgdggdgdgdggf«       t        j                  |d ¬«      }	d}
t        |	|
d¬«       y )N)r   r   r   r   r   r   rK   rL   rL   rL   rL   rM   r™   rc   rš   )rL   r™   r½   )rL   rK   rK   r    )rL   rM   r   )r   rL   )r   r   )rK   r   )r   r   r   r   r   r¢   rn   rL   r   r™   )Úmoderm  Tré   )	r+   rb  r,   re   Úmasked_wherer
   r-   rÓ  r   )r/   Úa1Úa2Úa3Úa4Úma1Úma2Úma3Úma4Úa1_resrð   s              r3   Ú	test_modezTestMode.test_mode  sE  € Ú,ˆÜ�Z‰Z˜˜EÓ"ˆÜ�X‰X’mÓ$ˆÜ�Z‰Z˜˜EÓ"ˆÜ�o‰oœbŸh™h r›l¨QÑ.°Ó3ˆÜ�o‰o˜b 1™f bÓ)ˆÜ�o‰o˜b 1™f bÓ)ˆÜ�o‰oœbŸh™h r›l¨QÑ.°Ó3ˆÜ”V—[‘[ ¨$Ô/°Ô7Ü”V—[‘[ ¨!Ô,¨eÔ4Ü”V—[‘[ ¨4Ô0°%Ô8Ü”V—[‘[ ¨$Ô/°Ô7Ü”V—[‘[ ¨4Ô0°%Ô8Ü”V—[‘[ ¨$Ô/°Ô7Ü”V—[‘[ ¨4Ô0°%Ô8Ü”V—[‘[ ¨!Ô,²¨}º{¸mÐ.LÔMÜ”V—[‘[ ¨1Ô-²°Â¸}Ð/MÔNÜ”V—[‘[ ¨"Ô-°!°°a°S¸!¸°ÀÀÀQÀCÈÈ¸}Ð/MÔNÜ”V—[‘[ ¨2Ô.°1°#°q°c¸1¸#°À!ÀÀaÀSÈ!ÈÀÐ0NÔOÜ”V—[‘[ ¨1Ô-°!°A°°¸1¸Q¸%¸Ð/AÔBÜ”V—[‘[ ¨2Ô.°1°#°q°c¸1¸#°À!ÀÀaÀSÈ!ÈÀÐ0NÔOä—‘˜R dÔ+ˆð 'ˆ
Ü˜F J°4Ö8r5   c                 óÞ   — t        j                  d«      }|d d…d d …fxx   dz  cc<   |d d …d d…fxx   dz  cc<   |j                  «       }t        j                  |d «       t        ||«       y )N)r_   r_   r(   r   )r+   r²  Úcopyr-   rÓ  r
   )r/   ÚimÚcps      r3   Útest_mode_modifies_inputz!TestMode.test_mode_modifies_input!  s_   € ô �X‰X�jÓ!ˆØ
ˆ3ˆBˆ3’ˆ6‹
�a‰‹
Ø
Š1ˆcˆrˆcˆ6‹
�a‰‹
Ø�W‰W‹YˆÜ�‰�B˜ÔÜ�R˜Õr5   N)r6   r7   r8   rÞ  rã  r9   r5   r3   rÑ  rÑ    s   „ ò9ó:r5   rÑ  c                   ó   — e Zd Zd„ Zd„ Zd„ Zy)ÚTestPercentilec                 ó:   — g d¢| _         g d¢| _        g d¢| _        y )N)rL   rM   r™   rW   éýÿÿÿéûÿÿÿrc   )rL   éúÿÿÿéþÿÿÿr�   rš   rM   rK   r   )r    rM   r™   rW   rç  rè  ré  r   )rÕ  rÖ  r×  ©r/   s    r3   Úsetup_methodzTestPercentile.setup_method-  s   € Ú*ˆŒÚ,ˆŒÚ1ˆ�r5   c                 óô   — t        j                  d«      dz  }t        t        j                  |d«      d«       t        t        j                  |d«      d«       t        t        j                  |d«      d«       y )	Nr�   r�   r   r¬   r_   r¡   r(   g      ü?)r+   r­   r
   r-   Úscoreatpercentiler¥   s     r3   Útest_percentilezTestPercentile.test_percentile2  sY   € Ü�I‰I�a‹L˜3ÑˆÜ”V×-Ñ-¨a°Ó3°RÔ8Ü”V×-Ñ-¨a°Ó5°sÔ;Ü”V×-Ñ-¨a°Ó4°dÕ;r5   c                 ó�   — t        j                  g d¢g d¢g d¢g d¢g d¢g«      }t        t        j                  |d«      g d¢«       y )N©r   r   r   )rM   rM   rL   r(   )re   r,   r
   r-   rî  r¥   s     r3   Útest_2DzTestPercentile.test_2D8  s<   € Ü�H‰H’iÚÚÚÚð	!ó "ˆô
 	”V×-Ñ-¨a°Ó4²iÕ@r5   N)r6   r7   r8   rì  rï  rò  r9   r5   r3   rå  rå  ,  s   „ ò2ò
<óAr5   rå  c                   ód   — e Zd ZdZ ej
                  ddddej                  g«      Zd„ Z	d„ Z
d„ Zy	)
ÚTestVariabilityz[  Comparison numbers are found using R v.1.5.1
         note that length(testcase) = 4
    r   rK   rL   rM   c                 óZ  — t        j                  | j                  «      }t        |d«       | j                  j	                  «       }t        t        j                  | j                  d¬«      t        j                  ||dz
  z  «      z  t        j                  | j                  d¬«      «       y )Ngã�ŽËé§ä?r   r˜  rK   )r-   Úsemr¯  r   rm  r   r+   r³   )r/   rÂ   rÿ   s      r3   Útest_semzTestVariability.test_semH  sq   € ä�J‰J�t—}‘}Ó%ˆÜ˜A˜|Ô,Ø�M‰M×ÑÓ!ˆÜœŸ
™
 4§=¡=°qÔ9¼B¿G¹GÀAÀqÈÁsÁGÓ<LÑLÜŸ
™
 4§=¡=°qÔ9õ	;r5   c                 ó°   — t        j                  | j                  | j                  «      }g d¢}t        ||j                  |j
                  dk(     d¬«       y )N)çÙòOT\wõ¿çþí¿Å%ŸÜ¿rµ  çÙòOT\wõ?Fr`  ©Údecimal)r-   Úzmapr¯  r   r0   rb   )r/   rÂ   Údesired_unmaskedvalss      r3   Ú	test_zmapzTestVariability.test_zmapP  sG   € ô �K‰K˜Ÿ™ t§}¡}Ó5ˆò!DÐä!Ð"6Ø"#§&¡&¨¯©°5©Ñ"9À2ö	Gr5   c                 ó®   — t        j                  | j                  «      }t        j                  ddddt
        j                  g«      }t        ||d¬«       y )Nrù  rú  rµ  rû  r`  rü  )r-   Úzscorer¯  re   rë   r+   r   r   )r/   rÂ   r1   s      r3   Útest_zscorezTestVariability.test_zscoreY  sH   € ô �M‰M˜$Ÿ-™-Ó(ˆÜ—.‘.Ð"2Ð4EØ"2°OÄRÇVÁVð"Mó Nˆä˜G Q°Ö3r5   N)r6   r7   r8   Ú__doc__re   rë   r+   r   r¯  r÷  r   r  r9   r5   r3   rô  rô  A  s8   „ ñð ˆr�~‰~˜q  1 Q r§v¡vÐ.Ó/€Hò;òGó4r5   rô  c                   ó   — e Zd Zd„ Zd„ Zd„ Zy)ÚTestMiscc                 óš  — dgdz  dgdz  z   dgdz  z   dgdz  z   dgdz  z   d	gdz  z   dgdgdz  z   dgd
z  z   dgdz  z   d	gdz  z   g}ddgz  ddgz  z   ddgz  z   ddgz  z   ddgz  z   ddgz  z   dgddgz  z   d
dgz  z   ddgz  z   ddgz  z   g}t        t        j                  t        j                  |Ž j
                  d
«      |d
«       y )Nr™   rc   r	  rš   r¾   r�   rL   rK   rW   rM   r#  gŸ«­Ø_v	@gù1æ®%äá?g"ýöuàœ¡?g”‡…ZÓ¼ù?gÿ!ýöu @g³q¬‹&@g;pÎˆÒÞ$@gTR' ‰p@gúíëÀ9#ö?gðHPüx?g“©‚QIç?)r   r+   r+  r-   Úobrientransformr)  )r/   Úargsr  s      r3   Útest_obrientransformzTestMisc.test_obrientransforme  s  € Ø��A‘�q�c˜"‘f‘˜a˜S ™UÑ" A 3 q¡5Ñ(¨!¨¨Q©Ñ.°¨t°A©vÑ5Ø��Q�C˜‘E‘	˜1˜#˜a™%‘   A¡Ñ% r d¨2¡gÑ-ð/ˆà�V�H‘*˜R  ™[Ñ(¨¨F¨8©Ñ3°A°v°h±JÑ>¸qÀ&À¹zÑIÈ!ÈWÈIÉ+ÑUØ�)˜A˜v˜h™JÑ& q¨&¨¡zÑ1°!°V°H±*Ñ<¸RÀÀ¹[ÑHðJˆäœBŸH™H¤V×%;Ñ%;¸TÐ%B×%DÑ%DÀaÓHØ" Aõ	'r5   c                 ó.  — t         t         dddddddddddgg d¢dddd	d
ddddt         ddt         gt         d	dddt         d	ddddddgg}t        j                  |«      j                  }|j                  \  }}}}t	        t        j                  t        j                  ||«      d«      d«       t	        t        j                  t        j                  ||d«      d«      d«       t	        t        j                  t        j                  ||d«      d«      d«       y )NrM   rK   r#  r$  r™   r   rL   r%  rc   r&  r¾   r	  r'  )ç¡Ö4ï8EÇ?gs×òAÏî?Úg)g¯%äƒžÍÂ?g·Ñ Þ	æ?Úl)r  gê46<ã?)	r   re   rë   r)  r   r+   r+  r-   Úks_2samp)r/   rB   ÚwinterÚspringÚsummerÚfalls         r3   Útest_ks_2sampzTestMisc.test_ks_2sampm  sû   € Ü”#�q˜!˜R  Q¨¨1¨a°°A°qÐ9Ú4Ø��A�q˜"˜a  A q¬#¨q°!´SÐ9Ü�1�b˜!˜R¤ a¨¨A¨q°!°Q¸Ð:ð<ˆô �N‰N˜1Ó×ÑˆØ)*¯©Ñ&ˆ�˜ äœBŸH™H¤V§_¡_°V¸VÓ%DÀaÓHØ,ô	.äœBŸH™H¤V§_¡_°V¸VÀSÓ%IÈ1ÓMØ,ô	.äœBŸH™H¤V§_¡_°V¸VÀSÓ%IÈ1ÓMØ,õ	.r5   c                 óº  — g d¢g d¢g d¢f}t        j                  |Ž }t        |d   dd«       t        |d   dd	«       t        t        dd
ddddddd
ddgg d¢dd
dd	dddddt        ddt        gt        d	dddt        d	ddd
dddgg}t	        j
                  |«      }t        j                  |Ž }t        |d   dd«       t        |d   dd«       d}t        ||d¬«       y )N)
g      "@ç      #@ç      @ç      @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   gû:pÎˆò$@rM   r   gí(ÎQGÇu?rc   rK   r#  r$  r™   rL   r%  r&  r¾   r	  r'  grŠŽäò @g-Cëâ6â?©Ú	statisticrè   Tré   )r-   Úfriedmanchisquarer   r   re   rë   r   )r/   r	  r  rB   rð   s        r3   Útest_friedmanchisqzTestMisc.test_friedmanchisq|  sû   € â9Ú9Ú9ð;ˆô ×)Ñ)¨4Ð0ˆÜ˜F 1™I w°Ô2Ü˜F 1™I x°Ô3ä”#�q˜!˜R  Q¨¨1¨a°°A°qÐ9Ú4Ø��A�q˜"˜a  A q¬¨a°´3Ð7Ü�1�b˜!˜R¤ Q¨¨1¨a°°A°qÐ9ð;ˆô �N‰N˜1ÓˆÜ×)Ñ)¨1Ð-ˆÜ˜F 1™I v¨qÔ1Ü˜F 1™I v¨qÔ1ð -ˆ
Ü˜F J°4Ö8r5   N)r6   r7   r8   r
  r  r  r9   r5   r3   r  r  b  s   „ ò'ò.ó9r5   r  c                  ó*  — t        j                  ddd«      } dt        j                  ddd«      z  dz   }|t        j                  t        j                  ddd«      «      z  }t        j                  | |«      }t
        j                  }t        t        ||«      «       d}t        ||d¬«       d	t        |«      v sJ ‚t        |j                  d
«       t        |j                  d«       t        |j                  d«       t        |j                  d«       y )Nr   r_   r(  rW   rX   )ÚslopeÚ	interceptÚrvaluerè   ÚstderrTré   Úintercept_stderrgÍ5E%É?gÎâ ’+l$@g·4_Q c?g8Äèê¿Á?)r+   rú   Úsinr-   Ú
linregressr   ÚLinregressResultr   Ú
isinstancer   Údirr   r  r   r"  r#  )rB   rÂ   r  Úlrrð   s        r3   Útest_regress_simpler*  “  sã   € ä
�‰�A�s˜CÓ €AØŒb�k‰k˜!˜S #Ó&Ñ&¨Ñ+€AØŒ�‰”—‘˜A˜r 3Ó'Ó	(Ñ(€Aä×Ñ˜q !Ó$€Fô 
×	#Ñ	#€BÜŒJ�v˜rÓ"Ô#ØE€JÜ˜ 
¨tÕ4Ø¤ V£Ñ,Ð,Ð,ô ˜Ÿ™Ð&9Ô:Ü˜×(Ñ(Ð*<Ô=Ü˜Ÿ™Ð';Ô<Ü˜×/Ñ/Ð1DÕEr5   c                  óè   — t        j                  d«      } t         j                  j                  d«      }d}t        t        |¬«      5  t        j                  | |«       d d d «       y # 1 sw Y   y xY w)NrW   zBCannot calculate a linear regression if all x values are identicalrø   )r+   r²  rû   r  rý   r-   r%  )rB   rÂ   Úmsgs      r3   Útest_linregress_identical_xr-  ©  sW   € Ü
�‰�‹€AÜ
�	‰	×Ñ˜Ó€AØ
N€CÜ	”z¨Ô	-ñ  Ü×Ñ˜!˜QÔ÷ ÷  ñ  ús   ÁA(Á(A1c                   ó$   — e Zd Zd„ Zd„ Zd„ Zd„ Zy)ÚTestTheilslopesc                 óB  — t        j                  g d¢«      \  }}}}t        |d«       t        |d«       t        j                  g d¢d¬«      \  }}}}t        |d«       t        |d«       t        j                  j                  g d¢g d¢¬«      }t        j                  |«      \  }}}}t        |d	«       t        |d
«       t        j                  |d¬«      \  }}}}t        |d	«       t        |d«       g d¢}g d¢}t        j                  ||d«      \  }}}}t        |d«       t        |d«       t        |dd¬«       t        |dd¬«       t        j                  ||dd¬«      \  }}}}t        |d«       t        |d«       t        |dd¬«       t        |dd¬«       y )N)r   r   r   r�   Újointr  r¬   )r   r   r_   r   )FFTFra   rd   gUUUUUUå?)r   rK   rL   rM   rW   r`  r&  )r¾   é   é   rX   r  é7   éN   gìQ¸…ë±?rM   r"   g…ëQ¸…@rK   rü  g®Gáz®@r   )r-   Útheilslopesr   r+   re   r,   )r/   r  r   ÚlowerÚupperrÂ   rB   s          r3   Útest_theilslopesz TestTheilslopes.test_theilslopes²  s  € ä)/×);Ñ);ºIÓ)FÑ&ˆˆy˜% Ü˜E 3Ô'Ü˜I sÔ+ä)/×);Ñ);ºIØCJô*LÑ&ˆˆy˜% ä˜E 3Ô'Ü˜I sÔ+ô �E‰E�K‰KšÒ-HˆKÓIˆÜ)/×);Ñ);¸AÓ)>Ñ&ˆˆy˜% Ü˜E 4Ô(Ü˜I tÔ,ä)/×);Ñ);¸AØCJô*LÑ&ˆˆy˜% ä˜E 4Ô(Ü˜I sÔ+ò %ˆÚ'ˆÜ)/×);Ñ);¸A¸qÀ$Ó)GÑ&ˆˆy˜% Ü˜E 1Ô%Ü˜I sÔ+Ü˜E 4°Õ3Ü˜E 4°Õ3ä)/×);Ñ);¸A¸qÀ$ØCJô*LÑ&ˆˆy˜% ä˜E 1Ô%Ü˜I sÔ+Ü˜E 4°Õ3Ü˜E 4°Ö3r5   c                 óè  — d}t        j                  t        |¬«      5  t        j                  ddgddg«      }t        j                  t        j                  |«      «      sJ ‚	 d d d «       t        «       5 }|j                  t        d«       t        j                  g d¢g d¢«      }t        |ddt
        j                  t
        j                  f«       d d d «       y # 1 sw Y   ŒyxY w# 1 sw Y   y xY w)NzFAll `x` coordinates.*|Mean of empty slice.|invalid value encountered.*rø   r   r   zinvalid value encountered...©r   r   r   )r   r   r   )rˆ   Úwarnsr±   r-   r6  r+   ÚallÚisnanr   Úfilterr   r   )r/   r,  rï   Úsups       r3   Útest_theilslopes_warningsz)TestTheilslopes.test_theilslopes_warningsÙ  s¾   € àVˆÜ�\‰\œ.°Ô4ñ 	)Ü×$Ñ$ a¨ V¨a°¨VÓ4ˆCÜ—6‘6œ"Ÿ(™( 3›-Ô(Ð(Ñ(÷	)ô Ó ð 	9 CØ�J‰J”~Ð'EÔFÜ×$Ñ$¢Y²	Ó:ˆCÜ˜C ! Q¬¯©´·±Ð!7Ô8÷	9ð 	9÷	)ð 	)ú÷	9ð 	9ús   žACÁ6AC(ÃC%Ã(C1c                 ó&  — g d¢}g d¢}t        j                  ||«      \  }}}}t        j                  ||«      }t        ||j                  «       t        ||j                  «       t        ||j
                  «       t        ||j                  «       y)zv
        Simple test to ensure tuple backwards-compatibility of the returned
        TheilslopesResult object
        ©r   rK   rM   ©rM   rc   r�   N)r-   r6  r
   r  r   Ú	low_slopeÚ
high_slope)r/   rÂ   rB   r  r   rE  rF  r  s           r3   Ú'test_theilslopes_namedtuple_consistencyz7TestTheilslopes.test_theilslopes_namedtuple_consistencyå  s{   € ò
 ˆÚˆÜ28×2DÑ2DÀQÈÓ2JÑ/ˆˆy˜) ZÜ×#Ñ# A qÓ)ˆô 	�U˜FŸL™LÔ)Ü�Y × 0Ñ 0Ô1Ü�Y × 0Ñ 0Ô1Ü�Z ×!2Ñ!2Õ3r5   c                 ó  — t         j                  j                  d«      }|j                  dddt         j                  ¬«      }t        j                  ||«      }t         j                  j                  |j                  d«       y )Nl    UÕ5<H r   éÿ   rW   )r…  r<   r   )
r+   rû   r‡  ÚintegersÚuint8r   r6  Útestingr   r  )r/   rŒ  rÂ   rï   s       r3   Útest_gh19678_uint8z"TestTheilslopes.test_gh19678_uint8õ  s]   € ô �i‰i×#Ñ#Ð$4Ó5ˆØ�L‰L˜˜C b´·±ˆLÓ9ˆÜ×Ñ  1Ó%ˆÜ
�
‰
×"Ñ" 3§9¡9¨aÕ0r5   N)r6   r7   r8   r9  rA  rG  rM  r9   r5   r3   r/  r/  ±  s   „ ò$4òN	9ò4ó 1r5   r/  c                  óV  — dt        j                  d«      z  dz   } t        t        j                  | «      d«       t        t        j                  | d¬«      d«       dt        j                  d«      z  }d|z  dz
  } t        t        j                  | |«      d	«       t        t        j                  | |d¬«      d	«       d
| d d t        t        j                  | |«      d	«       t        j                  d«      }dd|z  z   t
        j                  j                  dd¬«      z   } t        j                  || «      \  }}}}}t        j                  | |«      \  }}t        ||d¬«       t        ||d¬«       t        j                  | |d¬«      \  }}t        ||d¬«       t        ||d¬«       y )NrK   rW   r�   )r   r�   Úseparater  r™   r    )r  rÌ  r  rM   gffffffÀr×   éç   ©r…  Úrandom_stater÷   r=   )
r+   r­   r
   r-   Úsiegelslopesr   ÚnormÚrvsr%  r   )rÂ   rB   Ú	slope_olsÚintercept_olsÚ_r  r   s          r3   Útest_siegelslopesrY  þ  sg  € à	ŒB�I‰I�b‹MÑ˜CÑ€AÜ”×$Ñ$ QÓ'¨Ô4Ü”×$Ñ$ Q¨zÔ:¸JÔGà	ŒB�I‰I�b‹MÑ€AØ	ˆA‰�‰€AÜ”×$Ñ$ Q¨Ó*¨KÔ8Ü”×$Ñ$ Q¨°*Ô=¸{ÔKð €A€b€q€EÜ”×$Ñ$ Q¨Ó*¨KÔ8ô 	�	‰	�"‹€AØˆs�1‰u‰”u—z‘z—~‘~¨2¸C�~Ó@Ñ@€AÜ(-×(8Ñ(8¸¸AÓ(>Ñ%€Iˆ}˜a  Aä×*Ñ*¨1¨aÓ0Ñ€Eˆ9Ü�E˜9¨3Õ/Ü�I˜}°3Õ7ä×*Ñ*¨1¨a¸
ÔCÑ€Eˆ9Ü�E˜9¨3Õ/Ü�I˜}°3Ö7r5   c                  óÊ   — g d¢} g d¢}t        j                  | |«      \  }}t        j                  | |«      }t        ||j                  «       t        ||j                  «       y)zl
    Simple test to ensure tuple backwards-compatibility of the returned
    SiegelslopesResult object.
    rC  rD  N)r-   rS  r
   r  r   )rÂ   rB   r  r   r  s        r3   Ú(test_siegelslopes_namedtuple_consistencyr[    sW   € ò
 	€AÚ€AÜ×*Ñ*¨1¨aÓ0Ñ€Eˆ9Ü× Ñ   AÓ&€Fô �˜Ÿ™Ô%Ü�˜F×,Ñ,Õ-r5   c            
      óî  — t         j                  j                  d«      } | j                  d¬«      }t        j                  |«      \  }}d„ }t        d«      D ]1  }t        t        j                   ||d d …|f   «      «      ||   «       Œ3 t        j                  t        |j                  d   «      D �cg c]  } ||d d …|f   «      ‘Œ c}«      }t        t        j                  |«      |«       y c c}w )Nl   ,ûÎ'  )r_   rM   r„  c                 ó$  — t        | «      }t        j                  |«      }| | d d …t        j                  f   z
  }||d d …t        j                  f   z
  }t        j                  t        j
                  ||ft        ¬«      d¬«      }||   ||   z  S )Nr’   r   )Úk)Úlenr+   r­   ÚnewaxisÚtriuÚonesrÀ  )Úyirÿ   rB   ÚdyÚdxrb   s         r3   Údijkz&test_sen_seasonal_slopes.<locals>.dijk0  sx   € Ü�‹GˆÜ�I‰I�a‹LˆØ�"’QœŸ
™
�]Ñ#Ñ#ˆØ�’1”b—j‘j�=Ñ!Ñ!ˆÜ�w‰w”r—w‘w  1˜v¬TÔ2°aÔ8ˆØ�$‰x˜˜4™Ñ Ð r5   rM   r   )
r+   rû   r‡  r-   Úsen_seasonal_slopesrí   r   Úmedianr  ro  )rŒ  rB   Úintra_slopeÚinter_sloperf  rÅ  Ú
all_slopess          r3   Útest_sen_seasonal_slopesrl  *  sÈ   € Ü
�)‰)×
Ñ
Ð 3Ó
4€CØ�
‰
˜ˆ
Ó!€AÜ%×9Ñ9¸!Ó<Ñ€K�ò!ô �1‹Xò BˆÜœŸ	™	¡$ qª¨A¨¡w£-Ó0°+¸a±.ÕAðBô —‘¼¸a¿g¹gÀa¹jÓ8IÖ J°1¡ aª¨1¨¡g¥Ò JÓK€JÜ”B—I‘I˜jÓ)¨;Õ7ùò !Ks   Â5C2c                  ó®   — t        j                  t        j                  d«      dd«      } t	        | j
                  t        j                  g d¢«      «       y )NrL   r   )rO   r�   r  )r-   Úplotting_positionsr+   r­   r   r0   r,   )Úposs    r3   Útest_plotting_positionsrp  ?  s6   € ä
×
#Ñ
#¤B§I¡I¨a£L°!°QÓ
7€CÜ˜cŸh™h¬¯©Ò1BÓ(CÕDr5   c                   óÆ   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	e
j                  j                  dd	d
g«      e
j                  j                  dddg«      d„ «       «       Zd„ Zy)ÚTestNormalitytestsc                 ól  — t        j                  d«      dz  }t        t        j                  |«      t        j                  |«      «       t        t        j                  |«      t        j                  |«      «       t        t        j                  |«      t        j                  |«      «       t
        j                  t
        j                  t
        j                  g}t        j                  t        j                  t        j                  g}g d¢}t        ||«      D ]‰  \  }}t        j                  t        t        ¬«      5   ||«      }t        j                  |j                  «      sJ ‚t        j                  |j                  «      sJ ‚	 d d d «       t!        t"        ||«       Œ‹ y # 1 sw Y   ŒxY w)N©rê  rn   r   r   rK   rL   rê  rn   r   r   rK   rL   rê  rn   r   r   rK   rL   rê  rn   r   r   rK   rL   rK   rJ   rø   )r+   r,   r   r-   Ú
normaltestr   ÚskewtestÚkurtosistestÚziprˆ   r<  r   r   r>  r  rè   r  rý   )r/   rB   ÚfuncsÚmfuncsÚfuncÚmfuncrï   s          r3   Útest_vs_nonmaskedz$TestNormalitytests.test_vs_nonmaskedH  s<  € Ü�H‰HÐ+Ó,¨aÑ/ˆÜ!¤&×"3Ñ"3°AÓ"6Ü"'×"2Ñ"2°1Ó"5ô	7ä!¤&§/¡/°!Ó"4Ü"'§.¡.°Ó"3ô	5ä!¤&×"5Ñ"5°aÓ"8Ü"'×"4Ñ"4°QÓ"7ô	9ô ×!Ñ!¤5§>¡>´5×3EÑ3EÐFˆÜ×#Ñ#¤V§_¡_´f×6IÑ6IÐJˆÚˆÜ˜u fÓ-ò 	0‰KˆD�%Ü—‘Ô0Ô8MÔNñ ,Ù˜1“g�Ü—x‘x §¡Ô.Ð.Ð.Ü—x‘x §
¡
Ô+Ð+Ñ+÷,ô œ* e¨QÕ/ñ	0÷,ð ,ús   ÅAF*Æ*F3	c                 ól  — t        j                  d«      dz  }t        t        j                  |d ¬«      t        j                  |«      «       t        t        j
                  |d ¬«      t        j
                  |«      «       t        t        j                  |d ¬«      t        j                  |«      «       y )Nrt  rK   r    )r+   r,   r   r-   ru  rv  rw  r¥   s     r3   Útest_axis_Nonez!TestNormalitytests.test_axis_None[  s{   € ä�H‰HÐ&Ó'¨Ñ*ˆÜœ×)Ñ)¨!°$Ô7¼×9JÑ9JÈ1Ó9MÔNÜœŸ™¨°Ô5´v·±ÀqÓ7IÔJÜœ×+Ñ+¨A°DÔ9Ü×+Ñ+¨AÓ.õ	0r5   c                 ó*  — t        j                  d«      dz  }t         j                  j                  t         j                  t         j                  |df   t         j                  ddg|j
                  z  df   ¬«      }t        t        j                  |«      t        j                  |«      «       t        t        j                  |«      t        j                  |«      «       t        t        j                  |«      t        j                  |«      «       y )Nrt  rK   rW   TFra   )r+   r,   re   Úr_Úinfr…  r   r-   ru  r   rv  rw  )r/   rB   Úxms      r3   Útest_maskedarray_inputz)TestNormalitytests.test_maskedarray_inputc  s½   € ä�H‰HÐ+Ó,¨aÑ/ˆÜ�U‰U�[‰[œŸ™œrŸv™v q¨"˜}Ñ-Ü Ÿe™e D¨5¨'°A·F±FÑ*:¸DÐ$@ÑAð ó Cˆäœ×)Ñ)¨"Ó-¬u×/?Ñ/?ÀÓ/BÔCÜœŸ™¨Ó+¬U¯^©^¸AÓ->Ô?Üœ×+Ñ+¨BÓ/´×1CÑ1CÀAÓ1FÕGr5   c                 óZ  — t        j                  d«      dz  }t        j                  |gdz  «      j                  }t        j
                  t        j                  t        j                  fD ]>  } ||«      } ||«      }t        |d   |d   gdz  «       t        |d   |d   gdz  «       Œ@ y )Nrt  rK   r   r   )	r+   r,   Úvstackr)  r-   ru  rv  rw  r   )r/   rB   Úx_2dr{  Úres_1dÚres_2ds         r3   Útest_nd_inputz TestNormalitytests.test_nd_inputl  s˜   € Ü�H‰HÐ+Ó,¨aÑ/ˆÜ�y‰y˜!˜˜q™Ó!×#Ñ#ˆÜ×&Ñ&¬¯©¼×9LÑ9LÐMò 	8ˆDÙ˜!“WˆFÙ˜$“ZˆFÜ˜F 1™I¨¨q©	 {°Q¡Ô7Ü˜F 1™I¨¨q©	 {°Q¡Õ7ñ		8r5   c                 ó~   — t        j                  d«      dz  }t        j                  |«      }d}t	        ||d¬«       y ©Nrt  rK   r  Tré   )r+   r,   r-   ru  r   ©r/   rB   rï   rð   s       r3   Ú!test_normaltest_result_attributesz4TestNormalitytests.test_normaltest_result_attributesu  s7   € Ü�H‰HÐ+Ó,¨aÑ/ˆÜ×Ñ Ó"ˆØ,ˆ
Ü˜C °Ö5r5   c                 ó~   — t        j                  d«      dz  }t        j                  |«      }d}t	        ||d¬«       y rŒ  )r+   r,   r-   rw  r   r�  s       r3   Ú#test_kurtosistest_result_attributesz6TestNormalitytests.test_kurtosistest_result_attributes{  s7   € Ü�H‰HÐ+Ó,¨aÑ/ˆÜ×!Ñ! !Ó$ˆØ,ˆ
Ü˜C °Ö5r5   c           
      óì   — g d¢}t        j                  t        |«      D ��cg c]  \  }}t        j                  ||«      ‘Œ c}}«      }t	        t        j                  |«      d   dk  d«       y c c}}w )N)
é€   r   é:   rš   r   é)   r#  r   r   é§   r   g{®Gáz„?T)r+   ÚhstackrÂ  Úfullr
   r-   rw  )r/   ÚcountsrÅ  ÚcrB   s        r3   Útest_regression_9033z'TestNormalitytests.test_regression_9033�  sZ   € ò 7ˆÜ�I‰I´¸6Ó1B×C©¨¨A”r—w‘w˜q !•}ÓCÓDˆÜ”V×(Ñ(¨Ó+¨AÑ.°Ñ5°tÕ<ùó Ds   ¢ A0
Útestrv  rw  rõ   rö   ró   c                 ó  — t         j                  j                  dddd¬«      }t        t         |«      }t        t        |«      } |||¬«      \  }} |||¬«      \  }}	t        ||d¬«       t        |	|d¬«       t        j                  |d	d
 t        j                  j                  |t        j                  |«      ¬«      } ||j                  «       |¬«      \  }} |||¬«      \  }}	t        ||d¬«       t        |	|d¬«       y )NrW   r'   rY   é{   ©ÚlocÚscaler…  rR  rô   gê-�™—q=©Úatolr   r™   ra   )r   rT  rU  Úgetattrr-   r   r+   r   re   r¿   r>  r5  )
r/   r›  rõ   rB   Ú
stats_testÚmstats_testÚz_exÚp_exr  r¸   s
             r3   Útest_alternativez#TestNormalitytests.test_alternativeˆ  sÞ   € ô �J‰J�N‰N˜r¨°2ÀCˆNÓHˆäœU DÓ)ˆ
Üœf dÓ+ˆá ¨{Ô;‰
ˆˆdÙ˜1¨+Ô6‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÕ,ô —‘ˆˆ!ˆAˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÙ §¡£¸KÔH‰
ˆˆdÙ˜1¨+Ô6‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÖ,r5   c                 óf  — t         j                  j                  dd¬«      }d}t        j                  t
        |¬«      5  t        j                  |d¬«       d d d «       t        j                  t
        |¬«      5  t        j                  |d¬«       d d d «       y # 1 sw Y   ŒExY w# 1 sw Y   y xY w)NrX   r�  rQ  z`alternative` must be...rø   rª   rô   )	r   rT  rU  rˆ   r	   rý   r-   rv  rw  )r/   rB   r,  s      r3   Útest_bad_alternativez'TestNormalitytests.test_bad_alternative�  sˆ   € Ü�J‰J�N‰N °ˆNÓ5ˆØ)ˆä�]‰]œ:¨SÔ1ñ 	4Ü�O‰O˜A¨7Õ3÷	4ô �]‰]œ:¨SÔ1ñ 	8Ü×Ñ ¨wÕ7÷	8ð 	8÷	4ð 	4ú÷	8ð 	8ús   ¿BÁ:B'ÂB$Â'B0N)r6   r7   r8   r}  r  r„  rŠ  rŽ  r�  rš  rˆ   r‰   rU  r¨  rª  r9   r5   r3   rr  rr  E  su   „ ò0ò&0òHò8ò6ò6ò=ð ‡[�[×Ñ˜V j°.Ð%AÓBØ‡[�[×Ñ˜]¨V°YÐ,?Ó@ñ-ó Aó Cð-ó&8r5   rr  c                   ó   — e Zd Zd„ Zy)ÚTestFOnewayc                 óì   — t        j                  ddgt         j                  ¬«      }t        j                  ddgt         j                  ¬«      }t        j                  ||«      }d}t        ||d¬«       y )	Ni�  i  r’   i  i  r  Tré   )r+   r,   Úuint16r-   Úf_onewayr   )r/   rS   r  rï   rð   s        r3   Útest_result_attributesz"TestFOneway.test_result_attributes©  sS   € Ü�H‰H�c˜3�Z¤r§y¡yÔ1ˆÜ�H‰H�c˜3�Z¤r§y¡yÔ1ˆÜ�o‰o˜a Ó#ˆØ,ˆ
Ü˜C °Ö5r5   N©r6   r7   r8   r°  r9   r5   r3   r¬  r¬  ¨  s   „ ó6r5   r¬  c                   óh   — e Zd Z ej                  g d¢«      Z ej                  g d¢«      Zd„ Zd„ Zy)ÚTestMannwhitneyu)ôr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r    r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   )žr   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r    r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   c                 óx   — t        j                  | j                  | j                  «      }d}t	        ||d¬«       y )Nr  Tré   )r-   ÚmannwhitneyurB   rÂ   r   )r/   rï   rð   s      r3   r°  z'TestMannwhitneyu.test_result_attributesÓ  s-   € Ü×!Ñ! $§&¡&¨$¯&©&Ó1ˆØ,ˆ
Ü˜C °Ö5r5   c                 ó"  — t        j                  | j                  | j                  «      }t	        j                  | j                  | j                  «      }|j
                  |j
                  k(  sJ ‚t        |j                  |j                  «       y r¦  )r-   rµ  rB   rÂ   r   r  r   rè   )r/   Úres1Úres2s      r3   Útest_against_statsz#TestMannwhitneyu.test_against_statsØ  s`   € ô ×"Ñ" 4§6¡6¨4¯6©6Ó2ˆÜ×!Ñ! $§&¡&¨$¯&©&Ó1ˆØ�~‰~ §¡Ò/Ð/Ð/Ü˜Ÿ™ T§[¡[Õ1r5   N)	r6   r7   r8   r+   r,   rB   rÂ   r°  r¹  r9   r5   r3   r³  r³  ±  s8   „ àˆ�‰ò *ó 	+€Að& 	ˆ�‰ò "ó 	#€Aò6ó
2r5   r³  c                   ó   — e Zd Zd„ Zy)ÚTestKruskalc                 ó`   — g d¢}g d¢}t        j                  ||«      }d}t        ||d¬«       y )N)r   rL   r™   rš   r¾   )rK   rM   rc   r�   rW   r  Tré   )r-   Úkruskalr   rî   s        r3   r°  z"TestKruskal.test_result_attributesã  s-   € ÚˆÚˆä�n‰n˜Q Ó"ˆØ,ˆ
Ü˜C °Ö5r5   Nr±  r9   r5   r3   r»  r»  â  s   „ ó6r5   r»  c                   ó~   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	e
j                  j                  dd	d
g«      d„ «       Zy)ÚTestTtest_relc                 ó  — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   |d d …df   «      }t        j
                  |d d …df   |d d …df   «      }t        ||«       t	        j
                  |d d …df   |d d …df   d ¬«      }t        j
                  |d d …df   |d d …df   d ¬«      }t        ||«       t	        j
                  |d d …d d…f   |d d …dd …f   d¬«      }t        j
                  |d d …d d…f   |d d …dd …f   d¬«      }t        ||«       t        j
                  |d d …d d…f   |d d …dd …f   «      }t        ||«       y )	Né‡Ö rX   rM   ©r   r   r   rK   r   r   r    rK   )r+   rû   r4  Úrandnr   Ú	ttest_relr-   r   )r/   Úoutcomer·  r¸  Úres3s        r3   r}  zTestTtest_rel.test_vs_nonmaskedî  s[  € Ü
�	‰	�‰�wÔÜ—)‘)—/‘/ " aÓ(ª<Ñ7ˆô �‰˜w¢q¨! t™}¨g²a¸°d©mÓ<ˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}Ó=ˆÜ˜˜dÔ#ô �‰˜w¢q¨! t™}¨g²a¸°d©mÀ$ÔGˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}À4ÔHˆÜ˜˜dÔ#Ü�‰˜w¢q¨"¨1¨" u™~¨w²q¸!¹"°u©~ÀAÔFˆÜ×Ñ ª¨2¨A¨2¨¡°º¸1¹2¸±ÀQÔGˆÜ˜˜dÔ#ô ×Ñ ª¨2¨A¨2¨¡°º¸1¹2¸±Ó?ˆÜ˜˜dÕ#r5   c                 ól  — t         j                  j                  d«       t        j                  t         j                  j                  dd«      g d¢g d¢g¬«      }t        «       5 }|j                  t        d«       |d d …df   |d d …d	f   ft         j                  t         j                  gd
dgffD ]k  }t        j                  |Ž \  }}t        |t         j                  t         j                  f«       t        |t         j                  t         j                  f«       Œm 	 d d d «       y # 1 sw Y   y xY w©NrÁ  rL   rK   rñ  r;  ra   ú%invalid value encountered in absoluter   r   r   r   )r+   rû   r4  re   r¿   rÃ  r   r?  r±   r   r-   rÄ  r   ©r/   rÅ  r@  ÚpairÚtr¸   s         r3   Útest_fully_maskedzTestTtest_rel.test_fully_masked  sé   € Ü
�	‰	�‰�wÔÜ—/‘/¤"§)¡)§/¡/°!°QÓ"7Ú(1²9Ð'=ô?ˆäÓ ð 	8 CØ�J‰J”~Ð'NÔOØ!¢! Q $™-¨²°A°©Ð7ÜŸ6™6¤2§6¡6Ð*¨S°#¨JÐ7ð9ò 8�ä×'Ñ'¨Ð.‘��1Ü" 1¤r§v¡v¬r¯v©vÐ&6Ô7Ü" 1¤r§v¡v¬r¯v©vÐ&6Õ7ñ	8÷	8÷ 	8ñ 	8úó   Á%B;D*Ä*D3c                 óô   — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   |d d …df   «      }d}t        ||d¬	«       y ©
NrÁ  rX   rM   rÂ  r   r   r  Tré   )r+   rû   r4  rÃ  r-   rÄ  r   ©r/   rÅ  rï   rð   s       r3   r°  z$TestTtest_rel.test_result_attributes  ó^   € Ü
�	‰	�‰�wÔÜ—)‘)—/‘/ " aÓ(ª<Ñ7ˆä×Ñ˜w¢q¨! t™}¨g²a¸°d©mÓ<ˆØ,ˆ
Ü˜C °Ö5r5   c           
      óÊ  — t        t        t        j                  t	        j
                  d«      t	        j
                  d«      «       t	        j
                  d«      }t        t        t        j                  |j                  ddd«      |j                  ddd«      d¬«       t        t        t        j                  |j                  ddd«      |j                  ddd«      d¬«       y )	NrW   r	  rN   rK   rL   rM   r   r    )r  rý   r-   rÄ  r+   r­   rb  r¥   s     r3   Útest_invalid_input_sizez%TestTtest_rel.test_invalid_input_size  sž   € Ü”j¤&×"2Ñ"2Ü—i‘i “m¤R§Y¡Y¨r£]ô	4ä�I‰I�b‹MˆÜ”j¤&×"2Ñ"2Ø—i‘i  1 aÓ(¨!¯)©)°A°q¸!Ó*<À1õ	Fä”j¤&×"2Ñ"2Ø—i‘i  1 aÓ(¨!¯)©)°A°q¸!Ó*<À1ö	Fr5   c                 ó’   — t        j                  g g «      }t        t        j                  t        j
                  |«      «      «       y r¦  )r-   rÄ  r   r+   r=  r>  ©r/   r·  s     r3   Ú
test_emptyzTestTtest_rel.test_empty   ó-   € Ü×Ñ  BÓ'ˆÜ”—‘”r—x‘x “~Ó&Õ'r5   c                 ó8  — t        j                  g d¢g d¢«      \  }}t        t        j                  |«      |ft        j
                  df«       t        «       5 }|j                  t        d«       t        j                  g d¢g d¢«      \  }}t        |t        j                  t        j                  t        j                  g«      «       t        |t        j                  t        j                  t        j                  g«      «       d d d «       y # 1 sw Y   y xY w)Nr;  rñ  r   rÉ  )r-   Ú	ttest_indr
   r+   Úabsr‚  r   r?  r±   r   r,   r   ©r/   rÌ  r¸   r@  s       r3   Útest_zero_divisionz TestTtest_rel.test_zero_division$  s¸   € Ü×Ñ¢	ª9Ó5‰ˆˆ1Ü”b—f‘f˜Q“i �^¤b§f¡f¨a [Ô1äÓ ð 	> CØ�J‰J”~Ð'NÔOÜ×#Ñ#¢IªyÓ9‰DˆAˆqÜ˜q¤"§(¡(¬B¯F©F´B·F±FÐ+;Ó"<Ô=Ü˜q¤"§(¡(¬B¯F©F´B·F±FÐ+;Ó"<Ô=÷		>÷ 	>ñ 	>ús   ÁB.DÄDc                 ó    — d}t        j                  t        |¬«      5  t        j                  g d¢g d¢d¬«       d d d «       y # 1 sw Y   y xY w©Nú4alternative must be 'less', 'greater' or 'two-sided'rø   ©r   rK   rL   )rM   r™   rc   Úfoorô   ©rˆ   r	   rý   r-   rÚ  ©r/   r,  s     r3   rª  z"TestTtest_rel.test_bad_alternative.  ó@   € ØEˆÜ�]‰]œ:¨SÔ1ñ 	FÜ×ÑšYª	¸uÕE÷	F÷ 	Fñ 	Fúó   žAÁArõ   rö   ró   c                 ó,  — t         j                  j                  dddd¬«      }t         j                  j                  dddd¬«      }t        j                  |||¬«      \  }}t	        j                  |||¬«      \  }}t        ||d	¬
«       t        ||d	¬
«       t        j                  |dd t        j                  |dd t        j                  j                  |t        j                  |«      ¬«      }t        j                  j                  |t        j                  |«      ¬«      }t	        j                  |||¬«      \  }}t        j                  |j                  «       |j                  «       |¬«      \  }}t        ||d	¬
«       t        ||d	¬
«       y )NrW   r™   r  é*   rž  r�   rK   rô   rP   r=   r   ra   )r   rT  rU  rÄ  r-   r   r+   r   re   r¿   r>  r5  ©r/   rõ   rB   rÂ   Út_exr§  rÌ  r¸   s           r3   r¨  zTestTtest_rel.test_alternative3  s3  € ä�J‰J�N‰N˜r¨°À"ˆNÓEˆÜ�J‰J�N‰N˜q¨°ÀˆNÓDˆä—_‘_ Q¨°{ÔC‰
ˆˆdÜ×Ñ  1°+Ô>‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÕ,ô —&‘&ˆˆ!ˆBˆÜ—&‘&ˆˆ!ˆBˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÜ×Ñ  1°+Ô>‰ˆˆ1Ü—_‘_ Q§\¡\£^°Q·\±\³^Ø1<ô>‰
ˆˆdä˜˜4 eÕ,Ü˜˜4 eÖ,r5   N)r6   r7   r8   r}  rÍ  r°  rÔ  r×  rÝ  rª  rˆ   r‰   rU  r¨  r9   r5   r3   r¿  r¿  í  sR   „ ò$ò*
8ò6òFò(ò>òFð
 ‡[�[×Ñ˜]¨V°YÐ,?Ó@ñ-ó Añ-r5   r¿  c                   óx   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Ze	j                  j                  ddd	g«      d
„ «       Zy)ÚTestTtest_indc                 ól  — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   |d d …df   «      }t        j
                  |d d …df   |d d …df   «      }t        ||«       t	        j
                  |d d …df   |d d …df   d ¬«      }t        j
                  |d d …df   |d d …df   d ¬«      }t        ||«       t	        j
                  |d d …d d…f   |d d …dd …f   d¬«      }t        j
                  |d d …d d…f   |d d …dd …f   d¬«      }t        ||«       t        j
                  |d d …d d…f   |d d …dd …f   «      }t        ||«       t	        j
                  |d d …df   |d d …df   d	¬
«      }t        j
                  |d d …df   |d d …df   d	¬
«      }t        ||«       t	        j
                  |d d …df   |d d …df   d¬
«      }t        j
                  |d d …df   |d d …df   d¬
«      }t        ||«       y )NrÁ  rX   rM   rÂ  r   r   r    rK   T©Ú	equal_varF)r+   rû   r4  rÃ  r   rÚ  r-   r   )r/   rÅ  r·  r¸  rÆ  Úres4Úres5s          r3   r}  zTestTtest_ind.test_vs_nonmaskedJ  sý  € Ü
�	‰	�‰�wÔÜ—)‘)—/‘/ " aÓ(ª<Ñ7ˆô �‰˜w¢q¨! t™}¨g²a¸°d©mÓ<ˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}Ó=ˆÜ˜˜dÔ#ô �‰˜w¢q¨! t™}¨g²a¸°d©mÀ$ÔGˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}À4ÔHˆÜ˜˜dÔ#Ü�‰˜w¢q¨"¨1¨" u™~¨w²q¸!¹"°u©~ÀAÔFˆÜ×Ñ ª¨2¨A¨2¨¡°º¸1¹2¸±ÀQÔGˆÜ˜˜dÔ#ô ×Ñ ª¨2¨A¨2¨¡°º¸1¹2¸±Ó?ˆÜ˜˜dÔ#ô �‰˜w¢q¨! t™}¨g²a¸°d©mÀtÔLˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}ÈÔMˆÜ˜˜dÔ#Ü�‰˜w¢q¨! t™}¨g²a¸°d©mÀuÔMˆÜ×Ñ ª¨1¨¡¨w²q¸!°t©}ÈÔNˆÜ˜˜dÕ#r5   c                 ól  — t         j                  j                  d«       t        j                  t         j                  j                  dd«      g d¢g d¢g¬«      }t        «       5 }|j                  t        d«       |d d …df   |d d …d	f   ft         j                  t         j                  gd
dgffD ]k  }t        j                  |Ž \  }}t        |t         j                  t         j                  f«       t        |t         j                  t         j                  f«       Œm 	 d d d «       y # 1 sw Y   y xY wrÈ  )r+   rû   r4  re   r¿   rÃ  r   r?  r±   r   r-   rÚ  r   rÊ  s         r3   rÍ  zTestTtest_ind.test_fully_maskedg  sæ   € Ü
�	‰	�‰�wÔÜ—/‘/¤"§)¡)§/¡/°!°QÓ"7ºyÊ)Ð>TÔUˆÜÓ ð 	8 CØ�J‰J”~Ð'NÔOØ!¢! Q $™-¨²°A°©Ð7ÜŸ6™6¤2§6¡6Ð*¨S°#¨JÐ7ð9ò 8�ä×'Ñ'¨Ð.‘��1Ü" 1¤r§v¡v¬r¯v©vÐ&6Ô7Ü" 1¤r§v¡v¬r¯v©vÐ&6Õ7ñ	8÷	8÷ 	8ñ 	8úrÎ  c                 óô   — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   |d d …df   «      }d}t        ||d¬	«       y rÐ  )r+   rû   r4  rÃ  r-   rÚ  r   rÑ  s       r3   r°  z$TestTtest_ind.test_result_attributesr  rÒ  r5   c                 ó’   — t        j                  g g «      }t        t        j                  t        j
                  |«      «      «       y r¦  )r-   rÚ  r   r+   r=  r>  rÖ  s     r3   r×  zTestTtest_ind.test_emptyz  rØ  r5   c                 ó  — t        j                  g d¢g d¢«      \  }}t        t        j                  |«      |ft        j
                  df«       t        «       5 }|j                  t        d«       t        j                  g d¢g d¢«      \  }}t        |t        j                  t        j                  f«       t        |t        j                  t        j                  f«       d d d «       t        j                  g d¢g d¢d¬«      \  }}t        t        j                  |«      |ft        j
                  df«       t        t        j                  g d¢g d¢d¬«      t        j                  t        j                  f«       y # 1 sw Y   ŒžxY w)Nr;  rñ  r   rÉ  Frî  )r-   rÚ  r
   r+   rÛ  r‚  r   r?  r±   r   r   rÜ  s       r3   rÝ  z TestTtest_ind.test_zero_division~  s  € Ü×Ñ¢	ª9Ó5‰ˆˆ1Ü”b—f‘f˜Q“i �^¤b§f¡f¨a [Ô1äÓ ð 	4 CØ�J‰J”~Ð'NÔOÜ×#Ñ#¢IªyÓ9‰DˆAˆqÜ˜q¤2§6¡6¬2¯6©6Ð"2Ô3Ü˜q¤2§6¡6¬2¯6©6Ð"2Ô3÷		4ô ×Ñ¢	ª9ÀÔF‰ˆˆ1Ü”b—f‘f˜Q“i �^¤b§f¡f¨a [Ô1Üœ6×+Ñ+ªI²yØ6;ô=Ü?A¿v¹vÄrÇvÁvÐ>Nõ	P÷	4ð 	4ús   ÁBE>Å>Fc                 ó    — d}t        j                  t        |¬«      5  t        j                  g d¢g d¢d¬«       d d d «       y # 1 sw Y   y xY wrß  rã  rä  s     r3   rª  z"TestTtest_ind.test_bad_alternative�  rå  ræ  rõ   rö   ró   c                 ó,  — t         j                  j                  dddd¬«      }t         j                  j                  dddd¬«      }t        j                  |||¬«      \  }}t	        j                  |||¬«      \  }}t        ||d¬	«       t        ||d¬	«       t        j                  |d
d t        j                  |dd t        j                  j                  |t        j                  |«      ¬«      }t        j                  j                  |t        j                  |«      ¬«      }t        j                  |j                  «       |j                  «       |¬«      \  }}t	        j                  |||¬«      \  }}t        ||d¬	«       t        ||d¬	«       y )NrW   rK   r_   r�  rž  r�   rô   rP   r=   r   r]   r^   ra   )r   rT  rU  rÚ  r-   r   r+   r   re   r¿   r>  r5  ré  s           r3   r¨  zTestTtest_ind.test_alternative’  s3  € ä�J‰J�N‰N˜r¨°À3ˆNÓGˆÜ�J‰J�N‰N˜q¨°À#ˆNÓFˆä—_‘_ Q¨°{ÔC‰
ˆˆdÜ×Ñ  1°+Ô>‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÕ,ô —&‘&ˆˆ!ˆBˆÜ—6‘6ˆˆ"ˆRˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÜ—_‘_ Q§\¡\£^°Q·\±\³^Ø1<ô>‰
ˆˆdä×Ñ  1°+Ô>‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÖ,r5   N©r6   r7   r8   r}  rÍ  r°  r×  rÝ  rª  rˆ   r‰   rU  r¨  r9   r5   r3   rì  rì  I  sM   „ ò$ò:	8ò6ò(òPòFð
 ‡[�[×Ñ˜]¨V°YÐ,?Ó@ñ-ó Añ-r5   rì  c                   óx   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Ze	j                  j                  ddd	g«      d
„ «       Zy)ÚTestTtest_1sampc                 ó  — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   d«      }t        j
                  |d d …df   d«      }t        ||«       y )NrÁ  rX   rM   rÂ  r   r   )r+   rû   r4  rÃ  r   Úttest_1sampr-   r   )r/   rÅ  r·  r¸  s       r3   r}  z!TestTtest_1samp.test_vs_nonmasked©  sk   € Ü
�	‰	�‰�wÔÜ—)‘)—/‘/ " aÓ(ª<Ñ7ˆô × Ñ  ª¨A¨¡°Ó2ˆÜ×!Ñ! 'ª!¨Q¨$¡-°Ó3ˆÜ˜˜dÕ#r5   c                 ó
  — t         j                  j                  d«       t        j                  t         j                  j                  d«      g d¢¬«      }t         j                  t         j                  f}t        «       5 }|j                  t        d«       t         j                  t         j                  fdf|dffD ]/  }t        j                  |Ž \  }}t        ||«       t        ||«       Œ1 	 d d d «       y # 1 sw Y   y xY w)NrÁ  rL   rñ  ra   rÉ  r¬   )r+   rû   r4  re   r¿   rÃ  r   r   r?  r±   r-   rü  r   )r/   rÅ  r  r@  rË  rÌ  r¸   s          r3   rÍ  z!TestTtest_1samp.test_fully_masked²  s¿   € Ü
�	‰	�‰�wÔÜ—/‘/¤"§)¡)§/¡/°!Ó"4º9ÔEˆÜ—F‘FœBŸF™FÐ#ˆÜÓ ð 	0 CØ�J‰J”~Ð'NÔOÜŸ6™6¤2§6¡6Ð*¨CÐ0°7¸C°.ÐAò 0�Ü×)Ñ)¨4Ð0‘��1Ü" 1 hÔ/Ü" 1 hÕ/ñ0÷	0÷ 	0ñ 	0ús   Â A/C9Ã9Dc                 óæ   — t         j                  j                  d«       t         j                  j                  dd«      g d¢z   }t	        j
                  |d d …df   d«      }d}t        ||d¬	«       y rÐ  )r+   rû   r4  rÃ  r-   rü  r   rÑ  s       r3   r°  z&TestTtest_1samp.test_result_attributes½  sV   € Ü
�	‰	�‰�wÔÜ—)‘)—/‘/ " aÓ(ª<Ñ7ˆä× Ñ  ª¨A¨¡°Ó2ˆØ,ˆ
Ü˜C °Ö5r5   c                 ó’   — t        j                  g d«      }t        t        j                  t        j
                  |«      «      «       y )Nr   )r-   rü  r   r+   r=  r>  rÖ  s     r3   r×  zTestTtest_1samp.test_emptyÅ  s-   € Ü×!Ñ! " aÓ(ˆÜ”—‘”r—x‘x “~Ó&Õ'r5   c                 óÌ  — t        j                  g d¢d«      \  }}t        t        j                  |«      |ft        j
                  df«       t        «       5 }|j                  t        d«       t        j                  g d¢d«      \  }}t        t        j                  |«      «       t        |t        j                  t        j                  f«       d d d «       y # 1 sw Y   y xY w)Nr;  r   r   rÉ  )r-   rü  r
   r+   rÛ  r‚  r   r?  r±   r   r>  r   r   rÜ  s       r3   rÝ  z"TestTtest_1samp.test_zero_divisionÉ  sŸ   € Ü×!Ñ!¢)¨QÓ/‰ˆˆ1Ü”b—f‘f˜Q“i �^¤b§f¡f¨a [Ô1äÓ ð 	4 CØ�J‰J”~Ð'NÔOÜ×%Ñ%¢i°Ó3‰DˆAˆqÜ”B—H‘H˜Q“KÔ Ü˜q¤2§6¡6¬2¯6©6Ð"2Ô3÷		4÷ 	4ñ 	4ús   ÁA:CÃC#c                 óœ   — d}t        j                  t        |¬«      5  t        j                  g d¢dd¬«       d d d «       y # 1 sw Y   y xY w)Nrà  rø   rá  rM   râ  rô   )rˆ   r	   rý   r-   rü  rä  s     r3   rª  z$TestTtest_1samp.test_bad_alternativeÓ  s@   € ØEˆÜ�]‰]œ:¨SÔ1ñ 	@Ü×Ñšy¨!¸Õ?÷	@÷ 	@ñ 	@ús   žAÁArõ   rö   ró   c                 ó<  — t         j                  j                  dddd¬«      }t        j                  |d|¬«      \  }}t	        j                  |d|¬«      \  }}t        ||d¬	«       t        ||d¬	«       t        j                  |d
d t        j                  j                  |t        j                  |«      ¬«      }t        j                  |j                  «       d|¬«      \  }}t	        j                  |d|¬«      \  }}t        ||d¬	«       t        ||d¬	«       y )NrW   rK   r_   r�  rž  r¾   rô   rP   r=   r   ra   )r   rT  rU  rü  r-   r   r+   r   re   r¿   r>  r5  )r/   rõ   rB   rê  r§  rÌ  r¸   s          r3   r¨  z TestTtest_1samp.test_alternativeØ  sè   € ä�J‰J�N‰N˜r¨°À3ˆNÓGˆä×&Ñ& q¨!¸ÔE‰
ˆˆdÜ×!Ñ! ! Q°KÔ@‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÕ,ô —&‘&ˆˆ!ˆBˆÜ�E‰E×Ñ˜q¤r§x¡x°£{ÐÓ3ˆÜ×&Ñ& q§|¡|£~°qØ3>ô@‰
ˆˆdä×!Ñ! ! Q°KÔ@‰ˆˆ1Ü˜˜4 eÕ,Ü˜˜4 eÖ,r5   Nrø  r9   r5   r3   rú  rú  ¨  sL   „ ò$ò	0ò6ò(ò4ò@ð
 ‡[�[×Ñ˜]¨V°YÐ,?Ó@ñ-ó Añ-r5   rú  c                   ó   — e Zd ZdZd„ Zy)ÚTestDescribez…
    Tests for mstats.describe.

    Note that there are also tests for `mstats.describe` in the
    class TestCompareWithStats.
    c                 óÒ  — t         j                  j                  g d¢g d¢gg d¢g d¢g¬«      }t        j                  |d¬«      }t        |j                  dd	g«       |j                  \  }}t        |d
dg«       t        |d	dg«       t        |j                  ddg«       t        |j                  ddg«       t        |j                  ddg«       t        |j                  ddg«       y )N)r   r   rK   rL   rM   r¾   )r™   r™   r   r¾   rL   rL   )r   r   r   r   r   r   )r   r   r   r   r   r   ra   r   r    r™   rM   r   rL   r   r"   r   r¬   gÍÌÌÌÌÌô¿ç       À)r+   re   r¿   r-   Údescriber
   ÚnobsÚminmaxÚmeanÚvarianceÚskewnessr   rÎ  )r/   rS   r  ÚaminÚamaxs        r3   Útest_basic_with_axisz!TestDescribe.test_basic_with_axisò  sÆ   € ä�E‰E×ÑÒ 2Ú 2ð 4â%7Ú%7ð%9ð ó :ˆô —‘ ¨Ô+ˆÜ�V—[‘[ 1 a &Ô)Ø—]‘]‰
ˆˆdÜ�T˜A˜q˜6Ô"Ü�T˜A˜q˜6Ô"Ü�V—[‘[ 3¨ *Ô-Ü�V—_‘_ s¨C jÔ1Ü�V—_‘_ s¨C jÔ1Ü˜Ÿ™¨$°¨Õ6r5   N)r6   r7   r8   r  r  r9   r5   r3   r  r  ë  s   „ ñó7r5   r  c                   óî   — e Zd ZdZd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	d„ Z
d	„ Zd
„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z d„ Z!d „ Z"d!„ Z#d"„ Z$d#„ Z%d$„ Z&d%„ Z'd&„ Z(y')(ÚTestCompareWithStatsa  
    Class to compare mstats results with stats results.

    It is in general assumed that scipy.stats is at a more mature stage than
    stats.mstats.  If a routine in mstats results in similar results like in
    scipy.stats, this is considered also as a proper validation of scipy.mstats
    routine.

    Different sample sizes are used for testing, as some problems between stats
    and mstats are dependent on sample size.

    Author: Alexander Loew

    NOTE that some tests fail. This might be caused by
    a) actual differences or bugs between stats and mstats
    b) numerical inaccuracies
    c) different definitions of routine interfaces

    These failures need to be checked. Current workaround is to have disabled these
    tests, but issuing reports on scipy-dev

    c                 ó
   — g d¢S )z9 Returns list of sample sizes to be used for comparison. )r  r_   rW   r™   r9   rë  s    r3   Úget_nzTestCompareWithStats.get_n  s   € â!Ð!r5   c                 ó  — t         j                  j                  d«       t         j                  j                  |«      }|t         j                  j                  |«      z   }t        j                  t        |«      dz   d«      }t        j                  t        |«      dz   d«      }||dt        |«       ||dt        |«       |dkD  }t         j                  j                  ||¬«      }t         j                  j                  ||¬«      }||||fS )NrÁ  r™   g €à7yÃACr   g €úÊsù?Cra   )r+   rû   r4  rÃ  r—  r_  re   r,   )r/   rÿ   rB   rÂ   rƒ  Úymrb   s          r3   Úgenerate_xy_samplez'TestCompareWithStats.generate_xy_sample  sÎ   € ô 	�	‰	�‰�wÔÜ�I‰I�O‰O˜AÓˆØ”—	‘	—‘ Ó"Ñ"ˆÜ�W‰W”S˜“V˜a‘Z Ó&ˆÜ�W‰W”S˜“V˜a‘Z Ó&ˆØˆˆ1ŒS�‹VˆØˆˆ1ŒS�‹VˆØ�D‰yˆÜ�U‰U�[‰[˜ $ˆ[Ó'ˆÜ�U‰U�[‰[˜ $ˆ[Ó'ˆØ�!�R˜ˆ|Ðr5   c                 óº  — t        j                  ||ft         j                  «      }t        j                  ||ft         j                  «      }t        j                  |dz   |ft         j                  «      }t        j                  |dz   |ft         j                  «      }t        |«      D ]&  }| j	                  |«      \  |d d …|f<   |d d …|f<   }}	Œ( |d| |d|…d d …f<   |d| |d|…d d …f<   t         j
                  j                  |t        j                  |«      ¬«      }t         j
                  j                  |t        j                  |«      ¬«      }||||fS )Nr™   r   ra   )r+   r—  r   rí   r  re   r,   r>  )
r/   rÿ   ÚnxrB   rÂ   rƒ  r  rÅ  re  rd  s
             r3   Úgenerate_xy_sample2Dz)TestCompareWithStats.generate_xy_sample2D.  s"  € Ü�G‰G�Q˜�GœRŸV™VÓ$ˆÜ�G‰G�Q˜�GœRŸV™VÓ$ˆÜ�W‰W�a˜‘c˜2�Y¤§¡Ó'ˆÜ�W‰W�a˜‘c˜2�Y¤§¡Ó'ˆä�r“ò 	BˆAØ'+×'>Ñ'>¸qÓ'AÑ$ˆAŠa�ˆd‰G�Q’q˜!�t‘W˜b¡"ð	Bð �q˜�Vˆˆ1ˆQˆ3’ˆ6‰
Ø�q˜�Vˆˆ1ˆQˆ3’ˆ6‰
Ü�U‰U�[‰[˜¤"§(¡(¨2£,ˆ[Ó/ˆÜ�U‰U�[‰[˜¤"§(¡(¨2£,ˆ[Ó/ˆØ�!�R˜ˆ|Ðr5   c                 ó*  — | j                  «       D ]€  }| j                  |«      \  }}}}t        j                  ||«      }t        j                  j                  ||«      }t        t        j                  |«      t        j                  |«      «       Œ‚ y r¦  )r  r  r   r%  r-   r   r+   r“   )r/   rÿ   rB   rÂ   rƒ  r  Úresult1Úresult2s           r3   Útest_linregressz$TestCompareWithStats.test_linregress=  su   € Ø—‘“ò 	FˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ×&Ñ& q¨!Ó,ˆGÜ—l‘l×-Ñ-¨b°"Ó5ˆGÜœBŸJ™J wÓ/´·±¸GÓ1DÕEñ		Fr5   c                 ó
  — | j                  «       D ]p  }| j                  |«      \  }}}}t        j                  ||«      \  }}t        j                  j                  ||«      \  }}	t        ||d¬«       t        ||	d¬«       Œr y )Nrk  rü  )r  r  r   r²   r-   r   ©
r/   rÿ   rB   rÂ   rƒ  r  r·   r¸   ÚrmÚpms
             r3   r»   z"TestCompareWithStats.test_pearsonrD  su   € Ø—‘“ò 	3ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—>‘> ! QÓ'‰DˆAˆqÜ—\‘\×*Ñ*¨2¨rÓ2‰FˆB�ä  2¨rÕ2Ü  2¨rÖ2ñ	3r5   c                 ó  — | j                  «       D ]n  }| j                  |«      \  }}}}t        j                  ||«      \  }}t        j                  j                  ||«      \  }}	t        ||d«       t        ||	d«       Œp y ©Nrk  )r  r  r   rê   r-   r   r  s
             r3   rñ   z#TestCompareWithStats.test_spearmanrM  su   € Ø—‘“ò 	+ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—?‘? 1 aÓ(‰DˆAˆqÜ—\‘\×+Ñ+¨B°Ó3‰FˆB�Ü  2 rÔ*Ü  2 rÕ*ñ	+r5   c                 óp   — t        j                  d«      }t        t        t        j
                  ||d«       y )Nrc   F)r+   r­   r  rý   r-   rê   r¥   s     r3   Ú!test_spearmanr_backcompat_usetiesz6TestCompareWithStats.test_spearmanr_backcompat_usetiesU  s'   € ô �I‰I�a‹LˆÜ”j¤&×"2Ñ"2°A°q¸%Õ@r5   c                 óª  — | j                  «       D ]À  }| j                  |«      \  }}}}t        j                  t	        |«      «      }t        j
                  j                  t	        |«      «      }t        ||d¬«       t        j                  t	        |«      «      }t        j
                  j                  t	        |«      «      }t        ||d¬«       ŒÂ y )Ng‚vIhÂ%<=r=   )r  r  r   r?   rÛ  r-   r   ©r/   rÿ   rB   rÂ   rƒ  r  r·   r   s           r3   Ú
test_gmeanzTestCompareWithStats.test_gmean[  s›   € Ø—‘“ò 	/ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—‘œC ›FÓ#ˆAÜ—‘×#Ñ#¤C¨£GÓ,ˆBÜ˜A˜r¨Õ.ä—‘œC ›FÓ#ˆAÜ—‘×#Ñ#¤C¨£GÓ,ˆBÜ˜A˜r¨Ö.ñ	/r5   c                 ó¦  — | j                  «       D ]¾  }| j                  |«      \  }}}}t        j                  t	        |«      «      }t        j
                  j                  t	        |«      «      }t        ||d«       t        j                  t	        |«      «      }t        j
                  j                  t	        |«      «      }t        ||d«       ŒÀ y ©NrW   )r  r  r   rE   rÛ  r-   r   r'  s           r3   Ú
test_hmeanzTestCompareWithStats.test_hmeanf  s›   € Ø—‘“ò 		+ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bä—‘œC ›FÓ#ˆAÜ—‘×#Ñ#¤C¨£GÓ,ˆBÜ  2 rÔ*ä—‘œC ›FÓ#ˆAÜ—‘×#Ñ#¤C¨£GÓ,ˆBÜ  2 rÕ*ñ		+r5   c                 ó^  — | j                  «       D ]š  }| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        ||d«       t        j                  |«      }t        j                  j                  |«      }t        ||d«       Œœ y r*  )r  r  r   r¾  r-   r   r'  s           r3   Ú	test_skewzTestCompareWithStats.test_skewr  s‰   € Ø—‘“ò 		+ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bä—
‘
˜1“ˆAÜ—‘×"Ñ" 2Ó&ˆBÜ  2 rÔ*ä—
‘
˜1“ˆAÜ—‘×"Ñ" 2Ó&ˆBÜ  2 rÕ*ñ		+r5   c                 ó^  — | j                  «       D ]š  }| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        ||d«       t        j                  |«      }t        j                  j                  |«      }t        ||d«       Œœ y r*  )r  r  r   r®  r-   r   r'  s           r3   r³  z TestCompareWithStats.test_moment~  s‰   € Ø—‘“ò 		+ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bä—‘˜Q“ˆAÜ—‘×$Ñ$ RÓ(ˆBÜ  2 rÔ*ä—‘˜Q“ˆAÜ—‘×$Ñ$ RÓ(ˆBÜ  2 rÕ*ñ		+r5   c                 óš  — | j                  «       D �]7  }| j                  |«      \  }}}}||j                  «       z
  |j                  «       z  }||j                  «       z
  |j                  «       z  }t	        t        j                  |«      |d¬«       t	        t        j                  |«      |d¬«       t	        t        j                  |«      t
        j                  j                  |dt        |«       «      d¬«       t	        t        j                  |«      t
        j                  j                  |dt        |«       «      d¬«       �Œ: y )Nç»½×Ùß|Û=r=   r   )	r  r  r
  r’  r   r   r  r-   r_  )r/   rÿ   rB   rÂ   rƒ  r  ÚzxÚzys           r3   r  z TestCompareWithStats.test_zscoreŠ  sñ   € Ø—‘“ó 	(ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bð �a—f‘f“h‘, !§%¡%£'Ñ)ˆBØ�a—f‘f“h‘, !§%¡%£'Ñ)ˆBô œEŸL™L¨›O¨R°eÕ<ÜœEŸL™L¨›O¨R°eÕ<ô œEŸL™L¨›O¬U¯\©\×-@Ñ-@ÀÀAÄcÈ!ÃfÀÓ-NØ!&õ(äœEŸL™L¨›O¬U¯\©\×-@Ñ-@ÀÀAÄcÈ!ÃfÀÓ-NØ!&÷(ñ	(r5   c                 ó^  — | j                  «       D ]š  }| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        ||d«       t        j                  |«      }t        j                  j                  |«      }t        ||d«       Œœ y r*  )r  r  r   rÎ  r-   r   r'  s           r3   rÏ  z"TestCompareWithStats.test_kurtosisœ  s‹   € Ø—‘“ò 	+ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—‘˜qÓ!ˆAÜ—‘×&Ñ& rÓ*ˆBÜ  2 rÔ*ä—‘˜qÓ!ˆAÜ—‘×&Ñ& rÓ*ˆBÜ  2 rÕ*ñ	+r5   c           	      ó®  — t        j                  d«      j                  dd«      }t         j                  j	                  |«      }t        j                  |d¬«      }t
        j                  j                  |d¬«      }t        |dd¬«       t        |dd¬«       | j                  «       D �])  }| j                  |«      \  }}}}	t        t
        j                  j                  |d d	¬
«      t        j                  |d d	¬
«      d¬«       t        t
        j                  j                  |	d d	¬
«      t        j                  |d d	¬
«      d¬«       t        t
        j                  j                  |d d¬
«      t        j                  |d d¬
«      d¬«       t        t
        j                  j                  |	d d¬
«      t        j                  |d d¬
«      d¬«       �Œ, y )NrX   r™   rM   r   r˜  gß'Üež @gñhãˆµøä>r¡  r   )r)   r™  é   rü  )r+   r­   rb  re   r,   r   rö  r-   r   r  r  r   )
r/   rS   Úamr·   r   rÿ   rB   rÂ   rƒ  r  s
             r3   r÷  zTestCompareWithStats.test_sem§  si  € ä�I‰I�b‹M×!Ñ! ! QÓ'ˆÜ�U‰U�[‰[˜‹^ˆÜ�I‰I�a˜aÔ ˆÜ�\‰\×Ñ˜b qÐÓ)ˆä˜˜:¨DÕ1Ü˜˜J¨TÕ2à—‘“ó 		MˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§¡× 0Ñ 0°¸$ÀQÐ 0Ó GÜ %§	¡	¨!°$¸QÔ ?ÈõMä¤§¡× 0Ñ 0°¸$ÀQÐ 0Ó GÜ %§	¡	¨!°$¸QÔ ?ÈõMä¤§¡× 0Ñ 0°¸$ÀQÐ 0Ó GÜ %§	¡	¨!°$¸QÔ ?ÈõMä¤§¡× 0Ñ 0°¸$ÀQÐ 0Ó GÜ %§	¡	¨!°$¸QÔ ?È÷Mñ		Mr5   c           	      ó^  — | j                  «       D ]š  }| j                  |«      \  }}}}t        j                  |d¬«      }t        j                  j                  |d¬«      }t        d«      D ]<  }t        t        j                  ||   «      t        j                  ||   «      d¬«       Œ> Œœ y )Nr   r˜  rc   r`  rü  )	r  r  r   r  r-   rí   r   r+   r“   )	r/   rÿ   rB   rÂ   rƒ  r  r·   r   Úiis	            r3   Útest_describez"TestCompareWithStats.test_describe¼  s•   € Ø—‘“ò 	0ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—‘˜q qÔ)ˆAÜ—‘×&Ñ& r°Ð&Ó2ˆBÜ˜A“hò 0�Ü#¤B§J¡J¨q°©uÓ$5Ü$&§J¡J¨r°"©vÓ$6Ø,.ö0ñ0ñ		0r5   c                 ót   — t        j                  t        j                  d«      «      }d}t	        ||d¬«       y )Nr™   )r  r	  r
  r  r  rÎ  Tré   )r-   r  r+   r­   r   )r/   r¨  rð   s      r3   Útest_describe_result_attributesz4TestCompareWithStats.test_describe_result_attributesÆ  s+   € Ü—‘¤§¡¨1£Ó.ˆð"ˆ
ä˜F J°4Ö8r5   c                 óÚ   — | j                  «       D ]X  }| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        ||«       ŒZ y r¦  )r  r  r   r¤   r-   r   r'  s           r3   Útest_rankdataz"TestCompareWithStats.test_rankdataÌ  s[   € Ø—‘“ò 	#ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—‘˜qÓ!ˆAÜ—‘×&Ñ& qÓ)ˆBÜ˜A˜rÕ"ñ		#r5   c                 óN  — | j                  «       D ]’  }| j                  |«      \  }}}}t        t        j                  |«      t        j
                  j	                  |«      d«       t        t        j                  |«      t        j
                  j	                  |«      d«       Œ” y r#  )r  r  r   r   Útmeanr-   ©r/   rÿ   rB   rÂ   rƒ  r  s         r3   Ú
test_tmeanzTestCompareWithStats.test_tmeanÓ  sw   € Ø—‘“ò 	KˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§¡¨A£¬u¯|©|×/AÑ/AÀ"Ó/EÀrÔJÜ¤§¡¨A£¬u¯|©|×/AÑ/AÀ"Ó/EÀrÕJñ	Kr5   c                 ó^  — | j                  «       D �]  }| j                  |«      \  }}}}t        t        j                  |d«      t        j
                  j	                  |d«      d«       t        t        j                  |d«      t        j
                  j	                  |d«      d«       t        t        j                  |d¬«      t        j
                  j	                  |d¬«      d«       t        t        j                  |d¬«      t        j
                  j	                  |d¬«      d«       �Œ y )Nr   rW   r    )Ú
upperlimit)r  r  r   r   Útmaxr-   r@  s         r3   Ú	test_tmaxzTestCompareWithStats.test_tmaxÙ  sç   € Ø—‘“ó 
	JˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§
¡
¨1¨RÓ 0Ü %§¡× 1Ñ 1°"°RÓ 8¸"ô>ä¤§
¡
¨1¨RÓ 0Ü %§¡× 1Ñ 1°"°RÓ 8¸"ô>ô  ¤§
¡
¨1¸Ô <Ü %§¡× 1Ñ 1°"ÀÐ 1Ó DÀbôJä¤§
¡
¨1¸Ô <Ü %§¡× 1Ñ 1°"ÀÐ 1Ó DÀböJñ
	Jr5   c                 óR  — | j                  «       D �]  }| j                  |«      \  }}}}t        t        j                  |«      t        j
                  j	                  |«      «       t        t        j                  |«      t        j
                  j	                  |«      «       t        t        j                  |d¬«      t        j
                  j	                  |d¬«      d«       t        t        j                  |d¬«      t        j
                  j	                  |d¬«      d«       �Œ y )Nr«   )Ú
lowerlimitrW   )r  r  r
   r   Útminr-   r   r@  s         r3   Ú	test_tminzTestCompareWithStats.test_tminæ  sÑ   € Ø—‘“ó 	KˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜœŸ™ A›¬¯©×(9Ñ(9¸"Ó(=Ô>ÜœŸ™ A›¬¯©×(9Ñ(9¸"Ó(=Ô>ä¤§
¡
¨1¸Ô =Ü %§¡× 1Ñ 1°"ÀÐ 1Ó EÀrôKä¤§
¡
¨1¸Ô =Ü %§¡× 1Ñ 1°"ÀÐ 1Ó EÀröKñ	Kr5   c           	      óú   — | j                  «       D ]h  }| j                  |«      \  }}}}t        j                  ||«      }t        j                  j                  ||«      }t        ||dt        |«       d¬«       Œj y )Nr   r0  r¡  )r  r  r   rþ  r-   r   r_  )r/   rÿ   rB   rÂ   rƒ  r  r  Úzms           r3   r   zTestCompareWithStats.test_zmapñ  sk   € Ø—‘“ò 	9ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ—
‘
˜1˜aÓ ˆAÜ—‘×"Ñ" 2 rÓ*ˆBÜ˜A˜r !¤C¨£F˜|°%Ö8ñ		9r5   c                 óR  — | j                  «       D ]”  }| j                  |«      \  }}}}t        t        j                  |«      t        j
                  j	                  |«      d¬«       t        t        j                  |«      t        j
                  j	                  |«      d¬«       Œ– y ©Nr`  rü  )r  r  r   r   r¶  r-   r@  s         r3   r·  z#TestCompareWithStats.test_variationø  s}   € Ø—‘“ò 	,ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§¡°Ó 2´E·L±L×4JÑ4JÈ2Ó4NØ(*õ,ä¤§¡°Ó 2´E·L±L×4JÑ4JÈ2Ó4NØ(*ö,ñ		,r5   c                 óR  — | j                  «       D ]”  }| j                  |«      \  }}}}t        t        j                  |«      t        j
                  j	                  |«      d¬«       t        t        j                  |«      t        j
                  j	                  |«      d¬«       Œ– y rM  )r  r  r   r   Útvarr-   r@  s         r3   Ú	test_tvarzTestCompareWithStats.test_tvar   s{   € Ø—‘“ò 	,ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§
¡
¨1£¬u¯|©|×/@Ñ/@ÀÓ/DØ(*õ,ä¤§
¡
¨1£¬u¯|©|×/@Ñ/@ÀÓ/DØ(*ö,ñ		,r5   c                 ó  — t        j                  d«      }t        j                  |d«      }t        j                  j                  |d«      }t        t        j                  |«      |j                  |j                      «       y )NrX   r÷   )	r+   r­   r   rl  r-   r   rˆ  r0   rb   )r/   rS   r  Úbms       r3   Útest_trimbothz"TestCompareWithStats.test_trimboth  sV   € Ü�I‰I�b‹MˆÜ�N‰N˜1˜cÓ"ˆÜ�\‰\×"Ñ" 1 cÓ*ˆÜœŸ™ ›
 B§G¡G¨R¯W©W¨HÑ$5Õ6r5   c                 óÖ  — | j                  «       D ]Ö  }| j                  |«      \  }}}}t        t        j                  |«      t        j
                  j	                  |«      d¬«       t        t        j                  |«      t        j
                  j	                  |«      d¬«       t        t        j                  |d¬«      t        j
                  j	                  |d¬«      d¬«       ŒØ y )Nrk  rü  )r  r   )r]  )r  r  r   r   Útsemr-   r@  s         r3   Ú	test_tsemzTestCompareWithStats.test_tsem  sª   € Ø—‘“ò 	,ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ¤§
¡
¨1£¬u¯|©|×/@Ñ/@ÀÓ/DØ(*õ,ä¤§
¡
¨1£¬u¯|©|×/@Ñ/@ÀÓ/DØ(*õ,ä¤§
¡
¨1°YÔ ?Ü %§¡× 1Ñ 1°"¸YÐ 1Ó GØ(*ö,ñ	,r5   c                 óæ   — | j                  «       D ]^  }|dkD  sŒ	| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        ||«       Œ` y )Nr�   )r  r  r   rv  r-   r   r'  s           r3   Útest_skewtestz"TestCompareWithStats.test_skewtest  sb   € à—‘“ò 	'ˆAØ�1‹uØ#×6Ñ6°qÓ9‘��1�b˜"Ü—N‘N 1Ó%�Ü—\‘\×*Ñ*¨2Ó.�Ü  2Õ&ñ	'r5   c                 ó~   — t        j                  d«      dz  }t        j                  |«      }d}t	        ||d¬«       y rŒ  )r+   r,   r-   rv  r   r�  s       r3   Útest_skewtest_result_attributesz4TestCompareWithStats.test_skewtest_result_attributes"  s5   € Ü�H‰HÐ+Ó,¨aÑ/ˆÜ�o‰o˜aÓ ˆØ,ˆ
Ü˜C °Ö5r5   c                 ó  — t         j                  j                  d«      dz  }t        j                  |«      }t        j                  j                  |«      }t        t        j                  |«      t        j                  |«      «       y )N)rX   rK   g      4@)r+   rû   r   rv  r-   r   r“   )r/   rB   r·   r   s       r3   Útest_skewtest_2D_notmaskedz/TestCompareWithStats.test_skewtest_2D_notmasked(  sW   € ä�I‰I×Ñ˜WÓ%¨Ñ+ˆÜ�N‰N˜1ÓˆÜ�\‰\×"Ñ" 1Ó%ˆÜœŸ
™
 1›¤r§z¡z°"£~Õ6r5   c                 ó<  — d}| j                  «       D ]‡  }|dkD  sŒ	| j                  ||«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        |d   d   |d   d   d¬«       t        |d   d   |d   d   d¬«       Œ‰ y )NrK   r�   r   rP   r=   r   )r  r  r   rv  r-   r   )	r/   r  rÿ   rB   rÂ   rƒ  r  r·   r   s	            r3   Útest_skewtest_2D_WithMaskz.TestCompareWithStats.test_skewtest_2D_WithMask/  s—   € ØˆØ—‘“ò 	?ˆAØ�1‹uØ#×8Ñ8¸¸BÓ?‘��1�b˜"Ü—N‘N 1Ó%�Ü—\‘\×*Ñ*¨2Ó.�ä  !¡ Q¡¨¨A©¨q©¸Õ>Ü  !¡ Q¡¨¨A©¨q©¸Ö>ñ	?r5   c           	      ó   — t        j                  d¬«      5  t        «       5 }|j                  t        d«       |j                  t        d«       | j                  «       D ]„  }|dkD  sŒ	| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        t        j                  |«      t        j                  |«      «       Œ† 	 d d d «       d d d «       y # 1 sw Y   ŒxY w# 1 sw Y   y xY w)NrŸ  )Úoverz(`kurtosistest` p-value may be inaccuratez!kurtosistest only valid for n>=20r�   )r+   rq   r   r?  ÚUserWarningr  r  r   ru  r-   r   r“   )	r/   r@  rÿ   rB   rÂ   rƒ  r  r·   r   s	            r3   Útest_normaltestz$TestCompareWithStats.test_normaltest:  sÚ   € Ü�[‰[˜gÔ&ñ 	CÔ(9Ó(;ð 	C¸sØ�J‰J”{Ð$NÔOØ�J‰J”{Ð$GÔHØ—Z‘Z“\ò C�Ø�q“5Ø#'×#:Ñ#:¸1Ó#=‘L�A�q˜"˜bÜ×(Ñ(¨Ó+�AÜŸ™×0Ñ0°Ó4�BÜ#¤B§J¡J¨q£M´2·:±:¸b³>ÕBñC÷	C÷ 	C÷ 	Cñ 	Cú÷ 	Cð 	Cús*   —D¢AC8Á(A>C8Ã'DÃ8D	Ã=DÄDc                 ó¤  — t        j                  g d¢«      j                  d«      }t        j                  g d¢«      j                  d«      }|dk(  }t         j                  j	                  ||¬«      }|j                  «       |j                  «       }}t        j                  |d¬«      \  }}||dkD     ||dkD     f}	t        j                  j                  |«      }
t        |	|
«       t        ||«       t        ||«       t        j                  j                  g «      \  }}t        |t        j                  d	t         j                  ¬
«      «       y )N)r   r   rK   rK   rL   rL   rL   rM   rM   rM   rM   r  )r   r   rK   rK   rL   rL   rL   rM   rM   rM   rM   r™   r™   r™   r™   r  ra   T)Úreturn_countsr   r   r’   )r+   r“   rx   re   r,   rà  Úuniquer   r-   Úfind_repeatsr
   Úintp)r/   rB   Útmprb   rƒ  Úx_origÚxm_origre  Úunique_countsr·   r   rX  r˜  s                r3   Útest_find_repeatsz&TestCompareWithStats.test_find_repeatsE  s  € Ü�J‰JÒ8Ó9×@Ñ@ÀÓIˆÜ�j‰jÒFÓG×NÑNÈwÓWˆØ�r‘	ˆÜ�U‰U�[‰[˜ 4ˆ[Ó(ˆØŸ&™&›( B§G¡G£I�ˆä "§	¡	¨!¸4Ô @Ñˆ�Ø�= 1Ñ$Ñ% }°]ÀQÑ5FÑ'GÐGˆÜ�\‰\×&Ñ& rÓ*ˆä�Q˜ÔÜ�Q˜ÔÜ�R˜Ô!ô —L‘L×-Ñ-¨bÓ1‰	ˆˆ6Ü�VœRŸX™X a¬r¯w©wÔ7Õ8r5   c                 ó  — | j                  «       D ]v  }| j                  |«      \  }}}}t        j                  ||«      }t        j                  j                  ||«      }t        |d   |d   d¬«       t        |d   |d   d¬«       Œx y )Nr   rW   rü  r   rš   )r  r  r   r  r-   r   r'  s           r3   r  z$TestCompareWithStats.test_kendalltauY  s   € Ø—‘“ò 	8ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ× Ñ   AÓ&ˆAÜ—‘×(Ñ(¨¨RÓ0ˆBÜ  !¡ b¨¡e°RÕ8Ü  !¡ b¨¡e°QÖ7ñ	8r5   c           	      ó  — | j                  «       D ]n  }| j                  |«      \  }}}}t        j                  |«      }t        j                  j                  |«      }t        |j                  |dt        |«       «       Œp y )Nr   )r  r  r   r  r-   r   r)  r_  r'  s           r3   r
  z)TestCompareWithStats.test_obrientransforma  sk   € Ø—‘“ò 	3ˆAØ×2Ñ2°1Ó5‰LˆAˆq�"�bÜ×%Ñ% aÓ(ˆAÜ—‘×-Ñ-¨bÓ1ˆBÜ §¡ R¨¬#¨a«& \Õ2ñ		3r5   c           
      óÆ  — dD �]N  }t        «       5  dD �]1  }| j                  «       D �]  }| j                  |«      \  }}}}t        j                  |t        j
                  j                  ||¬«      }t        j                  j	                  |t        j
                  j                  ||¬«      }	t        t        j                  |«      t        j                  |	«      «       t        j                  |t        j
                  j                  ||¬«      }
t        t        j                  |«      t        j                  |
«      «       �Œ �Œ4 	 ddd«       �ŒQ y# 1 sw Y   �Œ]xY w)zFChecks that mstats.ks_1samp and stats.ks_1samp agree on masked arrays.©Úautor.  Úasymp©rö   ró   r0  ©rõ   rÓ  N)r   r  r  r   Úks_1samprT  Úcdfr-   r
   r+   r“   ©r/   rÓ  rõ   rÿ   rB   rÂ   rƒ  r  r·  r¸  rÆ  s              r3   Útest_ks_1sampz"TestCompareWithStats.test_ks_1samph  s  € à.ó 	IˆDÜ"Ó$ñ IØ#Có 
I�KØ!ŸZ™Z›\ó 	I˜Ø'+×'>Ñ'>¸qÓ'A™˜˜1˜b "Ü$Ÿ~™~¨a´·±·±Ø:EÈDô R˜ä$Ÿ|™|×4Ñ4°R¼¿¹¿¹ØALÐSWð  5ó  Y˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÔHÜ$Ÿ~™~¨b´%·*±*·.±.Ø:EÈDô R˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÖHò	Iñ
I÷Iñ Iñ	I÷Iñ Iús   ‘D8EÅE 	c           
      ó.  — dD �]  }t        «       5  dD ]æ  }| j                  «       D ]Ñ  }| j                  |«      \  }}}}t        j                  |d||¬«      }t        j
                  j	                  |d||¬«      }	t        t        j                  |«      t        j                  |	«      «       t        j                  |d||¬«      }
t        t        j                  |«      t        j                  |
«      «       ŒÓ Œè 	 ddd«       �Œ y# 1 sw Y   �ŒxY w)z]
        Checks that 1-sample mstats.kstest and stats.kstest agree on masked arrays.
        rp  rs  rT  rt  N)	r   r  r  r   Úkstestr-   r
   r+   r“   rw  s              r3   Útest_kstest_1sampz&TestCompareWithStats.test_kstest_1sampx  s  € ð /ó 	IˆDÜ"Ó$ñ IØ#Cò 
I�KØ!ŸZ™Z›\ò 	I˜Ø'+×'>Ñ'>¸qÓ'A™˜˜1˜b "Ü$Ÿ|™|¨A¨vØ8CÈ$ô P˜ä$Ÿ|™|×2Ñ2°2°vØ?JÐQUð  3ó  W˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÔHÜ$Ÿ|™|¨B°Ø8CÈ$ô P˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÕHñ	Iñ
I÷Iñ Iñ	I÷Iñ Iús   ‘C,D
Ä
D	c           
      óf  — dD �]  }t        «       5 }|dv rd}|j                  t        |«       dD ]æ  }| j                  «       D ]Ñ  }| j	                  |«      \  }}}}	t        j                  ||||¬«      }
t
        j                  j                  ||	||¬«      }t        t        j                  |
«      t        j                  |«      «       t        j                  ||||¬«      }t        t        j                  |
«      t        j                  |«      «       ŒÓ Œè 	 ddd«       �Œ! y# 1 sw Y   �Œ-xY w)zVChecks that mstats.ks_2samp and stats.ks_2samp agree on masked arrays.
        gh-8431rp  ©rq  r.  ú)ks_2samp: Exact calculation unsuccessful.rs  rt  N)r   r?  r±   r  r  r   r  r-   r
   r+   r“   ©r/   rÓ  r@  Úmessagerõ   rÿ   rB   rÂ   rƒ  r  r·  r¸  rÆ  s                r3   r  z"TestCompareWithStats.test_ks_2sampŠ  s$  € ð /ó 	IˆDÜ"Ó$ð I¨ØÐ,Ñ,ØI�GØ—J‘Jœ~¨wÔ7Ø#Cò 
I�KØ!ŸZ™Z›\ò 	I˜Ø'+×'>Ñ'>¸qÓ'A™˜˜1˜b "Ü$Ÿ~™~¨a°Ø:EÈDô R˜ä$Ÿ|™|×4Ñ4°R¸ØALÐSWð  5ó  Y˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÔHÜ$Ÿ~™~¨b°!Ø:EÈDô R˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÕHñ	Iñ
I÷	Iñ Iñ	I÷Iñ Iúó   ‘DD&Ä&D0	c           
      óf  — dD �]  }t        «       5 }|dv rd}|j                  t        |«       dD ]æ  }| j                  «       D ]Ñ  }| j	                  |«      \  }}}}	t        j                  ||||¬«      }
t
        j                  j                  ||	||¬«      }t        t        j                  |
«      t        j                  |«      «       t        j                  ||||¬«      }t        t        j                  |
«      t        j                  |«      «       ŒÓ Œè 	 ddd«       �Œ! y# 1 sw Y   �Œ-xY w)z]
        Checks that 2-sample mstats.kstest and stats.kstest agree on masked arrays.
        rp  r}  r~  rs  rt  N)r   r?  r±   r  r  r   rz  r-   r
   r+   r“   r  s                r3   Útest_kstest_2sampz&TestCompareWithStats.test_kstest_2sampž  s$  € ð /ó 	IˆDÜ"Ó$ð I¨ØÐ,Ñ,ØI�GØ—J‘Jœ~¨wÔ7Ø#Cò 
I�KØ!ŸZ™Z›\ò 	I˜Ø'+×'>Ñ'>¸qÓ'A™˜˜1˜b "Ü$Ÿ|™|¨A¨qØ8CÈ$ô P˜ä$Ÿ|™|×2Ñ2°2°rØ?JÐQUð  3ó  W˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÔHÜ$Ÿ|™|¨B°Ø8CÈ$ô P˜ä$¤R§Z¡Z°Ó%5´r·z±zÀ$Ó7GÕHñ	Iñ
I÷	Iñ Iñ	I÷Iñ Iúr�  N))r6   r7   r8   r  r  r  r  r  r»   rñ   r%  r(  r+  r-  r³  r  rÏ  r÷  r9  r;  r=  rA  rE  rI  r   r·  rP  rS  rV  rX  rZ  r\  r^  rb  rl  r  r
  rx  r{  r  rƒ  r9   r5   r3   r  r    sË   „ ñò,"òòòFò3ò+òAò	/ò
+ò
+ò
+ò(ò$	+òMò*0ò9ò#òKòJò	Kò9ò,ò,ò7ò	,ò'ò6ò7ò	?ò	Cò9ò(8ò3òIò Iò$Ió(Ir5   r  c                   ó  — e Zd Zej                  j                  dddddej                  dddddddddej                  g«      Zej                  j                  ddddej                  dddddddg«      ZdZ	d„ Z
d„ Zd	„ Zd
„ Zd„ Zd„ Zd„ Zy)ÚTestBrunnerMunzelr   rK   rM   rL   r™   rk  c                 ó  — t        j                  | j                  | j                  d¬«      \  }}t        j                  | j                  | j                  d¬«      \  }}t        j                  | j                  | j                  d¬«      \  }}t        j                  | j                  | j                  d¬«      \  }}t	        ||| j
                  ¬«       t	        ||| j
                  ¬«       t        ||k7  «       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       y )	Nrö   rô   ró   rü  ç|	&�ˆ	@ç|	&�ˆ	Àg°¿Ò�G³g?g@-p¸Lèï?)r-   ÚbrunnermunzelÚXÚYr   Úsignificantr   )	r/   Úu1r  Úu2r  Úu3r  Úu4Úp4s	            r3   Útest_brunnermunzel_one_sidedz.TestBrunnerMunzel.test_brunnermunzel_one_sided»  s6  € ä×%Ñ% d§f¡f¨d¯f©fÀ&ÔI‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÀ)ÔL‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÀ)ÔL‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÀ&ÔI‰ˆˆBä˜B ¨D×,<Ñ,<Õ=Ü˜B ¨D×,<Ñ,<Õ=Ü��b‘ÔÜ˜BÐ 2Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4õ	6ä˜BÐ 2Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4õ	6ä˜BÐ 5Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4ö	6r5   c                 ó€  — t        j                  | j                  | j                  d¬«      \  }}t        j                  | j                  | j                  d¬«      \  }}t	        ||| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       y )Nr0  rô   rü  r‡  rˆ  ç ÀÒ�G³w?©r-   r‰  rŠ  r‹  r   rŒ  ©r/   r�  r  rŽ  r  s        r3   Útest_brunnermunzel_two_sidedz.TestBrunnerMunzel.test_brunnermunzel_two_sidedÒ  s–   € ä×%Ñ% d§f¡f¨d¯f©fÀ+ÔN‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÀ+ÔN‰ˆˆBä˜B ¨D×,<Ñ,<Õ=Ü˜BÐ 2Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4õ	6ä˜BÐ 5Ø$(×$4Ñ$4ö	6r5   c                 óx  — t        j                  | j                  | j                  «      \  }}t        j                  | j                  | j                  «      \  }}t	        ||| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       y )Nrü  r‡  rˆ  r”  r•  r–  s        r3   Útest_brunnermunzel_defaultz,TestBrunnerMunzel.test_brunnermunzel_defaultß  s’   € ä×%Ñ% d§f¡f¨d¯f©fÓ5‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÓ5‰ˆˆBä˜B ¨D×,<Ñ,<Õ=Ü˜BÐ 2Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4õ	6ä˜BÐ 5Ø$(×$4Ñ$4ö	6r5   c                 ó’   — d}d}t        |dv«       t        t        t        j                  | j
                  | j                  ||«       y )Nrª   rÌ  r/  ©r   r  rý   r-   r‰  rŠ  r‹  ©r/   rõ   Údistributions      r3   Ú$test_brunnermunzel_alternative_errorz6TestBrunnerMunzel.test_brunnermunzel_alternative_errorì  sA   € ØˆØˆÜ�Ð#CÐCÔDÜ”jÜ×*Ñ*Ø—f‘fØ—f‘fØ!Ø"õ	$r5   c                 ó€  — t        j                  | j                  | j                  d¬«      \  }}t        j                  | j                  | j                  d¬«      \  }}t	        ||| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       t	        |d| j
                  ¬«       y )NÚnormal)r�  rü  r‡  rˆ  g œ"H°ë[?r•  r–  s        r3   Ú$test_brunnermunzel_distribution_normz6TestBrunnerMunzel.test_brunnermunzel_distribution_norm÷  s–   € Ü×%Ñ% d§f¡f¨d¯f©fÀ8ÔL‰ˆˆBÜ×%Ñ% d§f¡f¨d¯f©fÀ8ÔL‰ˆˆBÜ˜B ¨D×,<Ñ,<Õ=Ü˜BÐ 2Ø$(×$4Ñ$4õ	6ä˜BÐ 3Ø$(×$4Ñ$4õ	6ä˜BÐ 5Ø$(×$4Ñ$4ö	6r5   c                 ó’   — d}d}t        |dv«       t        t        t        j                  | j
                  | j                  ||«       y )Nr0  rª   )rÌ  r   r›  rœ  s      r3   Ú%test_brunnermunzel_distribution_errorz7TestBrunnerMunzel.test_brunnermunzel_distribution_error  s@   € Ø!ˆØˆÜ� ?Ð2Ô3Ü”jÜ×*Ñ*Ø—f‘fØ—f‘fØ!Ø"õ	$r5   c                 ó*  — t        j                  | j                  g «      \  }}t        j                  g | j                  «      \  }}t        j                  g g «      \  }}t	        t        j                  |«      «       t	        t        j                  |«      «       t	        t        j                  |«      «       t	        t        j                  |«      «       t	        t        j                  |«      «       t	        t        j                  |«      «       y r¦  )r-   r‰  rŠ  r‹  r   r+   r>  )r/   r�  r  rŽ  r  r�  r  s          r3   Útest_brunnermunzel_empty_imputz0TestBrunnerMunzel.test_brunnermunzel_empty_imput  s¬   € Ü×%Ñ% d§f¡f¨bÓ1‰ˆˆBÜ×%Ñ% b¨$¯&©&Ó1‰ˆˆBÜ×%Ñ% b¨"Ó-‰ˆˆBä”—‘˜“ÔÜ”—‘˜“ÔÜ”—‘˜“ÔÜ”—‘˜“ÔÜ”—‘˜“ÔÜ”—‘˜“Õr5   N)r6   r7   r8   r+   re   Úmasked_invalidr   rŠ  r‹  rŒ  r’  r—  r™  rž  r¡  r£  r¥  r9   r5   r3   r…  r…  ´  s�   „ à
�‰×Ñ˜a  A q¨!¨R¯V©V°Q¸Ø  A q¨!¨Q°°2·6±6ð;ó 	<€Aà
�‰×Ñ˜a  A q¨"¯&©&°!°Q¸¸1¸aÀÀAÐFÓG€AØ€Kò6ò.6ò6ò	$ò	6ò	$ó
r5   r…  )NNgH¯¼šò×z>)Cr  r®   rR  Únumpyr+   r   Únumpy.mare   r   r   Úscipy.stats.mstatsr   r-   ÚscipyÚcommon_testsr   rˆ   r	   r  Únumpy.ma.testutilsr
   r   r   r   r   r   r   Únumpy.testingr   Úscipy.statsr   r   Úscipy.conftestr   Úscipy.stats._axis_nan_policyr   r   r   rC   rF   rH   r�   r—   r¨   rX  r£  rÑ  rå  rô  r  r*  r-  r/  rY  r[  rl  rp  rr  r¬  r³  r»  r¿  rì  rú  r  r  r…  r9   r5   r3   ú<module>r±     sú  ðñó Û ã Ý Ý ß #ç #Ð #Ý Ý -Û Ý *÷E÷ Eñ Eõ ,ß 0Ý .ß R÷-ñ -ó*!ó!ð ÷93ð 93ó ð93ðx ÷!3ð !3ó ð!3÷H9ñ 9÷,y6ñ y6ðx	 ÷u:ð u:ó ðu:ðp ÷V/ð V/ó ðV/÷r&ñ &÷RAñ Að* ÷4ð 4ó ð4ð@ ÷-9ð -9ó ð-9ò`Fò, ÷J1ñ J1òZ8ò:.ò8ò*Eð ÷_8ð _8ó ð_8÷D6ñ 6÷.2ñ .2÷b6ñ 6÷Y-ñ Y-÷x\-ñ \-÷~@-ñ @-÷F7ñ 7ð0 ÷mIð mIó ðmI÷`cò cr5   