Ë
    âQ(hj8 ã                   ó–  — 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-y) é    )Ú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                   ó0   — e Zd Zd„ Zd„ Zd„ Zd„ Zd„ Zd„ Zy)ÚTestEppsSingletonc                 ó¶   — t        j                  g d¢«      }t        j                  g d¢«      }t        ||«      \  }}t        |dd¬«       t        |dd¬«       y )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        ú^/var/www/skyplay_api_hub/venv/lib/python3.12/site-packages/scipy/stats/tests/test_hypotests.pyÚtest_statistic_1z"TestEppsSingleton.test_statistic_1   sR   € ô
 �H‰Hò 7ó 8ˆä�H‰Hò 4ó 5ˆä# A qÓ)‰ˆˆ1Ü˜A˜u¨aÕ0Ü˜A˜w°Ö2ó    c                 ó®   — t        j                  d«      }t        j                  d«      }t        ||«      \  }}t        |dd¬«       t	        |dd¬«       y )	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%   sP   € ä�H‰Hð "ó #ˆä�H‰Hð  ó !ˆä# A qÓ)‰ˆˆ1Ü˜˜5 uÕ-Ü˜A˜w°Ö2r,   c                 ó¦  — t         j                  j                  d«       t        j                  d«      t        j                  d«      }}t	        t        |«      t        |«      «      \  }}t	        t        |«      t        |«      «      \  }}t	        ||«      \  }}t        ||cxk(  xr |k(  nc «       t        ||cxk(  xr
 |k(  «       y c «       y )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‰y˜‹}œbŸi™i¨›mˆ1ˆä%¤d¨1£g¬t°A«wÓ7‰ˆˆBÜ%¤e¨A£h´°a³Ó9‰ˆˆBÜ% a¨Ó+‰ˆˆBä��b–˜B”ÔÜ��b–˜B‘Õ‘Õr,   c                 óB  — dt        j                  d«      }}t        j                  t        t
        ¬«      5  t        ||«      }t        |j                  t         j                  «       t        |j                  t         j                  «       d d d «       y # 1 sw Y   y xY w)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:   sh   € àœRŸY™Y r›]ˆ1ˆÜ�\‰\Ô,Ô4IÔJñ 	-Ü& q¨!Ó,ˆCÜ˜Ÿ™¬¯©Ô/Ü˜Ÿ™¤R§V¡VÔ,÷	-÷ 	-ñ 	-ús   ·ABÂBc                 ó†   — dddddt         j                  ft        j                  d«      }}t        t        t
        ||«       y )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   s3   € à�1�a˜˜AœrŸv™vÐ&¬¯	©	°"«ˆ1ˆÜ”jÔ"6¸¸1Õ=r,   c                 óŒ   — t        j                  d«      t        j                  d«      }}t        ||«      }d}t        ||«       y )Né   r;   )rO   rQ   )r"   r?   r
   r   )r%   r&   r'   rS   Ú
attributess        r*   Ú
test_nameszTestEppsSingleton.test_namesG   s6   € Ü�y‰y˜‹}œbŸi™i¨›mˆ1ˆÜ" 1 aÓ(ˆØ,ˆ
Ü˜C Õ,r,   N)	Ú__name__Ú
__module__Ú__qualname__r+   r8   rH   rT   rZ   r^   © r,   r*   r   r      s    „ ò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                  j                  d	d
dgg«      d„ «       Zd„ Zd„ Zy)ÚTestCvmc                 ó<   — t        t        g d¢d«      g d¢d¬«       y )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   s   € ÜÜÒ=¸qÓAÚ(Øö	r,   c                 ó<   — t        t        g d¢d«      g d¢d¬«       y )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   s   € ÜÜÒ=¸rÓBÚ(Øö	r,   c                 ó<   — t        t        g d¢d«      g d¢d¬«       y )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^   s   € ÜÜÒ=¸tÓDÚ(Øö	r,   c                 ó:   — t        t        g d¢«      g d¢d¬«       y )N)çaÃÓ+e™?ç+û®þ·¢?çÉ&pën¾?ç+MJA·—ò?rf   rj   r6   rk   rl   s    r*   Útest_cdf_infzTestCvm.test_cdf_infd   s   € ÜÜÒ=Ó>Ú(Øö	r,   c                 ól   — t        t        ddgd«      ddg«       t        t        ddgd«      ddg«       y )	Ng²Xª(~$?gUUUUU5f@i  r   r   g†a†ah?g«ªªªªª"@é   )r   r   rl   s    r*   Útest_cdf_supportzTestCvm.test_cdf_supportj   s6   € ä”X˜z¨5Ð1°3Ó7¸!¸Q¸Ô@Ü”X˜°
Ð;¸RÓ@À1ÀaÀ&ÕIr,   c                 óN   — t        t        g d¢d«      t        g d¢«      d¬«       y )N)rt   ru   rv   rw   éd   é'  rj   r6   rk   rl   s    r*   Útest_cdf_large_nzTestCvm.test_cdf_large_no   s!   € äÜÒBÀEÓJÜÒBÓCØö	r,   c                 ó˜   — t        dt        dd«      cxk  xr dk  nc «       t        dt        d«      cxk  xr
 dk  «       y c «       y )NgwJëÿï?gÍÌÌÌÌÔt@rq   ç      ð?)r   r   rl   s    r*   Útest_large_xzTestCvm.test_large_xv   s9   € ô 	�œ( 5¨$Ó/Ö5°#Ô5Ô6Ü�œ( 5›/Ö/¨CÑ/Õ0Ñ/Õ0r,   c                 ó¼   — d}t        t        j                  |«      dz  d«      }t        t	        |j
                  |«      dkD  «       t        |j                  d«       y )Né   çš™™™™™é?Únormr�   r   )r   r"   Úonesr   r   rO   r   rQ   )r%   ÚnrS   s      r*   Ú
test_low_pzTestCvm.test_low_p   sG   € ð ˆÜœRŸW™W Q›Z¨™^¨VÓ4ˆÜ”˜Ÿ™¨Ó*¨SÑ0Ô1Ü�S—Z‘Z Õ#r,   r&   rb   ç      ø?c                 ó  — t        j                  t        t        ¬«      5  t	        |d«      }t        |j                  t        j                  «       t        |j                  t        j                  «       d d d «       y # 1 sw Y   y xY w)NrK   r†   )
rM   rN   r   r   r   r   rO   r"   rP   rQ   ©r%   r&   rS   s      r*   Útest_invalid_inputzTestCvm.test_invalid_input‡   sW   € ä�\‰\Ô,Ô4IÔJñ 	-Ü   FÓ+ˆCÜ˜Ÿ™¬¯©Ô/Ü˜Ÿ™¤R§V¡VÔ,÷	-÷ 	-ñ 	-ús    AA>Á>Bc                 óz  — 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¬«       y )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ˆÜ˜Ÿ™ x°dÕ;Ü˜Ÿ
™
 I°DÕ9ô Ò;¸VÀXÓNˆÜ˜Ÿ™ y°tÕ<Ü˜Ÿ
™
 K°dÕ;ô ÒCÀWÓ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«       y )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Ñ2×6Ñ6Ó7ˆÜ˜A˜wÓ'ˆÜ�b—l‘l B§I¡IÐ.°·±¸r¿y¹yÐ0IÔJä˜Aœ}×1Ñ1×5Ñ5°tÓ<ˆÜ˜A˜v tÓ,ˆÜ�b—l‘l B§I¡IÐ.°·±¸r¿y¹yÐ0IÕJr,   N)r_   r`   ra   rm   ro   rr   rx   r{   r   r‚   r‰   rM   ÚmarkÚparametrizer�   r˜   r    rb   r,   r*   rd   rd   N   sa   „ òòòòòJò
ò1ò$ð ‡[�[×Ñ˜S 2¨ u +Ó.ñ-ó /ð-ò<ó$Kr,   rd   c                   óÒ  — e Zd Zej                  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                  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                  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                  j                  d?e«      d@„ «       ZdA„ Zej                  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                  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                  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                  j                  g dm¢e#«      dn„ «       Z$do„ Z%ej                  j                  dpg dq¢«      dr„ «       Z&ys)tÚTestMannWhitneyUÚkwargs_updater&   r'   ©r&   r'   c                 óª  — t        j                  ddg«      }t        j                  ddg«      }t        ||¬«      }|j                  |«       t	        j
                  t        t        ¬«      5  t        di |¤Ž}t        |j                  t         j                  «       t        |j                  t         j                  «       d d d «       y # 1 sw Y   y xY w)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‰H�a˜�VÓˆÜ�H‰H�a˜�VÓˆÜ˜˜Q”ˆØ�‰�mÔ$Ü�\‰\Ô,Ô4IÔJñ 	-ÜÑ( Ñ(ˆCÜ˜Ÿ™¬¯©Ô/Ü˜Ÿ™¤R§V¡VÔ,÷	-÷ 	-ñ 	-ús   Á,AC	Ã	Cc                 ó   — t        j                  ddg«      }t        j                  ddg«      }t        t        d¬«      5  t	        ||d¬«       d d d «       t        t        d	¬«      5  t	        ||d¬
«       d d d «       t        t        d¬«      5  t	        ||d¬«       d d d «       t        t        d¬«      5  t	        ||d¬«       d d d «       y # 1 sw Y   Œ‚xY w# 1 sw Y   ŒfxY w# 1 sw Y   ŒJxY w# 1 sw Y   y xY w)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‰H�a˜�VÓˆÜ�H‰H�a˜�VÓˆÜœ:Ð-KÔLñ 	6Ü˜˜A¨fÕ5÷	6äœ:Ð-KÔLñ 	3Ü˜˜A¨6Õ2÷	3äœ:Ð-HÔIñ 	)Ü˜˜A CÕ(÷	)äœ:Ð-FÔGñ 	.Ü˜˜A fÕ-÷	.ð 	.÷	6ð 	6ú÷	3ð 	3ú÷	)ð 	)ú÷	.ð 	.ús0   Á CÁ(CÂC(Â8C4ÃCÃC%Ã(C1Ã4C=c                 ó"  — t         j                  j                  d«       d}t         j                  j                  |dz
  «      }t         j                  j                  |dz
  «      }t	        ||«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k(  sJ ‚|j
                  |j
                  k7  sJ ‚t         j                  j                  |dz
  «      }t         j                  j                  |dz   «      }t	        ||«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k(  sJ ‚|j
                  |j
                  k7  sJ ‚t	        ||«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k(  sJ ‚|j
                  |j
                  k7  sJ ‚t         j                  j                  |dz   «      }t         j                  j                  |dz   «      }t	        ||«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k7  sJ ‚|j
                  |j
                  k(  sJ ‚t         j                  j                  |dz
  «      }t         j                  j                  |dz
  «      }|d   |d<   t	        ||«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k7  sJ ‚|j
                  |j
                  k(  sJ ‚y )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‰I�N‰N˜1˜Q™3ÓˆÜ�I‰I�N‰N˜1˜Q™3ÓˆÜ˜A˜qÓ!ˆÜ! ! Q¨|Ô<ˆ
