§
    fŠtj‹{ ã                   óª  — d dl mZ d dlZd dlZd dlZd dl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mZmZmZmZmZ d dlmZmZmZ d d	lmZ d d
lm Z! d dl"m#Z#m$Z$m%Z%m&Z& d dl'm(Z(m)Z) d dl*m+Z+m,Z,  e#e¦  «         G d„ d¦  «        ¦   «         Z- e#ej        ¦  «         G d„ d¦  «        ¦   «         Z. e#ej        ¦  «         G d„ d¦  «        ¦   «         Z/ G d„ de¦  «        Z0 G d„ d¦  «        Z1 G d„ d¦  «        Z2 e#e¦  «         G d„ d¦  «        ¦   «         Z3 G d„ d¦  «        Z4 G d„ d¦  «        Z5 G d „ d!¦  «        Z6 G d"„ d#¦  «        Z7dS )$é    )ÚproductN)Úassert_Úassert_equalÚassert_allclose)Úraises)ÚstatsÚspecial)Údistributions)Úepps_singleton_2sampÚcramervonmisesÚ_cdf_cvmÚcramervonmises_2sampÚ_pval_cvm_2samp_exactÚbarnard_exactÚboschloo_exact)ÚmannwhitneyuÚ
_mwu_stateÚ_MWU)Ú_TestPythranFunc)Úarray_api_extra)Úmake_xp_test_caseÚxp_default_dtypeÚis_numpyÚeager_warns)Úxp_assert_equalÚxp_assert_close)ÚSmallSampleWarningÚtoo_small_1d_not_omitc                   ór   — e Zd Zej                             dg d¢¦  «        d„ ¦   «         Zd„ Zd„ Zd„ Z	d„ Z
dS )	ÚTestEppsSingletonÚdtype©NÚfloat32Úfloat64c                 óÖ  — t          |¦  «        r%|j        dk     r|dk    rt          j        d¦  «         |€t	          |¦  «        nt          ||¦  «        }|                     g d¢|¬¦  «        }|                     g d¢|¬¦  «        }t          ||¦  «        \  }}t          ||                     d|¬¦  «        d¬	¦  «         t          ||                     d
|¬¦  «        d¬	¦  «         d S )Nú2.0r#   ú!Pre-NEP 50 doesn't respect dtypes)
gffffffÖ¿gffffff@g®Gáz®û?ç\�Âõ(\ç?çffffffÖ?g…ëQ¸…@gq=
×£pÝ?g®Gázî¿g®Gáz®×¿g¤p=
×#(@©r!   )
gffffffò¿g333333Ã¿g×£p=
×@g      
@g®Gáz®@g)\�Âõ(@ç      @gö(\�Âõ@gÃõ(\�Â @ç333333!@gHáz®G.@g¸…ëQ¸ž?©Úatolg˜Q,·´r?ç-Cëâ6?)	r   Ú__version__ÚpytestÚskipr   ÚgetattrÚasarrayr   r   )Úselfr!   ÚxpÚxÚyÚwÚps          ú^/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/stats/tests/test_hypotests.pyÚtest_statistic_1z"TestEppsSingleton.test_statistic_1   s  € õ �B‰<Œ<ð 	=˜BœN¨UÒ2Ð2°uÀ	Ò7IÐ7IÝŒKÐ;Ñ<Ô<Ð<Ø(-¨Õ  Ñ$Ô$Ð$½7À2ÀuÑ;MÔ;MˆØ�JŠJð 9ð 9ð 9Ø@Eð ñ Gô Gˆà�JŠJð 3ð 3ð 3Ø:?ð ñ Aô Aˆå# A qÑ)Ô)‰ˆˆ1Ý˜˜2Ÿ:š: e°5˜:Ñ9Ô9ÀÐEÑEÔEÐEÝ˜˜2Ÿ:š: g°U˜:Ñ;Ô;À&ÐIÑIÔIÐIÐIÐIó    c                 ó  — |                      g d¢¦  «        }|                      g d¢¦  «        }t          ||¦  «        \  }}t          ||                      d¦  «        d¬¦  «         t          ||                      d¦  «        d¬¦  «         d S )N)r   é   é   r@   r@   r@   é   rA   rA   rA   é   é   rC   rC   rC   é   é
   rE   rE   rE   )rE   rB   r   rC   rE   rE   r   rC   rD   é   rE   rA   r?   rF   r   é   r?   rC   rG   rE   gÍÌÌÌÌÌ!@çü©ñÒMbP?r-   g&ª·¶J°?g-Cëâ6
?)r4   r   r   )r5   r6   r7   r8   r9   r:   s         r;   Útest_statistic_2z"TestEppsSingleton.test_statistic_2+   s¡   € à�JŠJð 7ð 7ð 7ñ 8ô 8ˆà�JŠJð 8ð 8ð 8ñ 9ô 9ˆå# A qÑ)Ô)‰ˆˆ1Ý˜˜2Ÿ:š: eÑ,Ô,°4Ð8Ñ8Ô8Ð8Ý˜˜2Ÿ:š: gÑ.Ô.°TÐ:Ñ:Ô:Ð:Ð:Ð:r=   c                 ó¤  — t          j        d¦  «        t          j        d¦  «        }}t          t          |¦  «        t          |¦  «        ¦  «        \  }}t          t	          |¦  «        t	          |¦  «        ¦  «        \  }}t          ||¦  «        \  }}t          ||cxk    o|k    nc ¦  «         t          ||cxk    o|k    nc ¦  «         d S )Né   é   )ÚnpÚaranger   ÚlistÚtupler   )	r5   r7   r8   Úw1Úp1Úw2Úp2Úw3Úp3s	            r;   Útest_epps_singleton_array_likez0TestEppsSingleton.test_epps_singleton_array_like5   s¶   € ÝŒy˜‰}Œ}�bœi¨™mœmˆ1ˆå%¥d¨1¡g¤g­t°A©w¬wÑ7Ô7‰ˆˆBÝ%¥e¨A¡h¤hµ°a±´Ñ9Ô9‰ˆˆBÝ% a¨Ñ+Ô+‰ˆˆBå��b��’�˜B’���ÑÔÐÝ��b��’�˜B’���ÑÔÐÐÐr=   c                 óœ  — |                      g d¢¦  «        |                     d¦  «        }}t          t          t          |¬¦  «        5  t          ||¦  «        }t          |j        |                      |j        ¦  «        ¦  «         t          |j	        |                      |j        ¦  «        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )N©r?   r@   rA   rB   rE   ©Úmatchr6   )
r4   rN   r   r   r   r   r   Ú	statisticÚnanÚpvalue©r5   r6   r7   r8   Úress        r;   Útest_epps_singleton_sizez*TestEppsSingleton.test_epps_singleton_size?   sï   € à�zŠz˜,˜,˜,Ñ'Ô'¨¯ª°2©¬ˆ1ˆÝÕ+Õ3HÈRÐPÑPÔPð 	<ð 	<Ý& q¨!Ñ,Ô,ˆCÝ˜CœM¨2¯:ª:°b´fÑ+=Ô+=Ñ>Ô>Ð>Ý˜CœJ¨¯
ª
°2´6Ñ(:Ô(:Ñ;Ô;Ð;ð	<ð 	<ð 	<ñ 	<ô 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<øøøð 	<ð 	<ð 	<ð 	<ð 	<ð 	<s   Á	A+CÃCÃCc                 óÄ  — t           j                             d¦  «        }|                     d¬¦  «        }|                     d¬¦  «        }t          j        g d¢¦  «        }t	          ||d¬¦  «        \  }}t           j        ||<   t           j        ||<   t           j        ||d         df<   t           j        ||d	         d	f<   t           j         ||d
         d
f<   |                     |¦  «        |                     |¦  «        }}t	          ||d¬¦  «        \  }}	t          ||                     |¦  «        ¦  «         t          |	|                     |¦  «        ¦  «         d S )Nl   2$>Ü. )rE   é   ©Úsize)rE   é   )r?   rB   é	   éÿÿÿÿ©Úaxisr   r?   r@   )rM   ÚrandomÚdefault_rngr4   r   r]   Úinfr   )
r5   r6   Úrngr7   r8   ÚiÚw_refÚp_refÚw_resÚp_ress
             r;   Útest_epps_singleton_nonfinitez/TestEppsSingleton.test_epps_singleton_nonfiniteG   s,  € ÝŒi×#Ò# NÑ3Ô3ˆØ�JŠJ˜HˆJÑ%Ô%ˆØ�JŠJ˜HˆJÑ%Ô%ˆÝŒJ�y�y�yÑ!Ô!ˆå+¨A¨q°rÐ:Ñ:Ô:‰ˆˆuÝ”6ˆˆa‰Ý”6ˆˆa‰å”Vˆˆ!ˆAŒ$�ˆ'‰
Ý”Vˆˆ!ˆAŒ$�ˆ'‰
Ý”f�Wˆˆ!ˆAŒ$�ˆ'‰
à�zŠz˜!‰}Œ}˜bŸjšj¨™mœmˆ1ˆÝ+¨A¨q°rÐ:Ñ:Ô:‰ˆˆuÝ˜˜rŸzšz¨%Ñ0Ô0Ñ1Ô1Ð1Ý˜˜rŸzšz¨%Ñ0Ô0Ñ1Ô1Ð1Ð1Ð1r=   N)Ú__name__Ú
__module__Ú__qualname__r1   ÚmarkÚparametrizer<   rI   rW   ra   rt   © r=   r;   r    r       s‡   € € € € € à„[×Ò˜WÐ&BÐ&BÐ&BÑCÔCðJð Jñ DÔCðJð ;ð ;ð ;ð ð  ð  ð<ð <ð <ð2ð 2ð 2ð 2ð 2r=   r    c            
       ó’  — e Zd Zej                             ddg d¢g d¢fdg d¢g d¢fdg d	¢g d¢fd
g d¢g d¢fg¦  «        d„ ¦   «         Zd„ Zd„ Zd„ Z	d„ Z
ej                             dd¬¦  «        ej                             dddgg¦  «        d„ ¦   «         ¦   «         Zej                             dg d¢¦  «        d„ ¦   «         Zd„ Zd
S )ÚTestCvmz	n, x, refrB   )gy;ÂiÁ‹ž?g#¾³^¥?gþE�>‘¿?gåD»
)î?)ç{®Gáz„?çš™™™™™©?ç      à?g+‡ÙÎ÷ï?rE   )g…8„*5›?g@¤ß¾œ£?g÷Hmâä¾?g%¯Î1 â?)r}   r~   r   g333333ï?éè  )g€}têÊg™?g˜`�º¢?gäIÒ5“o¾?gÜ×�sF”ò?N)çaÃÓ+e™?ç+û®þ·¢?çÉ&pën¾?ç+MJA·—ò?c                 ó’   — t          t          |                     |¦  «        |¦  «        |                     |¦  «        d¬¦  «         d S ©Nr/   r-   )r   r   r4   )r5   Únr7   Úrefr6   s        r;   Útest_cdf_refzTestCvm.test_cdf_ref`   s=   € õ 	� §¢¨A¡¤°Ñ2Ô2°B·J²J¸s±O´OÈ$ÐOÑOÔOÐOÐOÐOr=   c                 ó&  — t          t          |                     ddg¦  «        d¦  «        |                     ddg¦  «        ¦  «         t          t          |                     ddg¦  «        d¦  «        |                     ddg¦  «        ¦  «         d S )	Ng²Xª(~$?gUUUUU5f@i  ç        ç      ð?g†a†ah?g«ªªªªª"@é   )r   r   r4   )r5   r6   s     r;   Útest_cdf_supportzTestCvm.test_cdf_supporti   s�   € å� §¢¨Z¸Ð,?Ñ!@Ô!@À#ÑFÔFØŸ
š
 B¨ 8Ñ,Ô,ñ	.ô 	.ð 	.å� §¢¨_¸jÐ,IÑ!JÔ!JÈBÑOÔOØŸ
š
 B¨ 8Ñ,Ô,ñ	.ô 	.ð 	.ð 	.ð 	.r=   c                 óŽ   — |                      g d¢¦  «        }t          t          |d¦  «        t          |¦  «        d¬¦  «         d S )N)r�   r‚   rƒ   r„   éd   é'  r/   r-   )r4   r   r   ©r5   r6   r7   s      r;   Útest_cdf_large_nzTestCvm.test_cdf_large_np   sE   € à�JŠJÐ@Ð@Ð@ÑAÔAˆÝ�  EÑ*Ô*­H°Q©K¬K¸dÐCÑCÔCÐCÐCÐCr=   c                 ó¸   — |                      d|j        ¬¦  «        }dt          |d¦  «        cxk     rdk     sn J ‚dt          |¦  «        cxk     rdk     sn J ‚d S )NgÍÌÌÌÌÔt@r*   gwJëÿï?r€   rŒ   )r4   r$   r   r’   s      r;   Útest_large_xzTestCvm.test_large_xu   s{   € ð �JŠJ�u B¤JˆJÑ/Ô/ˆØ�( 1 dÑ+Ô+Ð1Ð1Ò1Ð1¨cÒ1Ð1Ð1Ð1Ð1Ð1Ø�( 1™+œ+Ð+Ð+Ò+Ð+¨Ò+Ð+Ð+Ð+Ð+Ð+Ð+Ð+r=   c                 óô   — d}t          |                     |¦  «        dz  t          j        ¦  «        }t	          |j        ||¬¦  «        dk    sJ ‚t          |j        |                     d¦  «        ¦  «         d S )Nrf   çš™™™™™é?)r6   rŒ   r‹   )	r   Úonesr	   Úndtrr   r\   r   r^   r4   )r5   r6   r‡   r`   s       r;   Ú
test_low_pzTestCvm.test_low_p   sl   € ð ˆÝ˜RŸWšW Q™ZœZ¨™^­W¬\Ñ:Ô:ˆÝ˜œ q¨RÐ0Ñ0Ô0°3Ò6Ð6Ð6Ð6Ý˜œ
 B§J¢J¨r¡N¤NÑ3Ô3Ð3Ð3Ð3r=   ú	jax.numpyúlazy -> no _axis_nan_policy©Úreasonr7   rz   ç      ø?c                 ó†  — t          j        t          t          ¬¦  «        5  t	          |                     |¦  «        t          j        ¦  «        }t          |j	        |                     |j
        ¦  «        ¦  «         t          |j        |                     |j
        ¦  «        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )N©r[   )r1   Úwarnsr   r   r   r4   r	   r™   r   r\   r]   r^   )r5   r7   r6   r`   s       r;   Útest_invalid_inputzTestCvm.test_invalid_input‡   sÚ   € õ Œ\Õ,Õ4IÐJÑJÔJð 	<ð 	<Ý  §¢¨A¡¤µ´Ñ=Ô=ˆCÝ˜CœM¨2¯:ª:°b´fÑ+=Ô+=Ñ>Ô>Ð>Ý˜CœJ¨¯
ª
°2´6Ñ(:Ô(:Ñ;Ô;Ð;ð	<ð 	<ð 	<ñ 	<ô 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<øøøð 	<ð 	<ð 	<ð 	<ð 	<ð 	<s   ¡BB6Â6B:Â=B:r!   r"   c                 ó°  ‡— t          ‰¦  «        r%‰j        dk     r|dk    rt          j        d¦  «         |€t	          ‰¦  «        nt          ‰|¦  «        }t          ‰                     g d¢|¬¦  «        t          j	        ¦  «        }t          |j        ‰                     d|¬¦  «        ¦  «         t          |j        ‰                     d|¬¦  «        ¦  «         t          ‰                     g d¢|¬¦  «        d„ ¦  «        }t          |j        ‰                     d	|¬¦  «        ¦  «         t          |j        ‰                     d
|¬¦  «        ¦  «         t          ‰                     g d¢|¬¦  «        ˆfd„¦  «        }t          |j        ‰                     d|¬¦  «        ¦  «         t          |j        ‰                     d|¬¦  «        ¦  «         d S )Nr&   r#   r'   )g333333û¿r@   r   gÍÌÌÌÌÌô?rB   çš™™™™™¹?ç333333ã?r*   gÓˆk&qÒ?gDC+·šÂ?c                 ó6   — t          j        | dz
  dz  ¦  «        S )NrA   rŸ   )r	   r™   )r7   s    r;   ú<lambda>z'TestCvm.test_values_R.<locals>.<lambda>Ÿ   s   € ¥w¤|°Q¸±U¸C±KÑ'@Ô'@€ r=   g�ºrW*î?g¤ju2´™`?)	r?   r@   rC   çffffffö?gìQ¸…ëÁ?rc   é   çÍÌÌÌÌÌì?ç      @c                 ó2   •— ‰                      |  ¦  «         S ©N)Úexpm1)r7   r6   s    €r;   r¨   z'TestCvm.test_values_R.<locals>.<lambda>¦   s   ø€ �r—x’x  ‘|”|�m€ r=   g  ¹Ð.óê?góûÎ`Ã(r?)r   r0   r1   r2   r   r3   r   r4   r	   r™   r   r\   r^   )r5   r!   r6   r`   s     ` r;   Útest_values_RzTestCvm.test_values_R�   sÊ  ø€ å�B‰<Œ<ð 	=˜BœN¨UÒ2Ð2°uÀ	Ò7IÐ7IÝŒKÐ;Ñ<Ô<Ð<Ø(-¨Õ  Ñ$Ô$Ð$½7À2ÀuÑ;MÔ;Mˆõ
 ˜RŸZšZÐ(FÐ(FÐ(FÈe˜ZÑTÔTÝ$œ\ñ+ô +ˆå˜œ r§z¢zÐ2BÈ% zÑ'PÔ'PÑQÔQÐQÝ˜œ
 B§J¢J°Àu JÑ$MÔ$MÑNÔNÐNõ ˜RŸZšZÐ(FÐ(FÐ(FÈe˜ZÑTÔTØ@Ð@ñBô Bˆå˜œ r§z¢zÐ2BÈ% zÑ'PÔ'PÑQÔQÐQÝ˜œ
 B§J¢JÐ/@È JÑ$NÔ$NÑOÔOÐOõ Ø�JŠJÐ=Ð=Ð=ÀUˆJÑKÔKØ#Ð#Ð#Ð#ñ%ô %ˆõ 	˜œ r§z¢zÐ2BÈ% zÑ'PÔ'PÑQÔQÐQÝ˜œ
 B§J¢JÐ/@È JÑ$NÔ$NÑOÔOÐOÐOÐOr=   c                 óh  — t          |¦  «        s\d}t          j        t          |¬¦  «        5  t	          |                     g d¢¦  «        d¦  «         d d d ¦  «         n# 1 swxY w Y   d S t          j        d¦  «        d}}t	          |t          j	        j
        ¦  «        }t	          |d¦  «        }t          |j        |j        f|j        |j        f¦  «         t	          |t          j        j
        |¦  «        }t	          |d|¦  «        }t          |j        |j        f|j        |j        f¦  «         d S )Nz7`cdf` must be a callable if `rvs` is a non-NumPy array.r¡   ©r?   r@   rA   ÚbetarC   )r©   çffffffæ?Úexpon)r   r1   r   Ú
