§
    µ³Ùgj8 ã                   óÔ  — d dl mZ d dlZd dlZd dlZd dlZd dlmZm	Z	m
Z
mZ d dlmZ d dl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#  G d„ d¦  «        Z$ G d„ d¦  «        Z% G d„ d¦  «        Z& G d„ de ¦  «        Z' G d„ d¦  «        Z( G d„ d¦  «        Z) G d„ d¦  «        Z* G d„ d¦  «        Z+ G d„ d¦  «        Z, G d„ d¦  «        Z-dS ) é    )ÚproductN)Úassert_Úassert_equalÚassert_allcloseÚassert_almost_equal)Úraises)Údistributions)Úepps_singleton_2sampÚcramervonmisesÚ_cdf_cvmÚcramervonmises_2sampÚ_pval_cvm_2samp_exactÚbarnard_exactÚboschloo_exact)ÚmannwhitneyuÚ
_mwu_stateÚ_MWUé   )Úcheck_named_results)Ú_TestPythranFunc)ÚSmallSampleWarningÚtoo_small_1d_not_omitc                   ó2   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ ZdS )ÚTestEppsSingletonc                 óÌ   — t          j        g d¢¦  «        }t          j        g d¢¦  «        }t          ||¦  «        \  }}t          |dd¬¦  «         t          |dd¬¦  «         d S )N)
gffffffÖ¿gffffff@g®Gáz®û?ç\�Âõ(\ç?çffffffÖ?g…ëQ¸…@gq=
×£pÝ?g®Gázî¿g®Gáz®×¿g¤p=
×#(@)
gffffffò¿g333333Ã¿g×£p=
×@g      
@g®Gáz®@g)\�Âõ(@g      @gö(\�Âõ@gÃõ(\�Â @ç333333!@gHáz®G.@r   ©Údecimalg˜Q,·´r?é   )ÚnpÚarrayr
   r   ©ÚselfÚxÚyÚwÚps        úV/granary/edev/quantum/lib/python3.11/site-packages/scipy/stats/tests/test_hypotests.pyÚtest_statistic_1z"TestEppsSingleton.test_statistic_1   s‡   € õ
 ŒHð 7ð 7ð 7ñ 8ô 8ˆåŒHð 4ð 4ð 4ñ 5ô 5ˆå# A qÑ)Ô)‰ˆˆ1Ý˜A˜u¨aÐ0Ñ0Ô0Ð0Ý˜A˜w°Ð2Ñ2Ô2Ð2Ð2Ð2ó    c                 óÄ   — t          j        d¦  «        }t          j        d¦  «        }t          ||¦  «        \  }}t          |dd¬¦  «         t	          |dd¬¦  «         d S )	N)r   r   é   r.   r.   r.   r!   r!   r!   r!   é   é   r0   r0   r0   é   é
   r2   r2   r2   )r2   r/   r   r0   r2   r2   r   r0   r1   é   r2   r!   r   r3   r   é   r   r0   r4   r2   gÍÌÌÌÌÌ!@çü©ñÒMbP?©Úatolg&ª·¶J°?r!   r   )r"   r#   r
   r   r   r$   s        r*   Útest_statistic_2z"TestEppsSingleton.test_statistic_2%   sq   € åŒHð "ñ #ô #ˆåŒHð  ñ !ô !ˆå# A qÑ)Ô)‰ˆˆ1Ý˜˜5 uÐ-Ñ-Ô-Ð-Ý˜A˜w°Ð2Ñ2Ô2Ð2Ð2Ð2r,   c                 óâ  — t           j                             d¦  «         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éÒ  é   é   )r"   ÚrandomÚseedÚaranger
   ÚlistÚtupler   )	r%   r&   r'   Úw1Úp1Úw2Úp2Úw3Úp3s	            r*   Útest_epps_singleton_array_likez0TestEppsSingleton.test_epps_singleton_array_like/   sÊ   € Ý
Œ	�Š�tÑÔÐÝŒ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                 ó@  — dt          j        d¦  «        }}t          j        t          t
          ¬¦  «        5  t          ||¦  «        }t          |j        t           j	        ¦  «         t          |j
        t           j	        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )N©r   r.   r!   r/   r2   ©Úmatch)r"   r?   ÚpytestÚwarnsr   r   r
   r   Ú	statisticÚnanÚpvalue©r%   r&   r'   Úress       r*   Útest_epps_singleton_sizez*TestEppsSingleton.test_epps_singleton_size:   sÅ   € à�RœY r™]œ]ˆ1ˆÝŒ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý& q¨!Ñ,Ô,ˆCÝ˜œ­¬Ñ/Ô/Ð/Ý˜œ¥R¤VÑ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð 	-ð 	-s   ·ABÂBÂBc                 óŠ   — dddddt           j        ft          j        d¦  «        }}t          t          t
          ||¦  «         d S )Nr   r.   r!   r/   r0   r2   )r"   Úinfr?   Úassert_raisesÚ
ValueErrorr
   ©r%   r&   r'   s      r*   Útest_epps_singleton_nonfinitez/TestEppsSingleton.test_epps_singleton_nonfiniteB   s=   € à�1�a˜˜A�rœvÐ&­¬	°"©¬ˆ1ˆÝ•jÕ"6¸¸1Ñ=Ô=Ð=Ð=Ð=r,   c                 óš   — t          j        d¦  «        t          j        d¦  «        }}t          ||¦  «        }d}t          ||¦  «         d S )Né   r;   )rO   rQ   )r"   r?   r
   r   )r%   r&   r'   rS   Ú