Ü˜Q ¨'Ô2ˆØ�{‰{˜eŸl™lÒ*Ð*Ð*Ø�{‰{˜j×/Ñ/Ò/Ð/Ð/ô �I‰I�N‰N˜1˜Q™3ÓˆÜ�I‰I�N‰N˜1˜Q™3ÓˆÜ˜A˜qÓ!ˆÜ! ! Q¨|Ô<ˆ
Ü˜Q ¨'Ô2ˆØ�{‰{˜eŸl™lÒ*Ð*Ð*Ø�{‰{˜j×/Ñ/Ò/Ð/Ð/ô ˜A˜qÓ!ˆÜ! ! Q¨|Ô<ˆ
Ü˜Q ¨'Ô2ˆØ�{‰{˜eŸl™lÒ*Ð*Ð*Ø�{‰{˜j×/Ñ/Ò/Ð/Ð/ô �I‰I�N‰N˜1˜Q™3ÓˆÜ�I‰I�N‰N˜1˜Q™3ÓˆÜ˜A˜qÓ!ˆÜ! ! Q¨|Ô<ˆ
Ü˜Q ¨'Ô2ˆØ�{‰{˜eŸl™lÒ*Ð*Ð*Ø�{‰{˜j×/Ñ/Ò/Ð/Ð/ô �I‰I�N‰N˜1˜Q™3ÓˆÜ�I‰I�N‰N˜1˜Q™3ÓˆØ�‰tˆˆ!‰Ü˜A˜qÓ!ˆÜ! ! Q¨|Ô<ˆ
Ü˜Q ¨'Ô2ˆØ�{‰{˜eŸl™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                 ó^   — t        | j                  | j                  fi |¤Ž}t        ||«       y ©N)r   r&   r'   r   ©r%   rÄ   rÅ   rS   s       r*   Ú
test_basiczTestMannWhitneyU.test_basic   s%   € ä˜4Ÿ6™6 4§6¡6Ñ2¨TÑ2ˆÜ˜˜XÕ&r,   T)r°   r®   )é   r¿   )rÊ   rÃ   )rÊ   rÁ   F)rÊ   gl,KÎNhä?)rÊ   gÊiÚ˜ØËå?)rÊ   gl,KÎNhÔ?c                 ób   — t        | j                  | j                  fddi|¤Ž}t        ||«       y )Nr´   r·   )r   r'   r&   r   rÈ   s       r*   Útest_continuityz TestMannWhitneyU.test_continuity2  s,   € ô ˜4Ÿ6™6 4§6¡6Ñ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                  |«       y )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¢   € ò ˆÜ�X‰X’oÓ&ˆÜ�X‰X’oÓ& tÑ+ˆÜ�h‰h’Ó'¨Ñ,ˆØ�‰W�b˜‘e˜R ™V R¨¨C©°°B±¸¸4¹Ð@ˆÜ˜1˜a b°Ô>ˆÚ/ˆ
ò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œ}|j                  «       D �]‘  \  }}|j                  «       D �]w  \  }}t        j                  dt        |«      «      }t        j                  j                  ||«       t        t        j                  j                  |¬«      |d¬«       t        j                  d||z  dz   «      }t        t        j                  j                  |¬«      t        j                  j                  |¬«      z   t        j                  j!                  |¬«      z
  d«       t        j                  j!                  |¬«      }t        ||d d d…   «       t        j                  j                  ||«       t        j                  j!                  |¬«      }	t        ||	«       �Œz �Œ” y )	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  sw  € ä”
˜C¤ a¨£Ô,à—x‘x D§H¡H°·±¸d¿h¹hÑGˆØ Ÿ™Ó(ó 	+‰HˆAˆuØŸ™›ó +‘��1ä—I‘I˜a¤ Q£Ó(�Ü—‘×'Ñ'¨¨1Ô-Ü¤
§¡× 0Ñ 0°1Ð 0Ó 5°q¸tÕDô —Y‘Y˜q ! A¡# a¡%Ó(�Ü¤
§¡× 0Ñ 0°2Ð 0Ó 6Ü",§,¡,§/¡/°B /Ó"7ñ!8ä",§,¡,×"2Ñ"2°RÐ"2Ó"8ñ!9à:;ô=ô
 !—l‘l×&Ñ&¨Ð&Ó,�Ü  S©¨2¨¡YÔ/ô —‘×'Ñ'¨¨1Ô-Ü!—|‘|×'Ñ'¨"Ð'Ó-�Ü  TÖ*ò)+ñ	+r,   c                 óÚ  — t         j                  j                  d«       t         j                  j                  d«      }t         j                  j                  d«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k(  sJ ‚t        j                  |j                  |j                  z
  «      dkD  sJ ‚t         j                  j                  d«      }t         j                  j                  d«      }t	        ||d¬«      }t	        ||d¬«      }|j
                  |j
                  k(  sJ ‚t        j                  |j                  |j                  z
  «      dk  sJ ‚y )	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‰I�N‰N˜1ÓˆÜ�I‰I�N‰N˜1ÓˆÜ˜A˜q¨Ô1ˆÜ˜A˜q¨Ô6ˆØ�~‰~ §¡Ò/Ð/Ð/Ü�v‰v�d—k‘k D§K¡KÑ/Ó0°4Ò7Ð7Ð7ô �I‰I�N‰N˜2ÓˆÜ�I‰I�N‰N˜2ÓˆÜ˜A˜q¨Ô1ˆÜ˜A˜q¨Ô6ˆØ�~‰~ §¡Ò/Ð/Ð/Ü�v‰v�d—k‘k D§K¡KÑ/Ó0°4Ò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kD  sJ ‚t        g d¢ddgd	d¬«      }t        |d
«       y )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¡  sw   € ô šY¨¨c¨
ÀØ$+ô-ˆäšY¨¨c¨
À	Ø$+ô-ˆä�U—\‘\ 5§<¡<Ô0Ø�|‰|˜cÒ!Ð!Ð!äš9 s¨C j¸kØ")ô+ˆä�S˜&Õ!r,   ©r   r   )r   ri   )r   g‹·éƒ¡Eï?)rÄ   Úresultc                 ó.   — t        t        di |¤Ž|«       y )N©r   r.   )r   r   )r%   rÄ   r  s      r*   Útest_scalar_dataz!TestMannWhitneyU.test_scalar_data½  s   € ô 	œÑ2¨TÑ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«       y )	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Â  sP   € ô
 	”\ ! Q¨wÔ7¸ÔBÜ”\ ! Q¨|Ô<¸hÔGô 	”\ ! Q¨|Ø16ô8Ø:=¼r¿v¹v¸õ	Hr,   r´   c                 óx  — t         j                  j                  d«       d}d\  }}t         j                  j                  |dd«      }t         j                  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 �cg c]  }t        |«      ‘Œ c}Ž D ]8  }||   }||   }t	        |||¬«      }|j                  |	|<   |j
                  |
|<   Œ: t         j                  j                  |j
                  |
«       t         j                  j                  |j                  |	«       y c c}w )Nr   éýÿÿÿ)r3   r2   r!   r4   r1   r   r�   )r´   r²   )r1   r!   r4   rÏ   )N.r³   )r"   r=   r>   r¹   r   rQ   ÚshaperO   ÚmoveaxisÚndimÚbroadcast_toÚzerosr   ÚrangeÚtestingr   )r%   r´   r²   r  rˆ   r&   r'   rS   r  Ú
statisticsÚpvaluesÚiÚindicesÚxiÚyiÚtemps                   r*   Útest_gh_12837_11113z$TestMannWhitneyU.test_gh_12837_11113Ñ  sî  € ô 	�	‰	�‰�qÔð ˆØ‰ˆˆ1Ü�I‰I�N‰N˜1˜a Ó#ˆÜ�I‰I�N‰N˜1˜a  AÓ&¨Ñ,ˆÜ˜1˜a¨°TÔ:ˆàˆØ�z‰z×Ñ 5Ò(Ð(Ð(Ø�}‰}×"Ñ" eÒ+Ð+Ð+ô �{‰{˜1˜d BÓ'¬¯©°Q¸¸bÓ)Aˆ1ˆàˆi‰LˆØ�v‰v˜Ÿ™ÒÐÐä�O‰O˜A˜u¨ t™|Ó,ˆÜ�O‰O˜A˜u¨ t™|Ó,ˆØ�w‰w�s˜ˆ|˜uÒ$Ð$Ð$Ø�w‰w�s˜ˆ|˜uÒ$Ð$Ð$ô —X‘X˜e“_ˆ
Ü—(‘(˜5“/ˆÜ°5Ö 9¨a¤ q¥Ò 9Ð:ò 	+ˆGØ�7‘ˆBØ�7‘ˆBÜ  B¨vÔ6ˆDØ"&§.¡.ˆJ�wÑØ#Ÿ{™{ˆG�GÒð	+ô 	�
‰
×Ñ §
¡
¨GÔ4Ü
�
‰
×Ñ §¡¨zÕ:ùò !:s   ÆH7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                  «       y )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–   € âˆÚ-ˆÜ˜A˜qÓ!ˆô �v‰vˆˆ!‰Ü˜A˜qÓ!ˆä�T—^‘^ T§^¡^Ô4Ü�T—[‘[ $§+¡+Ô.ô �v‰vˆˆ!‰Ü˜A˜qÓ!ˆÜ�T—^‘^¤R§V¡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¬«       y )Nr·   r³   çê-�™—q=r6   )r   r   rO   rQ   )r%   r&   r'   rO   rQ   rS   s         r*   Útest_gh_11355bzTestMannWhitneyU.test_gh_11355b   s2   € ô ˜1˜a¨Ô5ˆÜ˜Ÿ™ y°uÕ=Ü˜Ÿ
™
 F°Ö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                  |«       y )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´   r0  Ústatistic_expr&   r'   rS   s	            r*   Útest_gh_9184zTestMannWhitneyU.test_gh_91841  sE   € ð. ˆØHˆØ*ˆÜ˜1˜a°Ø'2¸6ôCˆä�S—]‘] MÔ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                  «       y rÇ   )r"   r#   rP   r   r   rO   rQ   )r%   ÚaÚbrS   s       r*   Útest_gh_4067zTestMannWhitneyU.test_gh_4067P  sŒ   € ä�H‰H”b—f‘fœbŸf™f¤b§f¡f¬b¯f©f´b·f±fÐ=Ó>ˆÜ�H‰H”b—f‘fœbŸf™f¤b§f¡f¬b¯f©f´b·f±fÐ=Ó>ˆÜ˜1˜aÓ ˆÜ�S—]‘]¤B§F¡FÔ+Ü�S—Z‘Z¤§¡Õ(r,   rŠ   r  )r!   ga×€}¢ã?)r!   r�   )rŠ   g¤Û?hÍå?)rŠ   r   )r.   gæ5&#\å?)r.   r   )r&   r'   r°   rÅ   c                 ó@   — t        ||d|d¬«      }t        ||d¬«       y )NTr·   r3  r.  ©Úrtol)r   r   )r%   r&   r'   r°   rÅ   rS   s         r*   Útest_gh_2118zTestMannWhitneyU.test_gh_2118h  s%   € ô ˜1˜a°À+Ø".ô0ˆä˜˜X¨EÖ2r,   c                 ó¼  — t         j                  j                  d«      }d\  }}|j                  |¬«      }|j                  |¬«      }t        t        dt        dd«      «       t        j                  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                  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 ‚y )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  sq  € ô
 �i‰i×#Ñ#Ð$7Ó8ˆØ‰ˆˆ1Ø�J‰J˜AˆJÓˆØ�J‰J˜AˆJÓˆä”