ValueErrorr   r4   rM   rN   r
   rµ   Úcdfr   r\   r^   r³   )r5   r6   Úmessager7   Úargsrˆ   r`   s          r;   Útest_str_cdfzTestCvm.test_str_cdfª   sK  € Ý˜‰|Œ|ð 	ØOˆGÝ”�z°Ð9Ñ9Ô9ð >ð >Ý˜rŸzšz¨)¨)¨)Ñ4Ô4°fÑ=Ô=Ð=ð>ð >ð >ñ >ô >ð >ð >ð >ð >ð >ð >øøøð >ð >ð >ð >àˆFå”)˜A‘,”, 
ˆ4ˆÝ˜Q¥Ô 3Ô 7Ñ8Ô8ˆÝ˜Q Ñ(Ô(ˆÝ�c”m S¤ZÐ0°3´=À#Ä*Ð2MÑNÔNÐNå˜Q¥Ô 2Ô 6¸Ñ=Ô=ˆÝ˜Q ¨Ñ-Ô-ˆÝ�c”m S¤ZÐ0°3´=À#Ä*Ð2MÑNÔNÐNÐNÐNs   ­&AÁA#Á&A#)ru   rv   rw   r1   rx   ry   r‰   rŽ   r“   r•   rš   Úskip_xp_backendsr£   r°   rº   rz   r=   r;   r|   r|   [   s�  € € € € € ð
 „[×Ò˜[Ø	
Ð0Ð0Ð0Ð2JÐ2JÐ2JÐKØ	Ð1Ð1Ð1Ð3KÐ3KÐ3KÐLØ	Ð3Ð3Ð3Ð5MÐ5MÐ5MÐNØ	Ð3Ð3Ð3Ð5MÐ5MÐ5MÐNð	+ñ ô ðPð Pñô ðPð.ð .ð .ðDð Dð Dð
-ð -ð -ð4ð 4ð 4ð „[×!Ò! +Ð6SÐ!ÑTÔTØ„[×Ò˜S 2¨ u +Ñ.Ô.ð<ð <ñ /Ô.ñ UÔTð<ð „[×Ò˜WÐ&BÐ&BÐ&BÑCÔCðPð Pñ DÔCðPð4Oð Oð Oð Oð Or=   r|   c                   ó¸  — e Zd Zej                             dd¬¦  «        d„ ¦   «         Zd„ Zd„ Zg d¢Z	g d¢Z
d	d
dœdgdd
dœdgdd
dœdgd	ddœdgdddœdgdddœdggZej                             de¦  «        ej                             dg d¢¦  «        d„ ¦   «         ¦   «         Zd	ddœdgdddœdgdddœdgd	ddœdgdddœd gdddœd!ggZej                             de¦  «        d"„ ¦   «         Zd#„ Zg d$¢g d%¢g d&¢d'œZg d(¢g d)¢g d*¢g d+¢d,œZg d-¢g d.¢g d/¢g d0¢g d1¢d2œZg d3¢g d4¢g d5¢g d6¢g d7¢g d8¢d9œZd:„ Zd;„ Zd<„ Zd	d
dœd=gdd
dœd>gdd
dœd?gd	ddœd@gdddœd>gdddœd@ggZej                             dAe¦  «        dB„ ¦   «         ZdC„ Zej                             dd¬¦  «        ej                             dDd
dg¦  «        dE„ ¦   «         ¦   «         ZdF„ Zej                             dd¬¦  «        dG„ ¦   «         Zg dH¢dIdJdKdLej        dIdMdNdOdOdPgdQdRfg dH¢dIdJdKdLej        ej        dMdNdOdOdPgdSdTfdNdMej        dOgdIdJdKdLej        dIdMdNdOdOdPgdUdVfdNdMej        dOgdIdJdKdLej        ej        dMdNdOdOdPgdWdXfdNej        ej        dOgdIdJdKdLej        ej        dMdNdOdOdPgdYdZfgZ ej                             d[e ¦  «        d\„ ¦   «         Z!g d]¢g d^¢g d_¢g d`¢g da¢g db¢g dc¢g dd¢g de¢g	Z"ej                             dfe"¦  «        dg„ ¦   «         Z#ej                             dd¬¦  «        dh„ ¦   «         Z$g di¢djdkgddlgg di¢djdkgddlgg di¢djdkgd	dmgg d'¢dMgddngg d'¢dMgddngg d'¢dMgd	dogdNdMgdNdMgddpgdNdMgdNdMgddpgdNdMgdNdMgd	dqgg	Z%ej                             g dr¢e%¦  «        ds„ ¦   «         Z&dt„ Z'ej                             dug dv¢¦  «        dw„ ¦   «         Z(g dx¢Z)g dy¢Z*dzZ+d{„ Z,d|„ Z-d}„ Z.d~„ Z/d„ Z0d€S )�ÚTestMannWhitneyUr›   rœ   r�   c                 ó|  — |                      ddg¦  «        }|                      ddg¦  «        }|                      g |j        ¬¦  «        }|                      |j        ¦  «        }t          j        t
          t          ¬¦  «        5  t          ||¦  «        }t          |j	        |¦  «         t          |j
        |¦  «         d d d ¦  «         n# 1 swxY w Y   t          j        t
          t          ¬¦  «        5  t          ||¦  «        }t          |j	        |¦  «         t          |j
        |¦  «         d d d ¦  «         n# 1 swxY w Y   t          j        t
          t          ¬¦  «        5  t          ||¦  «        }t          |j	        |¦  «         t          |j
        |¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr?   r@   rA   rB   r*   r¡   )r4   r!   r]   r1   r¢   r   r   r   r   r\   r^   )r5   r6   r7   r8   Úemptyr]   r`   s          r;   Ú
test_emptyzTestMannWhitneyU.test_emptyÃ   sF  € à�JŠJ˜˜1�vÑÔˆØ�JŠJ˜˜1�vÑÔˆØ—
’
˜2 Q¤W�
Ñ-Ô-ˆØ�jŠj˜œÑ Ô ˆåŒ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý˜q %Ñ(Ô(ˆCÝ˜CœM¨3Ñ/Ô/Ð/Ý˜CœJ¨Ñ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-õ
 Œ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý˜u aÑ(Ô(ˆCÝ˜CœM¨3Ñ/Ô/Ð/Ý˜CœJ¨Ñ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-õ
 Œ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý˜u eÑ,Ô,ˆCÝ˜CœM¨3Ñ/Ô/Ð/Ý˜CœJ¨Ñ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð 	-ð 	-s6   Â;CÃCÃCÃ7;D>Ä>EÅEÅ);F1Æ1F5Æ8F5c                 ód  — |                      ddg¦  «        }|                      ddg¦  «        }t          t          d¬¦  «        5  t          ||d¬¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          d	¬¦  «        5  t          ||d¬
¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          d¬¦  «        5  t          ||d¬¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          d¬¦  «        5  t          ||d¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr?   r@   rA   rB   z`use_continuity` must be oner¡   Úekki)Úuse_continuityz`alternative` must be one of©Úalternativez`axis` must be an integerrŸ   ri   z`method` must be one of©Úmethod)r4   Úassert_raisesr¶   r   )r5   r6   r7   r8   s       r;   Útest_input_validationz&TestMannWhitneyU.test_input_validationÙ   s3  € Ø�JŠJ˜˜1�vÑÔˆØ�JŠJ˜˜1�vÑÔˆÝ�:Ð-KÐLÑLÔLð 	6ð 	6Ý˜˜A¨fÐ5Ñ5Ô5Ð5ð	6ð 	6ð 	6ñ 	6ô 	6ð 	6ð 	6ð 	6ð 	6ð 	6ð 	6øøøð 	6ð 	6ð 	6ð 	6å�:Ð-KÐLÑLÔLð 	3ð 	3Ý˜˜A¨6Ð2Ñ2Ô2Ð2ð	3ð 	3ð 	3ñ 	3ô 	3ð 	3ð 	3ð 	3ð 	3ð 	3ð 	3øøøð 	3ð 	3ð 	3ð 	3å�:Ð-HÐIÑIÔIð 	)ð 	)Ý˜˜A CÐ(Ñ(Ô(Ð(ð	)ð 	)ð 	)ñ 	)ô 	)ð 	)ð 	)ð 	)ð 	)ð 	)ð 	)øøøð 	)ð 	)ð 	)ð 	)å�:Ð-FÐGÑGÔGð 	.ð 	.Ý˜˜A fÐ-Ñ-Ô-Ð-ð	.ð 	.ð 	.ñ 	.ô 	.ð 	.ð 	.ð 	.ð 	.ð 	.ð 	.ð 	.øøøð 	.ð 	.ð 	.ð 	.ð 	.ð 	.sH   ÁA$Á$A(Á+A(ÂB$Â$B(Â+B(ÃC$Ã$C(Ã+C(ÄD%Ä%D)Ä,D)c                 óÎ  — t           j                             d¦  «        }d}|                     |                     |dz
  ¦  «        ¦  «        }|                     |                     |dz
  ¦  «        ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚|                     |                     |dz
  ¦  «        ¦  «        }|                     |                     |dz   ¦  «        ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚|                     |                     |dz   ¦  «        ¦  «        }|                     |                     |dz   ¦  «        ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚|                     |                     |dz
  ¦  «        ¦  «        }|                     |                     |dz
  ¦  «        ¦  «        }t          j        |¦  «        d                              |d         ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚d S )Nl   nE QWS�frG   r?   Ú
asymptoticrÆ   ÚexactrA   )	rM   rk   rl   r4   r   r^   ÚxpxÚatÚset)	r5   r6   rn   r‡   r7   r8   ÚautorË   rÌ   s	            r;   Ú	test_autozTestMannWhitneyU.test_autoå   s4  € õ Œi×#Ò#Ð$6Ñ7Ô7ˆØˆð �JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆØ�JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/ð �JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆØ�JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/õ ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/ð �JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆØ�JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/ð �JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆØ�JŠJ�s—z’z ! A¡#‘”Ñ'Ô'ˆÝŒF�1‰IŒI�aŒL×Ò˜Q˜qœTÑ"Ô"ˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/Ð/Ð/r=   )gm9—âªAj@g¼è+H3Œ[@gi¢²>s@)g›#ëhA{@gð±lËûz@gcžÚDf@gÇ³*âú¹h@gZ¬£ÿõA@gI9^¿YQa@g­ò§Ü`@g©î£ÞÕžp@gÑ¥Š–:q@gõó&º‚‚@g§Zµ|@g’`À×ÖÉr@gM‡cä3g@ú	two-sidedrË   ©rÅ   rÇ   )ç      0@ç
+ž×÷å?Úless)rÔ   ç
+ž×÷Õ?Úgreater)rÔ   ç¼ŒÉ%c�æ?rÌ   )rÔ   g9§ƒ:¨ƒæ?)rÔ   g9§ƒ:¨ƒÖ?)rÔ   g*…:¨ƒ:æ?)ÚkwdsÚexpectedr!   r"   c                 óü  — t          |¦  «        r%|j        dk     r|dk    rt          j        d¦  «         |€t	          |¦  «        nt          ||¦  «        }|                     | j        |¬¦  «        |                     | j        |¬¦  «        }}t          ||fi |¤Ž}t          |j        |                     |d         |¬¦  «        ¦  «         t          |j        |                     |d         |¬¦  «        ¦  «         d S )Nr&   r#   z*Scalar dtypes only respected after NEP 50.r*   r   r?   )r   r0   r1   r2   r   r3   r4   r7   r8   r   r   r\   r^   )r5   rÚ   rÛ   r!   r6   r7   r8   r`   s           r;   Ú
test_basiczTestMannWhitneyU.test_basic;  sñ   € õ �B‰<Œ<ð 	F˜BœN¨UÒ2Ð2°uÀ	Ò7IÐ7IÝŒKÐDÑEÔEÐEØ(-¨Õ  Ñ$Ô$Ð$½7À2ÀuÑ;MÔ;MˆØ�zŠz˜$œ&¨ˆzÑ.Ô.°·
²
¸4¼6È°
Ñ0OÔ0Oˆ1ˆÝ˜1˜aÐ(Ð( 4Ð(Ð(ˆÝ˜œ r§z¢z°(¸1´+ÀU zÑ'KÔ'KÑLÔLÐLÝ˜œ
 B§J¢J¨x¸¬{À% JÑ$HÔ$HÑIÔIÐIÐIÐIr=   T)rÅ   rÃ   )ç      7@rÕ   )rÞ   rÙ   )rÞ   r×   F)rÞ   gl,KÎNhä?)rÞ   gÊiÚ˜ØËå?)rÞ   gl,KÎNhÔ?c                 óF  — |                      | j        ¦  «        |                      | j        ¦  «        }}t          ||fddi|¤Ž}t	          |j        |                      |d         ¦  «        ¦  «         t	          |j        |                      |d         ¦  «        ¦  «         d S )NrÇ   rË   r   r?   )r4   r7   r8   r   r   r\   r^   )r5   rÚ   rÛ   r6   r7   r8   r`   s          r;   Útest_continuityz TestMannWhitneyU.test_continuityS  s�   € ð �zŠz˜$œ&Ñ!Ô! 2§:¢:¨d¬fÑ#5Ô#5ˆ1ˆÝ˜1˜aÐ=Ð=¨Ð=¸Ð=Ð=ˆÝ˜œ r§z¢z°(¸1´+Ñ'>Ô'>Ñ?Ô?Ð?Ý˜œ
 B§J¢J¨x¸¬{Ñ$;Ô$;Ñ<Ô<Ð<Ð<Ð<r=   c           
      óü  — |                      g d¢¦  «        }|                      g d¢¦  «        }|                      g d¢¦  «        dz  }|                      g d¢¦  «        dz  }|                     |dz
  ||z
  ||z
  |||z   ||z   |dz   g¦  «        }t          ||dd¬¦  «        }g d	¢}g d
¢}	t          |j        |                      |¦  «        ¦  «         t          |j        |                      |	¦  «        ¦  «         d S )N©rŒ   ç       @ç      @ç      @)rŒ   rã   rä   rå   r+   )r‹   rŒ   r‹   rŒ   r‹   r}   )r‹   r‹   rŒ   r‹   r‹   rh   rË   )rj   rÇ   )rE   rg   ç      !@rG   r¬   rF   rD   )r?   g]�UÄÚì?g…æ[»é?giœ…\­æ?g°ZX<_Õã?gãžxº.á?gåß ê…
Ù?)r4   Ústackr   r   r\   r   r^   )
r5   r6   r7   Úy0ÚdyÚdy2r8   r`   Ú
U_expectedÚ
p_expecteds
             r;   Útest_tie_correctz!TestMannWhitneyU.test_tie_correctc  s  € ð �JŠJÐ'Ð'Ð'Ñ(Ô(ˆØ�ZŠZÐ,Ð,Ð,Ñ-Ô-ˆØ�ZŠZÐ,Ð,Ð,Ñ-Ô-¨dÑ2ˆØ�jŠjÐ-Ð-Ð-Ñ.Ô.¨tÑ3ˆØ�HŠH�b˜‘g˜r "™u b¨¡f¨b°"°S±&¸"¸R¹%ÀÀDÁÐIÑJÔJˆÝ˜1˜a b°Ð>Ñ>Ô>ˆØ/Ð/Ð/ˆ
ðIð Ið Iˆ
å˜œ r§z¢z°*Ñ'=Ô'=Ñ>Ô>Ð>Ý˜œ
 B§J¢J¨zÑ$:Ô$:Ñ;Ô;Ð;Ð;Ð;r=   )g      Ð?r   g      è?)r¥   çš™™™™™É?çš™™™™™Ù?r¦   )r~   r¥   rî   r)   r   gÍÌÌÌÌÌä?r²   )rî   rï   r¦   )gôýÔxé&±?g /Ý$Á?gJ+‡Ñ?rï   r¦   )çyé&1¬œ?çÉv¾Ÿ/­?gÉv¾Ÿ/½?rî   gj¼t“Ô?çÛù~j¼tÛ?çƒÀÊ¡Eâ?)	gyé&1¬Œ?gV-²�?rñ   r¥   gÙÎ÷SãÅ?g´Èv¾ŸÏ?gÁÊ¡E¶óÕ?g'1¬ZÜ?gmçû©ñÒá?rY   )gÇK7‰A`Å?gZd;ßOÕ?r   gòÒMbXå?)çªñÒMb¨?çR¸…ëQ¸?çR¸…ëQÈ?ççû©ñÒMÒ?rò   ró   )	g;ßO�—n’?ç;ßO�—n¢?ç“V-²?g      À?gJ+‡É?r÷   gôýÔxé&Ù?r   g�•C‹lã?)çü©ñÒMb€?çü©ñÒMb�?çü©ñÒMb ?gyé&1¬¬?rõ   ççû©ñÒMÂ?g‘í|?5^Ê?g˜nƒÀÊÑ?g\�Âõ(\×?ç!°rh‘íÜ?çð§ÆK7‰á?)çü©ñÒMbp?rú   rû   rð   çú~j¼t“¨?g333333³?gÑ"Ûù~j¼?ç×£p=
×Ã?gáz®GáÊ?çð§ÆK7‰Ñ?g®GázÖ?g‹lçû©ñÚ?r   gºI+‡â?©r?   r@   rA   rB   rC   )rý   r÷   g1¬ZdÛ?ró   )rø   rù   rý   ç1¬ZdË?g%�•C‹Ô?rò   ró   )
gú~j¼t“ˆ?gú~j¼t“˜?r  gsh‘í|?µ?gøSã¥›ÄÀ?rö   r  g+‡ÙÖ?rþ   rÿ   )g{®Gázt?r}   gÛù~j¼t“?gL7‰A`å ?rñ   gj¼t“¶?gP�—nƒÀ?gºI+‡Æ?gX9´ÈvÎ?g…ëQ¸…Ó?gü©ñÒMbØ?gsh‘í|?Ý?gÇK7‰A`á?)çü©ñÒMb`?r   g;ßO�—n‚?g¸…ëQ¸Ž?g9´Èv¾Ÿš?gË¡E¶óý¤?gTã¥›Ä °?gbX9´È¶?g°rh‘í|¿?g…ëQ¸Å?r  gôýÔxé&Ñ?gÉv¾Ÿ/Õ?gòÒMbXÙ?gÃõ(\�ÂÝ?g…ëQ¸á?)rH   r  r   rú   g9´Èv¾ŸŠ?g/Ý$�•?rü   rô   gL7‰A`å°?g
×£p=
·?g¸…ëQ¸¾?r  gžï§ÆK7É?g`åÐ"ÛùÎ?g7‰A`åÐÒ?r)   g“V-Ú?gj¼t“Þ?gË¡E¶óýà?)r?   r@   rA   rB   rC   rD   c           	      ó   — t          t          dt          dd¦  «        ¦  «         | j        | j        | j        | j        dœ}|                     ¦   «         D �]¨\  }}|                     ¦   «         D �]Œ\  }}t          j	        dt          |¦  «        ¦  «        }t          j                             ||¦  «         t          t          j                             |¬¦  «        |d¬¦  «         t          j	        d||z  dz   ¦  «        }t          t          j                             |¬¦  «        t          j                             |¬¦  «        z   t          j                             |¬¦  «        z
  d¦  «         t          j                             |¬¦  «        }t          ||d d d…         ¦  «         t          j                             ||¦  «         t          j                             |¬¦  «        }	t          ||	¦  «         �ŒŽ�Œªd S )	NÚsr   )rA   rB   rC   rD   )ÚkrH   r-   r?   rh   )Úsetattrr   r   Úpn3Úpn4Úpm5Úpm6ÚitemsrM   rN   Úlenr  Ú
set_shapesr   r·   ÚsfÚpmf)
r5   Úp_tablesr‡   ÚtableÚmr:   ÚuÚu2r  Úpmf2s
             r;   Útest_exact_distributionz(TestMannWhitneyU.test_exact_distribution“  sÄ  € å•
˜C¥ a¨¡¤Ñ,Ô,Ð,à”x D¤H°´¸d¼hÐGÐGˆØ ŸšÑ(Ô(ð 	+ñ 	+‰HˆAˆuØŸš™œð +ñ +‘��1å”I˜a¥ Q¡¤Ñ(Ô(�Ý”×'Ò'¨¨1Ñ-Ô-Ð-Ý¥
¤× 0Ò 0°1Ð 0Ñ 5Ô 5°q¸tÐDÑDÔDÐDõ ”Y˜q ! A¡# a¡%Ñ(Ô(�Ý¥
¤× 0Ò 0°2Ð 0Ñ 6Ô 6Ý",¤,§/¢/°B /Ñ"7Ô"7ñ!8å",¤,×"2Ò"2°RÐ"2Ñ"8Ô"8ñ!9à:;ñ=ô =ð =õ
 !”l×&Ò&¨Ð&Ñ,Ô,�Ý  S¨¨¨2¨¤YÑ/Ô/Ð/õ ”×'Ò'¨¨1Ñ-Ô-Ð-Ý!”|×'Ò'¨"Ð'Ñ-Ô-�Ý  TÑ*Ô*Ð*Ñ*ñ)+ð	+ð 	+r=   c                 óü  — t           j                             d¦  «        }|                     |                     d¦  «        ¦  «        }|                     |                     d¦  «        ¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|                     |j        |j        z
  ¦  «        dk    sJ ‚|                     |                     d¦  «        ¦  «        }|                     |                     d¦  «        ¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|                     |j        |j        z
  ¦  «        dk     sJ ‚d S )	Niÿ0  rC   rÌ   rÆ   rË   r}   é(   rH   )rM   rk   rl   r4   r   r\   Úabsr^   )r5   r6   rn   r7   r8   Úres1Úres2s          r;   Útest_asymptotic_behaviorz)TestMannWhitneyU.test_asymptotic_behavior¯  sS  € ÝŒi×#Ò# EÑ*Ô*ˆð �JŠJ�s—z’z !‘}”}Ñ%Ô%ˆØ�JŠJ�s—z’z !‘}”}Ñ%Ô%ˆÝ˜A˜q¨Ð1Ñ1Ô1ˆÝ˜A˜q¨Ð6Ñ6Ô6ˆØŒ~ ¤Ò/Ð/Ð/Ð/Ø�vŠv�d”k D¤KÑ/Ñ0Ô0°4Ò7Ð7Ð7Ð7ð �JŠJ�s—z’z "‘~”~Ñ&Ô&ˆØ�JŠJ�s—z’z "‘~”~Ñ&Ô&ˆÝ˜A˜q¨Ð1Ñ1Ô1ˆÝ˜A˜q¨Ð6Ñ6Ô6ˆØŒ~ ¤Ò/Ð/Ð/Ð/Ø�vŠv�d”k D¤KÑ/Ñ0Ô0°4Ò7Ð7Ð7Ð7Ð7Ð7r=   c                 óÂ  — |                      g d¢¦  «        |                      ddg¦  «        }}t          ||dd¬¦  «        }t          ||dd¬¦  «        }t          |j        |j        ¦  «         |j        dk    sJ ‚t          ||d	d¬¦  «        }t          |j        |                      d