attributess        r*   Ú
test_nameszTestEppsSingleton.test_namesG   sE   € ÝŒy˜‰}Œ}�bœi¨™mœmˆ1ˆÝ" 1 aÑ(Ô(ˆØ,ˆ
Ý˜C Ñ,Ô,Ð,Ð,Ð,r,   N)	Ú__name__Ú
__module__Ú__qualname__r+   r8   rH   rT   rZ   r^   © r,   r*   r   r      sn   € € € € € ð3ð 3ð 3ð3ð 3ð 3ð	 ð 	 ð 	 ð-ð -ð -ð>ð >ð >ð
-ð -ð -ð -ð -r,   r   c                   ó˜   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Z	d„ Z
ej                             d	d
dgg¦  «        d„ ¦   «         Zd„ Zd„ ZdS )ÚTestCvmc                 óN   — t          t          g d¢d¦  «        g d¢d¬¦  «         d S )N)gy;ÂiÁ‹ž?g#¾³^¥?gþE�>‘¿?gåD»
)î?r/   ©ç{®Gáz„?çš™™™™™©?ç      à?g+‡ÙÎ÷ï?ç-Cëâ6?r6   ©r   r   ©r%   s    r*   Ú
test_cdf_4zTestCvm.test_cdf_4R   sC   € ÝÝÐ=Ð=Ð=¸qÑAÔAØ(Ð(Ð(Øð	ñ 	ô 	ð 	ð 	ð 	r,   c                 óN   — t          t          g d¢d¦  «        g d¢d¬¦  «         d S )N)g…8„*5›?g@¤ß¾œ£?g÷Hmâä¾?g%¯Î1 â?r2   )rg   rh   ri   g333333ï?rj   r6   rk   rl   s    r*   Útest_cdf_10zTestCvm.test_cdf_10X   sC   € ÝÝÐ=Ð=Ð=¸rÑBÔBØ(Ð(Ð(Øð	ñ 	ô 	ð 	ð 	ð 	r,   c                 óN   — t          t          g d¢d¦  «        g d¢d¬¦  «         d S )N)g€}têÊg™?g˜`�º¢?gäIÒ5“o¾?gÜ×�sF”ò?éè  rf   rj   r6   rk   rl   s    r*   Útest_cdf_1000zTestCvm.test_cdf_1000^   sC   € ÝÝÐ=Ð=Ð=¸tÑDÔDØ(Ð(Ð(Øð	ñ 	ô 	ð 	ð 	ð 	r,   c                 óL   — t          t          g d¢¦  «        g d¢d¬¦  «         d S )N)çaÃÓ+e™?ç+û®þ·¢?çÉ&pën¾?ç+MJA·—ò?rf   rj   r6   rk   rl   s    r*   Útest_cdf_infzTestCvm.test_cdf_infd   sA   € ÝÝÐ=Ð=Ð=Ñ>Ô>Ø(Ð(Ð(Øð	ñ 	ô 	ð 	ð 	ð 	r,   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  r   r   g†a†ah?g«ªªªªª"@é   )r   r   rl   s    r*   Útest_cdf_supportzTestCvm.test_cdf_supportj   sN   € å•X˜z¨5Ð1°3Ñ7Ô7¸!¸Q¸Ñ@Ô@Ð@Ý•X˜°
Ð;¸RÑ@Ô@À1ÀaÀ&ÑIÔIÐIÐIÐIr,   c                 óh   — t          t          g d¢d¦  «        t          g d¢¦  «        d¬¦  «         d S )N)rt   ru   rv   rw   éd   é'  rj   r6   rk   rl   s    r*   Útest_cdf_large_nzTestCvm.test_cdf_large_no   sL   € åÝÐBÐBÐBÀEÑJÔJÝÐBÐBÐBÑCÔCØð	ñ 	ô 	ð 	ð 	ð 	r,   c                 ó°   — t          dt          dd¦  «        cxk     odk     nc ¦  «         t          dt          d¦  «        cxk     odk     nc ¦  «         d S )NgwJëÿï?gÍÌÌÌÌÔt@rq   ç      ð?)r   r   rl   s    r*   Útest_large_xzTestCvm.test_large_xv   sn   € õ 	��( 5¨$Ñ/Ô/Ð5Ð5Ò5Ð5°#Ò5Ð5Ð5Ð5Ñ6Ô6Ð6Ý��( 5™/œ/Ð/Ð/Ò/Ð/¨CÒ/Ð/Ð/Ð/Ñ0Ô0Ð0Ð0Ð0r,   c                 óÊ   — d}t          t          j        |¦  «        dz  d¦  «        }t          t	          |j        |¦  «        dk    ¦  «         t          |j        d¦  «         d S )Né   çš™™™™™é?Únormr�   r   )r   r"   Úonesr   r   rO   r   rQ   )r%   ÚnrS   s      r*   Ú
test_low_pzTestCvm.test_low_p   s[   € ð ˆÝ�RœW Q™ZœZ¨™^¨VÑ4Ô4ˆÝ•˜œ¨Ñ*Ô*¨SÒ0Ñ1Ô1Ð1Ý�S”Z Ñ#Ô#Ð#Ð#Ð#r,   r&   rb   ç      ø?c                 ó  — t          j        t          t          ¬¦  «        5  t	          |d¦  «        }t          |j        t          j        ¦  «         t          |j	        t          j        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )NrK   r†   )
rM   rN   r   r   r   r   rO   r"   rP   rQ   ©r%   r&   rS   s      r*   Útest_invalid_inputzTestCvm.test_invalid_input‡   s´   € åŒ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý   FÑ+Ô+ˆCÝ˜œ­¬Ñ/Ô/Ð/Ý˜œ¥R¤VÑ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð 	-ð 	-s   ¡AA=Á=BÂBc                 óˆ  — t          g d¢d¦  «        }t          |j        dd¬¦  «         t          |j        dd¬¦  «         t          g d¢dd¦  «        }t          |j        dd¬¦  «         t          |j        d	d¬¦  «         t          g d
¢d¦  «        }t          |j        dd¬¦  «         t          |j        dd¬¦  «         d S )N)g333333û¿r.   r   gÍÌÌÌÌÌô?r/   çš™™™™™¹?ç333333ã?r†   góZ	Ý%qÒ?ç�íµ ÷Æ°>r6   g®EÐ¶šÂ?)r!   rŠ   g!O!W*î?gz"ÇW´™`?)	r   r.   r0   çffffffö?gìQ¸…ëÁ?é   é   çÍÌÌÌÌÌì?ç      @ÚexpongÐeÅË.óê?gnz\Ã(r?)r   r   rO   rQ   )r%   rS   s     r*   Útest_values_RzTestCvm.test_values_RŽ   sÜ   € õ Ð;Ð;Ð;¸VÑDÔDˆÝ˜œ x°dÐ;Ñ;Ô;Ð;Ý˜œ
 I°DÐ9Ñ9Ô9Ð9õ Ð;Ð;Ð;¸VÀXÑNÔNˆÝ˜œ y°tÐ<Ñ<Ô<Ð<Ý˜œ
 K°dÐ;Ñ;Ô;Ð;õ ÐCÐCÐCÀWÑMÔMˆÝ˜œ y°tÐ<Ñ<Ô<Ð<Ý˜œ
 K°dÐ;Ñ;Ô;Ð;Ð;Ð;r,   c                 ó’  — 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 )Nr0   )r’   çffffffæ?r—   Úbeta)
r"   r?   r   r	   r—   Úcdfr   rO   rQ   r›   )r%   r&   ÚargsÚr1Úr2s        r*   Útest_callable_cdfzTestCvm.test_callable_cdf    s©   € Ý”)˜A‘,”, 
ˆ4ˆÝ˜A�}Ô2Ô6Ñ7Ô7ˆÝ˜A˜wÑ'Ô'ˆÝ�b”l B¤IÐ.°´¸r¼yÐ0IÑJÔJÐJå˜A�}Ô1Ô5°tÑ<Ô<ˆÝ˜A˜v tÑ,Ô,ˆÝ�b”l B¤IÐ.°´¸r¼yÐ0IÑJÔJÐJÐJÐJr,   N)r_   r`   ra   rm   ro   rr   rx   r{   r   r‚   r‰   rM   ÚmarkÚparametrizer�   r˜   r    rb   r,   r*   rd   rd   N   sæ   € € € € € ðð ð ðð ð ðð ð ðð ð ðJð Jð Jð
ð ð ð1ð 1ð 1ð$ð $ð $ð „[×Ò˜S 2¨ u +Ñ.Ô.ð-ð -ñ /Ô.ð-ð<ð <ð <ð$Kð Kð Kð Kð Kr,   rd   c                   ó€  — e Zd Zej                             ddg idg ig g dœg¦  «        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¦  «        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 d/¢g d0¢d1œZg d2¢g d3¢g d4¢g d5¢g d6¢g d7¢d8œZd9„ 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                             d?e¦  «        d@„ ¦   «         ZdA„ Zej                             dBddg¦  «        dC„ ¦   «         ZdD„ Zg d+¢dEdFdGdHej        dEdIdJdKdKdLgdMdNfg d+¢dEdFdGdHej        ej        dIdJdKdKdLgdOdPfdJdIej        dKgdEdFdGdHej        dEdIdJdKdKdLgdQdRfdJdIej        dKgdEdFdGdHej        ej        dIdJdKdKdLgdSdTfdJej        ej        dKgdEdFdGdHej        ej        dIdJdKdKdLgdUdVfgZej                             dWe¦  «        dX„ ¦   «         Zg dY¢g dZ¢g d[¢g d\¢g d]¢g d^¢g d_¢g d`¢g da¢g	Z ej                             dbe ¦  «        dc„ ¦   «         Z!dd„ Z"g d&¢dedfgddggg d&¢dedfgddggg d&¢dedfgd
dhgg d&¢dIgddigg d&¢dIgddigg d&¢dIgd
djgdJdIgdJdIgddkgdJdIgdJdIgddkgdJdIgdJdIgd
dlgg	Z#ej                             g dm¢e#¦  «        dn„ ¦   «         Z$do„ Z%ej                             dpg dq¢¦  «        dr„ ¦   «         Z&dsS )tÚTestMannWhitneyUÚkwargs_updater&   r'   ©r&   r'   c                 ó°  — t          j        ddg¦  «        }t          j        ddg¦  «        }t          ||¬¦  «        }|                     |¦  «         t	          j        t          t          ¬¦  «        5  t          di |¤Ž}t          |j
        t           j        ¦  «         t          |j        t           j        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr   r.   r!   r/   r¦   rK   rb   )r"   r#   ÚdictÚupdaterM   rN   r   r   r   r   rO   rP   rQ   )r%   r¥   r&   r'   ÚkwargsrS   s         r*   Ú