˜C¤ a¨£Ô,Ü�‰×ÑÔä× Ñ   A¨gÔ6ˆÜ—‘×+Ñ+×1Ñ1ˆØ�R‰yœC §¡¨q°©s°S·]±]Ñ/BÓCÀaÑGÒGÐGÐGÜ×Ñ˜1˜a¨Õ0Øœ
Ÿ™×3Ñ3×9Ñ9Ò9Ð9Ð9ô
 	�‰×ÑÔÜ×Ñ˜1˜a ™c¨'¸yÕIÜ—‘×+Ñ+×1Ñ1ˆØ�R‰y˜AŠ~Ðˆ~Ü×Ñ˜1˜Q™3 ¨'¸yÕIØœ
Ÿ™×3Ñ3×9Ñ9Ò9Ð9Ñ9r,   r°   )rÀ   rÂ   r¼   c                 óž  — t         j                  j                  d«      }|j                  d¬«      }|j                  d¬«      }t        j                  ||t        j
                  «       |d¬«      }t        j                  ||d|d¬«      }t        |j                  |j                  d¬	«       t        |j                  |j                  d¬	«       y )
Nr?  )r.   r0   r@  )r.   r1   r   )r´   r°   r²   r¸   çVçž¯Ò<r;  )	r"   r=   rB  rD  r   ÚPermutationMethodr   rO   rQ   )r%   r°   rG  r&   r'   rS   r
  s          r*   Útest_permutation_methodz(TestMannWhitneyU.test_permutation_method�  s    € ä�i‰i×#Ñ#Ð$7Ó8ˆØ�J‰J˜FˆJÓ#ˆØ�J‰J˜FˆJÓ#ˆÜ× Ñ   A¬e×.EÑ.EÓ.GØ-8¸qôBˆä×!Ñ! ! Q¨wØ.9ÀôCˆä˜Ÿ™ t§~¡~¸EÕBÜ˜Ÿ
™
 D§K¡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_11355r/  Ú
cases_9184r5  r9  Ú
cases_2118r=  rH  rL  rb   r,   r*   r¤   r¤   «   sL  „ ð ‡[�[×Ñ˜_°°R¨y¸3À¸)Ø57¸bÑ/Að/Có Dñ-óDð-ò
.ò10òj 	-€Aò	€Að& %0¸<ÑHØ)ð+à$*°lÑCØ)ð+à$-¸ÑFØ)ð+à$/¸7ÑCØ)ð+à$*°gÑ>Ø)ð+à$-¸ÑAØ)ð+ð,€Kð ‡[�[×ÑÐ1°;Ó?ñ'ó @ð'ð *5ÈÑMØ.ð0à)/À4ÑHØ.ð0à)2ÀdÑKØ.ð0à)4ÈÑNØ.ð0à)/À5ÑIØ.ð0à)2ÀeÑLØ.ð0ð1Ðð ‡[�[×ÑÐ1Ð3CÓDñ'ó Eð'ò0ò,  Ò$8Ú.ñ0€CâÒ"AÚ<ÚKñM€Cò )Ú7ÚKò2ò<ñ=€Cò +Ú?ò1ò9òGò$ñ%€Cò+ò88ò*"ð$ &1¸LÑIØðà%+°|ÑDØðà%.¸,ÑGØ(ð*à%0¸GÑDÀfÐMØ%+°wÑ?ÀÐJØ%.¸'ÑBÀFÐKðM€Lð ‡[�[×ÑÐ/°Ó>ñ<ó ?ð<òHð ‡[�[×Ñ˜X¨°gÐ'>Ó?ñ);ó @ð);òV*ò& !Ø˜˜1˜a §¡¨¨A¨q°!°Q¸Ð:Ø˜ð)ò !Ø˜˜1˜a §¡¨¯©°°A°q¸!¸QÐ?ØÐ)ð+ð ˜˜2Ÿ6™6 1Ð%Ø˜˜1˜a §¡¨¨A¨q°!°Q¸Ð:Ø˜/ð+ð ˜˜2Ÿ6™6 1Ð%Ø˜˜1˜a §¡¨¯©°°A°q¸!¸QÐ?Ø˜ð)ð ˜Ÿ™ §¡¨Ð*Ø˜˜1˜a §¡¨¯©°°A°q¸!¸QÐ?Ø˜/ð+ð,€Kð  ‡[�[×ÑÐ>ÀÓLñ8ó Mð8ò ?ÚBÚCÚ@ÚCÚEÚ:Ú=Ú?ðA€Jð ‡[�[×Ñð 6Ø7AóCñ0óCð0ò:)ò ˜s C˜j¨)Ð5IÐJÚ˜s C˜j¨&Ð2FÐGÚ˜s C˜j¨+°xÐ@Ú˜q˜c 9Ð.CÐDÚ˜q˜c 6Ð+@ÐAÚ˜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ÀJÓOñ3ó Pð3ò:ð: ‡[�[×Ñ˜]Ò,LÓMñ	=ó Nñ	=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                  j!                  dd«      d„ «       Zd„ Zy)ÚTestSomersDc                 óÜ  — | j                   | j                  z   | _        t        j                  d«      | j                   | j                  z   ft        j                  d«      | j                   | j                  z   fdœ| _        | j
                  D �cg c]  }| j
                  |   d   ‘Œ }}t        j                  t        j                  d¬«      | _
         | j                  |Ž | _        y c c}w )Nr2   r  r   r¼   r¯   )ÚALL_INTEGERÚ	ALL_FLOATÚdtypesr"   r?   Ú	argumentsÚ	functoolsÚpartialrD  ÚsomersdÚpartialfuncrÅ   )r%   ÚidxÚinput_arrays      r*   Úsetup_methodzTestSomersD.setup_method›  sÈ   € Ø×&Ñ&¨¯©Ñ7ˆŒÜ Ÿi™i¨›mØ"×.Ñ.°·±Ñ?ðAä Ÿi™i¨›mØ"×.Ñ.°·±Ñ?ðAñBˆŒð :>¿¹ÖH°#�t—~‘~ cÑ*¨1Ó-ÐHˆÐHô
 %×,Ñ,¬U¯]©]Ø9DôFˆÔà(˜×(Ñ(¨+Ð6ˆ�ùò Is   ÂC)c                 óÒ   —  | j                   |Ž }t        |j                  | j                  j                  d¬«       t        |j                  | j                  j                  d¬«       y )NrJ  r6   )r]  r   rO   rÅ   rQ   ©r%   r�   rS   s      r*   ÚpythranfunczTestSomersD.pythranfuncª  sH   € Øˆd×Ñ Ð%ˆÜ˜Ÿ™ t§}¡}×'>Ñ'>ÀUÕKÜ˜Ÿ
™
 D§M¡M×$8Ñ$8¸uÖEr,   c                 ó6  — g d¢g d¢g d¢g}t        j                  |«      }| j                  t         j                  «      }t        j                  |fi |¤Ž}t        |j                  |j                  d¬«       t        |j
                  |j
                  d¬«       y )N)rz   é   é   r3   r   )r3   rf  é   r2  r„   )r   r!   r.   r3   é   rJ  r6   )rD  r\  Úget_optional_argsr   rO   rQ   )r%   r   r	  Úoptional_argsr
  s        r*   Útest_pythranfunc_keywordsz%TestSomersD.test_pythranfunc_keywords¯  sm   € â#Ò%8Ò:JÐKˆÜ�}‰}˜UÓ#ˆà×.Ñ.¬u¯}©}Ó=ˆÜ�}‰}˜UÑ4 mÑ4ˆä˜Ÿ™¨¯©¸UÕ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                  ||«       y )N)r0   r.   r   r!   r1   r/   r3   r4   )r0   r.   r1   r!   r   r4   r3   r/   ©ç        r�   r   rJ  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Ò'}Ò'í?rn  rÏ   )g      ð¿r   )