¦  «        ¦  «         t          |j        |                      d¦  «        ¦  «         d S )N©rŒ   rã   rä   rŸ   ç      @rÖ   rÌ   rÓ   rØ   r   rÒ   rä   rŒ   )r4   r   r   r^   r\   )r5   r6   r7   r8   Úres_lÚres_gr`   s          r;   Útest_exact_U_equals_meanz)TestMannWhitneyU.test_exact_U_equals_meanÄ  sÔ   € ð �zŠz˜,˜,˜,Ñ'Ô'¨¯ª°S¸#°JÑ)?Ô)?ˆ1ˆÝ˜Q ¨v¸gÐFÑFÔFˆÝ˜Q ¨yÀÐIÑIÔIˆÝ˜œ e¤lÑ3Ô3Ð3ØŒ|˜cÒ!Ð!Ð!Ð!å˜1˜a¨[ÀÐIÑIÔIˆÝ˜œ r§z¢z°"¡~¤~Ñ6Ô6Ð6Ý˜œ
 B§J¢J¨r¡N¤NÑ3Ô3Ð3Ð3Ð3r=   ©r‹   rŒ   )r‹   r   )r‹   g‹·éƒ¡Eï?)r‹   r?   )rÚ   Úresultc                 ó:   — t          t          di |¤Ž|¦  «         d S )Nr?   r@   ©r?   r@   )r   r   )r5   rÚ   r(  s      r;   Útest_scalar_dataz!TestMannWhitneyU.test_scalar_dataß  s(   € õ 	�Ð2Ð2¨TÐ2Ð2°FÑ;Ô;Ð;Ð;Ð;r=   c                 óà   — t          t          ddd¬¦  «        d¦  «         t          t          ddd¬¦  «        d¦  «         t          t          dddd¬¦  «        dt          j        f¦  «         d S )	Nr?   rÌ   rÆ   )r   r?   rË   F)rÇ   rÃ   r   )r   r   rM   r]   ©r5   s    r;   Útest_equal_scalar_dataz'TestMannWhitneyU.test_equal_scalar_dataä  sˆ   € õ
 	•\ ! Q¨wÐ7Ñ7Ô7¸ÑBÔBÐBÝ•\ ! Q¨|Ð<Ñ<Ô<¸hÑGÔGÐGõ 	•\ ! Q¨|Ø16ð8ñ 8ô 8Ø:=½r¼v¸ñ	Hô 	Hð 	Hð 	Hð 	Hr=   rÇ   c                 óL  — t           j                             d¦  «        }d}d\  }}|                     |ddf¦  «        }|                     d|ddf¦  «        dz   }t          |                     |¦  «        |                     |¦  «        ||¬	¦  «        }	d
}
|	j        j        |
k    sJ ‚|	j        j        |
k    sJ ‚t          j        ||d¦  «        t          j        ||d¦  «        }}|d         }|j	        |j	        k    sJ ‚t          j
        ||
|fz   ¦  «        }t          j
        ||
|fz   ¦  «        }|j        d d…         |
k    sJ ‚|j        d d…         |
k    sJ ‚t          j        |
¦  «        }t          j        |
¦  «        }t          d„ |
D ¦   «         Ž D ]8}||         }||         }t          |||¬¦  «        }|j        ||<   |j        ||<   Œ9t          |	j        |                     |¦  «        d¬¦  «         t          |	j        |                     |¦  «        d¬¦  «         d S )Nl   2=U éýÿÿÿ)rF   rE   rA   rG   rD   r?   r¥   )rÇ   rj   )rD   rA   rG   rh   )N.c                 ó,   — g | ]}t          |¦  «        ‘ŒS rz   )Úrange)Ú.0ro   s     r;   ú
<listcomp>z8TestMannWhitneyU.test_gh_12837_11113.<locals>.<listcomp>  s   € Ð 9Ð 9Ð 9¨a¥ q¡¤Ð 9Ð 9Ð 9r=   rÆ   ç¼‰Ø—²Òœ<r-   )rM   rk   rl   r   r4   r^   Úshaper\   ÚmoveaxisÚndimÚbroadcast_toÚzerosr   r   )r5   rÇ   r6   rn   rj   r  r‡   r7   r8   r`   r6  Ú
statisticsÚpvaluesÚindicesÚxiÚyiÚtemps                    r;   Útest_gh_12837_11113z$TestMannWhitneyU.test_gh_12837_11113ó  s.  € õ Œi×#Ò# JÑ/Ô/ˆð ˆØ‰ˆˆ1Ø�JŠJ˜˜1˜a�yÑ!Ô!ˆØ�JŠJ˜˜1˜a �|Ñ$Ô$ sÑ*ˆÝ˜2Ÿ:š: a™=œ=¨"¯*ª*°Q©-¬-ÀÈTÐRÑRÔRˆàˆØŒzÔ 5Ò(Ð(Ð(Ð(ØŒ}Ô" eÒ+Ð+Ð+Ð+õ Œ{˜1˜d BÑ'Ô'­¬°Q¸¸bÑ)AÔ)Aˆ1ˆàˆiŒLˆØŒv˜œÒÐÐÐåŒO˜A˜u¨ t™|Ñ,Ô,ˆÝŒO˜A˜u¨ t™|Ñ,Ô,ˆØŒw�s˜�sŒ|˜uÒ$Ð$Ð$Ð$ØŒw�s˜�sŒ|˜uÒ$Ð$Ð$Ð$õ ”X˜e‘_”_ˆ
Ý”(˜5‘/”/ˆÝÐ 9Ð 9°5Ð 9Ñ 9Ô 9Ð:ð 	+ð 	+ˆGØ�7”ˆBØ�7”ˆBÝ  B¨vÐ6Ñ6Ô6ˆDØ"&¤.ˆJ�wÑØ#œ{ˆG�GÑÐå˜œ
 B§J¢J¨wÑ$7Ô$7¸eÐDÑDÔDÐDÝ˜œ r§z¢z°*Ñ'=Ô'=ÀEÐJÑJÔJÐJÐJÐJr=   c                 ó~  — g d¢}g d¢}t          |                     |¦  «        |                     |¦  «        ¦  «        }d|d<   t          j        |d<   t          |                     |¦  «        |                     |¦  «        ¦  «        }t	          |j        |j        ¦  «         t	          |j        |j        ¦  «         d S )NrY   )rA   rD   rF   rG   rg   rA   r@   r?   rB   rB   rC   rŒ   r   rB   )r   r4   rM   rm   r   r\   r^   )r5   r6   r7   r8   r  r  s         r;   Útest_gh_11355zTestMannWhitneyU.test_gh_11355   s    € àˆLˆLˆØ-Ð-Ð-ˆÝ˜BŸJšJ q™MœM¨2¯:ª:°a©=¬=Ñ9Ô9ˆð ˆˆ!‰ÝŒvˆˆ!‰Ý˜BŸJšJ q™MœM¨2¯:ª:°a©=¬=Ñ9Ô9ˆå˜œ¨¬Ñ7Ô7Ð7Ý˜œ T¤[Ñ1Ô1Ð1Ð1Ð1r=   c                 ó\  — g d¢}dddt           j        ddddddd	g}t          |                     |¦  «        |                     |¦  «        ¦  «        }t	          |j        |                     |j        ¦  «        ¦  «         t	          |j        |                     |j        ¦  «        ¦  «         d S )
Nrâ   rA   rD   rF   rg   r@   r?   rB   rC   )rM   r]   r   r4   r   r\   r^   )r5   r6   r7   r8   Úres3s        r;   Útest_gh11355_nanz!TestMannWhitneyU.test_gh11355_nan.  s–   € ð ÐÐˆØ��1•b”f˜a  A q¨!¨Q°Ð2ˆÝ˜BŸJšJ q™MœM¨2¯:ª:°a©=¬=Ñ9Ô9ˆÝ˜œ¨¯
ª
°2´6Ñ(:Ô(:Ñ;Ô;Ð;Ý˜œ R§Z¢Z°´Ñ%7Ô%7Ñ8Ô8Ð8Ð8Ð8r=   )rŒ   r@   rA   rB   rA   rD   rF   rG   r@   r?   rB   rC   ç      $@g+z’¿QœÀ?ræ   g}$kù\¶?g     €1@g!Ë›G*ã?rÔ   g›Á,ÖsÿÝ?ç     €8@gêFÏHïQé?)r7   r8   r\   r^   c                 ó  — t          |                     |¦  «        |                     |¦  «        d¬¦  «        }t          |j        |                     |¦  «        d¬¦  «         t          |j        |                     |¦  «        d¬¦  «         d S )NrË   rÆ   çê-�™—q=r-   )r   r4   r   r\   r^   )r5   r7   r8   r\   r^   r6   r`   s          r;   Útest_gh_11355bzTestMannWhitneyU.test_gh_11355bG  sx   € õ ˜2Ÿ:š: a™=œ=¨"¯*ª*°Q©-¬-ÀÐMÑMÔMˆÝ˜œ r§z¢z°)Ñ'<Ô'<À5ÐIÑIÔIÐIÝ˜œ
 B§J¢J¨vÑ$6Ô$6¸UÐCÑCÔCÐCÐCÐCr=   )TrÖ   rË   g¾ˆ²Ò&Óì?)TrØ   rË   g„Ïì‡ÓO¿?)TrÒ   rË   g„Ïì‡ÓOÏ?)FrÖ   rË   g9@VN!xì?)FrØ   rË   g9þM�õ>¼?)FrÒ   rË   g9þM�õ>Ì?)TrÖ   rÌ   g ?UÈV²ì?)TrØ   rÌ   gßºVJHÀ?)TrÒ   rÌ   gèÁVJHÐ?)rÃ   rÅ   rÇ   Ú
pvalue_expc                 ó.  — d}|                      g d¢¦  «        }|                      g d¢¦  «        }t          |||||¬¦  «        }	t          |	j        |                      |¦  «        ¦  «         t	          |	j        |                      |¦  «        ¦  «         d S )Ng     €A@)
r—   g�Âõ(\�ê?g=
×£p=þ?g¤p=
×£ð?g333333÷?g®Gázö?g�Âõ(\�þ?g=
×£p=ú?r(   g\�Âõ(\÷?)gffffffò?g)\�Âõ(ì?r«   g®Gáz®ç?g\�Âõ(\ó?©rÃ   rÅ   rÇ   )r4   r   r   r\   r   r^   )
r5   rÃ   rÅ   rÇ   rL  r6   Ústatistic_expr7   r8   r`   s
             r;   Útest_gh_9184zTestMannWhitneyU.test_gh_9184X  s�   € ð. ˆØ�JŠJÐSÐSÐSÑTÔTˆØ�JŠJÐ5Ð5Ð5Ñ6Ô6ˆÝ˜1˜a°Ø'2¸6ðCñ Cô Cˆå˜œ r§z¢z°-Ñ'@Ô'@ÑAÔAÐAÝ˜œ
 B§J¢J¨zÑ$:Ô$:Ñ;Ô;Ð;Ð;Ð;r=   c                 ó<  — |                      |j        ¦  «        }|                     |||||g¦  «        }|                     |||||g¦  «        }t          ||¦  «        }t	          |j        |                      |¦  «        ¦  «         t	          |j        |¦  «         d S r®   )r4   r]   rç   r   r   r\   r^   )r5   r6   r]   ÚaÚbr`   s         r;   Útest_gh_4067zTestMannWhitneyU.test_gh_4067w  s’   € ð �jŠj˜œÑ Ô ˆØ�HŠH�c˜3  S¨#Ð.Ñ/Ô/ˆØ�HŠH�c˜3  S¨#Ð.Ñ/Ô/ˆÝ˜1˜aÑ Ô ˆÝ˜œ r§z¢z°#¡¤Ñ7Ô7Ð7Ý˜œ
 CÑ(Ô(Ð(Ð(Ð(r=   r"  rŸ   r#  )rä   ga×€}¢ã?)rä   rŒ   )rŸ   g¤Û?hÍå?)rŸ   rŒ   )rã   gæ5&#\å?)rã   rŒ   )r7   r8   rÅ   rÛ   c                 ór  — t          |                     |¦  «        |                     |¦  «        d|d¬¦  «        }t          |¦  «        |j        k    rdnd}t	          |j        |                     |d         ¦  «        |¬¦  «         t	          |j        |                     |d         ¦  «        |¬¦  «         d S )	NTrË   rN  ç�íµ ÷Æ°>rJ  r   ©Úrtolr?   )r   r4   r   r#   r   r\   r^   )r5   r7   r8   rÅ   rÛ   r6   r`   rX  s           r;   Útest_gh_2118zTestMannWhitneyU.test_gh_2118‘  s«   € õ ˜2Ÿ:š: a™=œ=¨"¯*ª*°Q©-¬-ÈØ'2¸<ðIñ Iô Iˆå'¨Ñ+Ô+¨r¬zÒ9Ð9ˆtˆt¸uˆÝ˜œ r§z¢z°(¸1´+Ñ'>Ô'>ÀTÐJÑJÔJÐJÝ˜œ
 B§J¢J¨x¸¬{Ñ$;Ô$;À$ÐGÑGÔGÐGÐGÐGr=   c                 óp  — t           j                             d¦  «        }d\  }}|                     |¬¦  «        }|                     |¬¦  «        }t          t          dt          dd¦  «        ¦  «         t          j                             ¦   «          t          j	        ||d¬¦  «        }t          j        j
        j        }|d         t          |j        ||z  |j        z
  ¦  «        d	z   k    sJ ‚t          j	        ||d¬¦  «         |t          j        j
        j        k    sJ ‚t          j                             ¦   «          t          j	        |d|z  dd
¬¦  «         t          j        j
        j        }|d         d	k    sJ ‚t          j	        d|z  |dd
¬¦  «         |t          j        j
        j        k    sJ ‚d S )Nì   g>·mŽjÑK )rC   rc   rd   r  r   rÌ   rÆ   rh   r?   rØ   )rÇ   rÅ   )rM   rk   rl   r
  r   r   r  Úresetr   r   Úconfigurationsr6  Úminr\   )r5   rn   r  r‡   r7   r8   r`   r6  s           r;   Útest_gh19692_smaller_tablez+TestMannWhitneyU.test_gh19692_smaller_table›  sŒ  € õ
 Œi×#Ò#Ð$7Ñ8Ô8ˆØ‰ˆˆ1Ø�JŠJ˜AˆJÑÔˆØ�JŠJ˜AˆJÑÔˆå•
˜C¥ a¨¡¤Ñ,Ô,Ð,ÝŒ×ÒÑÔÐåÔ   A¨gÐ6Ñ6Ô6ˆÝ”Ô+Ô1ˆØ�RŒy�C ¤¨q°©s°S´]Ñ/BÑCÔCÀaÑGÒGÐGÐGÐGÝÔ˜1˜a¨Ð0Ñ0Ô0Ð0Ø�
œÔ3Ô9Ò9Ð9Ð9Ð9õ
 	Œ×ÒÑÔÐÝÔ˜1˜a ™c¨'¸yÐIÑIÔIÐIÝ”Ô+Ô1ˆØ�RŒy˜AŠ~ˆ~ˆ~ˆ~ÝÔ˜1˜Q™3 ¨'¸yÐIÑIÔIÐIØ�
œÔ3Ô9Ò9Ð9Ð9Ð9Ð9Ð9r=   rÅ   )rÖ   rØ   rÒ   c                 ó’  — t           j                             d¦  «        }|                     d¬¦  «        }|                     d¬¦  «        }t          j        ||t          j        ¦   «         |d¬¦  «        }t          j        ||d|d¬¦  «        }t          |j        |j        d¬	¦  «         t          |j        |j        d¬	¦  «         d S )
Nr[  )r@   rC   rd   )r@   rD   r?   )rÇ   rÅ   rj   rÌ   çVçž¯Ò<rW  )	rM   rk   rl   r   r   ÚPermutationMethodr   r\   r^   )r5   rÅ   rn   r7   r8   r`   r  s          r;   Útest_permutation_methodz(TestMannWhitneyU.test_permutation_method¸  sÉ   € åŒi×#Ò#Ð$7Ñ8Ô8ˆØ�JŠJ˜FˆJÑ#Ô#ˆØ�JŠJ˜FˆJÑ#Ô#ˆÝÔ   A­eÔ.EÑ.GÔ.GØ-8¸qðBñ Bô BˆåÔ! ! Q¨wØ.9ÀðCñ Cô Cˆå˜œ t¤~¸EÐBÑBÔBÐBÝ˜œ
 D¤K°eÐ<Ñ<Ô<Ð<Ð<Ð<r=   )g®Þ	ÂUå3@g^û�Â–‹3@g™òtÇ3@g½¢™]o·5@g·z³S4@gÔtÆ< 4@gM	‚«X3@gXÄmCók4@gLìo&£3@gfÉ˜ÊâÇ2@gžý{m;^3@ga A$|3@gMÑ¢Á_À4@gŽÌ#é3@g19T•_3@gg§„Cÿ2@gúqïR4@g˜Kk/±4@gè—Ü�Š3@gÌZ†û2@g„_
À²ƒ3@g«ÝÀ$`Ó3@g
ŒÃx„4@g3ƒ#Ð»5@g;VND÷ò1@g >®ƒåH2@g—µâr4@g³»ØyÜ¤2@gèà’?4@g=bòÉ‡Ž3@)gÔäìpG3@g„Lä:jþ0@gýî„¨Š2@g«¹JNnC1@gRÀ;`(3@g+Pi,ƒ2@g6:Z›2@gZ"3=ùÕ2@gF	uSa	3@gXâÒ€2@g¥ÙX3^83@g§2¡Õ1@g÷��J2@gZ)Eû‡ª2@g¦U8ø3@g°»°50@g�ˆ³Ô†3@gFz‰™33@gøèãz3@g¨¢²Ú´3@g›+¡†›„=c                 óx  — |                      | j        ¦  «        |                      | j        ¦  «        }}t          j        ||d¬¦  «        \  }}t          j        ||d¬¦  «        \  }}t          j        ||d¬¦  «        \  }}	t          j        ||d¬¦  «        \  }
}t          ||¦  «         t          |	|¦  «         ||	k    sJ ‚t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         t          |
|                      d¦  «        ¦  «         t          ||                      d¦  «        | j        ¬¦  «         |j        |j	        k    r| j        nd}t          |	|                      d	¦  «        |d
¬¦  «         d S )NrÖ   rÄ   rØ   ç      @ç     €Y@g=_A§ÿï?rW  çü©ñÒMb@?glP†z.?r5  ©rX  r.   )
r4   ÚXÚYr   r   r   r   rX  r!   r$   ©r5   r6   ri  rj  Úu1rR   r  rT   Úu3rV   Úu4Úp4rX  s                r;   Útest_mannwhitneyu_one_sidedz,TestMannWhitneyU.test_mannwhitneyu_one_sidedÖ  s—  € Ø�zŠz˜$œ&Ñ!Ô! 2§:¢:¨d¬fÑ#5Ô#5ˆ1ˆÝÔ# A q°fÐ=Ñ=Ô=‰ˆˆBÝÔ# A q°iÐ@Ñ@Ô@‰ˆˆBÝÔ# A q°iÐ@Ñ@Ô@‰ˆˆBÝÔ# A q°fÐ=Ñ=Ô=‰ˆˆBå˜˜BÑÔÐÝ˜˜BÑÔÐØ�RŠxˆxˆxˆxÝ˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJÐ'8Ñ9Ô9ÀÄ	ÐJÑJÔJÐJØœG r¤zÒ1Ð1ˆtŒyˆy°tˆÝ˜˜BŸJšJÐ':Ñ;Ô;À$ÈUÐSÑSÔSÐSÐSÐSr=   c                 ó   — |                      | j        ¦  «        |                      | j        ¦  «        }}t          j        ||d¬¦  «        \  }}t          j        ||d¬¦  «        \  }}t          ||¦  «         t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         |j        |j        k    r| j        nd}t          ||                      d¦  «        |d¬¦  «         d S )	NrÒ   rÄ   re  rf  rg  glP†z.?r5  rh  ©
r4   ri  rj  r   r   r   r!   r$   rX  r   ©	r5   r6   ri  rj  rl  rR   r  rT   rX  s	            r;   Útest_mannwhitneyu_two_sidedz,TestMannWhitneyU.test_mannwhitneyu_two_sidedè  sè   € Ø�zŠz˜$œ&Ñ!Ô! 2§:¢:¨d¬fÑ#5Ô#5ˆ1ˆÝÔ# A q°kÐBÑBÔB‰ˆˆBÝÔ# A q°kÐBÑBÔB‰ˆˆBå˜˜BÑÔÐÝ˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-ØœG r¤zÒ1Ð1ˆtŒyˆy°tˆÝ˜˜BŸJšJÐ'9Ñ:Ô:ÀÈEÐRÑRÔRÐRÐRÐRr=   c                 óx  — |                      | j        ¦  «        |                      | j        ¦  «        }}t          j        ||dd¬¦  «        \  }}t          j        ||dd¬¦  «        \  }}t          j        ||dd¬¦  «        \  }}	t          j        ||dd¬¦  «        \  }
}t          ||¦  «         t          |	|¦  «         ||	k    sJ ‚t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         t          |
|                      d¦  «        ¦  «         |j        |j        k    r| j        nd}t          ||                      d¦  «        |d	¬
¦  «         t          |	|                      d¦  «        |d	¬
¦  «         d S )NFrÖ   rÄ   rØ   re  rf  rg  gii5‡£ÿï?r5  rh  g¬2¥¥2?rr  rk  s                r;   Ú&test_mannwhitneyu_no_correct_one_sidedz7TestMannWhitneyU.test_mannwhitneyu_no_correct_one_sidedó  sŸ  € Ø�zŠz˜$œ&Ñ!Ô! 2§:¢:¨d¬fÑ#5Ô#5ˆ1ˆÝÔ# A q¨%¸VÐDÑDÔD‰ˆˆBÝÔ# A q¨%¸YÐGÑGÔG‰ˆˆBÝÔ# A q¨%¸YÐGÑGÔG‰ˆˆBÝÔ# A q¨%¸VÐDÑDÔD‰ˆˆBå˜˜BÑÔÐÝ˜˜BÑÔÐØ�RŠxˆxˆxˆxÝ˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-ØœG r¤zÒ1Ð1ˆtŒyˆy°tˆÝ˜˜BŸJšJÐ'8Ñ9Ô9ÀÈ5ÐQÑQÔQÐQÝ˜˜BŸJšJÐ';Ñ<Ô<À4ÈeÐTÑTÔTÐTÐTÐTr=   c                 ó  — |                      | j        ¦  «        |                      | j        ¦  «        }}t          j        ||dd¬¦  «        \  }}t          j        ||dd¬¦  «        \  }}t          ||¦  «         t          ||                      d¦  «        ¦  «         t          ||                      d¦  «        ¦  «         |j        |j        k    r| j        nd}t          ||                      d¦  «        |d¬	¦  «         d S )
NFrÒ   rÄ   re  rf  rg  g¬2¥¥2?r5  rh  rr  rs  s	            r;   Ú&test_mannwhitneyu_no_correct_two_sidedz7TestMannWhitneyU.test_mannwhitneyu_no_correct_two_sided  sì   € Ø�zŠz˜$œ&Ñ!Ô! 2§:¢:¨d¬fÑ#5Ô#5ˆ1ˆÝÔ# A q¨%¸[ÐIÑIÔI‰ˆˆBÝÔ# A q¨%¸[ÐIÑIÔI‰ˆˆBå˜˜BÑÔÐÝ˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-Ý˜˜BŸJšJ tÑ,Ô,Ñ-Ô-Ð-ØœG r¤zÒ1Ð1ˆtŒyˆy°tˆÝ˜˜BŸJšJÐ';Ñ<Ô<À4ÈeÐTÑTÔTÐTÐTÐTr=   c                 óÌ  — |                      g d¢¦  «        }|                      g d¢¦  «        }t          j        ||d¬¦  «        }t          |j        |                      d¦  «        ¦  «         t          |j        |                      d¦  «        ¦  «         t          j        ||d¬¦  «        }t          |j        |                      d¦  «        ¦  «         t          |j        |                      d¦  «        ¦  «         t          j        ||d	¬¦  «        }t          |j        |                      d¦  «        ¦  «         t          |j        |                      d
¦  «        ¦  «         d S )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Œ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rã   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rä   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   )žrŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rä   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rã   rŒ   rŒ   rã   rŒ   rŒ   rã   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rã   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rã   rã   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rŒ   rã   rŒ   rŒ   rã   rŒ   rŒ   rŒ   rŒ   rã   rŒ   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     •Ð@g €?Õ¹•ý>rØ   gÜ&)EÅÿï?rÒ   g�Ž?Õ¹•?)r4   r   r   r   r\   r^   r_   s        r;   Útest_mannwhitneyu_onesz'TestMannWhitneyU.test_mannwhitneyu_ones  sK  € à�JŠJð 0ð 0ð 0ñ 1ô 1ˆð& �JŠJð (ð (ð (ñ )ô )ˆõ Ô   A°6Ð:Ñ:Ô:ˆÝ˜œ r§z¢z°'Ñ':Ô':Ñ;Ô;Ð;Ý˜œ
 B§J¢JÐ/DÑ$EÔ$EÑFÔFÐFÝÔ   A°9Ð=Ñ=Ô=ˆÝ˜œ r§z¢z°'Ñ':Ô':Ñ;Ô;Ð;Ý˜œ
 B§J¢J¨Ñ$?Ô$?Ñ@Ô@Ð@ÝÔ   A°;Ð?Ñ?Ô?ˆÝ˜œ r§z¢z°'Ñ':Ô':Ñ;Ô;Ð;Ý˜œ
 B§J¢JÐ/@Ñ$AÔ$AÑBÔBÐBÐBÐBr=   N)1ru   rv   rw   r1   rx   r»   rÀ   rÉ   rÑ   r7   r8   Úcases_basicry   rÝ   Úcases_continuityrà   rí   r  r  r  r  r  r   r&  Úcases_scalarr+  r.  rA  rC  rF  rM   rm   Úcases_11355rK  Ú