test_emptyzTestMannWhitneyU.test_empty²   s  € õ ŒH�a˜�VÑÔˆÝŒH�a˜�VÑÔˆÝ˜˜Q�‘”ˆØ�Š�mÑ$Ô$Ð$ÝŒ\Õ,Õ4IÐJÑJÔJð 	-ð 	-ÝÐ(Ð( Ð(Ð(ˆCÝ˜œ­¬Ñ/Ô/Ð/Ý˜œ¥R¤VÑ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð 	-ð 	-s   Á3ACÃCÃCc                 ó`  — t          j        ddg¦  «        }t          j        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.   r!   r/   z`use_continuity` must be onerK   Úekki)Úuse_continuityz`alternative` must be one of©Úalternativez`axis` must be an integerrŠ   ©Úaxisz`method` must be one of©Úmethod)r"   r#   rW   rX   r   rY   s      r*   Útest_input_validationz&TestMannWhitneyU.test_input_validation¾   s/  € ÝŒH�a˜�VÑÔˆÝŒH�a˜�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}t           j                             |dz
  ¦  «        }t           j                             |dz
  ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚t           j                             |dz
  ¦  «        }t           j                             |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 ‚t           j                             |dz   ¦  «        }t           j                             |dz   ¦  «        }t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚t           j                             |dz
  ¦  «        }t           j                             |dz
  ¦  «        }|d         |d<   t	          ||¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚|j        |j        k    sJ ‚d S )Nr   r4   Ú
asymptoticr³   Úexactr!   )r"   r=   r>   Úrandr   rQ   )r%   rˆ   r&   r'   Úautor·   r¸   s          r*   Ú	test_autozTestMannWhitneyU.test_autoÊ   sÜ  € õ 	Œ	�Š�qÑÔÐØˆõ ŒI�NŠN˜1˜Q™3ÑÔˆÝŒI�NŠN˜1˜Q™3ÑÔˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/õ ŒI�NŠN˜1˜Q™3ÑÔˆÝŒI�NŠN˜1˜Q™3ÑÔˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/õ ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/õ ŒI�NŠN˜1˜Q™3ÑÔˆÝŒI�NŠN˜1˜Q™3ÑÔˆÝ˜A˜qÑ!Ô!ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆ
Ý˜Q ¨'Ð2Ñ2Ô2ˆØŒ{˜eœlÒ*Ð*Ð*Ð*ØŒ{˜jÔ/Ò/Ð/Ð/Ð/õ ŒI�NŠN˜1˜Q™3ÑÔˆÝŒI�NŠN˜1˜Q™3ÑÔˆØ�Œ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´   )é   ç
+ž×÷å?Úless)r¾   ç
+ž×÷Õ?Úgreater)r¾   ç¼ŒÉ%c�æ?r¸   )r¾   g9§ƒ:¨ƒæ?)r¾   g9§ƒ:¨ƒÖ?)r¾   g*…:¨ƒ:æ?)ÚkwdsÚexpectedc                 óV   — t          | j        | j        fi |¤Ž}t          ||¦  «         d S ©N)r   r&   r'   r   ©r%   rÄ   rÅ   rS   s       r*   Ú
test_basiczTestMannWhitneyU.test_basic   s3   € å˜4œ6 4¤6Ð2Ð2¨TÐ2Ð2ˆÝ˜˜XÑ&Ô&Ð&Ð&Ð&r,   T)r°   r®   )é   r¿   )rÊ   rÃ   )rÊ   rÁ   F)rÊ   gl,KÎNhä?)rÊ   gÊiÚ˜ØËå?)rÊ   gl,KÎNhÔ?c                 óZ   — t          | j        | j        fddi|¤Ž}t          ||¦  «         d S )Nr´   r·   )r   r'   r&   r   rÈ   s       r*   Útest_continuityz TestMannWhitneyU.test_continuity2  s:   € õ ˜4œ6 4¤6ÐGÐG°,ÐGÀ$ÐGÐGˆÝ˜˜XÑ&Ô&Ð&Ð&Ð&r,   c                 ó^  — g d¢}t          j        g d¢¦  «        }t          j        g d¢¦  «        dz  }t          j        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 )NrJ   ©r   r.   r!   r/   r0   )r   r   r   r   r   rg   )r   r   r   r   r   éÿÿÿÿr·   )r²   r´   )r2   é	   ç      !@r4   r–   r3   r1   )r   g]�UÄÚì?g…æ[»é?giœ…\­æ?g°ZX<_Õã?gãžxº.á?gåß ê…
Ù?)r"   r#   r   r   rO   r   rQ   )	r%   r&   Úy0ÚdyÚdy2r'   rS   Ú
U_expectedÚ
p_expecteds	            r*   Útest_tie_correctz!TestMannWhitneyU.test_tie_correct@  sÛ   € ð ˆLˆLˆÝŒX�o�o�oÑ&Ô&ˆÝŒX�o�o�oÑ&Ô& tÑ+ˆÝŒh���Ñ'Ô'¨Ñ,ˆØ�‰W�b˜‘e˜R ™V R¨¨C©°°B±¸¸4¹Ð@ˆÝ˜1˜a b°Ð>Ñ>Ô>ˆØ/Ð/Ð/ˆ
ðIð Ið Iˆ
å�S”] JÑ/Ô/Ð/Ý˜œ
 JÑ/Ô/Ð/Ð/Ð/r,   )g      Ð?ri   g      è?)r�   çš™™™™™É?çš™™™™™Ù?r�   )rh   r�   rØ   r   ri   gÍÌÌÌÌÌä?©r   r.   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çû©ñÒá?rJ   )gÇK7‰A`Å?gZd;ßOÕ?ri   gòÒMbXå?)çªñÒMb¨?çR¸…ëQ¸?çR¸…ëQÈ?ççû©ñÒMÒ?rÝ   rÞ   )	g;ßO�—n’?ç;ßO�—n¢?ç“V-²?g      À?gJ+‡É?râ   gôýÔxé&Ù?ri   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çû©ñÚ?ri   gºI+‡â?rÎ   )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?rg   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¸á?)r5   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.   r!   r/   r0   r1   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   )r!   r/   r0   r1   ©Úkr5   r6   r   rÏ   )Úsetattrr   r   Úpn3Úpn4Úpm5Úpm6Úitemsr"   r?   Úlenrò   Ú
set_shapesr   rœ   ÚsfÚpmf)
r%   Úp_tablesrˆ   ÚtableÚmr)   ÚuÚu2rþ   Úpmf2s
             r*   Útest_exact_distributionz(TestMannWhitneyU.test_exact_distributionp  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¦  «         t           j                             d¦  «        }t           j                             d¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚t          j        |j        |j        z
  ¦  «        dk    sJ ‚t           j                             d¦  «        }t           j                             d¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }|j        |j        k    sJ ‚t          j        |j        |j        z
  ¦  «        dk     sJ ‚d S )	Nr   r0   r¸   r³   r·   rg   é(   r5   )r"   r=   r>   r¹   r   rO   ÚabsrQ   )r%   r&   r'   Úres1Úres2s        r*   Útest_asymptotic_behaviorz)TestMannWhitneyU.test_asymptotic_behaviorŒ  s.  € Ý