rÐ   r3   r4   r1   r0   r!   r/   r.   r   r   )g}Ò'}Ò'í¿rn  )r„   r.   r   r„   r.   )r   r/   r3   r   r   )ç      à¿gÞ.Ê‚ƒÓ?)r.   r.   r.   )r.   r   r.   g      $@g      4@)rD  r\  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‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ò (ˆÚ'ˆà9ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ò "ˆÚ!ˆà:ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ô �I‰I�b‹MˆÜ�I‰I�b‹Mˆð *ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ô �I‰I�b‹MˆÜ�H‰HÒ3Ó4ˆà9ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ô �I‰I�b‹MˆÜ�I‰I�b‹M™$˜B˜$Ñˆð +ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ô �I‰I�b‹MˆÜ�H‰HÒ3Ó4ˆà;ˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ò ˆÚˆà:ˆÜ�m‰m˜B Ó#ˆÜ˜Ÿ™ x°¡{¸Õ?Ü˜Ÿ
™
 H¨Q¡K°eÕ<ô �m‰mšI¢yÓ1ˆÜ˜Ÿ™¤r§v¡vÔ.Ü˜Ÿ
™
¤B§F¡FÔ+ä�m‰mšI¢yÓ1ˆÜ˜Ÿ™¤r§v¡vÔ.Ü˜Ÿ
™
¤B§F¡FÔ+ä�m‰mšI¢yÓ1ˆÜ˜Ÿ™¤r§v¡vÔ.Ü˜Ÿ
™
¤B§F¡FÔ+ä�m‰m˜Q˜C ! Ó%ˆÜ˜Ÿ™¤r§v¡vÔ.Ü˜Ÿ
™
¤B§F¡FÔ+ô �m‰m˜B Ó#ˆÜ˜Ÿ™¤r§v¡vÔ.Ü˜Ÿ
™
¤B§F¡FÔ+ô �I‰I�c‹NˆÜ�I‰I�c‹NˆÜ”j¤%§-¡-°°AÕ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
«       y )N)r   r   r   r.   r.   r.   r.   r.   r!   r!   r   r.   r.   r.   r.   r.   r.   r.   r!   r!   r!   r!   r!   r!   )r   r   r   r   r   r   r   r   r   r   r.   r.   r.   r.   r.   r.   r.   r.   r.   r.   r.   r.   r.   r.   gCÑE]tÑ?gâ^ñ_ñÕ?gO((ÄÇ·?rJ  r6   rj   )r!   r.   ©r.   r!   )rD  r\  r   rO   rQ   r   r   r  )r%   r&   r'   Úd_crÚd_rcr)   rS   s          r*   Útest_asymmetryzTestSomersD.test_asymmetry)  s¥   € ò1ˆò1ˆð !ˆØ ˆØˆÜ�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ t°%Õ8Ü˜Ÿ
™
 A¨DÕ1Ü�S—Y‘Y—_‘_ fÔ-Ü�m‰m˜A˜qÓ!ˆÜ˜Ÿ™ t°%Õ8Ü˜Ÿ
™
 A¨EÕ2Ü�S—Y‘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                  |«       y )Nr4   r.   r1   r0   r!   r/   r   gHHHHHHØ?re  r   éU   r;   )gœÜMî&wã?r�   gtÑE]tá¿)r"   r#   ÚTr   rD  r\  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ˆà—‘ˆØˆÜœŸ™ eÓ,×6Ñ6¸Ô<ô —‘˜2˜q˜' B¨ 7¨Q°¨GÐ4Ó5ˆØ'‰ˆˆSÜœŸ™ eÓ,×6Ñ6¸Ô<ÜœŸ™ e§g¡gÓ.×8Ñ8¸#Ô>ô —‘˜2˜q˜' A r 7¨R°¨GÐ4Ó5ˆØˆÜœŸ™ e§g¡gÓ.×8Ñ8¸#Õ>r,   c                 ó¨  — d}d}t        j                  |«      }t         j                  j                  d«       t        j
                  j                  |t        j                  |«      |z  ¬«      j                  |«      }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                  «       y )Nr}   ©r/   r1   r   ©r)   r.   r   r±   gaƒ¸yò½¿rJ  r6   gP—j¹$Ä?)r"   Úprodr=   r>   rD  ÚmultinomialÚrvsr‡   Úreshaper\  Úinsertr  r   rO   rQ   )r%   ÚNr  rA  rò   rS   Ús2r
  Ús3r*  Ús4Úres4s               r*   Ú*test_contingency_table_with_zero_rows_colsz6TestSomersD.test_contingency_table_with_zero_rows_colsX  s�  € ð ˆØˆÜ�w‰w�u‹~ˆä
�	‰	�‰�qÔÜ×Ñ×!Ñ! !¤r§w¡w¨t£}°TÑ'9Ð!Ó:×BÑBÀ5ÓIˆÜ�m‰m˜AÓˆä�Y‰Y�q˜!œRŸX™X e¨A¡hÓ/°aÔ8ˆÜ�}‰}˜RÓ ˆä�Y‰Y�q˜!œRŸX™X e¨A¡hÓ/°aÔ8ˆÜ�}‰}˜RÓ ˆä�Y‰Y�r˜1œbŸh™h u¨Q¡x°¡zÓ2¸Ô;ˆÜ�}‰}˜RÓ ˆô 	˜Ÿ™Ð'9ÀÕFÜ˜Ÿ™ t§~¡~Ô6Ü˜Ÿ™ t§~¡~Ô6Ü˜Ÿ™ t§~¡~Ô6ä˜Ÿ
™
Ð$5¸EÕBÜ˜Ÿ
™
 D§K¡KÔ0Ü˜Ÿ
™
 D§K¡KÔ0Ü˜Ÿ
™
 D§K¡KÕ0r,   c                 ó   — d}d}t        j                  |«      }t         j                  j                  d«       t        j
                  j                  |t        j                  |«      |z  ¬«      j                  |«      }|dz
  }d}t        t        |¬«      5  t	        j                  |«       d d d «       |dz   }d	}t        t        |¬«      5  t	        j                  |«       d d d «       d
}t        t        |¬«      5  t	        j                  g g«       d d d «       t        t        |¬«      5  t	        j                  dgg«       d d d «       t        j                  d«      }t        t        |¬«      5  t	        j                  |«       d d d «       d|d<   t        t        |¬«      5  t	        j                  |«       d d d «       y # 1 sw Y   �ŒxY w# 1 sw Y   ŒòxY w# 1 sw Y   ŒÌxY w# 1 sw Y   Œ§xY w# 1 sw Y   ŒoxY w# 1 sw Y   y xY w)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>   rD  rƒ  r„  r‡   r…  rW   rX   r\  r  )	r%   r‡  r  rA  rò   Ús5ÚmessageÚs6Ús7s	            r*   Útest_invalid_contingency_tablesz+TestSomersD.test_invalid_contingency_tablesw  s«  € ØˆØˆÜ�w‰w�u‹~ˆä
�	‰	�‰�qÔä×Ñ×!Ñ! !¤r§w¡w¨t£}°TÑ'9Ð!Ó:×BÑBÀ5ÓIˆà�‰UˆØNˆÜœ:¨WÔ5ñ 	Ü�M‰M˜"Ô÷	ð �‰XˆØIˆÜœ:¨WÔ5ñ 	Ü�M‰M˜"Ô÷	ð,ˆäœ:¨WÔ5ñ 	 Ü�M‰M˜2˜$Ô÷	 ô œ:¨WÔ5ñ 	!Ü�M‰M˜A˜3˜%Ô ÷	!ô �X‰X�fÓˆÜœ:¨WÔ5ñ 	Ü�M‰M˜"Ô÷	ð ˆˆ4‰Üœ:¨WÔ5ñ 	Ü�M‰M˜"Ô÷	ð 	÷+	ñ 	ú÷
	ð 	ú÷
	 ð 	 ú÷	!ð 	!ú÷	ð 	ú÷	ð 	úsH   ÂGÃGÃ>G Ä.G,Å4G8Æ(HÇGÇGÇ G)Ç,G5Ç8HÈH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                  «       y )NrÚ   rÏ   gÍÌÌÌÌÌ @)r!   r.   r   r   rp  )r"   rV   rD  r\  r   rO   rQ   )r%   r&   rr  r'   Úy2rS   r
  s          r*   Útest_only_ranks_matterz"TestSomersD.test_only_ranks_matterš  sr   € âˆØ�#”r—v‘vÐˆÚˆØ�œŸ™�wÐˆÜ�m‰m˜A˜qÓ!ˆÜ�}‰}˜R Ó$ˆÜ�S—]‘] D§N¡NÔ3Ü�S—Z‘Z §¡Õ-r,   c                 óÖ   — t        j                  d«      }t        j                  d«      }t        j                  ||«      }t	        |j
                  t        j                  d«      «       y )Nr2   )r"   r?   rD  r\  r   r   ÚeyerR   s       r*   Útest_contingency_table_returnz)TestSomersD.test_contingency_table_return¥  sB   € ä�I‰I�b‹MˆÜ�I‰I�b‹MˆÜ�m‰m˜A˜qÓ!ˆÜ�S—Y‘Y¤§¡ r£
Õ+r,   c                 óL  — g d¢}g d¢}t        j                  ||d¬«      }|j                  dkD  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  «       |j                  «        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 «       y # 1 sw Y   y xY w)NrÎ   )r0   r1   r3   r4   r3   r¼   r¯   r   rÀ   r   r.   rÂ   z`alternative` must be...rK   ú	ekki-ekki)
rD  r\  rO   r   r   rQ   ÚreverserM   r   rX   )r%   rq  rr  rÅ   rS   s        r*   Útest_somersd_alternativez$TestSomersD.test_somersd_alternative¬  s£  € ò ˆÚˆô —=‘=  R°[ÔAˆØ×!Ñ! AÒ%Ð%Ð%ô �m‰m˜B °Ô7ˆÜ�S—]‘] H×$6Ñ$6Ô7Ü˜Ÿ
™
 A¨¯©¸1Ñ)<Ñ$=Ô>ô �m‰m˜B °	Ô:ˆÜ�S—]‘] H×$6Ñ$6Ô7Ü˜Ÿ
™
 H§O¡O°aÑ$7Ô8ð 	�
‰
Œô —=‘=  R°[ÔAˆØ×!Ñ! AÒ%Ð%Ð%ô �m‰m˜B °	Ô:ˆÜ�S—]‘] H×$6Ñ$6Ô7Ü˜Ÿ
™
 A¨¯©¸1Ñ)<Ñ$=Ô>ô �m‰m˜B °Ô7ˆÜ�S—]‘] H×$6Ñ$6Ô7Ü˜Ÿ