cases_9184rP  rT  Ú
cases_2118rY  r_  rc  ri  rj  rX  rp  rt  rv  rx  rz  rz   r=   r;   r½   r½   »   s_  € € € € € ð „[×!Ò! +Ð6SÐ!ÑTÔTð-ð -ñ UÔTð-ð*
.ð 
.ð 
.ð10ð 10ð 10ðj 	-Ð,Ð,€Að	ð 	ð 	€Að& %0¸<ÐHÐHØ*ð,à$*°lÐCÐCØ*ð,à$-¸ÐFÐFØ*ð,à$/¸7ÐCÐCØ*ð,à$*°gÐ>Ð>Ø*ð,à$-¸ÐAÐAØ*ð,ð-€Kð „[×ÒÐ1°;Ñ?Ô?Ø„[×Ò˜WÐ&BÐ&BÐ&BÑCÔCðJð Jñ DÔCñ @Ô?ðJð *5ÈÐMÐMØ/ð1à)/À4ÐHÐHØ/ð1à)2ÀdÐKÐKØ/ð1à)4ÈÐNÐNØ/ð1à)/À5ÐIÐIØ/ð1à)2ÀeÐLÐLØ/ð1ð2Ðð „[×ÒÐ1Ð3CÑDÔDð=ð =ñ EÔDð=ð<ð <ð <ð,  ÐÐÐ$8Ð$8Ð$8Ø.Ð.Ð.ð0ð 0€CàˆoˆoÐ"AÐ"AÐ"AØ<Ð<Ð<ØKÐKÐKðMð M€Cð )Ð(Ð(Ø7Ð7Ð7ØKÐKÐKð2ð 2ð 2ð<ð <ð <ð=ð =€Cð +Ð*Ð*Ø?Ð?Ð?ð1ð 1ð 1ð9ð 9ð 9ðGð Gð Gð$ð $ð $ð%ð %€Cð+ð +ð +ð88ð 8ð 8ð*4ð 4ð 4ð" &1¸LÐIÐIØðà%+°|ÐDÐDØð à%.¸,ÐGÐGØ)ð+à%0¸GÐDÐDÀgÐNØ%+°wÐ?Ð?ÀÐKØ%.¸'ÐBÐBÀGÐLðN€Lð „[×ÒÐ/°Ñ>Ô>ð<ð <ñ ?Ô>ð<ðHð Hð Hð „[×!Ò! +Ð6SÐ!ÑTÔTØ„[×Ò˜X¨°gÐ'>Ñ?Ô?ð)Kð )Kñ @Ô?ñ UÔTð)KðV2ð 2ð 2ð „[×!Ò! +Ð6SÐ!ÑTÔTð9ð 9ñ UÔTð9ð "�M�MØ˜˜1˜a ¤¨¨A¨q°!°Q¸Ð:Ø˜ð*ð "�M�MØ˜˜1˜a ¤¨¬°°A°q¸!¸QÐ?ØÐ)ð+ð ˜˜2œ6 1Ð%Ø˜˜1˜a ¤¨¨A¨q°!°Q¸Ð:Ø˜/ð+ð ˜˜2œ6 1Ð%Ø˜˜1˜a ¤¨¬°°A°q¸!¸QÐ?Ø˜ð*ð ˜œ ¤¨Ð*Ø˜˜1˜a ¤¨¬°°A°q¸!¸QÐ?Ø˜/ð+ð,€Kð  „[×ÒÐ>ÀÑLÔLðDð Dñ MÔLðDð ?Ð>Ð>ØBÐBÐBØCÐCÐCØ@Ð@Ð@ØCÐCÐCØEÐEÐEØ:Ð:Ð:Ø=Ð=Ð=Ø?Ð?Ð?ðA€Jð „[×Òð 6Ø7AñCô Cð<ð <ñCô Cð<ð: „[×!Ò! +Ð6SÐ!ÑTÔTð)ð )ñ UÔTð)ð  �<�< # s ¨YÐ8MÐNØ�<�< # s ¨VÐ5JÐKØ�<�< # s ¨[¸)ÐDØ�9�9˜q˜c 9Ð.CÐDØ�9�9˜q˜c 6Ð+@ÐAØ�9�9˜q˜c ;°	Ð:Ø�q�6˜A˜q˜6 9Ð.BÐCØ�q�6˜A˜q˜6 6Ð+?Ð@Ø�q�6˜A˜q˜6 ;°Ð9ð;€Jð „[×ÒÐBÐBÐBÀJÑOÔOðHð Hñ PÔOðHð:ð :ð :ð: „[×Ò˜]Ð,LÐ,LÐ,LÑMÔMð	=ð 	=ñ NÔMð	=ð	-ð 	-ð 	-€Að	Nð 	Nð 	N€Að €DðTð Tð Tð$	Sð 	Sð 	SðUð Uð Uð$	Uð 	Uð 	Uð+Cð +Cð +Cð +Cð +Cr=   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ej                             dd¦  «        d„ ¦   «         Zd„ ZdS )ÚTestSomersDc                 óV  ‡ — ‰ j         ‰ j        z   ‰ _        t          j        d¦  «        ‰ j         ‰ j        z   ft          j        d¦  «        ‰ j         ‰ j        z   fdœ‰ _        ˆ fd„‰ j        D ¦   «         }t          j        t          j	        d¬¦  «        ‰ _
         ‰ j
        |Ž ‰ _        d S )NrE   ©r   r?   c                 ó6   •— g | ]}‰j         |         d          ‘ŒS )r   )Ú	arguments)r3  Úidxr5   s     €r;   r4  z,TestSomersD.setup_method.<locals>.<listcomp>E  s%   ø€ ÐHÐHÐH°#�t”~ cÔ*¨1Ô-ÐHÐHÐHr=   rÒ   rÄ   )ÚALL_INTEGERÚ	ALL_FLOATÚdtypesrM   rN   r†  Ú	functoolsÚpartialr   ÚsomersdÚpartialfuncrÛ   )r5   Úinput_arrays   ` r;   Úsetup_methodzTestSomersD.setup_method?  s¹   ø€ ØÔ&¨¬Ñ7ˆŒÝ œi¨™mœmØ"Ô.°´Ñ?ðAå œi¨™mœmØ"Ô.°´Ñ?ðAðBð BˆŒð IÐHÐHÐH¸¼ÐHÑHÔHˆõ
 %Ô,­U¬]Ø9DðFñ Fô FˆÔà(˜Ô(¨+Ð6ˆŒˆˆr=   c                 óž   —  | j         |Ž }t          |j        | j        j        d¬¦  «         t          |j        | j        j        d¬¦  «         d S )Nra  r-   )rŽ  r   r\   rÛ   r^   )r5   r¹   r`   s      r;   ÚpythranfunczTestSomersD.pythranfuncN  sO   € ØˆdÔ Ð%ˆÝ˜œ t¤}Ô'>ÀUÐKÑKÔKÐKÝ˜œ
 D¤MÔ$8¸uÐEÑEÔEÐEÐEÐEr=   c                 ó  — g d¢g d¢g d¢g}t          j        |¦  «        }|                      t           j        ¦  «        }t          j        |fi |¤Ž}t          |j        |j        d¬¦  «         t          |j        |j        d¬¦  «         d S )N)r�   é   é   rF   r   )rF   r•  é   é#   rf   )r?   rA   r@   rF   é   ra  r-   )r   r�  Úget_optional_argsr   r\   r^   )r5   r  r  Úoptional_argsr  s        r;   Útest_pythranfunc_keywordsz%TestSomersD.test_pythranfunc_keywordsS  s•   € à#Ð#Ð#Ð%8Ð%8Ð%8Ð:JÐ:JÐ:JÐKˆÝŒ}˜UÑ#Ô#ˆà×.Ò.­u¬}Ñ=Ô=ˆÝŒ}˜UÐ4Ð4 mÐ4Ð4ˆå˜œ¨¬¸UÐCÑCÔCÐCÝ˜œ T¤[°uÐ=Ñ=Ô=Ð=Ð=Ð=r=   c                 óœ
  — g d¢}g d¢}d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         g d¢}g d	¢}d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         g d
¢}g d¢}d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         t          j        d¦  «        }t          j        d¦  «        }d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         t          j        d¦  «        }t          j        g d¢¦  «        }d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         t          j        d¦  «        }t          j        d¦  «        d d d…         }d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         t          j        d¦  «        }t          j        g d¢¦  «        }d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         g d¢}g d¢}d}t          j        ||¦  «        }t          |j        |d         d¬¦  «         t          |j        |d         d¬¦  «         t          j        g d¢g d¢¦  «        }t          |j        t
          j        ¦  «         t          |j        t
          j        ¦  «         t          j        g d¢g d¢¦  «        }t          |j        t
          j        ¦  «         t          |j        t
          j        ¦  «         t          j        g d¢g d¢¦  «        }t          |j        t
          j        ¦  «         t          |j        t
          j        ¦  «         t          j        dgdg¦  «        }t          |j        t
          j        ¦  «         t          |j        t
          j        ¦  «         t          j        g g ¦  «        }t          |j        t
          j        ¦  «         t          |j        t
          j        ¦  «         t          j        d¦  «        }t          j        d¦  «        }t          t          t           j        ||¦  «         d S )N)rC   r@   r?   rA   rD   rB   rF   rG   )rC   r@   rD   rA   r?   rG   rF   rB   r'  r   ra  r-   r?   )	r   rC   r@   r?   rA   rD   rB   rF   rG   )	rC   r@   r   rD   rA   r?   rG   rF   rB   )rC   r@   r?   rA   rD   rB   rF   )rC   r@   rD   rA   r?   rF   rB   )g+$I’$IÂ¿g=/3n£+ä?rE   ©rŒ   r   )
r   r@   r?   rA   rB   rD   rC   rF   rG   rg   )gsÒ'}Ò'í?r‹   rh   )g      ð¿r   )
rg   rF   rG   rD   rC   rA   rB   r@   r?   r   )g}Ò'}Ò'í¿r‹   )rf   r@   r?   rf   r@   )r?   rB   rF   r?   r   )ç      à¿gÞ.Ê‚ƒÓ?)r@   r@   r@   )r@   r   r@   rG  ç      4@)r   r�  r   r\   r^   rM   rN   Úarrayr]   rÈ   r¶   )r5   r7   r8   rÛ   r`   Úx1Úx2s          r;   Útest_like_kendalltauz TestSomersD.test_like_kendalltau^  sÆ  € ð %Ð$Ð$ˆØ$Ð$Ð$ˆà9ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<ð (Ð'Ð'ˆØ'Ð'Ð'ˆà9ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<ð "Ð!Ð!ˆØ!Ð!Ð!ˆà:ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<õ ŒI�b‰MŒMˆÝŒI�b‰MŒMˆð *ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<õ ŒI�b‰MŒMˆÝŒHÐ3Ð3Ð3Ñ4Ô4ˆà9ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<õ ŒI�b‰MŒMˆÝŒI�b‰MŒM˜$˜$˜B˜$Ôˆð +ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<õ ŒI�b‰MŒMˆÝŒHÐ3Ð3Ð3Ñ4Ô4ˆà;ˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<ð ÐÐˆØˆ_ˆ_ˆà:ˆÝŒm˜B Ñ#Ô#ˆÝ˜œ x°¤{¸Ð?Ñ?Ô?Ð?Ý˜œ
 H¨Q¤K°eÐ<Ñ<Ô<Ð<õ Œm˜I˜I˜I y y yÑ1Ô1ˆÝ˜œ¥r¤vÑ.Ô.Ð.Ý˜œ
¥B¤FÑ+Ô+Ð+åŒm˜I˜I˜I y y yÑ1Ô1ˆÝ˜œ¥r¤vÑ.Ô.Ð.Ý˜œ
¥B¤FÑ+Ô+Ð+åŒm˜I˜I˜I y y yÑ1Ô1ˆÝ˜œ¥r¤vÑ.Ô.Ð.Ý˜œ
¥B¤FÑ+Ô+Ð+åŒm˜Q˜C ! Ñ%Ô%ˆÝ˜œ¥r¤vÑ.Ô.Ð.Ý˜œ
¥B¤FÑ+Ô+Ð+õ Œm˜B Ñ#Ô#ˆÝ˜œ¥r¤vÑ.Ô.Ð.Ý˜œ
¥B¤FÑ+Ô+Ð+õ ŒI�c‰NŒNˆÝŒI�c‰NŒNˆÝ•j¥%¤-°°AÑ6Ô6Ð6Ð6Ð6r=   c                 ó–  — g d¢}g d¢}d}d}d}t          j        ||¦  «        }t          |j        |d¬¦  «         t          |j        |d¬¦  «         t          |j        j        d	¦  «         t          j        ||¦  «        }t          |j        |d¬¦  «         t          |j        |d¬¦  «         t          |j        j        d
¦  «         d S )N)r?   r?   r?   r@   r@   r@   r@   r@   rA   rA   r?   r@   r@   r@   r@   r@   r@   r@   rA   rA   rA   rA   rA   rA   )r?   r?   r?   r?   r?   r?   r?   r?   r?   r?   r@   r@   r@   r@   r@   r@   r@   r@   r@   r@   r@   r@   r@   r@   gCÑE]tÑ?gâ^ñ_ñÕ?gO((ÄÇ·?ra  r-   r/   )rA   r@   ©r@   rA   )r   r�  r   r\   r^   r   r  r6  )r5   r7   r8   Úd_crÚd_rcr:   r`   s          r;   Útest_asymmetryzTestSomersD.test_asymmetryÍ  sá   € ð1ð 1ð 1ˆð1ð 1ð 1ˆð !ˆØ ˆØˆÝŒm˜A˜qÑ!Ô!ˆÝ˜œ t°%Ð8Ñ8Ô8Ð8Ý˜œ
 A¨DÐ1Ñ1Ô1Ð1Ý�S”Y”_ fÑ-Ô-Ð-ÝŒm˜A˜qÑ!Ô!ˆÝ˜œ t°%Ð8Ñ8Ô8Ð8Ý˜œ
 A¨EÐ2Ñ2Ô2Ð2Ý�S”Y”_ fÑ-Ô-Ð-Ð-Ð-r=   c                 ó,  — t          j        ddgddgddgddgddgg¦  «        }|j        }d}t          t	          j        |¦  «        j        |¦  «         t          j        d	d
gdd
gd
dgg¦  «        }d\  }}t          t	          j        |¦  «        j        |¦  «         t          t	          j        |j        ¦  «        j        |¦  «         t          j        d	d
gd
dgdd
gg¦  «        }d}t          t	          j        |j        ¦  «        j        |¦  «         d S )NrG   r@   rD   rC   rA   rB   r?   gHHHHHHØ?r”  r   éU   rK   )gœÜMî&wã?rŒ   gtÑE]tá¿)rM   r   ÚTr   r   r�  r\   )r5   r  ÚdyxÚdxys       r;   Útest_somers_originalz TestSomersD.test_somers_originalæ  s  € õ
 ”˜1˜a˜& 1 a &¨1¨a¨&°1°a°&¸1¸a¸&ÐAÑBÔBˆà”ˆØˆÝ�œ eÑ,Ô,Ô6¸Ñ<Ô<Ð<õ ”˜2˜q˜' B¨ 7¨Q°¨GÐ4Ñ5Ô5ˆØ'‰ˆˆSÝ�œ eÑ,Ô,Ô6¸Ñ<Ô<Ð<Ý�œ e¤gÑ.Ô.Ô8¸#Ñ>Ô>Ð>õ ”˜2˜q˜' A r 7¨R°¨GÐ4Ñ5Ô5ˆØˆÝ�œ e¤gÑ.Ô.Ô8¸#Ñ>Ô>Ð>Ð>Ð>r=   c                 ób  — d}d}t          j        |¦  «        }t           j                             d¦  «        }t          j                             |t          j        |¦  «        |z  |¬¦  «                             |¦  «        }t	          j	        |¦  «        }t          j
        |dt          j        |d         ¦  «        d¬¦  «        }t	          j	        |¦  «        }t          j
        |dt          j        |d         ¦  «        d¬¦  «        }	t	          j	        |	¦  «        }
t          j
        |dt          j        |d         dz   ¦  «        d¬¦  «        }t	          j	        |¦  «        }t          |j        dd	¬
¦  «         t          |j        |j        ¦  «         t          |j        |
j        ¦  «         t          |j        |j        ¦  «         t          |j        dd	¬
¦  «         t          |j        |j        ¦  «         t          |j        |
j        ¦  «         t          |j        |j        ¦  «         d S )Nr�   ©rB   rD   r   ©r:   Úrandom_stater@   r?   ri   gaƒ¸yò½¿ra  r-   gP—j¹$Ä?)rM   Úprodrk   ÚRandomStater   ÚmultinomialÚrvsr˜   Úreshaper�  Úinsertr:  r   r\   r^   )r5   ÚNr6  re   rn   r  r`   Ús2r  Ús3rE  Ús4Úres4s                r;   Ú*test_contingency_table_with_zero_rows_colsz6TestSomersD.test_contingency_table_with_zero_rows_colsü  sÑ  € ð ˆØˆÝŒw�u‰~Œ~ˆåŒi×#Ò# AÑ&Ô&ˆÝÔ×!Ò! !¥r¤w¨t¡}¤}°TÑ'9Ø/2ð "ñ 4ô 4ß4;²G¸E±N´Nð 	
åŒm˜AÑÔˆåŒY�q˜!�RœX e¨A¤hÑ/Ô/°aÐ8Ñ8Ô8ˆÝŒ}˜RÑ Ô ˆåŒY�q˜!�RœX e¨A¤hÑ/Ô/°aÐ8Ñ8Ô8ˆÝŒ}˜RÑ Ô ˆåŒY�r˜1�bœh u¨Q¤x°¡zÑ2Ô2¸Ð;Ñ;Ô;ˆÝŒ}˜RÑ Ô ˆõ 	˜œÐ'9ÀÐFÑFÔFÐFÝ˜œ t¤~Ñ6Ô6Ð6Ý˜œ t¤~Ñ6Ô6Ð6Ý˜œ t¤~Ñ6Ô6Ð6å˜œ
Ð$5¸EÐBÑBÔBÐBÝ˜œ
 D¤KÑ0Ô0Ð0Ý˜œ
 D¤KÑ0Ô0Ð0Ý˜œ
 D¤KÑ0Ô0Ð0Ð0Ð0r=   c                 óz  — d}d}t          j        |¦  «        }t           j                             d¦  «        }t          j                             |t          j        |¦  «        |z  |¬¦  «                             |¦  «        }|dz
  }d}t          t          |¬¦  «        5  t	          j        |¦  «         d d d ¦  «         n# 1 swxY w Y   |dz   }d	}t          t          |¬¦  «        5  t	          j        |¦  «         d d d ¦  «         n# 1 swxY w Y   d
}t          t          |¬¦  «        5  t	          j        g g¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          |¬¦  «        5  t	          j        dgg¦  «         d d d ¦  «         n# 1 swxY w Y   t          j        d¦  «        }	t          t          |¬¦  «        5  t	          j        |	¦  «         d d d ¦  «         n# 1 swxY w Y   d|	d<   t          t          |¬¦  «        5  t	          j        |	¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr�   r°  r   r±  r@   z:All elements of the contingency table must be non-negativer¡   r}   z5All elements of the contingency table must be integerz?At least two elements of the contingency table must be nonzero.r?   )rA   rA   r„  )rM   r³  rk   rl   r   rµ  r¶  r˜   r·  rÈ   r¶   r�  r:  )
r5   r¹  r6  re   rn   r  Ús5r¸   Ús6Ús7s
             r;   Útest_invalid_contingency_tablesz+TestSomersD.test_invalid_contingency_tables  s—  € ØˆØˆÝŒw�u‰~Œ~ˆåŒi×#Ò# AÑ&Ô&ˆåÔ×!Ò! !¥r¤w¨t¡}¤}°TÑ'9Ø/2ð "ñ 4ô 4ß4;²G¸E±N´Nð 	