Œ	�Š�qÑÔÐõ ŒI�NŠN˜1ÑÔˆÝŒI�NŠN˜1ÑÔˆÝ˜A˜q¨Ð1Ñ1Ô1ˆÝ˜A˜q¨Ð6Ñ6Ô6ˆØŒ~ ¤Ò/Ð/Ð/Ð/ÝŒv�d”k D¤KÑ/Ñ0Ô0°4Ò7Ð7Ð7Ð7õ ŒI�NŠN˜2ÑÔˆÝŒI�NŠN˜2ÑÔˆÝ˜A˜q¨Ð1Ñ1Ô1ˆÝ˜A˜q¨Ð6Ñ6Ô6ˆØŒ~ ¤Ò/Ð/Ð/Ð/ÝŒv�d”k D¤KÑ/Ñ0Ô0°4Ò7Ð7Ð7Ð7Ð7Ð7r,   c                 óþ   — t          g d¢ddgdd¬¦  «        }t          g d¢ddgdd¬¦  «        }t          |j        |j        ¦  «         |j        dk    sJ ‚t          g d¢ddgd	d¬¦  «        }t          |d
¦  «         d S )NrÚ   rŠ   ç      @rÀ   r¸   r½   rÂ   ri   r¼   )r!   r   )r   r   rQ   )r%   Úres_lÚres_grS   s       r*   Útest_exact_U_equals_meanz)TestMannWhitneyU.test_exact_U_equals_mean¡  s°   € õ ˜Y˜Y˜Y¨¨c¨
ÀØ$+ð-ñ -ô -ˆå˜Y˜Y˜Y¨¨c¨
À	Ø$+ð-ñ -ô -ˆå�U”\ 5¤<Ñ0Ô0Ð0ØŒ|˜cÒ!Ð!Ð!Ð!å˜9˜9˜9 s¨C j¸kØ")ð+ñ +ô +ˆå�S˜&Ñ!Ô!Ð!Ð!Ð!r,   ©r   r   )r   ri   )r   g‹·éƒ¡Eï?)rÄ   Úresultc                 ó:   — t          t          di |¤Ž|¦  «         d S )Nr   r.   ©r   r.   )r   r   )r%   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³   )ri   r   r·   F)r´   r®   ri   )r   r   r"   rP   rl   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                 ó  — t           j                             d¦  «         d}d\  }}t           j                             |dd¦  «        }t           j                             d|dd¦  «        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                             |j        |
¦  «         t           j                             |j        |	¦  «         d S )Nr   éýÿÿÿ)r3   r2   r!   r4   r1   r   r�   )r´   r²   )r1   r!   r4   rÏ   )N.c                 ó,   — g | ]}t          |¦  «        ‘ŒS rb   )Úrange)Ú.0Úis     r*   ú
<listcomp>z8TestMannWhitneyU.test_gh_12837_11113.<locals>.<listcomp>ó  s   € Ð 9Ð 9Ð 9¨a¥ q¡¤Ð 9Ð 9Ð 9r,   r³   )r"   r=   r>   r¹   r   rQ   ÚshaperO   ÚmoveaxisÚndimÚbroadcast_toÚzerosr   Útestingr   )r%   r´   r²   r  rˆ   r&   r'   rS   r  Ú
statisticsÚpvaluesÚindicesÚxiÚyiÚtemps                  r*   Útest_gh_12837_11113z$TestMannWhitneyU.test_gh_12837_11113Ñ  s  € õ 	Œ	�Š�qÑÔÐð ˆØ‰ˆˆ1ÝŒI�NŠN˜1˜a Ñ#Ô#ˆÝŒI�NŠN˜1˜a  AÑ&Ô&¨Ñ,ˆÝ˜1˜a¨°TÐ:Ñ:Ô:ˆàˆØŒ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ÑÐå
Œ
×Ò ¤
¨GÑ4Ô4Ð4Ý
Œ
×Ò ¤¨zÑ:Ô:Ð:Ð:Ð:r,   c                 ó–  — g d¢}g d¢}t          ||¦  «        }t          j        |d<   t          ||¦  «        }t          |j        |j        ¦  «         t          |j        |j        ¦  «         t          j        |d<   t          ||¦  «        }t          |j        t          j        ¦  «         t          |j        t          j        ¦  «         d S )NrJ   )r!   r1   r3   r4   rÐ   r!   r.   r   r/   r/   r0   r/   )r   r"   rV   r   rO   rQ   rP   )r%   r&   r'   r	  r
  Úres3s         r*   Útest_gh_11355zTestMannWhitneyU.test_gh_11355ý  s³   € àˆLˆLˆØ-Ð-Ð-ˆÝ˜A˜qÑ!Ô!ˆõ Œvˆˆ!‰Ý˜A˜qÑ!Ô!ˆå�T”^ T¤^Ñ4Ô4Ð4Ý�T”[ $¤+Ñ.Ô.Ð.õ Œvˆˆ!‰Ý˜A˜qÑ!Ô!ˆÝ�T”^¥R¤VÑ,Ô,Ð,Ý�T”[¥"¤&Ñ)Ô)Ð)Ð)Ð)r,   r!   r1   r3   r4   r.   r   r/   r0   r2   g+z’¿QœÀ?rÑ   g}$kù\¶?g     €1@g!Ë›G*ã?r¾   g›Á,ÖsÿÝ?ç     €8@gêFÏHïQé?)r&   r'   rO   rQ   c                 ó†   — t          ||d¬¦  «        }t          |j        |d¬¦  «         t          |j        |d¬¦  «         d S )Nr·   r³   çê-�™—q=r6   )r   r   rO   rQ   )r%   r&   r'   rO   rQ   rS   s         r*   Útest_gh_11355bzTestMannWhitneyU.test_gh_11355b   sL   € õ ˜1˜a¨Ð5Ñ5Ô5ˆÝ˜œ y°uÐ=Ñ=Ô=Ð=Ý˜œ
 F°Ð7Ñ7Ô7Ð7Ð7Ð7r,   )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}d}d}t          |||||¬¦  «        }t          |j        |¦  «         t          |j        |¦  «         d S )Né#   )
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´   )r   r   rO   r   rQ   )	r%   r®   r°   r´   r3  Ústatistic_expr&   r'   rS   s	            r*   Útest_gh_9184zTestMannWhitneyU.test_gh_91841  s_   € ð. ˆØHˆØ*ˆÝ˜1˜a°Ø'2¸6ðCñ Cô Cˆå�S”] MÑ2Ô2Ð2Ý˜œ
 JÑ/Ô/Ð/Ð/Ð/r,   c                 óÎ  — t          j        t           j        t           j        t           j        t           j        t           j        g¦  «        }t          j        t           j        t           j        t           j        t           j        t           j        g¦  «        }t          ||¦  «        }t	          |j        t           j        ¦  «         t	          |j        t           j        ¦  «         d S rÇ   )r"   r#   rP   r   r   rO   rQ   )r%   ÚaÚbrS   s       r*   Útest_gh_4067zTestMannWhitneyU.test_gh_4067P  s‡   € åŒH•b”f�bœf¥b¤f­b¬fµb´fÐ=Ñ>Ô>ˆÝŒH•b”f�bœf¥b¤f­b¬fµb´fÐ=Ñ>Ô>ˆÝ˜1˜aÑ Ô ˆÝ�S”]¥B¤FÑ+Ô+Ð+Ý�S”Z¥¤Ñ(Ô(Ð(Ð(Ð(r,   rŠ   r  )r!   ga×€}¢ã?)r!   r�   )rŠ   g¤Û?hÍå?)rŠ   r   )r.   gæ5&#\å?)r.   r   )r&   r'   r°   rÅ   c                 óR   — t          ||d|d¬¦  «        }t          ||d¬¦  «         d S )NTr·   r6  r1  ©Úrtol)r   r   )r%   r&   r'   r°   rÅ   rS   s         r*   Útest_gh_2118zTestMannWhitneyU.test_gh_2118h  s>   € õ ˜1˜a°À+Ø".ð0ñ 0ô 0ˆå˜˜X¨EÐ2Ñ2Ô2Ð2Ð2Ð2r,   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 )r0   r“   ©Úsizerò   r   r¸   r³   rÏ   r   rÂ   )r´   r°   )r"   r=   Údefault_rngrõ   r   r   rò   ÚresetÚstatsr   Úconfigurationsr  ÚminrO   )r%   Úrngr  rˆ   r&   r'   rS   r  s           r*   Útest_gh19692_smaller_tablez+TestMannWhitneyU.test_gh19692_smaller_tablep  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 )