™
 H§O¡O°aÑ$7Ô8ä�]‰]œ:Ð-GÔHñ 	;Ü�M‰M˜"˜b¨kÕ:÷	;÷ 	;ñ 	;ús   Ç8HÈH#Ú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 ‚y )	Nr2   r   rÏ   r¼   r¯   r   rÀ   rÂ   )r"   r?   ÚfliprD  r\  rO   rQ   )r%   r�  rq  rr  Úexpected_statisticrS   s         r*   Ú test_somersd_perfect_correlationz,TestSomersD.test_somersd_perfect_correlationÕ  sß   € ô �Y‰Y�r‹]ˆÙ'‰R¬R¯W©W°R«[ˆÙ"6™Q¸BÐô �m‰m˜B °Ô<ˆØ�}‰}Ð 2Ò2Ð2Ð2Ø�z‰z˜QŠÐˆô �m‰m˜B °Ô7ˆØ�}‰}Ð 2Ò2Ð2Ð2Ø�z‰zÑ#7™a¸QÒ?Ð?Ð?ô �m‰m˜B °	Ô:ˆØ�}‰}Ð 2Ò2Ð2Ð2Ø�z‰zÑ#7™a¸QÒ?Ð?Ñ?r,   c                 óö   — ddg}d}t        j                  d«       t        j                  ||¬«      }t        j                  ||¬«      }d}t        j                  ||«      j
                  }t        ||d¬«       y )	Nr   r.   é@B i¡´_ ró   g H�„z	Y¿rJ  r6   )r=   r>   ÚchoicesrD  r\  rO   r   )r%   ÚclassesÚ	n_samplesr&   r'   Úval_sklearnÚ	val_scipys          r*   Ú!test_somersd_large_inputs_gh18132z-TestSomersD.test_somersd_large_inputs_gh18132î  si   € ð �a�&ˆØˆ	Ü�‰�GÔÜ�N‰N˜7 iÔ0ˆÜ�N‰N˜7 iÔ0ˆð ,ˆô —M‘M ! QÓ'×1Ñ1ˆ	Ü˜ Y°UÖ;r,   N)r_   r`   ra   r`  rc  rk  rs  rx  r~  rŒ  r’  r•  r˜  rœ  rM   r¡   r¢   r¡  r©  rb   r,   r*   rT  rT  š  sl   „ ò7òFò
	>òm7ò^.ò2?ò,1ò>!òF	.ò,ò';ðR ‡[�[×ÑÐ3°]ÓCñ@ó Dð@ó0<r,   rT  c                   óÔ  — e Zd ZdZej
                  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
                  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
                  j                  dddgddggd-fg«      d.„ «       Z
ej
                  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
                  j                  dd	dgdd	ggd1fdd2gd3dggd4fd5d6gd7dggd8fg«      ej
                  j                  d9d:d;g«      d<„ «       «       Zy=)>Ú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¾   re  )g»ÔÌÑÁù¿g¼ÂJ"¸ò½?)g_½Òc1÷¿gwÝ�Ù„¹Æ?r   r/   )gý7ŠËé§@g      Œ?r!   )gš~øtñ†ñ¿ç,À�?3Oà?r1   )gÝr?ô~Éø¿çCæ°Yð7É?c                 óf   — t        |«      }|j                  |j                  }}t        ||g|«       y)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
  s.   € ô2 ˜LÓ)ˆØŸM™M¨3¯:©:�6ˆ	Ü˜ FÐ+¨XÕ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                 ój   — t        |d¬«      }|j                  |j                  }}t        ||g|«       y)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'  s0   € ô2 ˜L°Ô7ˆØŸM™M¨3¯:©:�6ˆ	Ü˜ FÐ+¨XÕ6r,   c                 ó  — d}t        t        |¬«      5  t        ddgddggd¬«       d d d «       d	}t        t        |¬«      5  t        t        j                  d
«      j                  dd«      «       d d d «       d}t        t        |¬«      5  t        ddgddgg«       d d d «       d}t        t        |¬«      5  t        ddgddggd«       d d d «       y # 1 sw Y   Œ¯xY w# 1 sw Y   ŒqxY w# 1 sw Y   ŒPxY w# 1 sw Y   y xY w)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  s  € ð Fð 	ô œ:¨YÔ7ñ 	1Ü˜A˜q˜6 A q 6Ð*¨aÕ0÷	1ð Eˆ	Üœ:¨YÔ7ñ 	6Üœ"Ÿ)™) A›,×.Ñ.¨q°!Ó4Ô5÷	6ð Aˆ	Üœ:¨YÔ7ñ 	-Ü˜B ˜7 Q¨ FÐ+Ô,÷	-ð
ð 	ô œ:¨YÔ7ñ 	;Ü˜A˜q˜6 A q 6Ð*¨MÔ:÷	;ð 	;÷%	1ð 	1ú÷
	6ð 	6ú÷
	-ð 	-ú÷	;ð 	;úó/   ”CÁ/C"ÂC.Â:C:ÃCÃ"C+Ã.C7Ã:Dro  c                 ó†   — t        |«      }|j                  |j                  }}t        ||d   «       t        ||d   «       y ©Nr   r   ©r   rO   rQ   r   rµ  s         r*   Útest_edge_casesz TestBarnardExact.test_edge_cases^  s;   € ô ˜LÓ)ˆØŸM™M¨3¯:©:�6ˆ	Ü�V˜X a™[Ô)Ü�Y ¨¡Õ,r,   r�   c                 ó†   — t        |«      }|j                  |j                  }}t        ||d   «       t        ||d   «       y rÆ  rÇ  rµ  s         r*   Útest_row_or_col_zeroz%TestBarnardExact.test_row_or_col_zeroj  s;   € ô ˜LÓ)ˆØŸM™M¨3¯:©:�6ˆ	Ü�V˜X a™[Ô)Ü�Y ¨¡Õ,r,   )r¯  gEç\/?„?éÈ   é,  )gÐgÑìíQ5Àrn  é   r<   i¥  )gü&çìX}>Àrn  r°   rÂ   rÀ   c                 óÊ   — |\  }}|dk(  r"t        j                  |«      dd…ddd…f   }| }t        ||¬«      }|j                  |j                  }}t        ||g||gd¬«       y)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  sp   € ð, -5Ñ)ˆÐ)à˜)Ò#ÜŸ8™8 LÓ1²!±T°r°T°'Ñ:ˆLØ*˜NˆMä˜L°kÔBˆØŸM™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«    su  „ Ù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ð ‡[�[×ÑØà�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ò;ð4 ‡[�[×ÑØà�!ˆf�q˜!�fÐ˜xÐ(ð	
óñ-óð-ð ‡[�[×ÑØà�!ˆf�q˜"�gÐ  b§f¡f Ð.Ø�!ˆf�r˜1�gÐ  b§f¡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                  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                  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                  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                  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                  j                  d5d6«      d7„ «       Zy8)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\   re  )ç_¬VÃÑ¶?gÂñ¥…Ö­?)ç‘uÝ Ø%Å?gc'�íµ?r   r/   ©r   r   r!   )ri   g      Ö?r„   )çß+Âf°?gXc}Áv¡?é   é%   )gZÑ‹D¸É?g‰¦¢gi]Ã?c                 ó‚   — t        |d¬«      }|j                  |j                  }}t        ||g|| j                  ¬«       y)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Ÿ  s6   € ô2 ˜\°vÔ>ˆØŸM™M¨3¯:©:�6ˆ	Ü˜ FÐ+¨X¸D¿I¹IÖ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                 ó‚   — t        |d¬«      }|j                  |j                  }}t        ||g|| j                  ¬«       y)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¼  s6   € ô6 ˜\°yÔAˆØŸM™M¨3¯:©:�6ˆ	Ü˜ FÐ+¨X¸D¿I¹IÖFr,   )rá  gôqQSé,5?)rÖ  gG’Þ?/?”?)rØ  gKE`ÕÇ?)r×  gÓhr1Ö½?)râ  grÒfbÛŒ?)ri   g      æ?)rÚ  gP:pRÁv±?c                 ó„   — t        |dd¬«      }|j                  |j                  }}t        ||g|| j                  ¬«       y)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Û  s8   € ô0 ˜\°{ÀbÔIˆàŸM™M¨3¯:©:�6ˆ	Ü˜ FÐ+¨X¸D¿I¹IÖFr,   c                 ó  — d}t        t        |¬«      5  t        ddgddggd¬«       d d d «       d	}t        t        |¬«      5  t        t        j                  d
«      j                  dd«      «       d d d «       d}t        t        |¬«      5  t        ddgddgg«       d d d «       d}t        t        |¬«      5  t        ddgddggd«       d d d «       y # 1 sw Y   Œ¯xY w# 1 sw Y   ŒqxY w# 1 sw Y   ŒPxY w# 1 sw Y   y xY w)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ø  s  € ð Fð 	ô œ:¨YÔ7ñ 	2Ü˜Q ˜F Q¨ FÐ+¨qÕ1÷	2ð Eˆ	Üœ:¨YÔ7ñ 	7Üœ2Ÿ9™9 Q›<×/Ñ/°°1Ó5Ô6÷	7ð Aˆ	Üœ:¨YÔ7ñ 	.Ü˜R ˜G a¨ VÐ,Ô-÷	.ð
%ð 	ô œ:¨YÔ7ñ 	<Ü˜Q ˜F Q¨ FÐ+¨]Ô;÷	<ð 	<÷%	2ð 	2ú÷
	7ð 	7ú÷
	.ð 	.ú÷	<ð 	<úrÄ  c                 ó†   — t        |«      }|j                  |j                  }}t        ||d   «       t        ||d   «       y rÆ  )r   rO   rQ   r   rµ  s         r*   rÊ  z&TestBoschlooExact.test_row_or_col_zero  s;   € ô ˜\Ó*ˆØŸM™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kD  sJ ‚t        |d¬«      j                  }|d	k(  sJ ‚y )