ð �‰UˆØNˆÝ�:¨WÐ5Ñ5Ô5ð 	ð 	ÝŒM˜"ÑÔÐð	ð 	ð 	ñ 	ô 	ð 	ð 	ð 	ð 	ð 	ð 	øøøð 	ð 	ð 	ð 	ð �‰XˆØIˆÝ�:¨WÐ5Ñ5Ô5ð 	ð 	ÝŒM˜"ÑÔÐð	ð 	ð 	ñ 	ô 	ð 	ð 	ð 	ð 	ð 	ð 	øøøð 	ð 	ð 	ð 	ð,ˆå�:¨WÐ5Ñ5Ô5ð 	 ð 	 ÝŒM˜2˜$ÑÔÐð	 ð 	 ð 	 ñ 	 ô 	 ð 	 ð 	 ð 	 ð 	 ð 	 ð 	 øøøð 	 ð 	 ð 	 ð 	 õ �:¨WÐ5Ñ5Ô5ð 	!ð 	!ÝŒM˜A˜3˜%Ñ Ô Ð ð	!ð 	!ð 	!ñ 	!ô 	!ð 	!ð 	!ð 	!ð 	!ð 	!ð 	!øøøð 	!ð 	!ð 	!ð 	!õ ŒX�fÑÔˆÝ�:¨WÐ5Ñ5Ô5ð 	ð 	ÝŒM˜"ÑÔÐð	ð 	ð 	ñ 	ô 	ð 	ð 	ð 	ð 	ð 	ð 	øøøð 	ð 	ð 	ð 	ð ˆˆ4‰Ý�:¨WÐ5Ñ5Ô5ð 	ð 	ÝŒM˜"ÑÔÐð	ð 	ð 	ñ 	ô 	ð 	ð 	ð 	ð 	ð 	ð 	ð 	øøøð 	ð 	ð 	ð 	ð 	ð 	sl   ÂC Ã CÃCÃ(D	Ä	DÄDÄ,EÅEÅEÅ/FÆFÆFÇG(Ç(G,Ç/G,ÈH0È0H4È7H4c                 ó  — g d¢}ddt           j        g}g d¢}ddt           j         g}t          j        ||¦  «        }t          j        ||¦  «        }t	          |j        |j        ¦  «         t	          |j        |j        ¦  «         d S )Nr²   rh   gÍÌÌÌÌÌ @)rA   r@   r?   r   rž  )rM   rm   r   r�  r   r\   r^   )r5   r7   r¢  r8   Úy2r`   r  s          r;   Útest_only_ranks_matterz"TestSomersD.test_only_ranks_matter@  s‚   € àˆIˆIˆØ�#•r”vÐˆØˆIˆIˆØ��œ�wÐˆÝŒm˜A˜qÑ!Ô!ˆÝŒ}˜R Ñ$Ô$ˆÝ�S”] D¤NÑ3Ô3Ð3Ý�S”Z ¤Ñ-Ô-Ð-Ð-Ð-r=   c                 óÎ   — t          j        d¦  «        }t          j        d¦  «        }t          j        ||¦  «        }t	          |j        t          j        d¦  «        ¦  «         d S )NrE   )rM   rN   r   r�  r   r  Úeye©r5   r7   r8   r`   s       r;   Útest_contingency_table_returnz)TestSomersD.test_contingency_table_returnK  sM   € åŒI�b‰MŒMˆÝŒI�b‰MŒMˆÝŒm˜A˜qÑ!Ô!ˆÝ�S”Y¥¤ r¡
¤
Ñ+Ô+Ð+Ð+Ð+r=   c                 óà  — g d¢}g d¢}t          j        ||d¬¦  «        }|j        dk    sJ ‚t          j        ||d¬¦  «        }t          |j        |j        ¦  «         t	          |j        d|j        dz  z
  ¦  «         t          j        ||d	¬¦  «        }t          |j        |j        ¦  «         t	          |j        |j        dz  ¦  «         |                     ¦   «          t          j        ||d¬¦  «        }|j        dk     sJ ‚t          j        ||d	¬¦  «        }t          |j        |j        ¦  «         t	          |j        d|j        dz  z
  ¦  «         t          j        ||d¬¦  «        }t          |j        |j        ¦  «         t	          |j        |j        dz  ¦  «         t          j        t          d
¬¦  «        5  t          j        ||d¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr  )rC   rD   rF   rG   rF   rÒ   rÄ   r   rÖ   r?   r@   rØ   z`alternative` must be...r¡   ú	ekki-ekki)
r   r�  r\   r   r   r^   Úreverser1   r   r¶   )r5   r¡  r¢  rÛ   r`   s        r;   Útest_somersd_alternativez$TestSomersD.test_somersd_alternativeR  s0  € ð ˆ_ˆ_ˆØˆ_ˆ_ˆõ ”=  R°[ÐAÑAÔAˆØÔ! AÒ%Ð%Ð%Ð%õ Œm˜B °Ð7Ñ7Ô7ˆÝ�S”] HÔ$6Ñ7Ô7Ð7Ý˜œ
 A¨¬¸1Ñ)<Ñ$=Ñ>Ô>Ð>õ Œm˜B °	Ð:Ñ:Ô:ˆÝ�S”] HÔ$6Ñ7Ô7Ð7Ý˜œ
 H¤O°aÑ$7Ñ8Ô8Ð8ð 	�
Š
‰Œˆõ ”=  R°[ÐAÑAÔAˆØÔ! AÒ%Ð%Ð%Ð%õ Œm˜B °	Ð:Ñ:Ô:ˆÝ�S”] HÔ$6Ñ7Ô7Ð7Ý˜œ
 A¨¬¸1Ñ)<Ñ$=Ñ>Ô>Ð>õ Œm˜B °Ð7Ñ7Ô7ˆÝ�S”] HÔ$6Ñ7Ô7Ð7Ý˜œ
 H¤O°aÑ$7Ñ8Ô8Ð8åŒ]�:Ð-GÐHÑHÔHð 	;ð 	;ÝŒM˜"˜b¨kÐ:Ñ:Ô:Ð:ð	;ð 	;ð 	;ñ 	;ô 	;ð 	;ð 	;ð 	;ð 	;ð 	;ð 	;ð 	;øøøð 	;ð 	;ð 	;ð 	;ð 	;ð 	;s   Æ>G#Ç#G'Ç*G'Úpositive_correlation)FTc                 ó   — t          j        d¦  «        }|r|nt          j        |¦  «        }|rdnd}t          j        ||d¬¦  «        }|j        |k    sJ ‚|j        dk    sJ ‚t          j        ||d¬¦  «        }|j        |k    sJ ‚|j        |rdndk    sJ ‚t          j        ||d¬¦  «        }|j        |k    sJ ‚|j        |rdndk    sJ ‚d S )	NrE   r?   rh   rÒ   rÄ   r   rÖ   rØ   )rM   rN   Úflipr   r�  r\   r^   )r5   rÏ  r¡  r¢  Úexpected_statisticr`   s         r;   Ú test_somersd_perfect_correlationz,TestSomersD.test_somersd_perfect_correlation{  s	  € õ ŒY�r‰]Œ]ˆØ'Ð8ˆRˆR­R¬W°R©[¬[ˆØ"6Ð>˜Q˜Q¸BÐõ Œm˜B °Ð<Ñ<Ô<ˆØŒ}Ð 2Ò2Ð2Ð2Ð2ØŒz˜QŠˆˆˆõ Œm˜B °Ð7Ñ7Ô7ˆØŒ}Ð 2Ò2Ð2Ð2Ð2ØŒzÐ#7Ð>˜a˜a¸QÒ?Ð?Ð?Ð?õ Œm˜B °	Ð:Ñ:Ô:ˆØŒ}Ð 2Ò2Ð2Ð2Ð2ØŒzÐ#7Ð>˜a˜a¸QÒ?Ð?Ð?Ð?Ð?Ð?r=   c                 ó  — ddg}d}t           j                             d¦  «        }|                     ||¦  «        }|                     ||¦  «        }d}t	          j        ||¦  «        j        }t          ||d¬¦  «         d S )Nr?   r@   é@B l   ÿ.E5 g 0ôÜuD?ra  r-   )rM   rk   rl   Úchoicer   r�  r\   r   )r5   ÚclassesÚ	n_samplesrn   r7   r8   Úval_sklearnÚ	val_scipys           r;   Ú!test_somersd_large_inputs_gh18132z-TestSomersD.test_somersd_large_inputs_gh18132”  s…   € ð �a�&ˆØˆ	ÝŒi×#Ò# JÑ/Ô/ˆØ�JŠJ�w 	Ñ*Ô*ˆØ�JŠJ�w 	Ñ*Ô*ˆð +ˆõ ”M ! QÑ'Ô'Ô1ˆ	Ý˜ Y°UÐ;Ñ;Ô;Ð;Ð;Ð;r=   N)ru   rv   rw   r�  r’  r›  r£  r¨  r®  r¾  rÃ  rÆ  rÊ  rÎ  r1   rx   ry   rÓ  rÛ  rz   r=   r;   r‚  r‚  >  s  € € € € € ð7ð 7ð 7ðFð Fð Fð
	>ð 	>ð 	>ðm7ð m7ð m7ð^.ð .ð .ð2?ð ?ð ?ð,1ð 1ð 1ð@"ð "ð "ðH	.ð 	.ð 	.ð,ð ,ð ,ð';ð ';ð ';ðR „[×ÒÐ3°]ÑCÔCð@ð @ñ DÔCð@ð0<ð <ð <ð <ð <r=   r‚  c                   óÚ  — e Zd ZdZej                             dddgddggdfdd	gd
dggdfd	dgdd	ggdfddgddggdfddgddggdfddgddggdfddgddggdfddgddggdfddgdd	ggdfdd	gddggdfd	dgdd	ggdfg¦  «        d „ ¦   «         Zej                             dddgddggd!fdd	gd
dggd"fd	dgdd	ggd#fddgddggd$fddgddggd%fddgddggd&fddgddggd'fddgddggd(fddgdd	ggd)fdd	gddggd*fd	dgdd	ggd#fg¦  «        d+„ ¦   «         Zd,„ Z	ej                             dddgddggd-fg¦  «        d.„ ¦   «         Z
ej                             dddgddggd/ej        ffddgddggd/ej        ffg¦  «        d0„ ¦   «         Zej                             dd	dgdd	ggd1fdd2gd3dggd4fd5d6gd7dggd8fg¦  «        ej                             d9d:d;g¦  «        d<„ ¦   «         ¦   «         Zd=S )>ÚTestBarnardExactz8Some tests to show that barnard_exact() works correctly.úinput_sample,expectedé+   r  rE   é'   )gÏXyq@g{2s&QÇ7?r�   r@   r€   rC   )gøl¯€lü¿gEA]0×KÁ?rF   rG   )ç*õ)1Ö%ÀgØ_  ¬“?r?   )g_½Òc1÷?gèß=ç Ä?é   é   )g5PyQ ý¿g�Q@Ì2ý°?é   r”  )g»ÔÌÑÁù¿g¼ÂJ"¸ò½?)g_½Òc1÷¿gwÝ�Ù„¹Æ?r   rB   )gý7ŠËé§@g      Œ?rA   )gš~øtñ†ñ¿ç,À�?3Oà?rD   )gÝr?ô~Éø¿çCæ°Yð7É?c                 ód   — t          |¦  «        }|j        |j        }}t          ||g|¦  «         dS )zþThe expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-6 :
        ```R
        library(Barnard)
        options(digits=10)
        barnard.test(43, 40, 10, 39, dp=1e-6, pooled=TRUE)
        ```
        N©r   r\   r^   r   ©r5   Úinput_samplerÛ   r`   r\   r^   s         r;   Útest_precisezTestBarnardExact.test_precise°  s9   € õ2 ˜LÑ)Ô)ˆØœM¨3¬:�6ˆ	Ý˜ FÐ+¨XÑ6Ô6Ð6Ð6Ð6r=   )gµ’¹7ç\@gAÞÄƒø2?)gXS‘ç;ò¿gÙŽþ hî?)g>!ÉŽ¢Àg6¦  ¨”?)gS¡y@õù?gúÒ^¶ÛFÃ?)g‘ª-º©˜ÿ¿g¼X«I#°?)g°a£ÑÐ�û¿goÇÜðÁ?)gbø]?ü¿gFug¹H	Ð?)g§6Ò­è@g      €?)gi(	Ž˜ó¿rå  )gÓNXèz¶û¿ræ  c                 óh   — t          |d¬¦  «        }|j        |j        }}t          ||g|¦  «         dS )zÿThe expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-6 :
        ```R
        library(Barnard)
        options(digits=10)
        barnard.test(43, 40, 10, 39, dp=1e-6, pooled=FALSE)
        ```
        F)ÚpooledNrè  ré  s         r;   Útest_pooled_paramz"TestBarnardExact.test_pooled_paramÍ  s>   € õ2 ˜L°Ð7Ñ7Ô7ˆØœM¨3¬:�6ˆ	Ý˜ FÐ+¨XÑ6Ô6Ð6Ð6Ð6r=   c                 óv  — d}t          t          |¬¦  «        5  t          ddgddggd¬¦  «         d d d ¦  «         n# 1 swxY w Y   d	}t          t          |¬¦  «        5  t          t          j        d
¦  «                             dd¦  «        ¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          ddgddgg¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          ddgddggd¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nú7Number of points `n` must be strictly positive, found 0r¡   r?   r@   rA   rB   r   ©r‡   ú,The input `table` must be of shape \(2, 2\).rD   ú*All values in `table` must be nonnegative.rh   zI`alternative` should be one of {'two-sided', 'less', 'greater'}, found .*únot-correct)rÈ   r¶   r   rM   rN   r·  ©r5   Ú	error_msgs     r;   Útest_raiseszTestBarnardExact.test_raisesê  sX  € ð Fð 	õ �:¨YÐ7Ñ7Ô7ð 	1ð 	1Ý˜A˜q˜6 A q 6Ð*¨aÐ0Ñ0Ô0Ð0ð	1ð 	1ð 	1ñ 	1ô 	1ð 	1ð 	1ð 	1ð 	1ð 	1ð 	1øøøð 	1ð 	1ð 	1ð 	1ð Eˆ	Ý�:¨YÐ7Ñ7Ô7ð 	6ð 	6Ý�"œ) A™,œ,×.Ò.¨q°!Ñ4Ô4Ñ5Ô5Ð5ð	6ð 	6ð 	6ñ 	6ô 	6ð 	6ð 	6ð 	6ð 	6ð 	6ð 	6øøøð 	6ð 	6ð 	6ð 	6ð Aˆ	Ý�:¨YÐ7Ñ7Ô7ð 	-ð 	-Ý˜B ˜7 Q¨ FÐ+Ñ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð
ð 	õ �:¨YÐ7Ñ7Ô7ð 	;ð 	;Ý˜A˜q˜6 A q 6Ð*¨MÑ:Ô:Ð:ð	;ð 	;ð 	;ñ 	;ô 	;ð 	;ð 	;ð 	;ð 	;ð 	;ð 	;ð 	;øøøð 	;ð 	;ð 	;ð 	;ð 	;ð 	;óE   ™=½AÁAÁ 6B"Â"B&Â)B&ÃC'Ã'C+Ã.C+Ä
D.Ä.D2Ä5D2r�  c                 ó˜   — t          |¦  «        }|j        |j        }}t          ||d         ¦  «         t          ||d         ¦  «         d S ©Nr   r?   ©r   r\   r^   r   ré  s         r;   Útest_edge_casesz TestBarnardExact.test_edge_cases  sL   € õ ˜LÑ)Ô)ˆØœM¨3¬:�6ˆ	Ý�V˜X aœ[Ñ)Ô)Ð)Ý�Y ¨¤Ñ,Ô,Ð,Ð,Ð,r=   rŒ   c                 ó˜   — t          |¦  «        }|j        |j        }}t          ||d         ¦  «         t          ||d         ¦  «         d S rú  rû  ré  s         r;   Útest_row_or_col_zeroz%TestBarnardExact.test_row_or_col_zero  sL   € õ ˜LÑ)Ô)ˆØœM¨3¬:�6ˆ	Ý�V˜X aœ[Ñ)Ô)Ð)Ý�Y ¨¤Ñ,Ô,Ð,Ð,Ð,r=   )rá  gEç\/?„?éÈ   é,  )gÐgÑìíQ5Àr‹   é   rL   i¥  )gü&çìX}>Àr‹   rÅ   rØ   rÖ   c                 óÎ   — |\  }}|dk    r$t          j        |¦  «        dd…ddd…f         }| }t          ||¬¦  «        }|j        |j        }}t          ||g||gd¬¦  «         dS )aˆ  
        "The expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-6 :
        ```R
        library(Barnard)
        options(digits=10)
        a = barnard.test(2, 7, 8, 2, dp=1e-6, pooled=TRUE)
        a$p.value[1]
        ```
        In this test, we are using the "one-sided" return value `a$p.value[1]`
        to test our pvalue.
        rØ   Nrh   rÄ   çH¯¼šò×z>r-   )rM   r   r   r\   r^   r   )	r5   rê  rÛ   rÅ   Úexpected_statÚless_pvalue_expectr`   r\   r^   s	            r;   Útest_less_greaterz"TestBarnardExact.test_less_greater  s•   € ð, -5Ñ)ˆÐ)à˜)Ò#Ð#Ýœ8 LÑ1Ô1°!°!°!°T°T°r°T°'Ô:ˆLØ*˜NˆMå˜L°kÐBÑBÔBˆØœM¨3¬:�6ˆ	ÝØ˜Ð -Ð1CÐ!DÈ4ð	