NrB  )r.   r0   rC  )r.   r1   r   )r´   r°   r²   r¸   çVçž¯Ò<r>  )	r"   r=   rE  rG  r   ÚPermutationMethodr   rO   rQ   )r%   r°   rJ  r&   r'   rS   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,   N)'r_   r`   ra   rM   r¡   r¢   r«   rµ   r»   r&   r'   Úcases_basicrÉ   Úcases_continuityrÌ   r×   rö   r÷   rø   rù   r  r  r  Úcases_scalarr  r  r+  r.  r"   rV   Úcases_11355r2  Ú
cases_9184r8  r<  Ú
cases_2118r@  rK  rO  rb   r,   r*   r¤   r¤   «   sM  € € € € € ð „[×Ò˜_°°R¨y¸3À¸)Ø57¸bÐ/AÐ/Að/Cñ Dô Dð-ð -ñDô Dð-ð
.ð 
.ð 
.ð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°;Ñ?Ô?ð'ð 'ñ @Ô?ð'ð *5ÈÐMÐMØ.ð0à)/À4ÐHÐHØ.ð0à)2ÀdÐKÐKØ.ð0à)4ÈÐNÐNØ.ð0à)/À5ÐIÐIØ.ð0à)2ÀeÐLÐLØ.ð0ð1Ðð „[×ÒÐ1Ð3CÑDÔDð'ð 'ñ EÔDð'ð0ð 0ð 0ð,  ÐÐÐ$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ð*"ð "ð "ð$ &1¸LÐIÐIØðà%+°|ÐDÐDØðà%.¸,ÐGÐGØ(ð*à%0¸GÐDÐDÀfÐMØ%+°wÐ?Ð?ÀÐJØ%.¸'ÐBÐBÀFÐKðM€Lð „[×ÒÐ/°Ñ>Ô>ð<ð <ñ ?Ô>ð<ðHð Hð Hð „[×Ò˜X¨°gÐ'>Ñ?Ô?ð);ð );ñ @Ô?ð);ðV*ð *ð *ð& !�L�LØ˜˜1˜a ¤¨¨A¨q°!°Q¸Ð:Ø˜ð)ð !�L�LØ˜˜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ð8ð 8ñ MÔLð8ð ?Ð>Ð>ØBÐBÐBØCÐCÐCØ@Ð@Ð@ØCÐCÐCØEÐEÐEØ:Ð:Ð:Ø=Ð=Ð=Ø?Ð?Ð?ðA€Jð „[×Òð 6Ø7AñCô Cð0ð 0ñCô Cð0ð:)ð )ð )ð �9�9˜s C˜j¨)Ð5IÐJØ�9�9˜s C˜j¨&Ð2FÐGØ�9�9˜s C˜j¨+°xÐ@Ø�9�9˜q˜c 9Ð.CÐDØ�9�9˜q˜c 6Ð+@ÐAØ�9�9˜q˜c ;°Ð9Ø�q�6˜A˜q˜6 9Ð.AÐBØ�q�6˜A˜q˜6 6Ð+>Ð?Ø�q�6˜A˜q˜6 ;°Ð7ð9€Jð „[×ÒÐBÐBÐBÀJÑOÔOð3ð 3ñ PÔOð3ð:ð :ð :ð: „[×Ò˜]Ð,LÐ,LÐ,LÑMÔMð	=ð 	=ñ NÔMð	=ð 	=ð 	=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 )Nr2   r  c                 ó6   •— g | ]}‰j         |         d          ‘ŒS )r   )Ú	arguments)r  Úidxr%   s     €r*   r  z,TestSomersD.setup_method.<locals>.<listcomp>¡  s%   ø€ ÐHÐHÐH°#�t”~ cÔ*¨1Ô-ÐHÐHÐHr,   r¼   r¯   )ÚALL_INTEGERÚ	ALL_FLOATÚdtypesr"   r?   rZ  Ú	functoolsÚpartialrG  ÚsomersdÚpartialfuncrÅ   )r%   Ú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 )NrM  r6   )rb  r   rO   rÅ   rQ   ©r%   r�   rS   s      r*   ÚpythranfunczTestSomersD.pythranfuncª  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)rz   é   é   r3   r   )r3   rj  é   r5  r„   )r   r!   r.   r3   é   rM  r6   )rG  ra  Úget_optional_argsr   rO   rQ   )r%   r   r	  Úoptional_argsr
  s        r*   Útest_pythranfunc_keywordsz%TestSomersD.test_pythranfunc_keywords¯  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)r0   r.   r   r!   r1   r/   r3   r4   )r0   r.   r1   r!   r   r4   r3   r/   ©ç        r�   r   rM  r6   r   )	r   r0   r.   r   r!   r1   r/   r3   r4   )	r0   r.   r   r1   r!   r   r4   r3   r/   )r0   r.   r   r!   r1   r/   r3   )r0   r.   r1   r!   r   r3   r/   )g+$I’$IÂ¿g=/3n£+ä?r2   ©r�   r   )
r   r.   r   r!   r/   r1   r0   r3   r4   rÐ   )gsÒ'}Ò'í?rr  rÏ   )g      ð¿r   )
rÐ   r3   r4   r1   r0   r!   r/   r.   r   r   )g}Ò'}Ò'í¿rr  )r„   r.   r   r„   r.   )r   r/   r3   r   r   )ç      à¿gÞ.Ê‚ƒÓ?)r.   r.   r.   )r.   r   r.   g      $@g      4@)rG  ra  r   rO   rQ   r"   r?   r#   rP   rW   rX   )r%   r&   r'   rÅ   rS   Ú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.   r!   r!   r   r.   r.   r.   r.   r.   r.   r.   r!   r!   r!   r!   r!   r!   )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((ÄÇ·?rM  r6   rj   )r!   r.   ©r.   r!   )rG  ra  r   rO   rQ   r   r   r  )r%   r&   r'   Úd_crÚd_rcr)   rS   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 )Nr4   r.   r1   r0   r!   r/   r   gHHHHHHØ?ri  r   éU   r;   )gœÜMî&wã?r�   gtÑE]tá¿)r"   r#   ÚTr   rG  ra  rO   )r%   r   ÚdyxÚdxys       r*   Útest_somers_originalz TestSomersD.test_somers_originalB  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                 ó`  — 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}   ©r/   r1   r   ©r)   r.   r   r±   gaƒ¸yò½¿rM  r6   gP—j¹$Ä?)r"   Úprodr=   r>   rG  ÚmultinomialÚrvsr‡   Úreshapera  Úinsertr#  r   rO   rQ   )r%   ÚNr  rD  rò   rS   Ús2r
  Ús3r-  Ús4Úres4s               r*   Ú*test_contingency_table_with_zero_rows_colsz6TestSomersD.test_contingency_table_with_zero_rows_colsX  sÇ  € ð ˆØˆÝŒw�u‰~Œ~ˆå