Nr   r”   r„   rÀ   r¯   rÂ   r.   r¼   r�   )r   rQ   rF  )r%   ÚtblÚplÚpgÚpts        r*   Útest_two_sided_gt_1z%TestBoschlooExact.test_two_sided_gt_1  sp   € ð �1ˆv˜˜B�xÐ ˆÜ˜C¨VÔ4×;Ñ;ˆÜ˜C¨YÔ7×>Ñ>ˆØ”�R˜“‰}˜qÒ Ð Ð Ü˜C¨[Ô9×@Ñ@ˆØ�SŠyÐ‰yr,   r°   )rÀ   rÂ   c                 óŽ   — ddgddgg}t        ||¬«      j                  }t        j                  ||¬«      d   }t	        ||«       y )Nr.   r3   r4   r¯   r   )r   rO   rD  Úfisher_exactr   )r%   r°   rë  Úboschloo_statÚfisher_ps        r*   Útest_against_fisher_exactz+TestBoschlooExact.test_against_fisher_exact)  sI   € ð �1ˆv˜˜1�vÐˆÜ& s¸ÔD×NÑNˆÜ×%Ñ% c°{ÔCÀAÑFˆÜ˜ xÕ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à€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ð ‡[�[×ÑØà�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ð ‡[�[×ÑØà�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ò <ð4 ‡[�[×ÑØà�!ˆf�q˜"�gÐ §¡¨¯©Ð 0Ð1Ø�!ˆf�r˜1�gÐ §¡¨¯©Ð 0Ð1ð	
óñ-óð-òð ‡[�[×Ñ˜]Ð,?Ó@ñ1ó Añ1r,   rÕ  c                   óB  — e Zd Zej                  j                  dg  ej                  d«      f ej                  d«      dgfg«      d„ «       Zd„ Z	d„ Z
d„ Zej                  j                  dg d	¢«      d
„ «       Zej                  j                  d„ «       Zd„ Zd„ Zd„ Zy)ÚTestCvm_2sampr�   r0   r   c                 ó  — t        j                  t        t        ¬«      5  t	        |Ž }t        |j                  t        j                  «       t        |j                  t        j                  «       d d d «       y # 1 sw Y   y xY w)NrK   )
rM   rN   r   r   r   r   rO   r"   rP   rQ   rb  s      r*   Útest_too_small_inputz"TestCvm_2samp.test_too_small_input4  sW   € ô �\‰\Ô,Ô4IÔJñ 	-Ü&¨Ð-ˆCÜ˜Ÿ™¬¯©Ô/Ü˜Ÿ™¤R§V¡VÔ,÷	-÷ 	-ñ 	-ús    AA:Á:Bc                 ó¬   — t        j                  d«      }d}t        j                  t        |¬«      5  t        ||d«       d d d «       y # 1 sw Y   y xY w)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<  sC   € Ü�I‰I�a‹LˆØ?ˆÜ�]‰]œ:¨SÔ1ñ 	.Ü   A uÔ-÷	.÷ 	.ñ 	.ús   ³A
Á
Ac                 ó   — g d¢}g d¢}t        ||«      }t        t        j                  |«      t        j                  |«      «      }t        |j                  |j
                  f|j                  |j
                  f«       y )N)r.   r!   r/   r3   r1   )rØ   rš   r„   rg  )r   r"   r#   r   rO   rQ   ©r%   r&   r'   rž   rŸ   s        r*   Útest_list_inputzTestCvm_2samp.test_list_inputB  sX   € ÚˆÚˆÜ! ! QÓ'ˆÜ!¤"§(¡(¨1£+¬r¯x©x¸«{Ó;ˆÜ�b—l‘l B§I¡IÐ.°·±¸r¿y¹yÐ0IÕJr,   c                 óŒ   — g d¢}g d¢}t        ||«      }t        |j                  dd¬«       t        |j                  dd¬«       y )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  s<   € ò =ˆòˆä   AÓ&ˆÜ˜Ÿ™ U°Õ6Ü˜Ÿ™ $¨TÖ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                 ó2   — t        t        |||«      |«       y rÇ   )r   r   )r%   rO   r  rˆ   Úpvals        r*   Útest_exact_pvaluezTestCvm_2samp.test_exact_pvalueS  s   € ô 	Ô*¨9°a¸Ó;¸TÕBr,   c                 óŽ  — t         j                  j                  d«       t        j                  j                  d¬«      }t        j                  j                  d¬«      }t        ||«      }t        d|j                  cxk  xr dk  nc «       t        ||dz   «      }t        d|j                  cxk  xr
 dk  «       y c «       y )Ni  r£  r@  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ˆÜ×Ñ×"Ñ"¨Ð"Ó/ˆÜ   AÓ&ˆÜ��A—H‘HÖ ˜qÔ Ô!Ü   A c¡EÓ*ˆÜ��A—H‘HÖ ˜qÑ Õ!Ñ Õ!r,   c                 óz  — t         j                  j                  d«       t         j                  j                  d«      }t         j                  j                  d«      }t	        ||d¬«      }t	        ||d¬«      }t        |j                  |j                  «       t        |j                  |j                  d¬«       y )	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  sx   € Ü
�	‰	�‰�qÔÜ�I‰I�N‰N˜1ÓˆÜ�I‰I�N‰N˜1ÓˆÜ! ! Q¨wÔ7ˆÜ! ! Q¨|Ô<ˆÜ�R—\‘\ 2§<¡<Ô0Ü˜Ÿ	™	 2§9¡9°4Ö8r,   c                 óP  — t        j                  d«      }g d¢}t        ||d¬«      }t        ||d¬«      }t        |j                  |j                  «       t        j                  d«      }t        ||d¬«      }t        ||d¬«      }t        |j                  |j                  «       y )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‰I�b‹MˆÚˆÜ! ! Q¨wÔ7ˆÜ! ! Q¨vÔ6ˆÜ�R—Y‘Y §	¡	Ô*ä�I‰I�b‹MˆÜ! ! Q¨|Ô<ˆÜ! ! Q¨vÔ6ˆÜ�R—Y‘Y §	¡	Õ*r,   c                 óò   — t        j                  d«      }t        ||«      }t        |j                  |j
                  fd«       t        |d d |d d «      }t        |j                  |j
                  fd«       y )Nr°  rm  r/   )r"   r?   r   r   rO   rQ   rŒ   s      r*   Útest_same_inputzTestCvm_2samp.test_same_input  sc   € ô �I‰I�b‹MˆÜ" 1 aÓ(ˆÜ�c—m‘m S§Z¡ZÐ0°*Ô=ä" 1 R a 5¨!¨B¨Q¨%Ó0ˆÜ�c—m‘m S§Z¡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  s¹   „ Ø‡[�[×Ñ˜V r¨9¨2¯9©9°Q«<Ð&8Ø'0 r§y¡y°£|°a°SÐ&9ð&;ó <ñ-ó<ð-ò.òKò3ð ‡[�[×ÑÐ4ò5ó6ñ
Có6ð
Cð ‡[�[×Ññ	"ó ð	"ò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                  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                  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                  j                  dd«      d„ «       Ze	j                  j                  d g d!¢«      d"„ «       Zd#„ Zy$)%Ú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                 ó"  — t        j                  |j                  dd«      j                  «       dd t        ¬«      j                  d«      }t        j                  |Ž }|j                  «       }|D ]�  \  }}}	}
}}t        |«      dz
  t        |«      dz
  }}t        |j                  ||f   |	|¬«       t        |j                  ||f   |
|¬«       t        |j                  ||f   ||¬«       t        |j                  ||f   d	k  |dk(  «       ŒŸ y)
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…  rD  Ú	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‘Z × 6Ñ 6°u¸cÓ B× HÑ HÓ JÈ1È2Ð NÜ&+ô-ß-4©W°V«_ð 	ä—O‘O TÐ*ˆ	Ø×,Ñ,Ó.ˆà",ò 	GÑˆAˆq�!�Q˜˜3Ü�q“6˜A‘:œs 1›v¨™zˆqˆAÜ˜DŸH™H Q¨ T™N¨A°DÕ9Ü˜I×/Ñ/°°1°Ñ5°q¸tÕDÜ˜DŸI™I a¨ d™O¨Q°TÕ:Ü˜Y×-Ñ-¨a°¨dÑ3°sÑ:¸SÀA¹XÕFñ	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
        r.  rÏ  r  zunequal size samplec                 óô  — t        j                  |j                  «       t        ¬«      j	                  d«      }t        j                  |Ž }|j                  «       }|D ]™  \  }}}	}
}}t        |«      dz
  t        |«      dz
  }}t        |j                  ||f   |	|¬«       t        |j                  ||f   |
|¬«       t        |j                  ||f   ||¬«       t        |j                  ||f   ||¬«       Œ› y)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…  rD  r)  r*  r+  r   r,  rO   r-  rQ   )r%   r.  r/  r7   r0  r1  r2  r#  r3  r4  rò   r5  r)   s                r*   Útest_compare_matlabz TestTukeyHSD.test_compare_matlabï  så   € ô —Z‘Z × 4Ñ 4Ó 6Ü&+ô-ß-4©W°V«_ð 	ä—O‘O TÐ*ˆ	Ø×,Ñ,Ó.ˆà *ò 	BÑˆAˆq�!�Q˜˜1Ü�q“6˜A‘:œs 1›v¨™zˆqˆAÜ˜DŸH™H Q¨ T™N¨A°DÕ9Ü˜I×/Ñ/°°1°Ñ5°q¸tÕDÜ˜DŸI™I a¨ d™O¨Q°TÕ:Ü˜I×,Ñ,¨Q°¨TÑ2°A¸DÖAñ	Br,   c                 ó4  — d}t        j                  |j                  dd«      j                  «       dd t        ¬«      j                  d«      }g d¢g d	¢g d
¢f}t        j                  |Ž }|j                  «       }|D ]™  \  }}}}	}
}t        |«      dz
  t        |«      dz
  }}t        |j                  ||f   |	d¬«       t        |j                  ||f   |d¬«       t        |j                  ||f   |
d¬«       t        |j                  ||f   |d¬«       Œ› y)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"  r9  )é   r;   é6   re  éF   é4   é3   r<  éC   rz   rf  é   é   rB  é   é)   r\   é,   )rg  rÍ  rB  rh  r„   rg  r2  r;   é$   é*   r<  rC  r¾   r®  r<   rÍ  r®  rB  )rG  rÍ  rÛ  rg  r2   r­  r<   r°  r<  r\   rÍ  rÛ  rh  r”   r°  r°  r¾   r<   r   rÏ  r6   r‘   gñhãˆµøä>r$  )r%   Ústr_resr0  r.  r1  r2  r#  r3  rò   r4  r5  r)   s               r*   Útest_compare_rzTestTukeyHSD.test_compare_r	  s  € ðˆô —Z‘Z §¡°°sÓ ;× AÑ AÓ CÀAÀBÐ GÜ&+ô-ß-4©W°V«_ð 	ò5ò5ò5ð	6ˆô —O‘O TÐ*ˆ	Ø×,Ñ,Ó.ˆà *ò 	BÑˆAˆq�!�Q˜˜1Ü�q“6˜A‘:œs 1›v¨™zˆqˆAä˜DŸH™H Q¨ T™N¨A°DÕ9Ü˜I×/Ñ/°°1°Ñ5°q¸tÕDÜ˜DŸI™I a¨ d™O¨Q°TÕ:Ü˜I×,Ñ,¨Q°¨TÑ2°A¸DÖAñ	Br,   c                 ó”  — g d¢}g d¢}g d¢}g d¢}t        j                  ||||«      }|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 ]I  \  }	}
t        |j                  |	|
f   ||	|
f   d¬«       t        |j                  |	|
f   ||	|
f   d¬«       ŒK y)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 @rL  r   )rM  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  ru  rg   r6   N)rD  r)  r*  r"   r%  r   r,  r-  )r%   Úgroup1Úgroup2Úgroup3Úgroup4rS   r2  ÚlowerÚupperr#  r3  s              r*   Útest_engineering_stat_handbookz+TestTukeyHSD.test_engineering_stat_handbook2  sÑ   € ò
 +ˆÚ*ˆÚ+ˆÚ*ˆÜ�o‰o˜f f¨f°fÓ=ˆØ×&Ñ&Ó(ˆÜ—
‘
ÚÚ ÚÚð	ó ˆô
 —
‘
ÚÚ!ÚÚð	ó ˆð ?ò 	E‰FˆQ�Ü˜DŸH™H Q¨ T™N¨E°!°Q°$©K¸dÕCÜ˜DŸI™I a¨ d™O¨U°1°a°4©[¸tÖDñ	Er,   c                 óB  — t         j                  j                  d«       t         j                  j                  dd«      }t	        j
                  |Ž }|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«       y )Nr:   r!   r}   ©r   r   r   r   )r"   r=   r>   r¹   rD  r)  r*  r   r,  r-  r{  ÚdiagonalrO   rQ   )r%   r.  rS   r2  s       r*   Útest_rand_symmzTestTukeyHSD.test_rand_symmL  s÷   € ä