ñ 	
ô 	
ð 	
ð 	
ð 	
r=   N)ru   rv   rw   Ú__doc__r1   rx   ry   rë  rî  r÷  rü  rM   r]   rþ  r  rz   r=   r;   rÝ  rÝ  ­  sÕ  € € € € € ØBÐBà„[×ÒØà�2ˆh˜˜R˜Ð!Ð#CÐDØ�Aˆh˜˜q˜	Ð"Ð$EÐFØ�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�r˜2�hÐÐ!AÐBØ�"ˆg˜˜B�xÐ Ð"CÐDØ�"ˆg˜˜B�xÐ Ð"CÐDØ�1ˆg˜˜A�wÐÐ!BÐCØ�!ˆf�q˜!�fÐÐ?Ð@Ø�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�q˜!�fÐÐ@ÐAð	
ñô ð 7ð 7ñ!ô ð 7ð „[×ÒØà�2ˆh˜˜R˜Ð!Ð#CÐDØ�Aˆh˜˜q˜	Ð"Ð$EÐFØ�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�r˜2�hÐÐ!AÐBØ�"ˆg˜˜B�xÐ Ð"CÐDØ�"ˆg˜˜B�xÐ Ð"CÐDØ�1ˆg˜˜A�wÐÐ!BÐCØ�!ˆf�q˜!�fÐÐ?Ð@Ø�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�q˜!�fÐÐ@ÐAØ�!ˆf�q˜!�fÐÐ@ÐAð	
ñô ð 7ð 7ñ!ô ð 7ð;ð ;ð ;ð4 „[×ÒØà�!ˆf�q˜!�fÐ˜xÐ(ð	
ñô ð-ð -ñô ð-ð „[×ÒØà�!ˆf�q˜"�gÐ  b¤f Ð.Ø�!ˆf�r˜1�gÐ  b¤f Ð.ð	
ñô ð-ð -ñô ð-ð „[×ÒØà�!ˆf�q˜!�fÐÐ@ÐAØ�#ˆh˜˜a˜Ð!Ð#:Ð;Ø�2ˆh˜˜q˜	Ð"Ð$;Ð<ð	
ñô ð „[×Ò˜]¨Y¸Ð,?Ñ@Ô@ð
ð 
ñ AÔ@ñô ð
ð 
ð 
r=   rÝ  c                   óÖ  — e Zd ZdZdZej                             dddgddggdfdd	gd
d
ggdfddgddggdfd
dgd
d	ggdfddgd	dggdfdd	gddggdfddgddggdfddgddggdfd
dgddggdfg	¦  «        d„ ¦   «         Zej                             dddgd
dggdfddgddggd fdd	gd
d
ggd!fdd"gddggd#fddgddggd$fddgd	dggd%fdd	gddggdfddgd&dggdfddgddggd fddgddggd'fd
dgddggd(fg¦  «        d)„ ¦   «         Z	ej                             dddgd
dggd*fddgddggd+fdd	gd
d
ggd,fddgddggd-fddgd	dggd.fdd	gddggd/fddgddggd+fddgddggd0fg¦  «        d1„ ¦   «         Z
d2„ Zej                             dddgdd
ggej        ej        ffddgd
dggej        ej        ffg¦  «        d3„ ¦   «         Zd4„ Zej                             d5d6¦  «        d7„ ¦   «         Zd8S )9ÚTestBoschlooExactz9Some tests to show that boschloo_exact() works correctly.r  rÞ  r@   rF   rG   )ç<vB\÷’?gèÁ–„/?„?rC   r?   rE   )g¬“ŽÍéMï?gßÏA>î?rä  rã  r”  )ç_¬VÃÑ¶?gÂñ¥…Ö­?)ç‘uÝ Ø%Å?gc'�íµ?r   rB   ©r?   r?   rA   )r   g      Ö?rf   )çß+Âf°?gXc}Áv¡?é   é%   )gZÑ‹D¸É?g‰¦¢gi]Ã?c                 óv   — t          |d¬¦  «        }|j        |j        }}t          ||g|| j        ¬¦  «         dS )a‰  The expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-8 :
        ```R
        library(Exact)
        options(digits=10)
        data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE)
        a = exact.test(data, method="Boschloo", alternative="less",
                       tsmethod="central", np.interval=TRUE, beta=1e-8)
        ```
        rÖ   rÄ   r-   N©r   r\   r^   r   ÚATOLré  s         r;   Ú	test_lesszTestBoschlooExact.test_lessE  sE   € õ2 ˜\°vÐ>Ñ>Ô>ˆØœM¨3¬:�6ˆ	Ý˜ FÐ+¨X¸D¼IÐFÑFÔFÐFÐFÐFr=   rß  r  rà  )çòk¾\Ø2?g–ìà0í,%?)gKïvî÷ï?gâN3î½ï?)r  g‡'&5Õ·?râ  )gÏw@_„ï?g™•7Ñøï?)gµ‹i¦{ï?g€ÁÉ‘)zî?)ç˜ïÖ…a˜?g³1×ëÛ|?rD   )gYðì<;ªï?g¨ýÖN”Dï?)gê­e×ì?g™Gþ` ë?c                 óv   — t          |d¬¦  «        }|j        |j        }}t          ||g|| j        ¬¦  «         dS )aŒ  The expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-8 :
        ```R
        library(Exact)
        options(digits=10)
        data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE)
        a = exact.test(data, method="Boschloo", alternative="greater",
                       tsmethod="central", np.interval=TRUE, beta=1e-8)
        ```
        rØ   rÄ   r-   Nr  ré  s         r;   Útest_greaterzTestBoschlooExact.test_greaterb  sE   € õ6 ˜\°yÐAÑAÔAˆØœM¨3¬:�6ˆ	Ý˜ FÐ+¨X¸D¼IÐFÑFÔFÐFÐFÐFr=   )r  gôqQSé,5?)r
  gG’Þ?/?”?)r  gKE`ÕÇ?)r  gÓhr1Ö½?)r  grÒfbÛŒ?)r   g      æ?)r  gP:pRÁv±?c                 óx   — t          |dd¬¦  «        }|j        |j        }}t          ||g|| j        ¬¦  «         dS )aŽ  The expected values have been generated by R, using a resolution
        for the nuisance parameter of 1e-8 :
        ```R
        library(Exact)
        options(digits=10)
        data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE)
        a = exact.test(data, method="Boschloo", alternative="two.sided",
                       tsmethod="central", np.interval=TRUE, beta=1e-8)
        ```
        rÒ   é@   )rÅ   r‡   r-   Nr  ré  s         r;   Útest_two_sidedz TestBoschlooExact.test_two_sided�  sG   € õ0 ˜\°{ÀbÐIÑIÔIˆàœM¨3¬:�6ˆ	Ý˜ FÐ+¨X¸D¼IÐFÑFÔFÐFÐFÐFr=   c                 óv  — d}t          t          |¬¦  «        5  t          ddgddggd¬¦  «         d d d ¦  «         n# 1 swxY w Y   d	}t          t          |¬¦  «        5  t          t          j        d
¦  «                             dd¦  «        ¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          ddgddgg¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          ddgddggd¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nrð  r¡   r?   r@   rA   rB   r   rñ  rò  rD   ró  rh   zK`alternative` should be one of \('two-sided', 'less', 'greater'\), found .*rô  )rÈ   r¶   r   rM   rN   r·  rõ  s     r;   r÷  zTestBoschlooExact.test_raisesž  sX  € ð Fð 	õ �:¨YÐ7Ñ7Ô7ð 	2ð 	2Ý˜Q ˜F Q¨ FÐ+¨qÐ1Ñ1Ô1Ð1ð	2ð 	2ð 	2ñ 	2ô 	2ð 	2ð 	2ð 	2ð 	2ð 	2ð 	2øøøð 	2ð 	2ð 	2ð 	2ð Eˆ	Ý�:¨YÐ7Ñ7Ô7ð 	7ð 	7Ý�2œ9 Q™<œ<×/Ò/°°1Ñ5Ô5Ñ6Ô6Ð6ð	7ð 	7ð 	7ñ 	7ô 	7ð 	7ð 	7ð 	7ð 	7ð 	7ð 	7øøøð 	7ð 	7ð 	7ð 	7ð Aˆ	Ý�:¨YÐ7Ñ7Ô7ð 	.ð 	.Ý˜R ˜G a¨ VÐ,Ñ-Ô-Ð-ð	.ð 	.ð 	.ñ 	.ô 	.ð 	.ð 	.ð 	.ð 	.ð 	.ð 	.øøøð 	.ð 	.ð 	.ð 	.ð
%ð 	õ �:¨YÐ7Ñ7Ô7ð 	<ð 	<Ý˜Q ˜F Q¨ FÐ+¨]Ñ;Ô;Ð;ð	<ð 	<ð 	<ñ 	<ô 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<øøøð 	<ð 	<ð 	<ð 	<ð 	<ð 	<rø  c                 ó˜   — t          |¦  «        }|j        |j        }}t          ||d         ¦  «         t          ||d         ¦  «         d S rú  )r   r\   r^   r   ré  s         r;   rþ  z&TestBoschlooExact.test_row_or_col_zero¸  sL   € õ ˜\Ñ*Ô*ˆØœM¨3¬:�6ˆ	Ý�V˜X aœ[Ñ)Ô)Ð)Ý�Y ¨¤Ñ,Ô,Ð,Ð,Ð,r=   c                 óÜ   — ddgddgg}t          |d¬¦  «        j        }t          |d¬¦  «        j        }dt          ||¦  «        z  dk    sJ ‚t          |d¬¦  «        j        }|d	k    sJ ‚d S )
Nr?   rª   rf   rÖ   rÄ   rØ   r@   rÒ   rŒ   )r   r^   r^  )r5   ÚtblÚplÚpgÚpts        r;   Útest_two_sided_gt_1z%TestBoschlooExact.test_two_sided_gt_1Å  sƒ   € ð �1ˆv˜˜B�xÐ ˆÝ˜C¨VÐ4Ñ4Ô4Ô;ˆÝ˜C¨YÐ7Ñ7Ô7Ô>ˆØ•�R˜‘”‰}˜qÒ Ð Ð Ð Ý˜C¨[Ð9Ñ9Ô9Ô@ˆØ�SŠyˆyˆyˆyˆyˆyr=   rÅ   )rÖ   rØ   c                 óš   — ddgddgg}t          ||¬¦  «        j        }t          j        ||¬¦  «        d         }t	          ||¦  «         d S )Nr@   rF   rG   rÄ   r?   )r   r\   r   Úfisher_exactr   )r5   rÅ   r  Úboschloo_statÚfisher_ps        r;   Útest_against_fisher_exactz+TestBoschlooExact.test_against_fisher_exactÏ  s[   € ð �1ˆv˜˜1�vÐˆÝ& s¸ÐDÑDÔDÔNˆÝÔ% c°{ÐCÑCÔCÀAÔFˆÝ˜ xÑ0Ô0Ð0Ð0Ð0r=   N)ru   rv   rw   r  r  r1   rx   ry   r  r  r  r÷  rM   r]   rþ  r#  r(  rz   r=   r;   r	  r	  @  só  € € € € € ØCÐCà€Dà„[×ÒØà�!ˆf�q˜!�fÐÐ8Ð9Ø�!ˆf�r˜2�hÐÐ!7Ð8Ø�"ˆg˜˜B�xÐ Ð":Ð;Ø�1ˆg˜˜A�wÐÐ!8Ð9Ø�!ˆf�q˜!�fÐ˜vÐ&Ø�!ˆf�q˜!�fÐ˜~Ð.Ø�!ˆf�q˜!�fÐÐ8Ð9Ø�"ˆg˜˜1�vÐÐ 8Ð9Ø�2ˆh˜˜R˜Ð!Ð#9Ð:ð
	
ñô ðGð Gñô ðGð „[×ÒØà�2ˆh˜˜R˜Ð!Ð#?Ð@Ø�!ˆf�q˜!�fÐÐ5Ð6Ø�!ˆf�r˜2�hÐÐ!8Ð9Ø�"ˆg˜˜B�xÐ Ð"8Ð9Ø�"ˆg˜˜B�xÐ Ð"7Ð8Ø�!ˆf�q˜!�fÐÐ8Ð9Ø�!ˆf�q˜!�fÐ˜vÐ&Ø�!ˆf�q˜!�fÐ˜vÐ&Ø�!ˆf�q˜!�fÐÐ5Ð6Ø�"ˆg˜˜1�vÐÐ 6Ð7Ø�2ˆh˜˜R˜Ð!Ð#9Ð:ð	
ñô ð Gð Gñ!ô ð Gð „[×ÒØà�2ˆh˜˜R˜Ð!Ð#?Ð@Ø�!ˆf�q˜!�fÐÐ7Ð8Ø�!ˆf�r˜2�hÐÐ!7Ð8Ø�"ˆg˜˜B�xÐ Ð"8Ð9Ø�!ˆf�q˜!�fÐÐ7Ð8Ø�!ˆf�q˜!�fÐ˜}Ð-Ø�!ˆf�q˜!�fÐÐ7Ð8Ø�"ˆg˜˜1�vÐÐ 8Ð9ð		
ñô ðGð Gñô ðGð <ð <ð <ð4 „[×ÒØà�!ˆf�q˜"�gÐ ¤¨¬Ð 0Ð1Ø�!ˆf�r˜1�gÐ ¤¨¬Ð 0Ð1ð	
ñô ð-ð -ñô ð-ðð ð ð „[×Ò˜]Ð,?Ñ@Ô@ð1ð 1ñ AÔ@ð1ð 1ð 1r=   r	  c                   óÀ  — e Zd Zej                             dg  ej        d¦  «        f ej        d¦  «        dgfg¦  «        ej                             dd¬¦  «        d„ ¦   «         ¦   «         Z	d„ Z
d	„ Zej                             d
g d¢¦  «        d„ ¦   «         Zej                             dg d¢¦  «        d„ ¦   «         Zej        j        d„ ¦   «         Zd„ Zd„ Zd„ ZdS )ÚTestCvm_2sampr¹   rC   r?   r›   zlazy -> no axis_nan_policyr�   c                 óT  ‡— ˆfd„|D ¦   «         }t          t          t          ‰¬¦  «        5  t          |Ž }t	          |j        ‰                     ‰j        ¦  «        ¦  «         t	          |j        ‰                     ‰j        ¦  «        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nc              3   ó`   •K  — | ](}‰                      |t          ‰¦  «        ¬ ¦  «        V — Œ)dS )r*   N)r4   r   )r3  Úargr6   s     €r;   ú	<genexpr>z5TestCvm_2samp.test_too_small_input.<locals>.<genexpr>ß  s<   øè è € ÐLÐLÀ�—
’
˜3Õ&6°rÑ&:Ô&:�
Ñ;Ô;ÐLÐLÐLÐLÐLÐLr=   rZ   )	r   r   r   r   r   r\   r4   r]   r^   )r5   r¹   r6   r`   s     ` r;   Útest_too_small_inputz"TestCvm_2samp.test_too_small_inputÛ  sä   ø€ ð MÐLÐLÐLÀtÐLÑLÔLˆÝÕ+Õ3HÈRÐPÑPÔPð 	<ð 	<Ý&¨Ð-ˆCÝ˜CœM¨2¯:ª:°b´fÑ+=Ô+=Ñ>Ô>Ð>Ý˜CœJ¨¯
ª
°2´6Ñ(:Ô(:Ñ;Ô;Ð;ð	<ð 	<ð 	<ñ 	<ô 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<ð 	<øøøð 	<ð 	<ð 	<ð 	<ð 	<ð 	<s   ¬A$BÂB!Â$B!c                 ó¾   — |                      d¦  «        }d}t          j        t          |¬¦  «        5  t	          ||d¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )NrC   z/method must be either auto, exact or asymptoticr¡   Úxyz)rN   r1   r   r¶   r   )r5   r6   r8   Úmsgs       r;   r£   z TestCvm_2samp.test_invalid_inputå  s¡   € Ø�IŠI�a‰LŒLˆØ?ˆÝŒ]�:¨SÐ1Ñ1Ô1ð 	.ð 	.Ý   A uÑ-Ô-Ð-ð	.ð 	.ð 	.ñ 	.ô 	.ð 	.ð 	.ð 	.ð 	.ð 	.ð 	.ð 	.øøøð 	.ð 	.ð 	.ð 	.ð 	.ð 	.s   ³AÁAÁAc                 óî   — g d¢}g d¢}t          ||¦  «        }t          t          j        |¦  «        t          j        |¦  «        ¦  «        }t          |j        |j        f|j        |j        f¦  «         d S )N)r@   rA   rB   rF   rD   )rî   r´   rf   r–  )r   rM   r   r   r\   r^   )r5   r7   r8   Úr1Úr2s        r;   Útest_list_inputzTestCvm_2samp.test_list_inputë  sl   € ØˆOˆOˆØÐÐˆÝ! ! QÑ'Ô'ˆÝ!¥"¤(¨1¡+¤+­r¬x¸©{¬{Ñ;Ô;ˆÝ�b”l B¤IÐ.°´¸r¼yÐ0IÑJÔJÐJÐJÐJr=   r!   r"   c                 óä  — t          |¦  «        r%|j        dk     r|dk    rt          j        d¦  «         |€t	          |¦  «        nt          ||¦  «        }|                     g d¢|¬¦  «        }|                     g d¢|¬¦  «        }t          ||¦  «        }t          |j	        |                     d|¬¦  «        d¬	¦  «         t          |j
        |                     d
|¬¦  «        d¬	¦  «         d S )Nr&   r#   r'   )	gffffff@gÍÌÌÌÌÌ @r,   gffffff!@çš™™™™™"@gÍÌÌÌÌÌ#@g333333$@g333333%@gffffff&@r*   )gÍÌÌÌÌÌ@gÍÌÌÌÌÌ@gš™™™™™@ç      @ç333333@gffffff @g333333"@gš™™™™™#@gš™™™™™%@gš™™™™™&@g      '@gš™™™™™(@g      )@gÍÌÌÌÌÌ*@g333333-@gøSã¥›ÄÐ?rH   r-   g
×£p=
Ç?r}   )r   r0   r1   r2   r   r3   r4   r   r   r\   r^   )r5   r!   r6   r7   r8   Úrs         r;   Útest_example_conoverz"TestCvm_2samp.test_example_conoverò  s  € õ �B‰<Œ<ð 	=˜BœN¨UÒ2Ð2°uÀ	Ò7IÐ7IÝŒKÐ;Ñ<Ô<Ð<Ø(-¨Õ  Ñ$Ô$Ð$½7À2ÀuÑ;MÔ;MˆØ�JŠJÐGÐGÐGÈuˆJÑUÔUˆØ�JŠJð Bð Bð BØINð ñ Pô Pˆå   AÑ&Ô&ˆÝ˜œ R§Z¢Z°¸U ZÑ%CÔ%CÈ$ÐOÑOÔOÐOÝ˜œ "§*¢*¨S¸ *Ñ">Ô">ÀTÐJÑJÔJÐJÐJÐJr=   zstatistic, m, n, pval))iÆ  rC   rD   g¾cj`ï˜º?)ii  rF   rF   gtÑE]t±?)i@  rB   rD   g8�8�ƒ?)iä  rD   rF   g¶ëéXwS?c                 óD   — t          t          |||¦  «        |¦  «         d S r®   )r   r   )r5   r\   r  r‡   Úpvals        r;   Útest_exact_pvaluezTestCvm_2samp.test_exact_pvalue   s'   € õ 	Õ*¨9°a¸Ñ;Ô;¸TÑBÔBÐBÐBÐBr=   c                 óº  — t           j                             d¦  «        }t          j                             d|¬¦  «        }t          j                             d|¬¦  «        }|                     |¦  «        |                     |¦  «        }}t          ||¦  «        }d|j        cxk     rdk     sn J ‚t          ||dz   ¦  «        }d|j        cxk     rdk     sn J ‚d S )Ni  rÕ  )re   r²  i » r   r?   r¥   )	rM   rk   rl   r
   Únormr¶  r4   r   r^   )r5   r6   rn   r7   r8   r;  s         r;   Útest_large_samplezTestCvm_2samp.test_large_sample  sâ   € õ Œi×#Ò# DÑ)Ô)ˆÝÔ×"Ò"¨¸cÐ"ÑBÔBˆÝÔ×"Ò"¨¸SÐ"ÑAÔAˆØ�zŠz˜!‰}Œ}˜bŸjšj¨™mœmˆ1ˆÝ   AÑ&Ô&ˆØ�1”8ÐÐÒÐ˜aÒÐÐÐÐÐÝ   A c¡EÑ*Ô*ˆØ�1”8ÐÐÒÐ˜aÒÐÐÐÐÐÐÐr=   c                 ó˜  — t           j                             d¦  «        }|                     |                     d¦  «        ¦  «        }|                     |                     d¦  «        ¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }t          |j        |j        ¦  «         t          |j        |j        d¬¦  «         d S )	Nr   rF   rG   rÌ   rÆ   rË   r}   r-   )	rM   rk   r´  r4   r   r   r\   r   r^   )r5   r6   rn   r7   r8   r4  r5  s          r;   Útest_exact_vs_asymptoticz&TestCvm_2samp.test_exact_vs_asymptotic  s¨   € ÝŒi×#Ò# AÑ&Ô&ˆØ�JŠJ�s—z’z !‘}”}Ñ%Ô%ˆØ�JŠJ�s—z’z !‘}”}Ñ%Ô%ˆÝ! ! Q¨wÐ7Ñ7Ô7ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆÝ˜œ b¤lÑ3Ô3Ð3Ý˜œ	 2¤9°4Ð8Ñ8Ô8Ð8Ð8Ð8r=   c                 ó€  — |                      d¦  «        }|                     g d¢¦  «        }t          ||d¬¦  «        }t          ||d¬¦  «        }t          |j        |j        ¦  «         |                      d¦  «        }t          ||d¬¦  «        }t          ||d¬¦  «        }t          |j        |j        ¦  «         d S )NrŸ  )r   gÍÌÌÌÌÌ@g333333*@rÌ   rÆ   rÐ   g      5@rË   )rN   r4   r   r   r^   )r5   r6   r7   r8   r4  r5  s         r;   Útest_method_autozTestCvm_2samp.test_method_auto!  s¶   € Ø�IŠI�c‰NŒNˆØ�JŠJÐ'Ð'Ð'Ñ(Ô(ˆÝ! ! Q¨wÐ7Ñ7Ô7ˆÝ! ! Q¨vÐ6Ñ6Ô6ˆÝ˜œ	 2¤9Ñ-Ô-Ð-à�IŠI�c‰NŒNˆÝ! ! Q¨|Ð<Ñ<Ô<ˆÝ! ! Q¨vÐ6Ñ6Ô6ˆÝ˜œ	 2¤9Ñ-Ô-Ð-Ð-Ð-r=   c                 óÐ  — |                      d¦  «        }t          ||¦  «        }t          |j        |                     d¦  «        ¦  «         t          |j        |                     d¦  «        ¦  «         t          |d d…         |d d…         ¦  «        }t          |j        |                     d¦  «        ¦  «         t          |j        |                     d¦  «        ¦  «         d S )Nrâ  r‹   rŒ   rB   )rN   r   r   r\   r4   r^   )r5   r6   r7   r`   s       r;   Útest_same_inputzTestCvm_2samp.test_same_input-  s·   € ð �IŠI�b‰MŒMˆÝ" 1 aÑ(Ô(ˆÝ˜œ r§z¢z°"¡~¤~Ñ6Ô6Ð6Ý˜œ
 B§J¢J¨r¡N¤NÑ3Ô3Ð3å" 1 R a R¤5¨!¨B¨Q¨B¬%Ñ0Ô0ˆÝ˜œ r§z¢z°"¡~¤~Ñ6Ô6Ð6Ý˜œ
 B§J¢J¨r¡N¤NÑ3Ô3Ð3Ð3Ð3r=   N)ru   rv   rw   r1   rx   ry   rM   rN   r»   r/  r£   r6  r<  r?  ÚxslowrB  rD  rF  rH  rz   r=   r;   r*  r*  Ù  s}  € € € € € à„[×Ò˜V r¨9¨2¬9°Q©<¬<Ð&8Ø'0 r¤y°¡|¤|°a°SÐ&9ð&;ñ <ô <à„[×!Ò! +Ð6RÐ!ÑSÔSð<ð <ñ TÔSñ<ô <ð<ð.ð .ð .ðKð Kð Kð „[×Ò˜WÐ&BÐ&BÐ&BÑCÔCðKð Kñ DÔCðKð „[×ÒÐ4ð5ð 5ð 5ñ6ô 6ð
Cð Cñ6ô 6ð
Cð „[Ôð
 ð 
 ñ Ôð
 ð9ð 9ð 9ð
.ð 
.ð 
.ð4ð 4ð 4ð 4ð 4r=   r*  c                   óð  — e Zd Zg d¢g d¢g d¢fZg d¢g d¢g d¢fZg d¢g d¢g d¢fZdZdZd	Ze	j
                             d
eedfeedfeedffg d¢¬¦  «        d„ ¦   «         ZdZdZe	j
                             d
eedfeedffddg¬¦  «        d„ ¦   «         Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Ze	j
                             dd¦  «        d „ ¦   «         Ze	j
                             d!g d"¢¦  «        d#„ ¦   «         Zd$„ Zd%S )&ÚTestTukeyHSD)rH  ç     €7@çffffff:@çš™™™™;@çfffffæ=@)çffffff<@çš™™™™A@ç     €=@çš™™™™@@çš™™™™>@)gš™™™™:@gÍÌÌÌÌL<@gÍÌÌÌÌL8@g333333:@gÍÌÌÌÌÌ;@)rH  rL  gHáz®G:@rM  rN  rO  rT  rT  )rH  rL  rM  )
rP  rQ  rR  rS  rT  rP  rQ  rR  rS  rT  aK  
    Comparison LowerCL Difference UpperCL Significance
    2 - 3	0.6908830568	4.34	7.989116943	    1
    2 - 1	0.9508830568	4.6 	8.249116943 	1
    3 - 2	-7.989116943	-4.34	-0.6908830568	1
    3 - 1	-3.389116943	0.26	3.909116943	    0
    1 - 2	-8.249116943	-4.6	-0.9508830568	1
    1 - 3	-3.909116943	-0.26	3.389116943	    0
    aS  
    Comparison LowerCL Difference UpperCL Significance
    2 - 1	0.2679292645	3.645	7.022070736	    1
    2 - 3	0.5934764007	4.34	8.086523599	    1
    1 - 2	-7.022070736	-3.645	-0.2679292645	1
    1 - 3	-2.682070736	0.695	4.072070736	    0
    3 - 2	-8.086523599	-4.34	-0.5934764007	1
    3 - 1	-4.072070736	-0.695	2.682070736	    0
    aS  
    Comparison LowerCL Difference UpperCL Significance
    2 - 3	1.561605075	    4.34	7.118394925	    1
    2 - 1	2.740784879	    6.08	9.419215121	    1
    3 - 2	-7.118394925	-4.34	-1.561605075	1
    3 - 1	-1.964526566	1.74	5.444526566	    0
    1 - 2	-9.419215121	-6.08	-2.740784879	1
    1 - 3	-5.444526566	-1.74	1.964526566	    0
    zdata,res_expect_str,atolr/   g»½×Ùß|Û=)úequal size sampleúunequal sample sizezextreme sample size differences©Úidsc                 óf  — t          j        |                     dd¦  «                             ¦   «         dd…         t          ¬¦  «                             d¦  «        }t          j        |Ž }|                     ¦   «         }|D ]¯\  }}}	}
}}t          |¦  «        dz
  t          |¦  «        dz
  }}t          |j        ||f         |	|¬¦  «         t          |j        ||f         |
|¬¦  «         t          |j        ||f         ||¬¦  «         t          |j        ||f         d	k    |dk    ¦  «         Œ°dS )
a£  
        SAS code used to generate results for each sample:
        DATA ACHE;
        INPUT BRAND RELIEF;
        CARDS;
        1 24.5
        ...
        3 27.8
        ;
        ods graphics on;   ODS RTF;ODS LISTING CLOSE;
           PROC ANOVA DATA=ACHE;
           CLASS BRAND;
           MODEL RELIEF=BRAND;
           MEANS BRAND/TUKEY CLDIFF;
           TITLE 'COMPARE RELIEF ACROSS MEDICINES  - ANOVA EXAMPLE';
           ods output  CLDiffs =tc;
        proc print data=tc;
            format LowerCL 17.16 UpperCL 17.16 Difference 17.16;
            title "Output with many digits";
        RUN;
        QUIT;
        ODS RTF close;
        ODS LISTING;
        ú - ú rC   Nr*   )rD   rD   r?   r-   r~   ©rM   r4   ÚreplaceÚsplitÚfloatr·  r   Ú	tukey_hsdÚconfidence_intervalÚintr   Úlowr\   Úhighr^   )r5   ÚdataÚres_expect_strr.   Ú