Œ	�Š�qÑÔÐÝÔ×!Ò! !¥r¤w¨t¡}¤}°TÑ'9Ð!Ñ:Ô:×BÒBÀ5ÑIÔIˆÝŒ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                 óx  — 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-negativerK   rg   z5All elements of the contingency table must be integerz?At least two elements of the contingency table must be nonzero.r   )r!   r!   r  )r"   r†  r=   r>   rG  r‡  rˆ  r‡   r‰  rW   rX   ra  r#  )	r%   r‹  r  rD  rò   Ús5ÚmessageÚs6Ús7s	            r*   Útest_invalid_contingency_tablesz+TestSomersD.test_invalid_contingency_tablesw  s‹  € ØˆØˆÝŒw�u‰~Œ~ˆå
Œ	�Š�qÑÔÐåÔ×!Ò! !¥r¤w¨t¡}¤}°TÑ'9Ð!Ñ:Ô:×BÒBÀ5ÑIÔIˆà�‰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   ÂB?Â?CÃCÃ'DÄDÄDÄ+EÅEÅEÅ.FÆFÆFÇG'Ç'G+Ç.G+ÈH/È/H3È6H3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Ú   rÏ   gÍÌÌÌÌÌ @)r!   r.   r   r   rt  )r"   rV   rG  ra  r   rO   rQ   )r%   r&   rv  r'   Úy2rS   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 )Nr2   )r"   r?   rG  ra  r   r   ÚeyerR   s       r*   Útest_contingency_table_returnz)TestSomersD.test_contingency_table_return¥  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Î   )r0   r1   r3   r4   r3   r¼   r¯   r   rÀ   r   r.   rÂ   z`alternative` must be...rK   ú	ekki-ekki)
rG  ra  rO   r   r   rQ   ÚreverserM   r   rX   )r%   ru  rv  rÅ   rS   s        r*   Útest_somersd_alternativez$TestSomersD.test_somersd_alternative¬  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 )	Nr2   r   rÏ   r¼   r¯   r   rÀ   rÂ   )r"   r?   ÚfliprG  ra  rO   rQ   )r%   r¡  ru  rv  Úexpected_statisticrS   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¦  «         t          j        ||¬¦  «        }t          j        ||¬¦  «        }d}t          j        ||¦  «        j        }t          ||d¬¦  «         d S )	Nr   r.   é@B i¡´_ ró   g H�„z	Y¿rM  r6   )r=   r>   ÚchoicesrG  ra  rO   r   )r%   ÚclassesÚ	n_samplesr&   r'   Úval_sklearnÚ	val_scipys          r*   Ú!test_somersd_large_inputs_gh18132z-TestSomersD.test_somersd_large_inputs_gh18132î  s‚   € ð �a�&ˆØˆ	ÝŒ�GÑÔÐÝŒN˜7 iÐ0Ñ0Ô0ˆÝŒN˜7 iÐ0Ñ0Ô0ˆð ,ˆõ ”M ! QÑ'Ô'Ô1ˆ	Ý˜ Y°UÐ;Ñ;Ô;Ð;Ð;Ð;r,   N)r_   r`   ra   rd  rg  ro  rw  r|  r‚  r�  r–  r™  rœ  r   rM   r¡   r¢   r¥  r­  rb   r,   r*   rW  rW  š  s  € € € € € ð7ð 7ð 7ðFð Fð Fð
	>ð 	>ð 	>ðm7ð m7ð m7ð^.ð .ð .ð2?ð ?ð ?ð,1ð 1ð 1ð>!ð !ð !ðF	.ð 	.ð 	.ð,ð ,ð ,ð';ð ';ð ';ðR „[×ÒÐ3°]ÑCÔCð@ð @ñ DÔCð@ð0<ð <ð <ð <ð <r,   rW  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  r2   é'   )gÏXyq@g{2s&QÇ7?r}   r.   rq   r0   )gøl¯€lü¿gEA]0×KÁ?r3   r4   )ç*õ)1Ö%ÀgØ_  ¬“?r   )g_½Òc1÷?gèß=ç Ä?é   r\   )g5PyQ ý¿g�Q@Ì2ý°?r¾   ri  )g»ÔÌÑÁù¿g¼ÂJ"¸ò½?)g_½Òc1÷¿gwÝ�Ù„¹Æ?r   r/   )gý7ŠËé§@g      Œ?r!   )gš~øtñ†ñ¿ç,À�?3Oà?r1   )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   rO   rQ   r   ©r%   Úinput_samplerÅ   rS   rO   rQ   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 0rK   r   r.   r!   r/   r   ©rˆ   ú,The input `table` must be of shape \(2, 2\).r1   ú*All values in `table` must be nonnegative.rÏ   zI`alternative` should be one of {'two-sided', 'less', 'greater'}, found .*únot-correct)rW   rX   r   r"   r?   r‰  ©r%   Ú	error_msgs     r*   Útest_raiseszTestBarnardExact.test_raisesD  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rs  c                 ó˜   — t          |¦  «        }|j        |j        }}t          ||d         ¦  «         t          ||d         ¦  «         d S ©Nr   r   ©r   rO   rQ   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_zeroj  sL   € õ ˜LÑ)Ô)ˆØœM¨3¬:�6ˆ	Ý�V˜X aœ[Ñ)Ô)Ð)Ý�Y ¨¤Ñ,Ô,Ð,Ð,Ð,r,   )r³  gEç\/?„?éÈ   é,  )gÐgÑìíQ5Àrr  é   r<   i¥  )gü&çìX}>Àrr  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Â   NrÏ   r¯   çH¯¼šò×z>r6   )r"   r#   r   rO   rQ   r   )	r%   rº  rÅ   r°   Úexpected_statÚless_pvalue_expectrS   rO   rQ   s	            r*   Útest_less_greaterz"TestBarnardExact.test_less_greaterw  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)r_   r`   ra   Ú__doc__rM   r¡   r¢   r»  r¾  rÇ  rÌ  r"   rP   rÎ  rÖ  rb   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.   r3   r4   )ç<vB\÷’?gèÁ–„/?„?r0   r   r2   )g¬“ŽÍéMï?gßÏA>î?r¾   r\   ri  )ç_¬VÃÑ¶?gÂñ¥…Ö­?)ç‘uÝ Ø%Å?gc'�íµ?r   r/   ©r   r   r!   )ri   g      Ö?r„   )çß+Â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¯   r6   N©r   rO   rQ   r   ÚATOLr¹  s         r*   Ú	test_lesszTestBoschlooExact.test_lessŸ  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×ëÛ|?r1   )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¯   r6   Nrâ  r¹  s         r*   Útest_greaterzTestBoschlooExact.test_greater¼  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ÛŒ?)ri   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ˆ   r6   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À  rK   r   r.   r!   r/   r   rÁ  rÂ  r1   rÃ  rÏ   zK`alternative` should be one of \('two-sided', 'less', 'greater'\), found .*rÄ  )rW   rX   r   r"   r?   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   rO   rQ   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”   r„   rÀ   r¯   rÂ   r.   r¼   r�   )r   rQ   rI  )r%   Ú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.   r3   r4   r¯   r   )r   rO   rG  Úfisher_exactr   )r%   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)r_   r`   ra   r×  rã  rM   r¡   r¢   rä  rè  rë  rÇ  r"   rP   rÎ  ró  rø  rb   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                   ó6  — e Zd Zej                             dg  ej        d¦  «        f ej        d¦  «        dgfg¦  «        d„ ¦   «         Zd„ Z	d„ Z