�	‰	�‰�tÔÜ�y‰y�~‰~˜a Ó%ˆÜ�o‰o˜tÐ$ˆØ×&Ñ&Ó(ˆä�T—X‘X §	¡	§¡˜|Ô,ô 	”R—[‘[ §¡Ó+¨T¯Y©Y°t©_Ô=Ü”R—[‘[ §¡Ó*¨D¯H©H°T©NÔ;ä�S—]‘] S§]¡]§_¡_Ð$4Ô5Ü”R—[‘[ §¡Ó/°Ô3ä�S—Z‘Z §¡§¡Ô.Ü”R—[‘[ §¡Ó,¨aÕ0r,   c                 ó¦   — t        t        d¬«      5  t        j                  g d¢dt        j
                  gg d¢«       d d d «       y # 1 sw Y   y xY w)Nz...must be finite.rK   rÚ   r.   )r1   r3   r!   )rW   rX   rD  r)  r"   rV   rl   s    r*   Útest_no_infzTestTukeyHSD.test_no_inf_  s:   € Üœ:Ð-AÔBñ 	?Ü�O‰OšI¨¬2¯6©6 {²IÔ>÷	?÷ 	?ñ 	?ús   ’,AÁAc                 ó’   — t        t        d¬«      5  t        j                  ddgddggddgg d¢«       d d d «       y # 1 sw Y   y xY w)Nz...must be one-dimensionalrK   r   r.   r!   r0   )r0   rÊ   r1   ©rW   rX   rD  r)  rl   s    r*   Ú
test_is_1dzTestTukeyHSD.test_is_1dc  sG   € Üœ:Ð-IÔJñ 	BÜ�O‰O˜a ˜V a¨ VÐ,¨q°!¨f²jÔA÷	B÷ 	Bñ 	Bús	   ’"=½Ac                 ó†   — t        t        d¬«      5  t        j                  g ddgg d¢«       d d d «       y # 1 sw Y   y xY w)Nz...must be greater than onerK   r.   r0   )r/   r0   r1   r]  rl   s    r*   Útest_no_emptyzTestTukeyHSD.test_no_emptyg  s6   € Üœ:Ð-JÔKñ 	3Ü�O‰O˜B  A ª	Ô2÷	3÷ 	3ñ 	3ús	   ’7·A Únargsr  c                 ó€   — t        t        d¬«      5  t        j                  g d¢g|z  Ž  d d d «       y # 1 sw Y   y xY w)Nz...more than 1 treatment.rK   ©rÊ   r3   r!   r]  )r%   ra  s     r*   Útest_not_enough_treatmentsz'TestTukeyHSD.test_not_enough_treatmentsk  s5   € äœ:Ð-HÔIñ 	5Ü�O‰Ošz˜l¨UÑ2Ñ4÷	5÷ 	5ñ 	5ús   ’4´=Úcl)rp  r   r   r.   c                 ó¬   — t        t        d¬«      5  t        j                  g d¢ddgddg«      }|j	                  |«       d d d «       y # 1 sw Y   y xY w)Nzmust be between 0 and 1rK   rc  r!   r/   rÐ   )rW   rX   rD  r)  r*  )r%   re  r  s      r*   Útest_conf_level_invalidz$TestTukeyHSD.test_conf_level_invalidp  sJ   € äœ:Ð-FÔGñ 	&Ü—‘¢
¨Q°¨F°Q¸°FÓ;ˆAØ×!Ñ! "Ô%÷	&÷ 	&ñ 	&ús   ’/A
Á
Ac                 ó  — t        j                  | j                  d d Ž }t        j                  | j                  d d Ž }t	        |j
                  |j
                  d   «       t	        |j
                  |j
                  d   «       y )Nr.   r  rN  )rD  r)  Údata_diff_sizeÚ	ttest_indr   rQ   )r%   r1  Ú	res_ttests      r*   Útest_2_args_ttestzTestTukeyHSD.test_2_args_ttestv  sn   € ä—O‘O T×%8Ñ%8¸¸!Ð%<Ð=ˆ	Ü—O‘O T×%8Ñ%8¸¸!Ð%<Ð=ˆ	Ü˜	×(Ñ(¨)×*:Ñ*:¸4Ñ*@ÔAÜ˜	×(Ñ(¨)×*:Ñ*:¸4Ñ*@ÕAr,   N)r_   r`   ra   Údata_same_sizeri  Úextreme_sizeÚsas_same_sizeÚsas_diff_sizeÚsas_extremerM   r¡   r¢   r7  Úmatlab_sm_sizÚmatlab_diff_szr:  rJ  rU  rY  r[  r^  r`  rd  rg  rl  rb   r,   r*   r  r  ‹  sy  „ â4Ú4Ú4ð6€Nò HÚ4Ú4ð6€Nò 'òâ2ð4€Lð
€Mð€Mð€Kð ‡[�[×ÑÐ7Ø-¨}¸dÐCØ-¨}¸dÐCØ+¨[¸%Ð@ð ò"Eð ó Fñ#GóFð#GðJ€Mð€Nð ‡[�[×ÑÐ7Ø-¨}¸eÐDØ-¨~¸tÐDðFà"5Ø"7ð"9ð ó :ñ
Bó:ð
Bò*'BòREò41ò&?òBò3ð ‡[�[×Ñ˜W fÓ-ñ5ó .ð5ð ‡[�[×Ñ˜T¢>Ó2ñ&ó 3ð&ó
Br,   r  c                   óâ   — e Zd Zej                  j                  dg d¢g d¢f«      d„ «       Zej                  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
y)ÚTestPoissonMeansTestzc1, n1, c2, n2, p_expect)r   r}   r!   r}   gþe÷äa¡¶?)r.   r}   r1   r}   gÞ	ŠcÆ?c                 ód   — t        j                  ||||«      }t        |j                  |d¬«       y )Nrj   r6   ©rD  Úpoisson_means_testr   rQ   )r%   Úc1Ún1Úc2Ún2Úp_expectrS   s          r*   Útest_paper_examplesz(TestPoissonMeansTest.test_paper_examples  s*   € ô ×&Ñ& r¨2¨r°2Ó6ˆÜ˜Ÿ
™
 H°4Ö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                 ól   — t        j                  ||||||¬«      }t        |j                  |dd¬«       y )N)r°   Údiffg�íµ ÷ÆÀ>g¼‰Ø—²Òœ<©r7   r<  rw  )	r%   ry  rz  r{  r|  r}  ÚaltÚdrS   s	            r*   Útest_fortran_authorsz)TestPoissonMeansTest.test_fortran_authorsˆ  s0   € ô$ ×&Ñ& r¨2¨r°2À3ÈQÔOˆÜ˜Ÿ
™
 H°4¸eÖDr,   c                 ót   — d\  }}d\  }}t        j                  ||||«      }t        |j                  d«       y )N)r~   r~   r   rw  ©r%   Úcount1Úcount2Únobs1Únobs2rS   s         r*   Útest_different_resultsz+TestPoissonMeansTest.test_different_results�  s:   € ð &‰ˆ�Ø#‰ˆˆuÜ×&Ñ& v¨u°f¸eÓDˆÜ˜Ÿ
™
 AÕ&r,   c                 ót   — d\  }}d\  }}t        j                  ||||«      }t        |j                  d«       y )NrW  rÙ  r   rw  r‡  s         r*   Útest_less_than_zero_lambda_hat2z4TestPoissonMeansTest.test_less_than_zero_lambda_hat2¦  s:   € ð ‰ˆ�Ø‰ˆˆuÜ×&Ñ& v¨u°f¸eÓDˆÜ˜Ÿ
™
 AÕ&r,   c                 ó  — d\  }}d\  }}d}t        t        |¬«      5  t        j                  d|||«       d d d «       t        t        |¬«      5  t        j                  ||d|«       d d d «       d}t        t        |¬«      5  t        j                  d|||«       d d d «       t        t        |¬«      5  t        j                  ||d|«       d d d «       d}t        t        |¬«      5  t        j                  |d||«       d d d «       t        t        |¬«      5  t        j                  |||d«       d d d «       d	}t        t        |¬«      5  t        j                  ||||d¬