res_expectÚ	res_tukeyÚconfro   ÚjÚlr  ÚhÚsigs                r;   Útest_compare_saszTestTukeyHSD.test_compare_sasf  s>  € õB ”Z × 6Ò 6°u¸cÑ BÔ B× HÒ HÑ JÔ JÈ1È2È2Ô NÝ&+ð-ñ -ô -ß-4ªW°V©_¬_ð 	å”O TÐ*ˆ	Ø×,Ò,Ñ.Ô.ˆà",ð 	Gð 	GÑˆAˆq�!�Q˜˜3Ý�q‘6”6˜A‘:�s 1™vœv¨™zˆqˆAÝ˜DœH Q¨ TœN¨A°DÐ9Ñ9Ô9Ð9Ý˜IÔ/°°1°Ô5°q¸tÐDÑDÔDÐDÝ˜DœI a¨ dœO¨Q°TÐ:Ñ:Ô:Ð:Ý˜YÔ-¨a°¨dÔ3°sÒ:¸SÀAºXÑFÔFÐFÐFð	Gð 	Gr=   zÐ
        1	2	-8.2491590248597	-4.6	-0.9508409751403	0.0144483269098
        1	3	-3.9091590248597	-0.26	3.3891590248597	0.9803107240900
        2	3	0.6908409751403	4.34	7.9891590248597	0.0203311368795
        zÛ
        1	2	-7.02207069748501	-3.645	-0.26792930251500 0.03371498443080
        1	3	-2.68207069748500	0.695	4.07207069748500 0.85572267328807
        2	3	0.59347644287720	4.34	8.08652355712281 0.02259047020620
        rJ  r  rU  zunequal size samplec                 ó"  — t          j        |                     ¦   «         t          ¬¦  «                             d¦  «        }t          j        |Ž }|                     ¦   «         }|D ]©\  }}}	}
}}t          |¦  «        dz
  t          |¦  «        dz
  }}t          |j
        ||f         |	|¬¦  «         t          |j        ||f         |