d„ Zej                             dg d	¢¦  «        d
„ ¦   «         Zej        j        d„ ¦   «         Zd„ Zd„ Zd„ ZdS )ÚTestCvm_2sampr�   r0   r   c                 ó  — t          j        t          t          ¬¦  «        5  t	          |Ž }t          |j        t          j        ¦  «         t          |j	        t          j        ¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )NrK   )
rM   rN   r   r   r   r   rO   r"   rP   rQ   rf  s      r*   Útest_too_small_inputz"TestCvm_2samp.test_too_small_input4  s±   € õ Œ\Õ,Õ4IÐJÑJÔJð 	-ð 	-Ý&¨Ð-ˆCÝ˜œ­¬Ñ/Ô/Ð/Ý˜œ¥R¤VÑ,Ô,Ð,ð	-ð 	-ð 	-ñ 	-ô 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-ð 	-øøøð 	-ð 	-ð 	-ð 	-ð 	-ð 	-s   ¡AA6Á6A:Á=A:c                 ó¼   — t          j        d¦  «        }d}t          j        t          |¬¦  «        5  t          ||d¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nr0   z/method must be either auto, exact or asymptoticrK   Úxyz)r"   r?   rM   r   rX   r   )r%   r'   Úmsgs      r*   r�   z TestCvm_2samp.test_invalid_input<  sŸ   € ÝŒ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.   r!   r/   r3   r1   )rØ   rš   r„   rk  )r   r"   r#   r   rO   rQ   ©r%   r&   r'   rž   rŸ   s        r*   Útest_list_inputzTestCvm_2samp.test_list_inputB  sl   € ØˆOˆOˆØÐÐˆÝ! ! QÑ'Ô'ˆÝ!¥"¤(¨1¡+¤+­r¬x¸©{¬{Ñ;Ô;ˆÝ�b”l B¤IÐ.°´¸r¼yÐ0IÑJÔJÐJÐJÐJr,   c                 ó’   — g d¢}g d¢}t          ||¦  «        }t          |j        dd¬¦  «         t          |j        dd¬¦  «         d S )N)	gffffff@gÍÌÌÌÌÌ @r   gffffff!@çš™™™™™"@gÍÌÌÌÌÌ#@g333333$@g333333%@gffffff&@)gÍÌÌÌÌÌ@gÍÌÌÌÌÌ@gš™™™™™@ç      @ç333333@gffffff @g333333"@gš™™™™™#@gš™™™™™%@gš™™™™™&@g      '@gš™™™™™(@g      )@gÍÌÌÌÌÌ*@g333333-@gøSã¥›ÄÐ?r5   r6   g
×£p=
Ç?rg   )r   r   rO   rQ   ©r%   r&   r'   Úrs       r*   Útest_example_conoverz"TestCvm_2samp.test_example_conoverI  sc   € ð =Ð<Ð<ˆðð ð ˆå   AÑ&Ô&ˆÝ˜œ U°Ð6Ñ6Ô6Ð6Ý˜œ $¨TÐ2Ñ2Ô2Ð2Ð2Ð2r,   zstatistic, m, n, pval))iÆ  r0   r1   g¾cj`ï˜º?)ii  r3   r3   gtÑE]t±?)i@  r/   r1   g8�8�ƒ?)iä  r1   r3   g¶ëéXwS?c                 óD   — t          t          |||¦  «        |¦  «         d S rÇ   )r   r   )r%   rO   r  rˆ   Úpvals        r*   Útest_exact_pvaluezTestCvm_2samp.test_exact_pvalueS  s'   € õ 	Õ*¨9°a¸Ñ;Ô;¸TÑBÔBÐBÐBÐBr,   c                 ó’  — t           j                             d¦  «         t          j                             d¬¦  «        }t          j                             d¬¦  «        }t          ||¦  «        }t          d|j        cxk     odk     nc ¦  «         t          ||dz   ¦  «        }t          d|j        cxk     odk     nc ¦  «         d S )Ni  r§  rC  i » r   r   r�   )	r"   r=   r>   r	   r†   rˆ  r   r   rQ   r  s       r*   Útest_large_samplezTestCvm_2samp.test_large_sample^  sÊ   € õ 	Œ	�Š�tÑÔÐÝÔ×"Ò"¨Ð"Ñ0Ô0ˆÝÔ×"Ò"¨Ð"Ñ/Ô/ˆÝ   AÑ&Ô&ˆÝ��A”HÐ Ð Ò Ð ˜qÒ Ð Ð Ð Ñ!Ô!Ð!Ý   A c¡EÑ*Ô*ˆÝ��A”HÐ Ð Ò Ð ˜qÒ Ð Ð Ð Ñ!Ô!Ð!Ð!Ð!r,   c                 ót  — t           j                             d¦  «         t           j                             d¦  «        }t           j                             d¦  «        }t	          ||d¬¦  «        }t	          ||d¬¦  «        }t          |j        |j        ¦  «         t          |j        |j        d¬¦  «         d S )	Nr   r3   r4   r¸   r³   r·   rg   r6   )	r"   r=   r>   r¹   r   r   rO   r   rQ   r  s        r*   Útest_exact_vs_asymptoticz&TestCvm_2samp.test_exact_vs_asymptoticj  s—   € Ý
Œ	�Š�qÑÔÐÝŒI�NŠN˜1ÑÔˆÝŒI�NŠN˜1ÑÔˆÝ! ! Q¨wÐ7Ñ7Ô7ˆÝ! ! Q¨|Ð<Ñ<Ô<ˆÝ�R”\ 2¤<Ñ0Ô0Ð0Ý˜œ	 2¤9°4Ð8Ñ8Ô8Ð8Ð8Ð8r,   c                 óV  — t          j        d¦  «        }g d¢}t          ||d¬¦  «        }t          ||d¬¦  «        }t          |j        |j        ¦  «         t          j        d¦  «        }t          ||d¬¦  «        }t          ||d¬¦  «        }t          |j        |j        ¦  «         d S )Nr\   )ri   gÍÌÌÌÌÌ@g333333*@r¸   r³   rº   rÑ  r·   )r"   r?   r   r   rQ   r  s        r*   Útest_method_autozTestCvm_2samp.test_method_autos  s¥   € ÝŒI�b‰MŒMˆØÐÐˆÝ! ! Q¨wÐ7Ñ7Ô7ˆÝ! ! Q¨vÐ6Ñ6Ô6ˆÝ�R”Y ¤	Ñ*Ô*Ð*åŒI�b‰MŒMˆÝ! ! Q¨|Ð<Ñ<Ô<ˆÝ! ! Q¨vÐ6Ñ6Ô6ˆÝ�R”Y ¤	Ñ*Ô*Ð*Ð*Ð*r,   c                 óþ   — t          j        d¦  «        }t          ||¦  «        }t          |j        |j        fd¦  «         t          |d d…         |d d…         ¦  «        }t          |j        |j        fd¦  «         d S )Nr´  rq  r/   )r"   r?   r   r   rO   rQ   rŒ   s      r*   Útest_same_inputzTestCvm_2samp.test_same_input  sw   € õ ŒI�b‰MŒMˆÝ" 1 aÑ(Ô(ˆÝ�c”m S¤ZÐ0°*Ñ=Ô=Ð=å" 1 R a R¤5¨!¨B¨Q¨B¬%Ñ0Ô0ˆÝ�c”m S¤ZÐ0°*Ñ=Ô=Ð=Ð=Ð=r,   N)r_   r`   ra   rM   r¡   r¢   r"   r?   rü  r�   r  r	  r  Úxslowr  r  r  r  rb   r,   r*   rú  rú  3  s1  € € € € € Ø„[×Ò˜V r¨9¨2¬9°Q©<¬<Ð&8Ø'0 r¤y°¡|¤|°a°SÐ&9ð&;ñ <ô <ð-ð -ñ<ô <ð-ð.ð .ð .ðKð Kð Kð3ð 3ð 3ð „[×ÒÐ4ð5ð 5ð 5ñ6ô 6ð
Cð Cñ6ô 6ð
Cð „[Ôð	"ð 	"ñ Ôð	"ð9ð 9ð 9ð
+ð 
+ð 
+ð	>ð 	>ð 	>ð 	>ð 	>r,   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e	j
                             dd¦  «        d„ ¦   «         Ze	j
                             d g d!¢¦  «        d"„ ¦   «         Zd#„ Zd$S )%ÚTestTukeyHSD)r/  ç     €7@çffffff:@çš™™™™;@çfffffæ=@)çffffff<@çš™™™™A@ç     €=@çš™™™™@@çš™™™™>@)gš™™™™:@gÍÌÌÌÌL<@gÍÌÌÌÌL8@g333333:@gÍÌÌÌÌÌ;@)r/  r  gHáz®G:@r  r  r  r   r   )r/  r  r  )
r  r  r  r  r   r  r  r  r  r   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,atolrj   g»½×Ùß|Û=)úequal size samplez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;
        ú - ú r0   N©Údtype)r1   r1   r   r6   rh   ©r"   ÚasarrayÚreplaceÚsplitÚfloatr‰  rG  Ú	tukey_hsdÚconfidence_intervalÚintr   ÚlowrO   ÚhighrQ   )r%   ÚdataÚres_expect_strr7   Ú
res_expectÚ	res_tukeyÚconfr  ÚjÚlrò   ÚhÚsigs                r*   Útest_compare_saszTestTukeyHSD.test_compare_sas¶  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
        r1  rÓ  r!  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&  ©r!   r1   r   r6   N)r"   r)  r+  r,  r‰  rG  r-  r.  r/  r   r0  rO   r1  rQ   )r%   r2  r3  r7   r4  r5  r6  r  r7  r8  rò   r9  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
        r$  r%  r0   Nr&  r=  )é   r;   é6   ri  éF   é4   é3   r@  éC   rz   rj  é   é   rF  é   é)   r\   é,   )rk  rÑ  rF  rl  r„   rk  r5  r;   é$   é*   r@  rG  r¾   r²  r<   rÑ  r²  rF  )rK  rÑ  rß  rk  r2   r±  r<   r´  r@  r\   rÑ  rß  rl  r”   r´  r´  r¾   r<   r   rÓ  r6   r‘   gñhãˆµøä>r(  )r%   Ústr_resr4  r2  r5  r6  r  r7  rò   r8  r9  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@g      @)gš™™™™™ @r  g333333@gffffff"@r  )g       @g      %@g333333 @rP  r  )rQ  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   r!   r  ry  rg   r6   N)rG  r-  r.  r"   r)  r   r0  r1  )r%   Úgroup1Úgroup2Úgroup3Úgroup4rS   r6  ÚlowerÚupperr  r7  s              r*   Útest_engineering_stat_handbookz+TestTukeyHSD.test_engineering_stat_handbook2  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¦  «         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 )Nr:   r!   r}   ©r   r   r   r   )r"   r=   r>   r¹   rG  r-  r.  r   r0  r1  r  ÚdiagonalrO   rQ   )r%   r2  rS   r6  s       r*   Útest_rand_symmzTestTukeyHSD.test_rand_symmL  s  € å
Œ	�Š�tÑÔÐÝŒy�~Š~˜a Ñ%Ô%ˆÝŒ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.rK   rÚ   r.   )r1   r3   r!   )rW   rX   rG  r-  r"   rV   rl   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-dimensionalrK   r   r.   r!   r0   )r0   rÊ   r1   ©rW   rX   rG  r-  rl   s    r*   Ú
test_is_1dzTestTukeyHSD.test_is_1dc  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 onerK   r.   r0   )r/   r0   r1   ra  rl   s    r*   Útest_no_emptyzTestTukeyHSD.test_no_emptyg  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Ú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.rK   ©rÊ   r3   r!   ra  )r%   re  s     r*   Útest_not_enough_treatmentsz'TestTukeyHSD.test_not_enough_treatmentsk  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)rt  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 1rK   rg  r!   r/   rÐ   )rW   rX   rG  r-  r.  )r%   ri  r  s      r*   Útest_conf_level_invalidz$TestTukeyHSD.test_conf_level_invalidp  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  rR  )rG  r-  Údata_diff_sizeÚ	ttest_indr   rQ   )r%   r5  Ú	res_ttests      r*   Útest_2_args_ttestzTestTukeyHSD.test_2_args_ttestv  sn   € å”O TÔ%8¸¸!¸Ô%<Ð=ˆ	Ý”O TÔ%8¸¸!¸Ô%<Ð=ˆ	Ý˜	Ô(¨)Ô*:¸4Ô*@ÑAÔAÐAÝ˜	Ô(¨)Ô*:¸4Ô*@ÑAÔAÐAÐAÐAr,   N)r_   r`   ra   Údata_same_sizerm  Úextreme_sizeÚsas_same_sizeÚsas_diff_sizeÚsas_extremerM   r¡   r¢   r;  Úmatlab_sm_sizÚmatlab_diff_szr>  rN  rY  r]  r_  rb  rd  rh  rk  rp  rb   r,   r*   r  r  ‹  sU  € € € € € à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ð „[×Ò˜W fÑ-Ô-ð5ð 5ñ .Ô-ð5ð „[×Ò˜T > > >Ñ2Ô2ð&ð &ñ 3Ô2ð&ð
Bð Bð 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}   r!   r}   gþe÷äa¡¶?)r.   r}   r1   r}   gÞ	ŠcÆ?c                 ób   — t          j        ||||¦  «        }t          |j        |d¬¦  «         d S )Nrj   r6   ©rG  Úpoisson_means_testr   rQ   )r%   Úc1Ún1Úc2Ún2Úp_expectrS   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\   r2   r\   r2   gŽŸ{}ÿÿï?r¼   r   )r2   r2   r2   r2   goPFªÿÿï?r¼   r   )é2   r´  r   r   góæ°ae¹?r¼   rh   )r!   r}   r\   rÐ  gò/V›-=¿?r¼   r   )r!   r„   r/   r\   g"Ü)‘ÅÝÙ?rÂ   r   )r/   r\   r!   r}   g_ù'Qm~€?rÂ   r   )r/   r\   r!   r2   g|¨Š»Ó?rÀ   r   )r   r   rƒ  r´  gï0�ëÝ·?rÀ   r   c                 ój   — t          j        ||||||¬¦  «        }t          |j        |dd¬¦  «         d S )N)r°   Údiffg�íµ ÷ÆÀ>g¼‰Ø—²Òœ<©r7   r?  r{  )	r%   r}  r~  r  r€  r�  ÚaltÚdrS   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{  ©r%   Úcount1Úcount2Únobs1Únobs2rS   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.rK   rš   z1`k1` and `k2` must be greater than or equal to 0.rÏ   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¯   )rW   Ú	TypeErrorrG  r|  rX   )r%   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)r_   r`   ra   rM   r¡   r¢   r‚  r‰  r�  r’  rµ   rb   r,   r*   ry  ry  ~  s  € € € € € Ø„[×ÒÐ7à Ð Ð Ø Ð Ð ð:ñ ô ð
9ð 9ñô ð
9ð „[×ÒÐ?ð
 	=Ð<Ð<Ø<Ð<Ð<Ø=Ð=Ð=Ø>Ð>Ð>Ø9Ð9Ð9Ø;Ð;Ð;ð 	6Ð5Ð5Ø6Ð6Ð6ðBñ ô ð"Eð Eñ#ô ð"Eð'ð 'ð 'ð'ð 'ð 'ð!Fð !Fð !Fð !Fð !Fr,   ry  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.   r3   rC  z,`x` and `y` must be exactly one-dimensional.rK   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...rL  r³   )	r"   r=   rE  rM   r   rX   rG  Úbws_testrP   )r%   rJ  r&   r'   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.   r!   r/   r1   r3   r4   )r0   rÐ   r2   r“   r„   r”   rj  r¼   r¯   gºI+‡@r5   r6   g¤f$/•Þg?)rG  r›  r   rO   r   rQ   rR   s       r*   Ú test_against_published_referencez,TestBWSTest.test_against_published_referenceî  s_   € ð "Ð!Ð!ˆØ&Ð&Ð&ˆÝŒn˜Q ¨{Ð;Ñ;Ô;ˆÝ˜œ u°4Ð8Ñ8Ô8Ð8Ý�S”Z Ñ)Ô)Ð)Ð)Ð)r,   )r°   rO   rQ   ))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š  rC  r¯   ç‚vIhÂ%<=r>  rg   r�   r†  ©r"   r=   rE  rG  r›  r   rO   rQ   ©r%   r°   rO   rQ   rJ  r&   r'   rS   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 rÐ   rC  r4   r¯   r¡  r>  rg   r�   r†  r¢  r£  s           r*   Útest_against_R_imbalancedz%TestBWSTest.test_against_R_imbalanced	  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.   r2   rC  r2   )Ún_resamplesrJ  r³   l   ÝVC£	ãA )