«       d d d «       d}t        t        |¬«      5  t        j                  ddddd¬«       d d d «       y # 1 sw Y   �ŒuxY w# 1 sw Y   �ŒPxY w# 1 sw Y   �Œ)xY w# 1 sw Y   �ŒxY w# 1 sw Y   ŒÜxY w# 1 sw Y   Œ¶xY w# 1 sw Y   ŒŒxY w# 1 sw Y   y xY w)NrW  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   Ú	TypeErrorrD  rx  rX   )r%   rˆ  r‰  rŠ  r‹  r�  s         r*   rµ   z*TestPoissonMeansTest.test_input_validation®  sê  € Ø‰ˆ�Ø‰ˆˆuð 4ˆÜœ9¨GÔ4ñ 	?Ü×$Ñ$ R¨°¸Ô>÷	?äœ9¨GÔ4ñ 	?Ü×$Ñ$ V¨U°B¸Ô>÷	?ð FˆÜœ:¨WÔ5ñ 	?Ü×$Ñ$ R¨°¸Ô>÷	?äœ:¨WÔ5ñ 	?Ü×$Ñ$ V¨U°B¸Ô>÷	?ð :ˆÜœ:¨WÔ5ñ 	@Ü×$Ñ$ V¨R°¸Ô?÷	@äœ:¨WÔ5ñ 	@Ü×$Ñ$ V¨U°F¸BÔ?÷	@ð =ˆÜœ:¨WÔ5ñ 	LÜ×$Ñ$ V¨U°F¸EÈÕK÷	Lð 3ˆÜœ:¨WÔ5ñ 	FÜ×$Ñ$ Q¨¨1¨a¸WÕE÷	Fð 	F÷5	?ñ 	?ú÷	?ñ 	?ú÷
	?ñ 	?ú÷	?ñ 	?ú÷
	@ð 	@ú÷	@ð 	@ú÷
	Lð 	Lú÷
	Fð 	Fús_   žF*ÁF7ÂGÂ6GÃ*GÄG*ÅG6ÆHÆ*F4Æ7GÇGÇGÇG'Ç*G3Ç6G?ÈHN)r_   r`   ra   rM   r¡   r¢   r~  r…  rŒ  rŽ  rµ   rb   r,   r*   ru  ru  ~  s‹   „ Ø‡[�[×ÑÐ7â Ú ð:ó ñ
9óð
9ð ‡[�[×ÑÐ?ò
 	=Ú<Ú=Ú>Ú9Ú;ò 	6Ú6ðBó ñ"Eó#ð"Eò'ò'ó!Fr,   ru  c                   ó´   — e Zd Zd„ Zd„ Zej                  j                  dg d¢«      d„ «       Zej                  j                  dg d¢«      d„ «       Z	d„ Z
d	„ Zy
)ÚTestBWSTestc                 ól  — t         j                  j                  d«      }|j                  d¬«      \  }}d}t        j                  t
        |¬«      5  t        j                  ||g||g«       d d d «       d}t        j                  t
        |¬«      5  t        j                  t         j                  g|«       d d d «       d}t        j                  t
        |¬«      5  t        j                  |g «       d d d «       d}t        j                  t
        |¬«      5  t        j                  ||d	¬
«       d d d «       d}t        j                  t
        |¬«      5  t        j                  ||d¬«       d d d «       y # 1 sw Y   �ŒxY w# 1 sw Y   ŒÏxY w# 1 sw Y   ŒŸxY w# 1 sw Y   ŒmxY w# 1 sw Y   y xY w)Nì   <oôvêT‘{ ©r.   r3   r@  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...rH  r³   )	r"   r=   rB  rM   r   rX   rD  Úbws_testrP   )r%   rG  r&   r'   r�  s        r*   Útest_bws_input_validationz%TestBWSTest.test_bws_input_validationÔ  sX  € Ü�i‰i×#Ñ#Ð$7Ó8ˆà�z‰z˜vˆzÓ&‰ˆˆ1à@ˆÜ�]‰]œ:¨WÔ5ñ 	+Ü�N‰N˜A˜q˜6 A q 6Ô*÷	+ð 7ˆÜ�]‰]œ:¨WÔ5ñ 	(Ü�N‰NœBŸF™F˜8 QÔ'÷	(ð 9ˆÜ�]‰]œ:¨WÔ5ñ 	"Ü�N‰N˜1˜bÔ!÷	"ð 3ˆÜ�]‰]œ:¨WÔ5ñ 	:Ü�N‰N˜1˜a¨[Õ9÷	:ð 6ˆÜ�]‰]œ:¨WÔ5ñ 	,Ü�N‰N˜1˜a¨Õ+÷	,ð 	,÷!	+ñ 	+ú÷	(ð 	(ú÷	"ð 	"ú÷	:ð 	:ú÷	,ð 	,ús<   ÁE9Â&FÃFÄFÅF*Å9FÆFÆFÆF'Æ*F3c                 ó    — g d¢}g d¢}t        j                  ||d¬«      }t        |j                  dd¬«       t	        |j
                  d«       y )	N)r   r.   r!   r/   r1   r3   r4   )r0   rÐ   r2   r“   r„   r”   rf  r¼   r¯   gºI+‡@r5   r6   g¤f$/•Þg?)rD  r—  r   rO   r   rQ   rR   s       r*   Ú test_against_published_referencez,TestBWSTest.test_against_published_referenceî  s>   € ò "ˆÚ&ˆÜ�n‰n˜Q ¨{Ô;ˆÜ˜Ÿ™ u°4Õ8Ü�S—Z‘Z Õ)r,   )r°   rO   rQ   ))r¼   g
»ù-ü?g4ˆñ³B/À?)rÀ   ç
»ù-ü¿g«î0&Àv­?)rÂ   r›  gÉñœý“(î?c                 óþ   — t         j                  j                  d«      }|j                  d¬«      \  }}t        j                  |||¬«      }t        |j                  |d¬«       t        |j                  |dd¬	«       y )
Nr•  r–  r@  r¯   ç‚vIhÂ%<=r;  rg   r�   r‚  ©r"   r=   rB  rD  r—  r   rO   rQ   ©r%   r°   rO   rQ   rG  r&   r'   rS   s           r*   Útest_against_RzTestBWSTest.test_against_Rø  sa   € ô �i‰i×#Ñ#Ð$7Ó8ˆØ�z‰z˜vˆzÓ&‰ˆˆ1Ü�n‰n˜Q ¨{Ô;ˆÜ˜Ÿ™ y°uÕ=Ü˜Ÿ
™
 F°¸DÖAr,   ))r¼   góD 5Hò?g«ÒdüÔ•Ò?)rÀ   ç`�Ñ‡©áï?gÈ²µ×šXë?)rÂ   r¡  gß4)¡”�Â?c                 ó  — t         j                  j                  d«      }|j                  d¬«      }|j                  d¬«      }t        j                  |||¬«      }t        |j                  |d¬«       t        |j                  |dd	¬
«       y )Nl   ‚.ÅÇs½Z rÐ   r@  r4   r¯   r�  r;  rg   r�   r‚  rž  rŸ  s           r*   Útest_against_R_imbalancedz%TestBWSTest.test_against_R_imbalanced	  sm   € ô �i‰i×#Ñ#Ð$7Ó8ˆØ�J‰J˜AˆJÓˆØ�J‰J˜AˆJÓˆÜ�n‰n˜Q ¨{Ô;ˆÜ˜Ÿ™ y°uÕ=Ü˜Ÿ
™
 F°¸DÖAr,   c                 ó  — t         j                  j                  d«      }|j                  d¬«      \  }}t         j                  j                  d«      }t        j                  d|¬«      }t        j
                  |||¬«      }t        |j                  «      dk(  sJ ‚t         j                  j                  d«      }t        j                  d|¬«      }t        j
                  |||¬«      }t        |j                  |j                  «       t         j                  j                  d«      }t        j                  d|¬«      }t        j
                  |||¬«      }t        j                  |j                  |j                  «      rJ ‚y )Nì   /°HøNÏ( )r.   r2   r@  r2   )Ún_resamplesrG  r³   l   ÝVC£	ãA )
r"   r=   rB  rD  rK  r—  rû   Únull_distributionr   Úallclose)r%   rG  r&   r'   r´   r	  r
  r*  s           r*   Útest_methodzTestBWSTest.test_method  s(  € ä�i‰i×#Ñ#Ð$7Ó8ˆØ�z‰z˜wˆzÓ'‰ˆˆ1ä�i‰i×#Ñ#Ð$7Ó8ˆÜ×(Ñ(°R¸SÔAˆÜ�~‰~˜a ¨6Ô2ˆä�4×)Ñ)Ó*¨bÒ0Ð0Ð0ä�i‰i×#Ñ#Ð$7Ó8ˆÜ×(Ñ(°R¸SÔAˆÜ�~‰~˜a ¨6Ô2ˆä˜×.Ñ.°×0FÑ0FÔGä�i‰i×#Ñ#Ð$7Ó8ˆÜ×(Ñ(°R¸SÔAˆÜ�~‰~˜a ¨6Ô2ˆä—;‘;˜t×5Ñ5°t×7MÑ7MÔNÐNÐNÐNr,   c                 óÀ  — t         j                  j                  d«      }|j                  d¬«      }|dz
  }t        j                  ||d¬«      }|j
                  dkD  sJ ‚t        |j                  dt        |j                  «      z  «       t        j                  ||d¬«      }|j
                  dkD  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«       y )	Nr¥  r0   r@  r   rÂ   r¯   r   rÀ   )
r"   r=   rB  rD  r—  rO   r   rQ   rû   r§  )r%   rG  r&   r'   rS   s        r*   Útest_directionszTestBWSTest.test_directions2  s  € ä�i‰i×#Ñ#Ð$7Ó8ˆØ�J‰J˜AˆJÓˆØ�‰Eˆä�n‰n˜Q ¨yÔ9ˆØ�}‰}˜qÒ Ð Ð Ü�S—Z‘Z ¤S¨×)>Ñ)>Ó%?Ñ!?Ô@ä�n‰n˜Q ¨vÔ6ˆØ�}‰}˜qÒ Ð Ð Ü�S—Z‘Z Ô#ä�n‰n˜Q ¨vÔ6ˆØ�}‰}˜qÒ Ð Ð Ü�S—Z‘Z ¤S¨×)>Ñ)>Ó%?Ñ!?Ô@ä�n‰n˜Q ¨yÔ9ˆØ�}‰}˜qÒ Ð Ð Ü�S—Z‘Z Õ#r,   N)r_   r`   ra   r˜  rš  rM   r¡   r¢   r   r£  r©  r«  rb   r,   r*   r“  r“  Ò  sw   „ ò,ò4*ð ‡[�[×ÑÐCòNóOñBó	OðBð ‡[�[×ÑÐCòNóOñBó	OðBòOó.$r,   r“  ).Ú	itertoolsr   Únumpyr"   r=   rZ  rM   Únumpy.testingr   r   r   r   r   rW   Úscipy.statsrD  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¤   rT  r«  rÕ  rö  r  ru  r“  rb   r,   r*   ú<module>rµ     sÐ   ðÝ ã Û Û Û ÷0ó 0å *å Ý %÷4÷ 4ñ 4÷ EÑ DÝ -Ý 2ß R÷4-ñ 4-÷nZKñ ZK÷zl=ñ l=ô^j<Ð"ô j<÷ZP
ñ P
÷fV1ñ V1÷rU>ñ U>÷ppBñ pB÷fQFñ QF÷ht$ò t$r,   