|¬¦  «         t          |j        ||f         ||¬¦  «         t          |j        ||f         ||¬¦  «         ŒªdS )an  
        vals = [24.5, 23.5,  26.4, 27.1, 29.9, 28.4, 34.2, 29.5, 32.2, 30.1,
         26.1, 28.3, 24.3, 26.2, 27.8]
        names = {'zero', 'zero', 'zero', 'zero', 'zero', 'one', 'one', 'one',
         'one', 'one', 'two', 'two', 'two', 'two', 'two'}
        [p,t,stats] = anova1(vals,names,"off");
        [c,m,h,nms] = multcompare(stats, "CType","hsd");
        r*   ©rA   rD   r?   r-   N)rM   r4   r^  r_  r·  r   r`  ra  rb  r   rc  r\   rd  r^   )r5   re  rf  r.   rg  rh  ri  ro   rj  rk  r  rl  r:   s                r;   Útest_compare_matlabz TestTukeyHSD.test_compare_matlabŸ  s   € õ ”Z × 4Ò 4Ñ 6Ô 6Ý&+ð-ñ -ô -ß-4ªW°V©_¬_ð 	å”O TÐ*ˆ	Ø×,Ò,Ñ.Ô.ˆà *ð 	Bð 	BÑˆAˆq�!�Q˜˜1Ý�q‘6”6˜A‘:�s 1™vœv¨™zˆqˆAÝ˜DœH Q¨ TœN¨A°DÐ9Ñ9Ô9Ð9Ý˜IÔ/°°1°Ô5°q¸tÐDÑDÔDÐDÝ˜DœI a¨ dœO¨Q°TÐ:Ñ:Ô:Ð:Ý˜IÔ,¨Q°¨TÔ2°A¸DÐAÑAÔAÐAÐAð	Bð 	Br=   c                 ót  — d}t          j        |                     dd¦  «                             ¦   «         dd…         t          ¬¦  «                             d¦  «        }g d¢g d	¢g d
¢f}t          j        |Ž }|                     ¦   «         }|D ]©\  }}}}	}
}t          |¦  «        dz
  t          |¦  «        dz
  }}t          |j        ||f         |	d¬¦  «         t          |j        ||f         |d¬¦  «         t          |j        ||f         |
d¬¦  «         t          |j        ||f         |d¬¦  «         ŒªdS )a+  
        Testing against results and p-values from R:
        from: https://www.rdocumentation.org/packages/stats/versions/3.6.2/
        topics/TukeyHSD
        > require(graphics)
        > summary(fm1 <- aov(breaks ~ tension, data = warpbreaks))
        > TukeyHSD(fm1, "tension", ordered = TRUE)
        > plot(TukeyHSD(fm1, "tension"))
        Tukey multiple comparisons of means
        95% family-wise confidence level
        factor levels have been ordered
        Fit: aov(formula = breaks ~ tension, data = warpbreaks)
        $tension
        zÞ
                diff        lwr      upr     p adj
        2 - 3  4.722222 -4.8376022 14.28205 0.4630831
        1 - 3 14.722222  5.1623978 24.28205 0.0014315
        1 - 2 10.000000  0.4401756 19.55982 0.0384598
        rZ  r[  rC   Nr*   rp  )é   rK   é6   r”  éF   é4   é3   rs  éC   r�   r•  é   é   ry  é   é)   rã  é,   )r–  r  ry  r˜  rf   r–  r—  rK   é$   é*   rs  rz  rä  rà  rL   r  rà  ry  )r~  r  r  r–  rE   rß  rL   râ  rs  rã  r  r  r˜  rª   râ  râ  rä  rL   r?   r  r-   rV  çñhãˆµøä>r\  )r5   Ústr_resrg  re  rh  ri  ro   rj  r  rk  rl  r:   s               r;   Útest_compare_rzTestTukeyHSD.test_compare_r¹  st  € ðˆõ ”Z §¢°°sÑ ;Ô ;× AÒ AÑ CÔ CÀAÀBÀBÔ GÝ&+ð-ñ -ô -ß-4ªW°V©_¬_ð 	ð5ð 5ð 5ð5ð 5ð 5ð5ð 5ð 5ð	6ˆõ ”O TÐ*ˆ	Ø×,Ò,Ñ.Ô.ˆà *ð 	Bð 	BÑˆAˆq�!�Q˜˜1Ý�q‘6”6˜A‘:�s 1™vœv¨™zˆqˆAå˜DœH Q¨ TœN¨A°DÐ9Ñ9Ô9Ð9Ý˜IÔ/°°1°Ô5°q¸tÐDÑDÔDÐDÝ˜DœI a¨ dœO¨Q°TÐ:Ñ:Ô:Ð:Ý˜IÔ,¨Q°¨TÔ2°A¸DÐAÑAÔAÐAÐAð	Bð 	Br=   c                 ó¨  — g d¢}g d¢}g d¢}g d¢}t          j        ||||¦  «        }|                     ¦   «         }t          j        g d¢g d¢g d¢g d¢g¦  «        }t          j        g d	¢g d
¢g d¢g d¢g¦  «        }dD ]S\  }	}
t          |j        |	|
f         ||	|
f         d¬¦  «         t          |j        |	|
f         ||	|
f         d¬¦  «         ŒTdS )zp
        Example sourced from:
        https://www.itl.nist.gov/div898/handbook/prc/section4/prc471.htm
        )çš™™™™™@gš™™™™™@ç333333@gffffff@rå   )gš™™™™™ @r:  g333333@gffffff"@r9  )g       @g      %@g333333 @r„  r8  )r…  gffffff@gffffff@gffffff@gÍÌÌÌÌÌ@)r   r   r   g      À)g�Âõ(\�Ò?r   gq=
×£pÀg¤p=
×£À?)g®Gázò?r   r   g
×£p=
ï?)r   r   r   r   )r   r   r   gáz®Gáþ?)gáz®Gá@r   g      ô?g=
×£p=@)g=
×£p=@r   r   gš™™™™™@)©r?   r   )r@   r   )r   rA   r*  r¥  r}   r-   N)r   r`  ra  rM   r4   r   rc  rd  )r5   Úgroup1Úgroup2Úgroup3Úgroup4r`   ri  ÚlowerÚupperro   rj  s              r;   Útest_engineering_stat_handbookz+TestTukeyHSD.test_engineering_stat_handbookâ  s:  € ð
 +Ð*Ð*ˆØ*Ð*Ð*ˆØ+Ð+Ð+ˆØ*Ð*Ð*ˆÝŒo˜f f¨f°fÑ=Ô=ˆØ×&Ò&Ñ(Ô(ˆÝ”
ØÐÐØ Ð Ð ØÐÐØˆLˆLð	ñ ô ˆõ
 ”
ØˆOˆOØ!Ð!Ð!ØÐÐØˆLˆLð	ñ ô ˆð ?ð 	Eð 	E‰FˆQ�Ý˜DœH Q¨ TœN¨E°!°Q°$¬K¸dÐCÑCÔCÐCÝ˜DœI a¨ dœO¨U°1°a°4¬[¸tÐDÑDÔDÐDÐDð	Eð 	Er=   c                 óÔ  — t           j                             d¦  «        }|                     d¦  «        }t          j        |Ž }|                     ¦   «         }t          |j        |j        j	         ¦  «         t          t          j
        |j        ¦  «        |j        d         ¦  «         t          t          j
        |j        ¦  «        |j        d         ¦  «         t          |j        |j        j	         ¦  «         t          t          j
        |j        ¦  «        d¦  «         t          |j        |j        j	        ¦  «         t          t          j
        |j        ¦  «        d¦  «         d S )Nl   ã]ÏA )rA   r�   ©r   r   r   r?   )rM   rk   rl   r   r`  ra  r   rc  rd  r«  Údiagonalr\   r^   )r5   rn   re  r`   ri  s        r;   Útest_rand_symmzTestTukeyHSD.test_rand_symmü  s  € åŒi×#Ò# JÑ/Ô/ˆØ�zŠz˜(Ñ#Ô#ˆÝŒo˜tÐ$ˆØ×&Ò&Ñ(Ô(ˆå�T”X ¤	¤˜|Ñ,Ô,Ð,õ 	•R”[ ¤Ñ+Ô+¨T¬Y°t¬_Ñ=Ô=Ð=Ý•R”[ ¤Ñ*Ô*¨D¬H°T¬NÑ;Ô;Ð;å�S”] S¤]¤_Ð$4Ñ5Ô5Ð5Ý•R”[ ¤Ñ/Ô/°Ñ3Ô3Ð3å�S”Z ¤¤Ñ.Ô.Ð.Ý•R”[ ¤Ñ,Ô,¨aÑ0Ô0Ð0Ð0Ð0r=   c                 ó°   — t          t          d¬¦  «        5  t          j        g d¢dt          j        gg d¢¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nz...must be finite.r¡   r²   r@   )rD   rF   rA   )rÈ   r¶   r   r`  rM   rm   r-  s    r;   Útest_no_infzTestTukeyHSD.test_no_inf  sœ   € Ý�:Ð-AÐBÑBÔBð 	?ð 	?ÝŒO˜I˜I˜I¨­2¬6 {°I°I°IÑ>Ô>Ð>ð	?ð 	?ð 	?ñ 	?ô 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?øøøð 	?ð 	?ð 	?ð 	?ð 	?ð 	?s   —'AÁAÁAc                 ó¤   — t          t          d¬¦  «        5  t          j        ddgddggddgg d¢¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nz...must be one-dimensionalr¡   r?   r@   rA   rC   )rC   é   rD   ©rÈ   r¶   r   r`  r-  s    r;   Ú
test_is_1dzTestTukeyHSD.test_is_1d  s·   € Ý�:Ð-IÐJÑJÔJð 	Bð 	BÝŒO˜a ˜V a¨ VÐ,¨q°!¨f°j°j°jÑAÔAÐAð	Bð 	Bð 	Bñ 	Bô 	Bð 	Bð 	Bð 	Bð 	Bð 	Bð 	Bð 	Bøøøð 	Bð 	Bð 	Bð 	Bð 	Bð 	Bs   —!AÁA	ÁA	c                 ó˜   — t          t          d¬¦  «        5  t          j        g ddgg d¢¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nz...must be greater than oner¡   r@   rC   )rB   rC   rD   r–  r-  s    r;   Útest_no_emptyzTestTukeyHSD.test_no_empty  s–   € Ý�:Ð-JÐKÑKÔKð 	3ð 	3ÝŒO˜B  A ¨	¨	¨	Ñ2Ô2Ð2ð	3ð 	3ð 	3ñ 	3ô 	3ð 	3ð 	3ð 	3ð 	3ð 	3ð 	3ð 	3øøøð 	3ð 	3ð 	3ð 	3ð 	3ð 	3s   —?¿AÁAc                 ó¤   — d}t          t          |¬¦  «        5  t          j        g d¢ddgddgd¬	¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )
Nz(Expected a boolean value for 'equal_var'r¡   r²   r@   rC   rD   rF   ÚFalse)Ú	equal_var)rÈ   Ú	TypeErrorr   r`  )r5   r2  s     r;   Útest_equal_var_input_validationz,TestTukeyHSD.test_equal_var_input_validation  s·   € Ø8ˆÝ�9¨CÐ0Ñ0Ô0ð 	Jð 	JÝŒO˜I˜I˜I¨¨1 v°°1¨vÀÐIÑIÔIÐIð	Jð 	Jð 	Jñ 	Jô 	Jð 	Jð 	Jð 	Jð 	Jð 	Jð 	Jð 	Jøøøð 	Jð 	Jð 	Jð 	Jð 	Jð 	Js   ™AÁA	ÁA	Únargsr„  c                 óŒ   — t          t          d¬¦  «        5  t          j        g d¢g|z  Ž  d d d ¦  «         d S # 1 swxY w Y   d S )Nz...more than 1 treatment.r¡   ©r•  rF   rA   r–  )r5   rŸ  s     r;   Útest_not_enough_treatmentsz'TestTukeyHSD.test_not_enough_treatments   s’   € å�:Ð-HÐIÑIÔIð 	5ð 	5ÝŒO˜z˜z˜z˜l¨UÑ2Ð4Ð4ð	5ð 	5ð 	5ñ 	5ô 	5ð 	5ð 	5ð 	5ð 	5ð 	5ð 	5ð 	5øøøð 	5ð 	5ð 	5ð 	5ð 	5ð 	5s   —9¹=Á =Úcl)rž  r   r?   r@   c                 óÆ   — t          t          d¬¦  «        5  t          j        g d¢ddgddg¦  «        }|                     |¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nzmust be between 0 and 1r¡   r¡  rA   rB   rg   )rÈ   r¶   r   r`  ra  )r5   r£  r;  s      r;   Útest_conf_level_invalidz$TestTukeyHSD.test_conf_level_invalid%  s­   € å�:Ð-FÐGÑGÔGð 	&ð 	&Ý” 
 
 
¨Q°¨F°Q¸°FÑ;Ô;ˆAØ×!Ò! "Ñ%Ô%Ð%ð	&ð 	&ð 	&ñ 	&ô 	&ð 	&ð 	&ð 	&ð 	&ð 	&ð 	&ð 	&øøøð 	&ð 	&ð 	&ð 	&ð 	&ð 	&s   —2AÁAÁAc                 óò   — t          j        | j        d d…         Ž }t          j        | j        d d…         Ž }t	          |j        |j        d         ¦  «         t	          |j        |j        d         ¦  «         d S )Nr@   r„  r†  )r   r`  Údata_diff_sizeÚ	ttest_indr   r^   )r5   rh  Ú	res_ttests      r;   Útest_2_args_ttestzTestTukeyHSD.test_2_args_ttest+  sn   € å”O TÔ%8¸¸!¸Ô%<Ð=ˆ	Ý”O TÔ%8¸¸!¸Ô%<Ð=ˆ	Ý˜	Ô(¨)Ô*:¸4Ô*@ÑAÔAÐAÝ˜	Ô(¨)Ô*:¸4Ô*@ÑAÔAÐAÐAÐAr=   N)ru   rv   rw   Údata_same_sizer§  Úextreme_sizeÚsas_same_sizeÚsas_diff_sizeÚsas_extremer1   rx   ry   rn  Úmatlab_sm_sizÚmatlab_diff_szrq  r‚  r�  r‘  r“  r—  r™  rž  r¢  r¥  rª  rz   r=   r;   rK  rK  ;  sg  € € € € € à4Ð4Ð4Ø4Ð4Ð4Ø4Ð4Ð4ð6€Nð HÐGÐGØ4Ð4Ð4Ø4Ð4Ð4ð6€Nð 'Ð&Ð&ðð ð à2Ð2Ð2ð4€Lð
€Mð€Mð€Kð „[×ÒÐ7Ø-¨}¸dÐCØ-¨}¸dÐCØ+¨[¸%Ð@ð ð"Eð "Eð "Eð ñ Fô Fð#Gð #GñFô Fð#GðJ€Mð€Nð „[×ÒÐ7Ø-¨}¸eÐDØ-¨~¸tÐDðFà"5Ø"7ð"9ð ñ :ô :ð
Bð Bñ:ô :ð
Bð*'Bð 'Bð 'BðREð Eð Eð41ð 1ð 1ð&?ð ?ð ?ðBð Bð Bð3ð 3ð 3ðJð Jð Jð
 „[×Ò˜W fÑ-Ô-ð5ð 5ñ .Ô-ð5ð „[×Ò˜T > > >Ñ2Ô2ð&ð &ñ 3Ô2ð&ð
Bð Bð Bð Bð Br=   rK  c                   ó,  — e Zd Zg d¢g d¢g d¢fZg d¢g d¢g d¢g d¢fZdZd	Zej        j	        ej         
                    d
eefeeffddg¬¦  «        d„ ¦   «         ¦   «         ZdZdZej         
                    d
eefeeffddg¬¦  «        d„ ¦   «         ZdS )ÚTestGamesHowell)g      8@rÞ   g      ?@ç     €I@)g      A@ç      2@rµ  g      :@)g      1@g      Q@g     €M@g      @)g      >@rÞ   r´  )g     @TÀg     ÀQ@g      ;Àg     €O@)g      E@g      &@g      =@g      3@g      I@)rÞ   g      6@rŸ  rµ  g      "@aK  
            Mean Diff      Lower Bound         Upper Bound         Sig
    0 - 1   8.25000000    -16.5492749527311    33.0492749527311    0.558733632413559
    0 - 2  -5.50000000    -63.6702454316458    52.6702454316458    0.941147750599221
    1 - 2  -13.7500000    -74.3174374251372    46.8174374251372    0.682983914946841
    a+  
             Mean Diff       Lower Bound        Upper Bound         Sig
    0 - 1	 28.16666667    -141.985416377670   198.318749711003	0.8727542747886180
    0 - 2	 4.466666667	-37.2830676783904   46.2164010117237	0.9752628408671710
    0 - 3	 16.26666667	-35.0933112382470   67.6266445715803	0.4262506151302880
    1 - 2	-23.70000000	-195.315617201249   147.915617201249	0.9148950609000590
    1 - 3	-11.90000000	-188.105478728519   164.305478728519	0.9861432250093960
    2 - 3	 11.80000000	-16.2894857524254	39.8894857524254    0.4755344436335670
    zdata, res_expect_strrU  rV  rW  c                 óV  — t          j        |                     dd¦  «                             ¦   «         dd…         t          ¬¦  «                             dd¦  «        }t          j        |dd	iŽ}|                     ¦   «         }|D ]£\  }}}}	}
}t          |¦  «        t          |¦  «        }}t          |j        ||f         |d
¬¦  «         t          |j        ||f         |d
¬¦  «         t          |j        ||f         |	d¬¦  «         t          |j        ||f         |
d¬¦  «         Œ¤dS )a^  
        DATA LIST LIST /Group (F1.0) Value (F8.2).
        BEGIN DATA
        0 24
        0 23
        0 31
        0 51
        1 34
        1 18
        1 18
        1 26
        2 17
        2 68
        2 59
        2 7
        END DATA.

        ONEWAY Value BY Group
            /MISSING ANALYSIS
            /POSTHOC=GH ALPHA(0.05).
        rZ  r[  rF   Nr*   rh   rD   rœ  Fg:Œ0âŽyE>r-   rV  r€  )rM   r4   r]  r^  r_  r·  r   r`  ra  rb  r   r\   r^   rc  rd  )r5   re  rf  rg  Ú	res_gamesri  ro   rj  r  rk  rl  r:   s               r;   Útest_compare_spssz!TestGamesHowell.test_compare_spssO  s:  € õ8 ”ZØ×"Ò" 5¨#Ñ.Ô.×4Ò4Ñ6Ô6°q°r°rÔ:Ýðñ ô ç š  Q™œð 	õ ”O TÐ;°UÐ;Ð;ˆ	Ø×,Ò,Ñ.Ô.ˆà *ð 	;ð 	;ÑˆAˆq�!�Q˜˜1Ý�q‘6”6�3˜q™6œ6ˆqˆAÝ˜IÔ/°°1°Ô5°q¸tÐDÑDÔDÐDÝ˜IÔ,¨Q°¨TÔ2°A¸DÐAÑAÔAÐAÝ˜DœH Q¨ TœN¨A°DÐ9Ñ9Ô9Ð9Ý˜DœI a¨ dœO¨Q°TÐ:Ñ:Ô:Ð:Ð:ð	;ð 	;r=   zÙ
                  q value             Pr(>|q|)
    1 - 0 == 0   -1.5467805948856344  0.55873362851759
    2 - 0 == 0    0.4726721776628535  0.94114775035993
    2 - 1 == 0    1.246837541297872   0.68298393799782
    aw  
                 q value             Pr(>|q|)
    1 - 0 == 0  -1.0589317485313876  0.87275427357438
    2 - 0 == 0  -0.5716222106144833  0.97526284087419
    3 - 0 == 0  -2.6209678382077000  0.42625067714691
    2 - 1 == 0   0.8971899898179028  0.91489506061850
    3 - 1 == 0   0.4579447210555352  0.98614322544695
    3 - 2 == 0  -2.198800177874794   0.47553444364614
    c                 óš  — t          j        |                     dd¦  «                             dd¦  «                             ¦   «         dd…         t          ¬¦  «                             dd¦  «        }t          j        |d	d
iŽ}|D ]E\  }}}}}t          |¦  «        t          |¦  «        }}t          |j
        ||f         |d¬¦  «         ŒFdS )a;  
        games-howell is provided by PMCMRplus package
        https://search.r-project.org/CRAN/refmans/PMCMRplus/html/gamesHowellTest.html
        > library("PMCMRplus")
        > options(digits=16)
        > table = data.frame(
            values = c(24., 23., 31., 51., 34., 18., 18., 26., 17., 68., 59.,  7.),
            groups = c("0", "0", "0", "0", "1", "1", "1", "1", "2", "2", "2", "2")
          )
        > table$groups = as.factor(table$groups)
        > fit <-aov(values ~ groups, table)
        > res = gamesHowellTest(fit)
        > summary(res)
        rZ  r[  z == rA   Nr*   rh   rC   rœ  Fr  r-   )rM   r4   r]  r^  r_  r·  r   r`  rb  r   r^   )	r5   re  rf  rg  r·  rj  ro   Ú_r:   s	            r;   r‚  zTestGamesHowell.test_compare_r‰  sÜ   € õ( ”ZØ×"Ò" 5¨#Ñ.Ô.ßŠW�V˜SÑ!Ô!§%¢%¡'¤'¨!¨"¨"ô.åðñ ô ÷ !š  Q™œð 	õ ”O TÐ;°UÐ;Ð;ˆ	ð (ð 	Bð 	B‰MˆAˆq�!�Q˜Ý�q‘6”6�3˜q™6œ6ˆqˆAÝ˜IÔ,¨Q°¨TÔ2°A¸DÐAÑAÔAÐAÐAð	Bð 	Br=   N)ru   rv   rw   r«  r§  Úspss_same_sizeÚspss_diff_sizer1   rx   rI  ry   r¸  Úr_same_sizeÚr_diff_sizer‚  rz   r=   r;   r³  r³  3  sC  € € € € € à*Ð*Ð*Ø*Ð*Ð*Ø*Ð*Ð*ð,€Nð &�o�oØ,Ð,Ð,Ø/Ð/Ð/Ø.Ð.Ð.ð0€Nð
€Nð€Nð „[ÔØ„[×ÒÐ3Ø,¨nÐ=Ø+¨^Ð<ð>à!4Ø!6ð!8ð ñ 9ô 9ð
!;ð !;ñ9ô 9ñ Ôð!;ðF€Kð€Kð „[×ÒÐ3Ø,¨kÐ:Ø+¨[Ð9ð;à!4Ø!6ð!8ð ñ 9ô 9ð
Bð Bñ9ô 9ð
Bð Bð Br=   r³  c                   óì   — e Zd Zej                             dg d¢g d¢f¦  «        d„ ¦   «         Zej                             dg d¢g d¢g d¢g d	¢g d
¢g d¢g d¢g d¢f¦  «        d„ ¦   «         Zd„ Zd„ Z	d„ Z
dS )ÚTestPoissonMeansTestzc1, n1, c2, n2, p_expect)r   r�   rA   r�   gþe÷äa¡¶?)r@   r�   rD   r�   gÞ	ŠcÆ?c                 ób   — t          j        ||||¦  «        }t          |j        |d¬¦  «         d S r†   ©r   Úpoisson_means_testr   r^   )r5   Úc1Ún1Úc2Ún2Úp_expectr`   s          r;   Útest_paper_examplesz(TestPoissonMeansTest.test_paper_examples­  s7   € õ Ô& r¨2¨r°2Ñ6Ô6ˆÝ˜œ
 H°4Ð8Ñ8Ô8Ð8Ð8Ð8r=   z c1, n1, c2, n2, p_expect, alt, d)rã  rE   rã  rE   gŽŸ{}ÿÿï?rÒ   r   )rE   rE   rE   rE   goPFªÿÿï?rÒ   r   )é2   râ  r?   r?   góæ°ae¹?rÒ   r~   )rA   r�   rã  r   gò/V›-=¿?rÒ   r   )rA   rf   rB   rã  g"Ü)‘ÅÝÙ?rØ   r   )rB   rã  rA   r�   g_ù'Qm~€?rØ   r   )rB   rã  rA   rE   g|¨Š»Ó?rÖ   r   )r?   r?   rÊ  râ  gï0�ëÝ·?rÖ   r   c                 ój   — t          j        ||||||¬¦  «        }t          |j        |dd¬¦  «         d S )N)rÅ   Údiffg�íµ ÷ÆÀ>r5  ©r.   rX  rÂ  )	r5   rÄ  rÅ  rÆ  rÇ  rÈ  ÚaltÚdr`   s	            r;   Útest_fortran_authorsz)TestPoissonMeansTest.test_fortran_authors¶  s@   € õ$ Ô& r¨2¨r°2À3ÈQÐOÑOÔOˆÝ˜œ
 H°4¸eÐDÑDÔDÐDÐDÐDr=   c                 ór   — d\  }}d\  }}t          j        ||||¦  «        }t          |j        d¦  «         d S )N)r‘   r‘   r?   rÂ  ©r5   Úcount1Úcount2Únobs1Únobs2r`   s         r;   Útest_different_resultsz+TestPoissonMeansTest.test_different_resultsË  sD   € ð &‰ˆ�Ø#‰ˆˆuÝÔ& v¨u°f¸eÑDÔDˆÝ˜œ
 AÑ&Ô&Ð&Ð&Ð&r=   c                 ór   — d\  }}d\  }}t          j        ||||¦  «        }t          |j        d¦  «         d S )Nr�  r  r?   rÂ  rÒ  s         r;   Útest_less_than_zero_lambda_hat2z4TestPoissonMeansTest.test_less_than_zero_lambda_hat2Ô  sD   € ð ‰ˆ�Ø‰ˆˆuÝÔ& v¨u°f¸eÑDÔDˆÝ˜œ
 AÑ&Ô&Ð&Ð&Ð&r=   c                 óˆ  — d\  }}d\  }}d}t          t          |¬¦  «        5  t          j        d|||¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          |¬¦  «        5  t          j        ||d|¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          j        d|||¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          |¬¦  «        5  t          j        ||d|¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          j        |d||¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          |¬¦  «        5  t          j        |||d¦  «         d d d ¦  «         n# 1 swxY w Y   d	}t          t          |¬¦  «        5  t          j        ||||d¬
¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          t          |¬¦  «        5  t          j        ddddd¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr�  r  z`k1` and `k2` must be integers.r¡   r´   z1`k1` and `k2` must be greater than or equal to 0.rh   z%`n1` and `n2` must be greater than 0.z(diff must be greater than or equal to 0.)rÌ  zAlternative must be one of ...r?   r@   ÚerrorrÄ   )rÈ   r�  r   rÃ  r¶   )r5   rÓ  rÔ  rÕ  rÖ  r¸   s         r;   rÉ   z*TestPoissonMeansTest.test_input_validationÜ  s�  € Ø‰ˆ�Ø‰ˆˆuð 4ˆÝ�9¨GÐ4Ñ4Ô4ð 	?ð 	?ÝÔ$ R¨°¸Ñ>Ô>Ð>ð	?ð 	?ð 	?ñ 	?ô 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?øøøð 	?ð 	?ð 	?ð 	?å�9¨GÐ4Ñ4Ô4ð 	?ð 	?ÝÔ$ V¨U°B¸Ñ>Ô>Ð>ð	?ð 	?ð 	?ñ 	?ô 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?øøøð 	?ð 	?ð 	?ð 	?ð FˆÝ�:¨WÐ5Ñ5Ô5ð 	?ð 	?ÝÔ$ R¨°¸Ñ>Ô>Ð>ð	?ð 	?ð 	?ñ 	?ô 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?øøøð 	?ð 	?ð 	?ð 	?å�:¨WÐ5Ñ5Ô5ð 	?ð 	?ÝÔ$ V¨U°B¸Ñ>Ô>Ð>ð	?ð 	?ð 	?ñ 	?ô 	?ð 	?ð 	?ð 	?ð 	?ð 	?ð 	?øøøð 	?ð 	?ð 	?ð 	?ð :ˆÝ�:¨WÐ5Ñ5Ô5ð 	@ð 	@ÝÔ$ V¨R°¸Ñ?Ô?Ð?ð	@ð 	@ð 	@ñ 	@ô 	@ð 	@ð 	@ð 	@ð 	@ð 	@ð 	@øøøð 	@ð 	@ð 	@ð 	@å�:¨WÐ5Ñ5Ô5ð 	@ð 	@ÝÔ$ V¨U°F¸BÑ?Ô?Ð?ð	@ð 	@ð 	@ñ 	@ô 	@ð 	@ð 	@ð 	@ð 	@ð 	@ð 	@øøøð 	@ð 	@ð 	@ð 	@ð =ˆÝ�:¨WÐ5Ñ5Ô5ð 	Lð 	LÝÔ$ V¨U°F¸EÈÐKÑKÔKÐKð	Lð 	Lð 	Lñ 	Lô 	Lð 	Lð 	Lð 	Lð 	Lð 	Lð 	Løøøð 	Lð 	Lð 	Lð 	Lð 3ˆÝ�:¨WÐ5Ñ5Ô5ð 	Fð 	FÝÔ$ Q¨¨1¨a¸WÐEÑEÔEÐEð	Fð 	Fð 	Fñ 	Fô 	Fð 	Fð 	Fð 	Fð 	Fð 	Fð 	Fð 	Føøøð 	Fð 	Fð 	Fð 	Fð 	Fð 	Fs�   £AÁAÁAÁ(BÂBÂBÂ/CÃCÃCÃ4DÄDÄDÄ;EÅE#Å&E#Æ F$Æ$F(Æ+F(ÇG-Ç-G1Ç4G1ÈH7È7H;È>H;N)ru   rv   rw   r1   rx   ry   rÉ  rÐ  r×  rÙ  rÉ   rz   r=   r;   rÀ  rÀ  ¬  s  € € € € € Ø„[×ÒÐ7à Ð Ð Ø Ð Ð ð:ñ ô ð
9ð 9ñô ð
9ð „[×ÒÐ?ð
 	=Ð<Ð<Ø<Ð<Ð<Ø=Ð=Ð=Ø>Ð>Ð>Ø9Ð9Ð9Ø;Ð;Ð;ð 	6Ð5Ð5Ø6Ð6Ð6ðBñ ô ð"Eð Eñ#ô ð"Eð'ð 'ð 'ð'ð 'ð 'ð!Fð !Fð !Fð !Fð !Fr=   rÀ  c                   ó¾   — e Zd Zd„ Zd„ Zej                             dg d¢¦  «        d„ ¦   «         Zej                             dg d¢¦  «        d„ ¦   «         Z	d„ Z
d	„ Zd
S )ÚTestBWSTestc                 ó‚  — t           j                             d¦  «        }|                     d¬¦  «        \  }}d}t          j        t
          |¬¦  «        5  t          j        ||g||g¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t          j        t           j        g|¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t          j        |g ¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t          j        ||d	¬
¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t          j        ||d¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nì   <oôvêT‘{ ©r@   rF   rd   z,`x` and `y` must be exactly one-dimensional.r¡   z"`x` and `y` must not contain NaNs.z$`x` and `y` must be of nonzero size.zalternative` must be one of...rÌ  rÄ   z!method` must be an instance of...r  rÆ   )	rM   rk   rl   r1   r   r¶   r   Úbws_testr]   )r5   rn   r7   r8   r¸   s        r;   Útest_bws_input_validationz%TestBWSTest.test_bws_input_validation  sá  € ÝŒi×#Ò#Ð$7Ñ8Ô8ˆà�zŠz˜vˆzÑ&Ô&‰ˆˆ1à@ˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	+ð 	+ÝŒN˜A˜q˜6 A q 6Ñ*Ô*Ð*ð	+ð 	+ð 	+ñ 	+ô 	+ð 	+ð 	+ð 	+ð 	+ð 	+ð 	+øøøð 	+ð 	+ð 	+ð 	+ð 7ˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	(ð 	(ÝŒN�BœF˜8 QÑ'Ô'Ð'ð	(ð 	(ð 	(ñ 	(ô 	(ð 	(ð 	(ð 	(ð 	(ð 	(ð 	(øøøð 	(ð 	(ð 	(ð 	(ð 9ˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	"ð 	"ÝŒN˜1˜bÑ!Ô!Ð!ð	"ð 	"ð 	"ñ 	"ô 	"ð 	"ð 	"ð 	"ð 	"ð 	"ð 	"øøøð 	"ð 	"ð 	"ð 	"ð 3ˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	:ð 	:ÝŒN˜1˜a¨[Ð9Ñ9Ô9Ð9ð	:ð 	:ð 	:ñ 	:ô 	:ð 	:ð 	:ð 	:ð 	:ð 	:ð 	:øøøð 	:ð 	:ð 	:ð 	:ð 6ˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	,ð 	,ÝŒN˜1˜a¨Ð+Ñ+Ô+Ð+ð	,ð 	,ð 	,ñ 	,ô 	,ð 	,ð 	,ð 	,ð 	,ð 	,ð 	,ð 	,øøøð 	,ð 	,ð 	,ð 	,ð 	,ð 	,sZ   ÁA<Á<B ÂB Â$!CÃCÃCÃ9DÄDÄ"DÅE'Å'E+Å.E+ÆF4Æ4F8Æ;F8c                 óœ   — g d¢}g d¢}t          j        ||d¬¦  «        }t          |j        dd¬¦  «         t	          |j        d¦  «         d S )	N)r?   r@   rA   rB   rD   rF   rG   )rC   rg   rE   rc   rf   rª   r•  rÒ   rÄ   gºI+‡@rH   r-   g¤f$/•Þg?)r   rá  r   r\   r   r^   rÉ  s       r;   Ú test_against_published_referencez,TestBWSTest.test_against_published_reference  s_   € ð "Ð!Ð!ˆØ&Ð&Ð&ˆÝŒn˜Q ¨{Ð;Ñ;Ô;ˆÝ˜œ u°4Ð8Ñ8Ô8Ð8Ý�S”Z Ñ)Ô)Ð)Ð)Ð)r=   )rÅ   r\   r^   ))rÒ   g
»ù-ü?g4ˆñ³B/À?)rÖ   ç
»ù-ü¿g«î0&Àv­?)rØ   rå  gÉñœý“(î?c                 ó  — t           j                             d¦  «        }|                     d¬¦  «        \  }}t          j        |||¬¦  «        }t          |j        |d¬¦  «         t          |j        |dd¬	¦  «         d S )
Nrß  rà  rd   rÄ   ç‚vIhÂ%<=rW  r}   r¥   rÍ  ©rM   rk   rl   r   rá  r   r\   r^   ©r5   rÅ   r\   r^   rn   r7   r8   r`   s           r;   Útest_against_RzTestBWSTest.test_against_R&  s}   € õ Œi×#Ò#Ð$7Ñ8Ô8ˆØ�zŠz˜vˆzÑ&Ô&‰ˆˆ1ÝŒn˜Q ¨{Ð;Ñ;Ô;ˆÝ˜œ y°uÐ=Ñ=Ô=Ð=Ý˜œ
 F°¸DÐAÑAÔAÐAÐAÐAr=   ))rÒ   góD 5Hò?g«ÒdüÔ•Ò?)rÖ   ç`�Ñ‡©áï?gÈ²µ×šXë?)rØ   rë  gß4)¡”�Â?c                 ó(  — t           j                             d¦  «        }|                     d¬¦  «        }|                     d¬¦  «        }t          j        |||¬¦  «        }t          |j        |d¬¦  «         t          |j        |dd	¬
¦  «         d S )Nl   ‚.ÅÇs½Z rg   rd   rG   rÄ   rç  rW  r}   r¥   rÍ  rè  ré  s           r;   Útest_against_R_imbalancedz%TestBWSTest.test_against_R_imbalanced7  sŒ   € õ Œi×#Ò#Ð$7Ñ8Ô8ˆØ�JŠJ˜AˆJÑÔˆØ�JŠJ˜AˆJÑÔˆÝŒn˜Q ¨{Ð;Ñ;Ô;ˆÝ˜œ y°uÐ=Ñ=Ô=Ð=Ý˜œ
 F°¸DÐAÑAÔAÐAÐAÐAr=   c                 óè  — t           j                             d¦  «        }|                     d¬¦  «        \  }}t           j                             d¦  «        }t          j        d|¬¦  «        }t          j        |||¬¦  «        }t          |j        ¦  «        dk    sJ ‚t           j                             d¦  «        }t          j        d|¬¦  «        }t          j        |||¬¦  «        }t          |j        |j        ¦  «         t           j                             d¦  «        }t          j        d|¬¦  «        }t          j        |||¬¦  «        }t          j	        |j        |j        ¦  «        rJ ‚d S )Nì   /°HøNÏ( )r@   rE   rd   rE   )Ún_resamplesrn   rÆ   l   ÝVC£	ãA )
rM   rk   rl   r   rb  rá  r  Únull_distributionr   Úallclose)r5   rn   r7   r8   rÇ   r  r  rE  s           r;   Útest_methodzTestBWSTest.test_methodI  sH  € åŒi×#Ò#Ð$7Ñ8Ô8ˆØ�zŠz˜wˆzÑ'Ô'‰ˆˆ1åŒi×#Ò#Ð$7Ñ8Ô8ˆÝÔ(°R¸SÐAÑAÔAˆÝŒ~˜a ¨6Ð2Ñ2Ô2ˆå�4Ô)Ñ*Ô*¨bÒ0Ð0Ð0Ð0åŒi×#Ò#Ð$7Ñ8Ô8ˆÝÔ(°R¸SÐAÑAÔAˆÝŒ~˜a ¨6Ð2Ñ2Ô2ˆå˜Ô.°Ô0FÑGÔGÐGåŒi×#Ò#Ð$7Ñ8Ô8ˆÝÔ(°R¸SÐAÑAÔAˆÝŒ~˜a ¨6Ð2Ñ2Ô2ˆå”;˜tÔ5°tÔ7MÑNÔNÐNÐNÐNÐNÐNr=   c                 ó–  — t           j                             d¦  «        }|                     d¬¦  «        }|dz
  }t          j        ||d¬¦  «        }|j        dk    sJ ‚t          |j        dt          |j	        ¦  «        z  ¦  «         t          j        ||d¬¦  «        }|j        dk    sJ ‚t          |j        d¦  «         t          j        ||d¬¦  «        }|j        dk     sJ ‚t          |j        dt          |j	        ¦  «        z  ¦  «         t          j        ||d¬¦  «        }|j        dk     sJ ‚t          |j        d¦  «         d S )	Nrï  rC   rd   r?   rØ   rÄ   r   rÖ   )
rM   rk   rl   r   rá  r\   r   r^   r  rñ  )r5   rn   r7   r8   r`   s        r;   Útest_directionszTestBWSTest.test_directions`  sF  € åŒi×#Ò#Ð$7Ñ8Ô8ˆØ�JŠJ˜AˆJÑÔˆØ�‰EˆåŒn˜Q ¨yÐ9Ñ9Ô9ˆØŒ}˜qÒ Ð Ð Ð Ý�S”Z ¥S¨Ô)>Ñ%?Ô%?Ñ!?Ñ@Ô@Ð@åŒn˜Q ¨vÐ6Ñ6Ô6ˆØŒ}˜qÒ Ð Ð Ð Ý�S”Z Ñ#Ô#Ð#åŒn˜Q ¨vÐ6Ñ6Ô6ˆØŒ}˜qÒ Ð Ð Ð Ý�S”Z ¥S¨Ô)>Ñ%?Ô%?Ñ!?Ñ@Ô@Ð@åŒn˜Q ¨yÐ9Ñ9Ô9ˆØŒ}˜qÒ Ð Ð Ð Ý�S”Z Ñ#Ô#Ð#Ð#Ð#r=   N)ru   rv   rw   râ  rä  r1   rx   ry   rê  rí  ró  rõ  rz   r=   r;   rÝ  rÝ     së   € € € € € ð,ð ,ð ,ð4*ð *ð *ð „[×ÒÐCðNð Nð NñOô OðBð Bñ	Oô OðBð „[×ÒÐCðNð Nð NñOô OðBð Bñ	Oô OðBðOð Oð Oð.$ð $ð $ð $ð $r=   rÝ  )8Ú	itertoolsr   ÚnumpyrM   r‹  r1   Únumpy.testingr   r   r   r   rÈ   Úscipyr   r	   Úscipy.statsr
   Úscipy.stats._hypotestsr   r   r   r   r   r   r   Úscipy.stats._mannwhitneyur   r   r   Úscipy._lib._testutilsr   Ú
scipy._libr   rÍ   Úscipy._lib._array_apir   r   r   r   Úscipy._lib._array_api_no_0dr   r   Úscipy.stats._axis_nan_policyr   r   r    r|   r½   r‚  rÝ  r	  r*  rK  r³  rÀ  rÝ  rz   r=   r;   ú<module>r     sN  ðØ Ð Ð Ð Ð Ð à Ð Ð Ð Ø Ð Ð Ð Ø €€€Ø @Ð @Ð @Ð @Ð @Ð @Ð @Ð @Ð @Ð @Ø *Ð *Ð *Ð *Ð *Ð *à  Ð  Ð  Ð  Ð  Ð  Ð  Ð  Ø %Ð %Ð %Ð %Ð %Ð %ð4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð 4ð EÐ DÐ DÐ DÐ DÐ DÐ DÐ DÐ DÐ DØ 2Ð 2Ð 2Ð 2Ð 2Ð 2Ø -Ð -Ð -Ð -Ð -Ð -ð0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0à HÐ HÐ HÐ HÐ HÐ HÐ HÐ HØ RÐ RÐ RÐ RÐ RÐ RÐ RÐ Rð ÐÐ'Ñ(Ô(ð?2ð ?2ð ?2ð ?2ð ?2ñ ?2ô ?2ñ )Ô(ð?2ðD Ð�5Ô'Ñ(Ô(ð\Oð \Oð \Oð \Oð \Oñ \Oô \Oñ )Ô(ð\Oð~ Ð�5Ô%Ñ&Ô&ð	Cð 	Cð 	Cð 	Cð 	Cñ 	Cô 	Cñ 'Ô&ð	CðDl<ð l<ð l<ð l<ð l<Ð"ñ l<ô l<ð l<ð^P
ð P
ð P
ð P
ð P
ñ P
ô P
ð P
ðfV1ð V1ð V1ð V1ð V1ñ V1ô V1ð V1ðr ÐÐ'Ñ(Ô(ð^4ð ^4ð ^4ð ^4ð ^4ñ ^4ô ^4ñ )Ô(ð^4ðBuBð uBð uBð uBð uBñ uBô uBð uBðpsBð sBð sBð sBð sBñ sBô sBð sBðrQFð QFð QFð QFð QFñ QFô QFð QFðht$ð t$ð t$ð t$ð t$ñ t$ô t$ð t$ð t$ð t$r=   