r"   r=   rE  rG  rN  r›  rû   Únull_distributionr   Úallclose)r%   rJ  r&   r'   r´   r	  r
  r-  s           r*   Útest_methodzTestBWSTest.test_method  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©  r0   rC  r   rÂ   r¯   r   rÀ   )
r"   r=   rE  rG  r›  rO   r   rQ   rû   r«  )r%   rJ  r&   r'   rS   s        r*   Útest_directionszTestBWSTest.test_directions2  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)r_   r`   ra   rœ  rž  rM   r¡   r¢   r¤  r§  r­  r¯  rb   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—  ).Ú	itertoolsr   Únumpyr"   r=   r_  rM   Únumpy.testingr   r   r   r   r   rW   Úscipy.statsrG  r	   Úscipy.stats._hypotestsr
   r   r   r   r   r   r   Úscipy.stats._mannwhitneyur   r   r   Úcommon_testsr   Úscipy._lib._testutilsr   Úscipy.stats._axis_nan_policyr   r   r   rd   r¤   rW  r¯  rÙ  rú  r  ry  r—  rb   r,   r*   ú<module>r¹     sr  ðØ Ð Ð Ð Ð Ð à Ð Ð Ð Ø €€€Ø Ð Ð Ð Ø €€€ð0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0ð 0à *Ð *Ð *Ð *Ð *Ð *à Ð Ð Ð Ð Ð Ø %Ð %Ð %Ð %Ð %Ð %ð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Ø RÐ RÐ RÐ RÐ RÐ RÐ RÐ Rð4-ð 4-ð 4-ð 4-ð 4-ñ 4-ô 4-ð 4-ðnZKð ZKð ZKð ZKð ZKñ ZKô ZKð ZKðzl=ð l=ð l=ð l=ð l=ñ l=ô l=ð l=ð^j<ð j<ð j<ð j<ð j<Ð"ñ j<ô j<ð j<ðZ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U>ð U>ð U>ð U>ð U>ñ U>ô U>ð U>ðppBð pBð pBð pBð pBñ pBô pBð pBðf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,   