o
    h                    @   s  d dl Z d dlZd dlZd dlmZ d dlmZmZm	Z	m
Z
mZmZmZmZ d dlZd dlmZ d dlZd dlmZ d dlmZ d dlmZmZmZ dd	lmZ d
dlmZmZ d dlm Z  d dl!m"Z" e#e"Z"zd dl$Z$de$j%d< d dl&m'Z( dZ)W n e*y   dZ)Y nw g dZ+g dZ,g dZ-g dZ.g dZ/g dZ0g dZ1g dZ2g dZ3g dZ4dd Z5G dd dZ6G d d! d!Z7G d"d# d#Z8G d$d% d%Z9G d&d' d'Z:G d(d) d)Z;G d*d+ d+Z<G d,d- d-Z=G d.d/ d/Z>G d0d1 d1Z?G d2d3 d3Z@d4d5 ZAG d6d7 d7ZBG d8d9 d9ZCG d:d; d;ZDG d<d= d=ZEG d>d? d?ZFG d@dA dAZGG dBdC dCZHG dDdE dEZIg dFZJG dGdH dHZKG dIdJ dJZLG dKdL dLZMG dMdN dNZNG dOdP dPZOG dQdR dRZPG dSdT dTZQG dUdV dVZRG dWdX dXZSG dYdZ dZZTdS )[    N)RandomState)assert_array_equalassert_almost_equalassert_array_lessassert_array_almost_equalassert_assert_allcloseassert_equalsuppress_warnings)raises)optimize)stats)
_abw_state_get_As_weibull_Avals_weibull   )check_named_results   )_get_wilcoxon_distr_get_wilcoxon_distr2)_binary_search_for_binom_tst)distcontAggbackendTF)
jt?Zd;?V-?      ?gX9v?-?x&1?+?+?r   )
r   r   r   r   灕Cl?r   r   r   r   V-?)
ʡE?g/$?r"   r    ףp=
?r!   r   r    r   r   )
gGz?r   r!   r   r%   r!   r   r   r   r   )
r   r   gCl?Gz?r   rh|?r   r   g\(\?r   )
羟/$?gh|?5?r(   r"   r#   r   r%   r   g1Zd?r   )
r&   Mb?r   gjt?r   r   g}?5^I?g)\(?r   r   )
r   r   r   r   r   r   r   r   r   r   )
r   r   r   r%   r   r)   r)   r   r    r$   )
r$   r%   r'   r!   r"   r"   r$   r   r)   r"   c                  O   s   t tjj| i |S N)nplogr   gammarvs)argskwargs r1   [/var/www/html/scripts/venv/lib/python3.10/site-packages/scipy/stats/tests/test_morestats.py_old_loggamma_rvs6      r3   c                   @   $   e Zd Zdd Zdd Zdd ZdS )TestBayes_mvsc                 C   st   g d}t |\}}}t|jd t|jddd t|jd t|jddd t|jd	d
d t|jddd d S )N   	            r<            "@g4`H#j@gO[%@+=rtol      $@gai	@gӇu8@&.>_紫@   decimalgYs?g1O@)r   	bayes_mvsr   	statisticr   minmaxselfdatameanvarstdr1   r1   r2   
test_basic;   s   
zTestBayes_mvs.test_basicc                 C      t ttjg  d S r*   )assert_raises
ValueErrorr   rK   rO   r1   r1   r2   test_empty_inputM      zTestBayes_mvs.test_empty_inputc                 C   s0   t d}d}t|}|D ]}t|| qd S )N   )rL   rM   )r+   aranger   rK   r   )rO   x
attributesresir1   r1   r2   test_result_attributesP   s   

z$TestBayes_mvs.test_result_attributesN)__name__
__module____qualname__rT   rY   ra   r1   r1   r1   r2   r6   :   s    r6   c                   @   ,   e Zd Zdd Zdd Zdd Zdd Zd	S )
TestMvsdistc                 C   s   g d}t |\}}}t| d t|dddd t| d t|ddd	d t| d
dd t|dddd d S )Nr7   r>   ?r?   r@   rA   rC   rD   rE   rF   rG   rH   rJ   )r   mvsdistr   rQ   r   intervalrN   r1   r1   r2   rT   Z   s   
zTestMvsdist.test_basicc                 C   rU   r*   rV   rW   r   rh   rX   r1   r1   r2   rY   i   rZ   zTestMvsdist.test_empty_inputc                 C      dg}t ttj| d S Nr   rj   rO   rP   r1   r1   r2   test_bad_argl      zTestMvsdist.test_bad_argc                 C   sj   t  ' t dt dd tg dD  dd tg dD  W d    d S 1 s.w   Y  d S )Nerrorc                 S      g | ]}|  qS r1   rQ   .0r]   r1   r1   r2   
<listcomp>v       z*TestMvsdist.test_warns.<locals>.<listcomp>r   r      c                 S   rq   r1   rr   rs   r1   r1   r2   ru   w   rv   r   r   rx         )warningscatch_warningssimplefilterRuntimeWarningr   rh   rX   r1   r1   r2   
test_warnsq   s
   
"zTestMvsdist.test_warnsN)rb   rc   rd   rT   rY   rn   r   r1   r1   r1   r2   rf   Y   s
    rf   c                   @   sL   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dS )TestShapiroc           
      C   s  g d}t |\}}t |}t|ddd t|jddd t|ddd t|jddd g d}t |\}}t |}t|ddd t|jddd t|dd	d t|jdd	d t jjd
d	ddd}t |\}}t |}t|ddd t|jddd t|dd	d t|jdd	d g d}d}d}	t |\}}t |}t||dd t|j|dd t||	d
d t|j|	d
d d S )N))\(?{Gz@q=
ףp@HzG$@@Q	@q=
ףp?Q@zG?      @Q?RQ?Q	@(\?zG?zG@(\#@333333@Gz?p=
ף@   ?r8   rH      Ɍ?)(\?=
ףp=?\(\@ffffff@\(\?(\?Gz@(\ףp=
@皙?)\(?RQr   p=
ף?@
ףp=
?{Gz?\(\@{Gz@Q@IddY?Dio?rx   r{   d   iNa locscalesizerandom_stateg   E?g   Oٴ?)gn?gjt?gffffff?gMb?gjt?g;On?gjt?gX9v?gS?g(\?gI+?gʡE?gx?g|?5^?r   gCl?g1Zd@gS@g(\	@gGz@g rh@g"~j@gQ@gtV@gʡE @ghyܝ?g0M?rz   )r   shapiror   rL   pvaluenormr.   )
rO   x1wpwshapiro_testx2x3x4
W_expected
p_expectedr1   r1   r2   rT   {   s<   



zTestShapiro.test_basicc                 C   s   g dg dg}t |\}}t |}t|ddd t|jddd t|ddd t|jddd g dg dg}t |\}}t |}t|d	dd t|jd	dd t|d
dd t|jd
dd d S )N)
r   r   r   r   r   r   r   r   r   r   )
r   r   r   r   r   r   r   r   r   r   r   r8   rH   r   )
r   r   r   r   r   r   r   r   r   r   )
r   r   r   r   r   r   r   r   r   r   r   r   rx   )r   r   r   rL   r   )rO   r   r   r   r   r   r1   r1   r2   test_2d   s$   

zTestShapiro.test_2dc                 C   s&   t ttjg  t ttjg g g g d S r*   rV   rW   r   r   rX   r1   r1   r2   rY         zTestShapiro.test_empty_inputc                 C   s4   t ttjddg t ttjtjg dggtd d S )Nr   r   dtype)rV   rW   r   r   r+   arrayobjectrX   r1   r1   r2   test_not_enough_values   s   "z"TestShapiro.test_not_enough_valuesc                 C   rk   rl   r   rO   r]   r1   r1   r2   rn      ro   zTestShapiro.test_bad_argc                 C   s`   t d}t j|d< t|\}}t|}t|t j t|jt j t|d t|jd d S )NrC   r9   r   )	r+   r\   nanr   r   r	   rL   r   r   )rO   r]   r   r   r   r1   r1   r2   test_nan_input   s   



zTestShapiro.test_nan_inputc                 C   s2   t g d\}}t |}d}t||dd d S )N)i i = i@ )gPq~?gF?h㈵>rA   )r   boxcoxr   r   )rO   	trans_valmaxlogr_   refr1   r1   r2   test_gh14462   s   
zTestShapiro.test_gh14462c                 C   sN   t g d}|jdksJ g d}t |}t|jd t|jddd d S )N)g9B.?        r   r   )g1gIۿg0!?gr??gel?ư>rA   )r   r   r   r   rL   )rO   r_   r]   r1   r1   r2   test_length_3_gh18322   s   
z!TestShapiro.test_length_3_gh18322N)rb   rc   rd   rT   r   rY   r   rn   r   r   r   r1   r1   r1   r2   r   z   s    +r   c                   @   s~   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zejdg ddd Zdd ZdS )TestAndersonc                 C   s   t d}|jdd}|jdd}t|\}}}t|d d | t|\}}}t||dd   td}d|d< t|\}}}t|d d S )	NI2   r   
   r   g	@)	r   standard_exponentialstandard_normalr   andersonr   r+   onesr   )rO   rsr   r   Acritsigvr1   r1   r2   test_normal   s   

zTestAnderson.test_normalc                 C   s   t d}|jdd}|jdd}t|d\}}}t||dd   tjdd t|d\}}}W d    n1 s<w   Y  t||d k d S )	Nr   r   r   exponr   ignore)allr   )	r   r   r   r   r   r   r+   errstater   )rO   r   r   r   r   r   r   r1   r1   r2   
test_expon  s   zTestAnderson.test_exponc              	   C   s   t d}d|d< t|d\}}}t|}tj|\}}tj|||}tj|||}	t 	d|d }
| t 
d|
 d ||	d d d    }t|| d S )Nr   r   r   gumbelr   r   r   )r+   r   r   r   lengumbel_lfitlogcdflogsfr\   rQ   r   )rO   r   a2r   r   nxbarsr   r   r`   expected_a2r1   r1   r2   test_gumbel  s   
*zTestAnderson.test_gumbelc                 C      t ttjdgdd d S Nr   plate_of_shrimpdist)rV   rW   r   r   rX   r1   r1   r2   rn   '  r4   zTestAnderson.test_bad_argc                 C   s0   t d}|jdd}t|}d}t|| d S )Nr   r   r   rL   critical_valuessignificance_level)r   r   r   r   r   )rO   r   r]   r_   r^   r1   r1   r2   ra   *  s
   
z#TestAnderson.test_result_attributesc           	      C   sF   t d}|jdd}t|d\}}}t|d\}}}t|| d S )Nr   r   r   r   r   )r   r   r   r   r   )	rO   r   r]   A1crit1sig1A2crit2sig2r1   r1   r2   test_gumbel_l1  s
   zTestAnderson.test_gumbel_lc           
      C   sp   t d}|jdd}td}d|d< t|d\}}}t|d\}}}	t||dd   t||d k d S )	Nr   r   r   r   r   gumbel_rr   r   )r   r   r+   r   r   r   r   r   )
rO   r   r   r   r   r   r   r   r  r  r1   r1   r2   test_gumbel_r;  s   
zTestAnderson.test_gumbel_rc           
      C   s   t g d}t|d}|jj\}}}t|||fddd t|jddd |j|jd k s0J d	| }t|d
dd t	d }t	d }||d d ||   }	t 
|j|	ksYJ t|j|	dd d S )N)
            r	        r
  u   r  weibull_min)g
ףp=
@gzGX@gQS@Mb`?rA   gp=
ף?MbP?r   r   g?r   皙?g?atol)r+   r   r   r   
fit_resultparamsr   rL   r   r   r   )
rO   r]   r_   mr   r   cAs40As45As_refr1   r1   r2   test_weibull_min_case_AJ  s   z$TestAnderson.test_weibull_min_case_Ac                 C   sR   t g d}d}tjt|d t|d W d    d S 1 s"w   Y  d S )N)J   9   0        r:   F      r    ;            iF  z/Maximum likelihood estimation has converged to matchr  )r+   r   pytestr   rW   r   r   )rO   r]   messager1   r1   r2   test_weibull_min_case_B]  s
   "z$TestAnderson.test_weibull_min_case_Bc              	   C   s   t g d }d}d}tjt|d}tjt|d}|) | t|d W d    n1 s0w   Y  W d    d S W d    d S 1 sHw   Y  d S )N)	r  K   r     k   r:   =   +   r  z:Critical values of the test statistic are given for the...z;An error occurred while fitting the Weibull distribution...r(  r  )	r+   r   r*  warnsUserWarningr   rW   r   r   )rO   r]   wmessageemessagewcontextecontextr1   r1   r2   test_weibull_warning_errore  s   Pz'TestAnderson.test_weibull_warning_errordistname)r   r   r   extreme1r   r  logisticr  c                 C   s\   t jd}|dv rdn|}tt|}t| }|j|d|d}t||}|jj	s,J d S )Nl   'Dah$>   r   r:  r     r   r   )
r+   randomdefault_rnggetattrr   r   r.   r   r  success)rO   r9  rngreal_distnamer   r  r]   r_   r1   r1   r2   test_anderson_fit_paramsp  s   
z%TestAnderson.test_anderson_fit_paramsc                 C   s:   d}t td| td  tj}t td| td  d S )Nr   r   r   )r	   r   r   r+   infrO   r  r1   r1   r2   test_anderson_weibull_As~  s   z%TestAnderson.test_anderson_weibull_AsN)rb   rc   rd   r   r   r   rn   ra   r  r  r  r,  r8  r*  markparametrizerD  rG  r1   r1   r1   r2   r      s     


r   c                   @   T   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd ZdS )TestAndersonKSampc                 C   s   g d}t g d}t g d}t g d}tj||||fdd\}}}t|dd tg d	|d
d d t|ddd d S )NgYC@g     D@gfffffE@g     @F@g     F@g      G@gG@g      M@gC@gfffffC@gC@g33333D@gfffffD@g33333sE@gfffffE@gfffffF@      A@g     A@g     C@g      D@     E@rP        F@g     F@rO  ffffffA@rS  g33333A@gB@gfffffB@gD@gffffffE@Fmidrankg@rx   gCl?gvq-?g,Ԛ?gl@gPk	@r   r{   rz   gJ4a?Mb0?r  r+   r   r   anderson_ksampr   r   r   rO   t1t2t3t4Tktmpr1   r1   r2   test_example1a  s   z TestAndersonKSamp.test_example1ac                 C   s   t g d}t g d}t g d}t g d}tj||||fdd\}}}t|dd tg d	|d
d d t|ddd d S )NrL  rM  rN  rR  TrT  gQ@rx   rV  r   r{   rz   r  rW  r  rX  rZ  r1   r1   r2   test_example1b  s   z TestAndersonKSamp.test_example1bc                 C      g d}g d}g d}g d}g d}g d}g d}g d}g d	}	g d
}
g d}g d}g d}g d}t j|||||||||	|
||||fdd\}}}t|dd tg d|dd d t|ddd d S )N   r[   )   r  !      i  rG   :   %   r   A   r9      i     $      v   
"         rv  C   r  >   r;      rt  Z   r   <      r0  1   rG      8      O   T   ,   r$  r  rr        i6  L   r'  r     rx        r!  e   r  r  r  r  r  r   r:   r!  r"  r  r#  r$  r%  7   i@  r  h         /            rh  r  i  W   r;   x   rG   rx  r  r  G   r  r"  *   r  r{   r:   r     rx   rG   r  r  rG   r     r{  r   r  4   _   a   3   r  rz      rv     D   M   P   r   r  j      R   6   ru     .   o   '   ?   rv     rv     r  r   r  f   H   ry  r  rx   r[         r  X   r  r{   r{   rp  ry        r     r  r=   rG   	g  r9   r:   i  i[  rx   r  r   i  r      r{   i  #   r:   i  rv  r   r;   b   r{   U   [   r1     rx   r  r     rG   r  r      rw  r$        r%  rG   r  B   r0  rt  FrT  gM
@rx   g^I+?g~:?g+e?gH.!@gcZB>h@r   r{   rz   	^)p?rW  r  r   rY  r   r   r   rO   r[  r\  r]  r^  t5t6t7t8t9t10t11t12t13t14r_  r`  ra  r1   r1   r2   test_example2a  s0   z TestAndersonKSamp.test_example2ac                 C   rd  )Nre  rj  rs  rz  r  r  r  r  r  r  r  r  r  r  TrT  g'1Z
@rx   r  r   r{   rz   r  rW  r  r  r  r1   r1   r2   test_example2b  s0   z TestAndersonKSamp.test_example2bc                 C   s~  t ddd}t }|jtdd tj||d gdd\}}}W d    n1 s*w   Y  t|dd	 t|d
 t }|jtdd t||d g\}}}W d    n1 s[w   Y  t|dd	 t|d
 t }|jtdd tj||d gdd\}}}W d    n1 sw   Y  t|dd t|d t }|jtdd t||d g\}}}W d    n1 sw   Y  t|dd t|d tj||d gdd\}}}t|dd t	|dddd t||d g\}}}t|dd t	|dddd tj||d gdd\}}}t|dd t	|dddd tj||d gdd\}}}t|d d t	|d!d"dd d S )#Nr   r   zp-value flooredr+  g     @D@FrT  g=
ףpD@rx   r  gGzD@zp-value capped      ?g?ܵrz         ?g?ܵ      @g!u?gףp=
׳?{Gzt?r   r  rB   r8   gBiq?g鷯?      '@g:M@g46<Rv?gMb@?g      +@g6>W[1@g1ZGUS?-C6?)
r+   linspacer
   filterr3  r   rY  r   r	   r   )rO   r   supr   _ra  r1   r1   r2   test_R_kSamples  sJ   !



z!TestAndersonKSamp.test_R_kSamplesc                 C   s   t ttjtd d S Nr{   rV   rW   r   rY  r+   r   rX   r1   r1   r2   test_not_enough_samples7  r4   z)TestAndersonKSamp.test_not_enough_samplesc                 C   s"   t ttjtdtdf d S r  r  rX   r1   r1   r2   test_no_distinct_observations:  s   z/TestAndersonKSamp.test_no_distinct_observationsc                 C   s   t ttjtdg f d S r  r  rX   r1   r1   r2   test_empty_sample>     z#TestAndersonKSamp.test_empty_samplec                 C   sH   g d}t g d}tj||fdd}d}t|| t|j|j d S )NrL  rM  FrT  r   )r+   r   r   rY  r   r	   r   r   )rO   r[  r\  r_   r^   r1   r1   r2   ra   A  s   
z(TestAndersonKSamp.test_result_attributesN)rb   rc   rd   rb  rc  r  r  r  r  r  r  ra   r1   r1   r1   r2   rK    s      PrK  c                   @   s   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	e
jdg dg dddfg dg dddfg dg dddfg dg dddfg dg dddfg dg dddfgdd Zdd ZdS ) 
TestAnsaric                 C   sn   g d}g d}t  }|td t||\}}W d    n1 s$w   Y  t|dd t|dd d S )Nr   r   rx   rx   rz   rx   r   r8   r   r8   r   rz   r   %Ties preclude use of exact statistic.g     7@r  g8#oG?)r
   r  r3  r   ansarir   )rO   r]   yr  Wpvalr1   r1   r2   
test_smallO  s   zTestAnsari.test_smallc                 C   sr   t d}t d}t }|td t||\}}W d    n1 s&w   Y  t|dd t|dd d S )Nr  r/  r   c   r  r  m   l   r  r  r  `   r  r  r/  q   t   r  n   r  r/  r  r  r  i   g   r  r  r  r   r  r  r  r  r   r  r  r  r  r  r  g     0g@r  g|:?)r+   r   r
   r  r3  r   r  r   )rO   ramsayparekhr  r  r  r1   r1   r2   test_approxX  s   

zTestAnsari.test_approxc                 C   s4   t g dg d\}}t|dd t|dd d S )Nr   r   rx   rz   )r[   r{   r  r<   r   r:   rC   r  g?r;   )r   r  r   )rO   r  r  r1   r1   r2   
test_exactf  s   zTestAnsari.test_exactc                 C   s(   t ttjg dg t ttjdgg  d S rl   rV   rW   r   r  rX   r1   r1   r2   rn   k  s   zTestAnsari.test_bad_argc                 C   s`   g d}g d}t  }|td t||}W d    n1 s"w   Y  d}t|| d S )Nr  r  r  rL   r   )r
   r  r3  r   r  r   )rO   r]   r  r  r_   r^   r1   r1   r2   ra   o  s   z!TestAnsari.test_result_attributesc                 C   sV   g d}g d}d}t t|d tj||dd W d    d S 1 s$w   Y  d S )Nr  r{   r8   r;   r<   z!'alternative' must be 'two-sided'r(  fooalternativer	  )rO   r   r   r)  r1   r1   r2   test_bad_alternativex  s   "zTestAnsari.test_bad_alternativec           
      C   s   g d}g d}t ||\}}t j||ddj}t j||ddj}|dks(J |dk s.J t|t|t|}t|| d| d	d
 t||d d	d
 t|d| |d  d	d
 t j||ddj}t j||ddj}	|dk srJ |	dksxJ d S )N)r   r{   r   r[   r  r  )r  g      !@      #@g      %@r  g      )@lessr  greaterffffff?皙?r   -q=r  r   )r   r  r   r   pmfr   r   )
rO   r   r   rL   r  pval_lpval_gprobpval_l_reversepval_g_reverser1   r1   r2   test_alternative_exact  s   z!TestAnsari.test_alternative_exactzx, y, alternative, expectedr  r  r  g@A?r  rw   )rz   r{   r8   r;   r<   gG$I?g%I$I?ry   )r8   r;   r<   c                 C   s$   t j|||dj}t||dd d S )Nr  r  r  )r   r  r   r   )rO   r]   r  r  expectedr  r1   r1   r2   test_alternative_exact_with_R  s   z(TestAnsari.test_alternative_exact_with_Rc                 C   s   t jjddddd}t jjddddd}t j||ddj}t j||d	dj}t|d
dd t|ddd t jjddddd}t jjddddd}t ||j}t j||ddj}t j||d	dj}t||d dd t|d|d  dd d S )Nr   r{   r   {   r=  r   r  r  r  r   r  r  r   r|        ?r   )r   r   r.   r  r   r   )rO   r   r   r  r  r  r1   r1   r2   test_alternative_approx  s   z"TestAnsari.test_alternative_approxN)rb   rc   rd   r  r  r  rn   ra   r  r  r*  rH  rI  r  r"  r1   r1   r1   r2   r  M  s(    		
r  c                   @   4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )TestBartlettc              
   C   B   t ttttttttt	g
}t
j| \}}t|dd t|dd d S )NgzU@ /4@r;   gG?)g1g2g3g4g5g6g7g8g9g10r   bartlettr   )rO   r/   Tr  r1   r1   r2   	test_data     zTestBartlett.test_datac                 C      t ttjdg d S rl   )rV   rW   r   r0  rX   r1   r1   r2   rn        zTestBartlett.test_bad_argc              
   C   4   t ttttttttt	g
}t
j| }d}t|| d S Nr
  )r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r   r0  r   rO   r/   r_   r^   r1   r1   r2   ra        
z#TestBartlett.test_result_attributesc                 C   s6   t ttttttttt	g f}t
tjtjftj|  d S r*   )r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r	   r+   r   r   r0  )rO   r/   r1   r1   r2   test_empty_arg  s   zTestBartlett.test_empty_argc                 C   *   t ddgddgg}tttjt| d S Nr   r   rx   rz   )r+   r   rV   rW   r   r0  r&  r   r1   r1   r2   test_1d_input     zTestBartlett.test_1d_inputN)rb   rc   rd   r2  rn   ra   r:  r=  r1   r1   r1   r2   r$    s    r$  c                   @   rJ  )
TestLevenec              
   C   r%  )NgYpK?r;   gI{]?)r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r   levener   )rO   r/   r  r  r1   r1   r2   r2    r3  zTestLevene.test_datac                 C   sF   t jtttdd\}}t jtttddd\}}t|| t|| d S )NrQ   centertrimmedr   rB  proportiontocut)r   r@  r&  r'  r(  r   )rO   W1pval1W2pval2r1   r1   r2   test_trimmed1  s   

zTestLevene.test_trimmed1c           
      C   s   g d}g d}t jd t j|}tj||ddd\}}tj||ddd\}}tj|dd |dd d	d
\}}	t|| t|| t||	 d S )Ng333333?       @      @      @      @      @      @g      Y@r   rM        @rN  g      @rO  g      @g      i@  rC        ?rD  r   r   rQ   rA  )r+   r>  seedpermutationr   r@  r   )
rO   r]   r  r   W0pval0rF  rG  rH  rI  r1   r1   r2   test_trimmed2  s   



$

zTestLevene.test_trimmed2c                 C   sn   t ddd}t jd t j|}|d }tj||dd\}}tj||dd\}}t|| t|| d S )	Nr   r   r"  rT  rx   rQ   rA  median)r+   r  r>  rV  rW  r   r@  r   )rO   r]   r   r  rF  rG  rH  rI  r1   r1   r2   test_equal_mean_median  s   
z!TestLevene.test_equal_mean_medianc                 C   &   t ddd}tttj||dd d S Nr   r   r"  皙?)portiontocut)r+   r  rV   	TypeErrorr   r@  r   r1   r1   r2   test_bad_keyword  r   zTestLevene.test_bad_keywordc                 C   r]  Nr   r   r"  trimrA  )r+   r  rV   rW   r   r@  r   r1   r1   r2   test_bad_center_value  r   z TestLevene.test_bad_center_valuec                 C   r4  rl   )rV   rW   r   r@  rX   r1   r1   r2   test_too_few_args#  s   zTestLevene.test_too_few_argsc              
   C   r6  r7  )r&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r   r@  r   r8  r1   r1   r2   ra   &  r9  z!TestLevene.test_result_attributesc                 C   r;  r<  )r+   r   rV   rW   r   r@  r&  r   r1   r1   r2   r=  -  r>  zTestLevene.test_1d_inputN)rb   rc   rd   r2  rJ  rZ  r\  rb  re  rf  ra   r=  r1   r1   r1   r2   r?    s    	
r?  c                   @   sd   e Zd ZdZedddZdd Zd	d
 Zdd Zdd Z	dd Z
ejjejdkdddd ZdS )TestBinomTestPzU
    Tests for stats.binomtest as a replacement for deprecated stats.binom_test.
    Nr  	two-sidedc                 C   s   t | t j} t| dkr| d | d  }| d } n t| dkr8| d } |d u s.|| k r2tdt |}ntdtj| |||d}|jS )Nr   r   r   zn must be >= xzIncorrect length for x.ra  r  )	r+   
atleast_1dastypeint_r   rW   r   	binomtestr   )r]   r   ra  r  resultr1   r1   r2   binom_test_func6  s   
zTestBinomTestP.binom_test_funcc                 C   sR   |  dd}t|dd |  dd}t|dd | j dd	gd
d}t|dd d S )Nr      g-"^?r  rq  i  g`w?i           ?ra  gz?ro  r   )rO   r  r1   r1   r2   r2  H  s   zTestBinomTestP.test_datac                 C   s   t t| jg d d S Nrw   rV   rW   ro  rX   r1   r1   r2   test_bad_len_xP     zTestBinomTestP.test_bad_len_xc                 C   s(   t t| jdg t t| jdgdd d S )Nr   r   r   rv  rX   r1   r1   r2   
test_bad_nT  s   zTestBinomTestP.test_bad_nc                 C   s   t t| jddgdd d S )Nr   rL  rs  rv  rX   r1   r1   r2   
test_bad_p[  s   
zTestBinomTestP.test_bad_pc                 C   sX   | j ddddd}t|d | j ddddd}t|d | j dddd	d}t|d
 d S )Nr     gUUUUUU?r  ri  g-$Ǻl?r  gs0,j.?rh  gJfJ"f?rt  rO   r_   r1   r1   r2   test_alternatives_  s   

z TestBinomTestP.test_alternativesl        z32-bit does not overflowreasonc                 C   s   t t| jddtjjd d S )NrO  r8   rs  )rV   OverflowErrorro  sys
float_infominrX   r1   r1   r2   test_boost_overflow_raisesi  s   z)TestBinomTestP.test_boost_overflow_raises)Nr  rh  )rb   rc   rd   __doc__staticmethodro  r2  rw  rz  r{  r~  r*  rH  skipifr  maxsizer  r1   r1   r1   r2   rg  2  s    
rg  c                   @   s  e Zd ZdZdd Zdd Zdd Zdd	 Zej	
d
g ddd Zej	
d
g ddd Zej	
dg ddd Zej	
dg ddd Zej	
dg dg dg dg dg dg dg d g d!g d"g d#g d$g d%g d&g d'g d(g d)g d*g d+g d,g d-g d.g d/g d0g d1gd2d3 Zd4d5 Zej	
d6d7d8gd9d: Zd;d< Zd=d> Zd?d@ ZdAdB ZdCS )DTestBinomTestzTests for stats.binomtest.c                 C   s   d}t ddd}t|jd|d t ddd}t|jd|d t d	dd}t|jd
|d t ddd}t|jddd d S )N绽|=iΙ i@o@Q?r   rA   iΙ gX7?i	ϙ g=?iϙ gybw?rE   r   rm  r   r   rO   rB   r_   r1   r1   r2   test_two_sided_pvalues1v  s   z%TestBinomTest.test_two_sided_pvalues1c                 C   s   d}t jdddd}t|jd|d t ddd}t|jd	|d t d
dd}t|jd|d t ddd}t|jd|d t ddd}t|jd|d t ddd}t|jd|d t ddd}t|jd|d d S )Nr  r9   r"  r  r   ra  gX-h?rA   rz   g㤶y?r  g	y}?r;   g+,?rx   r   r  g      ?r   r  g{Gz?g333333?g<,?r  r  r1   r1   r2   test_two_sided_pvalues2  s   z%TestBinomTest.test_two_sided_pvalues2c                 C   s   d}t ddd}t|jd|d t ddd	}t|jd|d t d
dd}t|jd|d t ddd}t|jd|d t dd
d}t|jd|d t ddd}t|jd|d d S )Nr  i  i  r  r   rA   rx   r  gAL &W?r=   r  g!YB?r[   r  gE]t?r;   r8   r  r  r  r1   r1   r2   test_edge_cases  s   zTestBinomTest.test_edge_casesc                    s,  d dd}t t    d }tj| }t jtj| |kdd}t fdd| t    } | t|tj| k }t	||d	d
 d}t t 
  d }tj| }t jtj| |kdd}t fdd|dt 
  }|d }t	||d	d
 d S )Nr   r  rx   r   r   axisc                    s   t j|   S r*   r   binomr  r   r  r1   r2   <lambda>  rZ   z>TestBinomTest.test_binary_srch_for_binom_tst.<locals>.<lambda>rE   rA   r;   c                    s   t j|  S r*   r  r  r  r1   r2   r       )r+   r\   ceilr   r  r  sumr   intr   floor)rO   kr`   dy1ixy2r1   r  r2   test_binary_srch_for_binom_tst  s(    z,TestBinomTest.test_binary_srch_for_binom_tstz"alternative, pval, ci_low, ci_high))r  g$W?r   g?)r  g}?gGI~?r   )rh  g}G?g>?6?g4M?c                 C   X   t jddd|d}t|j|dd t|jd |jdd	}t|j|jf||fdd d S )
Nr  r   r  r   ra  r  r  rA   皙?r  confidence_level	r   rm  r   r   r	   rL   proportion_cilowhighrO   r  r  ci_lowci_highr_   cir1   r1   r2   test_confidence_intervals1  
   z(TestBinomTest.test_confidence_intervals1))r  ga&+w?r   gny?)r  g5!Rx?g5%?r   )rh  g~th?gSuR&|?g|KnG?c                 C   r  )
Nrx   r   r  r  r   rA   gQ?r&   r  r  r  r1   r1   r2   test_confidence_intervals2  r  z(TestBinomTest.test_confidence_intervals2zalternative, pval, ci_high))r  gլ?g̓A?)r  r   r   )rh  gEow?g;Kڞj?c                 C   sN   t jddd|d}t|j|dd |jdd}t|jd	 t|j|dd d S )
Nr   r   r  ri  r   rA   r  r  r   )r   rm  r   r   r  r	   r  r  )rO   r  r  r  r_   r  r1   r1   r2   !test_confidence_interval_exact_k0  
   z/TestBinomTest.test_confidence_interval_exact_k0zalternative, pval, ci_low))r  r   r   )r  6>g@>}_?)rh  r  gcڒ ?c                 C   sN   t jddd|d}t|j|dd |jdd}t|jd t|j|dd d S )	Nr   r  ri  r   rA   r  r  r   )r   rm  r   r   r  r	   r  r  )rO   r  r  r  r_   r  r1   r1   r2   %test_confidence_interval_exact_k_is_n  r  z3TestBinomTest.test_confidence_interval_exact_k_is_nz+k, alternative, corr, conf, ci_low, ci_high)rx   rh  Tr  g`?g?)rx   rh  Tr&   grt!?g/?)rx   rh  Fr  gVK5?g5dM?)rx   rh  Fr&   gXCVu^?g*r[?)rx   r  Tr  r   gP-V?)rx   r  Tr&   r   g^'?)rx   r  Fr  r   gS5a?)rx   r  Fr&   r   g?)rx   r  Tr  g!Շ`ⰸ?r   )rx   r  Tr&   g$x"?r   )rx   r  Fr  gW~=?r   )rx   r  Fr&   g	L?r   )r   rh  Tr  r   gX?)r   rh  Fr  r   gN?)r   r  Tr  r   g/<7#9?)r   r  Fr  r   gAA?)r   r  Tr  r   r   )r   r  Fr  r   r   )r   rh  Tr  g7?r   )r   rh  Fr  g؍s?r   )r   r  Tr  r   r   )r   r  Fr  r   r   )r   r  Tr  ghadn?r   )r   r  Fr  gAB/?r   c           
      C   sL   t j|dd|d}|rd}nd}|j||d}	t|	j|	jf||fdd d S )	Nr   r_  r  wilsonccwilson)r  methodr   rA   )r   rm  r  r   r  r  )
rO   r  r  corrconfr  r  r_   r  r  r1   r1   r2   test_ci_wilson_method  s   z#TestBinomTest.test_ci_wilson_methodc                 C   s*   t ddd}t|jd t|jd d S )Nrz   r  r  r   )r   rm  r	   rL   r   r}  r1   r1   r2   &test_estimate_equals_hypothesized_prop#  s   z4TestBinomTest.test_estimate_equals_hypothesized_propzk, n)r   r   )r   r   c                 C   s@   t jtdd t|| W d    d S 1 sw   Y  d S )Nz must be an integer not less thanr(  r*  r   rW   r   rm  )rO   r  r   r1   r1   r2   test_invalid_k_n+  s
   "zTestBinomTest.test_invalid_k_nc                 C   sB   t jtdd tddd W d    d S 1 sw   Y  d S )Nzk must not be greater than nr(  r  r   r  r  rX   r1   r1   r2   test_invalid_k_too_big1  s
   "z$TestBinomTest.test_invalid_k_too_bigc                 C   P   t jdddd}tjtdd |jdd W d    d S 1 s!w   Y  d S )	Nrx   r   r_  r  zmust be in the intervalr(  r   r  r   rm  r*  r   rW   r  r}  r1   r1   r2   test_invalid_confidence_level6     "z+TestBinomTest.test_invalid_confidence_levelc                 C   r  )	Nrx   r   r_  r  zmethod must ber(  zplate of shrimpr  r  r}  r1   r1   r2   test_invalid_ci_method;  r  z$TestBinomTest.test_invalid_ci_methodc                 C   s"   t jdddd}t|j|j d S )Nrx   r   r_  r  )r   rm  r	   proportion_estimaterL   r}  r1   r1   r2   
test_alias@  s   zTestBinomTest.test_aliasN)rb   rc   rd   r  r  r  r  r  r*  rH  rI  r  r  r  r  r  r  r  r  r  r  r  r1   r1   r1   r2   r  o  st     


	
	


r  c                   @   sD   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dS )TestFlignerc                 C   s&   t d}tt||d dd d S )Nr{   r   )gnf	@gDq?r  )r+   r\   r   r   fligner)rO   r   r1   r1   r2   r2  G  s
   
zTestFligner.test_datac           	         sv   t jd  fdd}|t}|t}|t}tj|||dd\}}tj|||ddd\}}t|| t|| d S )	Nr   c                    s    t | d t|    S )Nr  )r+   asarrayrandnr   tolist)gr   r1   r2   _perturbS      z+TestFligner.test_trimmed1.<locals>._perturbrQ   rA  rC  r   rD  )	r+   r>  r   r&  r'  r(  r   r  r   )	rO   r  g1_g2_g3_Xsq1rG  Xsq2rI  r1   r  r2   rJ  N  s   

zTestFligner.test_trimmed1c                 C   sb   g d}g d}t j||ddd\}}t j|dd |dd dd	\}}t|| t|| d S )
NrK  rR  rC  rU  rD  r   r   rQ   rA  )r   r  r   )rO   r]   r  r  rG  r  rI  r1   r1   r2   rZ  a  s   

$
zTestFligner.test_trimmed2c                 C   r]  r^  )r+   r  rV   ra  r   r  r   r1   r1   r2   rb  {  r   zTestFligner.test_bad_keywordc                 C   r]  rc  )r+   r  rV   rW   r   r  r   r1   r1   r2   re    r   z!TestFligner.test_bad_center_valuec                 C   r4  rl   )rV   rW   r   r  rX   r1   r1   r2   test_bad_num_args  r5  zTestFligner.test_bad_num_argsc                 C   s.   t d}tt jt jft||d g  d S )Nr{   r   )r+   r\   r	   r   r   r  r   r1   r1   r2   r:    s   
$zTestFligner.test_empty_argN)
rb   rc   rd   r2  rJ  rZ  rb  re  r  r:  r1   r1   r1   r2   r  E  s    r  c                  c   s    g d} g d}t |D ]S\}}tj|}|d}|jdddd}|jdddd}t||D ]\}}	t|d ||	 D ]}
|| ||
< q>q1|| t|d	\}}||d
g| | R V  qd S )N))g]YCgQɣ?)gme5gϵ5?)g;jgH?)iFeiHcMi%:r   r   r  r{   )r  r  r   r   r   r  )		enumerater+   r>  r?  integersziprangeshufflesplit)expected_resultsseedssirV  rB  xytie_indnum_ties_per_indr`   r   jr]   r  r1   r1   r2   mood_cases_with_ties  s    

r  c                   @   s   e Zd Zejde dd Zejdg ddd Zdd	 Z	d
d Z
dd Zdd Zdd Zdd Zejdg ddd ZdS )TestMoodz$x,y,alternative,stat_expect,p_expectc                 C   s4   t j|||d\}}t||dd t||dd dS )am  
        Example code used to generate SAS output:
        DATA myData;
        INPUT X Y;
        CARDS;
        1 0
        1 1
        1 2
        1 3
        1 4
        2 0
        2 1
        2 4
        2 9
        2 16
        ods graphics on;
        proc npar1way mood data=myData ;
           class X;
            ods output  MoodTest=mt;
        proc contents data=mt;
        proc print data=mt;
          format     Prob1 17.16 Prob2 17.16 Statistic 17.16 Z 17.16 ;
            title "Mood Two-Sample Test";
        proc print data=myData;
            title "Data for above results";
          run;
        r  gؗҜ<r  Nr   moodr   )rO   r]   r  r  stat_expectp_expectrL   r   r1   r1   r2   test_against_SAS  s   zTestMood.test_against_SASzalternative, expected))rh  )3oXQ?gpL?)r  )r  gԣl?)r  )r  gpL?c                 C   s.   g d}g d}t j|||d}t|| d S )Nr  r   r  r  )rO   r  r  r]   r  r_   r1   r1   r2   test_against_SAS_2  s   	zTestMood.test_against_SAS_2c                 C   s`   t jd t jdd}t jdd}t||\}}t||\}}t||g| |g d S )NrT  r   r   r[   )r+   r>  rV  r  r   r  r   )rO   r   r   z1p1z2p2r1   r1   r2   test_mood_order_of_args  s   z TestMood.test_mood_order_of_argsc                 C   sN   g d}g d}t |}t |}d|_d|_ttj||d dddg d S )N)gg۟?gy<7xg1)E?g9cç?g
FAg<=ӟ	2?g{![?g_l?gڻ{ӿgGSsA0?g1?g 3ggL?g_lg,*.gː53?gvNqG?gW=?)gd"g `Y#?g^g?gNgBagp5_%?g(?gv]lοgg?ggr6ʺ?ḡj?gK!ٿg{g&?g*h}gZ?gyBU?g .ۋ4?g|w;?g&ck @gRj2g}+w(o?g@F?g}_9t?g =g	S?gx(gg7Y?g{ih?)r   r   )r[   r   r  g4g^}E	?)r+   r   shaper   r   r  )rO   r   r   r1   r1   r2   test_mood_with_axis_none  s   

z!TestMood.test_mood_with_axis_nonec              
   C   s   d}t jd t jd|}t jd|}t||\}}t|D ]}t|| || gt|d d |f |d d |f  q"| }| }tj||dd\}}t|D ]}t|| || gt||d d f ||d d f  qVd S )Nr{   rT  r   r[   r   r  )	r+   r>  rV  r  r   r  r  r   	transpose)rO   nyr   r   	z_vectestpval_vectestr  r`   r1   r1   r2   test_mood_2d  s"   ""zTestMood.test_mood_2dc              	   C   s>  d}t jd t jj| }t jj| }tdD ]}tj|||d\}}ddgddgddgf}t||| d  D ]b}t||| d  D ]U}	|dkr^|d d ||	f }
|d d ||	f }n)|dkru||d d |	f }
||d d |	f }n|||	d d f }
|||	d d f }t|||	f |||	f gt|
| qEq9qd S )N)r   r{   r8   rT  rx   r  r   r   r   )r+   r>  rV  r  r  r   r  r   )rO   r  r   r   r  r  r  axes_idxr`   r  slice1slice2r1   r1   r2   test_mood_3d  s2   

zTestMood.test_mood_3dc                 C   s   t ttjdgg  d S rl   )rV   rW   r   r  rX   r1   r1   r2   test_mood_bad_arg0     zTestMood.test_mood_bad_argc           	      C   s   t jd tjjddd}tjjddd}tj||dd\}}tj||dd\}}tj||d	d\}}||  kr?|ksBJ  J t|dd
d t||d  t|d|d   tj	t
dd tj||dd W d    d S 1 stw   Y  d S )Nr   rr  r   )r   r   g      ?rh  r  r  r  Hz>r  r   r   zalternative must be...r(  z	ekki-ekki)r+   r>  rV  r   r   r.   r  r   r*  r   rW   )	rO   r]   r  stat1r  stat2r  stat3p3r1   r1   r2   test_mood_alternative5  s   "zTestMood.test_mood_alternativer  )rh  r  r  c                 C   sF   t jd}|d}|d}tj|||d}t|j|jf| d S )Nl	   hTt fU6H~ )r   r   )r[   r   r  )	r+   r>  r?  r   r   r  r	   rL   r   )rO   r  rB  r   r   r_   r1   r1   r2   test_resultG  s
   

zTestMood.test_resultN)rb   rc   rd   r*  rH  rI  r  r  r  r  r  r  r  r  r  r  r1   r1   r1   r2   r    s"    
 
	
r  c                   @   sV   e Zd Zdd Zdd Zdd Zejje	 ddd	d
 Z
dd Zdd Zdd ZdS )TestProbplotc                 C   sl   t jjddd}t j|dd\}}g d}t|t| t|| t j|dd\}}g d}t|| d S )	Nr  90  r=  Fr   )g+/gY3gtg5Z0gFjg!2gsR{ܿg7Կg>hǿgcz隇gcz隇?g>h?g7?gsR{?g!2?gFj?g5Z0?gt?gY3?g+/?T)gƝV?g?gà?)r   r   r.   probplotr   r+   sort)rO   r]   osmosrosm_expectedr_   res_fitres_fit_expectedr1   r1   r2   rT   S  s   
zTestProbplot.test_basicc           
      C   s   t jjddd}t j|d dd\}}t j|ddd\}}t j|ddd\}}t|| t|| t|| t|| t j|ddd\}}	d S )Nr   i@ r=  Fsparamsr   r   r1   )r   r   r.   r  r   )
rO   r]   osm1osr1osm2osr2osm3osr3r  r  r1   r1   r2   test_sparams_keywordb  s   



z!TestProbplot.test_sparams_keywordc                 C   s   t jjddd}t j|dddd\}}t j|dt jdd\}}t|| t|| ttt j|dd	 ttt j|g d	 G d
d d}t j|ddd\}}t j|| dd\}}t|| t|| d S )Nr  r	  r=  Ft)rx   )r   r   r  zwrong-dist-namer   c                   @   s   e Zd ZdZdd ZdS )z3TestProbplot.test_dist_keyword.<locals>.custom_distz6Some class that looks just enough like a distribution.c                 S   s   t jj|ddS )Nr   )r   )r   r   ppf)rO   qr1   r1   r2   r  |  r  z7TestProbplot.test_dist_keyword.<locals>.custom_dist.ppfN)rb   rc   rd   r  r  r1   r1   r1   r2   custom_distz  s    r  )r   r  )r   r   )	r   r   r.   r  r  r   rV   rW   AttributeError)rO   r]   r  r  r  r  r  r1   r1   r2   test_dist_keywordp  s   


zTestProbplot.test_dist_keywordno matplotlibr  c           
      C   s  t  }|d tjjdddd}tj|t d\}}t   tj|d d\}}tj|dt d}t   tj|dd d}tt	|t	|  koXt	|  koXt	|  koXd	kn   t
|| t
|| t
|| t
|| t  }|d}	tj|d|	d t   d S )
Nr  rx   r   t r=  plotF)r   r$  r   )pltfigureadd_subplotr   r  r.   r  closer   r   r   )
rO   figr]   res1fitres1res2fitres2res3res4axr1   r1   r2   test_plot_kwarg  s$   
@




zTestProbplot.test_plot_kwargc                 C   r   r   )rV   rW   r   r  rX   r1   r1   r2   test_probplot_bad_args     z#TestProbplot.test_probplot_bad_argsc                 C   sZ   t tjg ddtg tg f t tjg ddtg tg ftjtjdff d S )NFr
  Tr   )r	   r   r  r+   r   r   rX   r1   r1   r2   
test_empty  s   zTestProbplot.test_emptyc                 C   sj   t jdd% ttjdgddt dgt dgft jt jdff W d    d S 1 s.w   Y  d S )Nr   )invalidr   Tr
  r   )r+   r   r	   r   r  r   r   rX   r1   r1   r2   test_array_of_size_one  s   "z#TestProbplot.test_array_of_size_oneN)rb   rc   rd   rT   r  r   r*  rH  r  have_matplotlibr1  r2  r4  r6  r1   r1   r1   r2   r  Q  s    
r  c                   @   s   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zdd Zejdg dg dg dgdd Zdd ZdS ) TestWilcoxonc                 C   sh   t ttjdgddg t ttjddgddgd t ttjddgddgdd t ttjdgd dd d S )Nr   r   dummyr  r   xyzmode)rV   rW   r   wilcoxonrX   r1   r1   r2   test_wilcoxon_bad_arg  s   z"TestWilcoxon.test_wilcoxon_bad_argc                 C   sR   t d}tttj||ddd tttj||ddd ttj||dddd d S )Nr  wilcoxapproxr;  prattzsplit)g     @Z@r   )r+   r\   rV   rW   r   r=  r	   r   r1   r1   r2   test_zero_diff  s   
zTestWilcoxon.test_zero_diffc                 C   sd   g d}g d}t  }|jtdd tj||ddd}W d    n1 s&w   Y  t|d d S )	Nr  )r   r   rx   r{   Sample size too smallr  rA  r@  zero_methodr<  )r   g|N?)r
   r  r3  r   r=  r   )rO   r]   r  r  r_   r1   r1   r2   
test_pratt  s   zTestWilcoxon.test_prattc                 C   s<   g d}t j|ddd}t j|ddd}t j|ddd}d S )N)r   r   rx   r   r   rx   r   r   r   r   r   rA  r@  rE  rB  r?  )r   r=  )rO   arrr  r1   r1   r2   test_wilcoxon_arg_type  s   z#TestWilcoxon.test_wilcoxon_arg_typec                 C   s2  g d}t dd}tdd t||D }t|j}tj||ddd\}}t|d	 t|d
 tj||ddd\}}t|d t|d tj||ddd\}}t|d t|d t	g d}t	g d}tj||ddd\}}t
|d t|ddd tj||ddd\}}t
|d t|ddd d S )N)	r   rz   r  r[   r<   rz   r{   r   r   r{   c                 S   s   g | ]	\}}|g| qS r1   r1   )rt   ur   r1   r1   r2   ru     s    z7TestWilcoxon.test_accuracy_wilcoxon.<locals>.<listcomp>rA  r@  r;  i  g.i?rB  i  gFYv]Uj?r?  iG  g%^ Ez?r  r  ri  r        y   r     r     p         w   r	  rR     r  r  s   rO  r     F
correctionr<  rt  g^҂<?r   rA   Tg:b+?)r  r+   concatenater  zerosr   r   r=  r   r   r	   )rO   freqnumsr]   r  r1  ra  r1   r1   r2   test_accuracy_wilcoxon  s*   








z#TestWilcoxon.test_accuracy_wilcoxonc                 C   s@   t g d}t g d}tj||ddd}d}t|| d S )NrL  rS  Fr@  rZ  r
  )r+   r   r   r=  r   )rO   r]   r  r_   r^   r1   r1   r2   test_wilcoxon_result_attributes  s
   z,TestWilcoxon.test_wilcoxon_result_attributesc                 C   s   t jd}|d|d}}tj||dd}tj|jd }t|j	| tj||dd}t
|dr7J t||}t
|drDJ d S )Nl   	m$S r[   r@  r;  r   exact
zstatistic)r+   r>  r?  r   r=  r   r  r   r   rc  hasattr)rO   rB  r]   r  r_   r   r1   r1   r2   test_wilcoxon_has_zstatistic  s   z)TestWilcoxon.test_wilcoxon_has_zstatisticc                 C   sn   t jdgd dd\}}d}t|d t||dd t jdgd d	dd
\}}d}t|d t||dd d S )Nr_  r   r@  r;  g4IťY?r   r   rA   TrZ  gc2_?)r   r=  r	   r   )rO   statra  
expected_pr1   r1   r2   test_wilcoxon_tie	  s   	

zTestWilcoxon.test_wilcoxon_tiec                 C   s  g d}g d}t  }|jtdd tj||ddd\}}W d    n1 s(w   Y  t|d t|d	d
d t  }|jtdd tj||dddd\}}W d    n1 sZw   Y  t|d t|dd
d t  }|jtdd tj||ddd\}}W d    n1 sw   Y  t|d t|dd
d t  }|jtdd tj||dddd\}}W d    n1 sw   Y  t|d t|dd
d d S )N)
}   rX  r     rj  rX  rj  ri  rj  r
  )
r  z   ri  r  rj  |   r      r
     rD  r  r  r@  r  r<  r%  g13}?r8   rH   T)r  r[  r<  g7a%?r  g?gڏoJT?)r
   r  r3  r   r=  r	   r   )rO   r]   r  r  r   ra  r1   r1   r2   test_onesided  s<   	



zTestWilcoxon.test_onesidedc                 C   sZ   t ddD ]%}t|}t|}t||d  d d t| tt|d t|| qd S )Nr   r  r   )r  r   r   r	   r   r  r   )rO   r   pmf1pmf2r1   r1   r2   test_exact_basicB  s   zTestWilcoxon.test_exact_basicc                 C   s  t g d}t g d}tj||ddd\}}t|ddd tj||d	dd\}}t|d
dd tj||ddd\}}t|ddd t ddd }t ddd}tj||ddd\}}t|ddd tj||d	dd\}}t|ddd tj||ddd\}}t|ddd d S )N)
g(\?g=
ףp=?g(\?gQ޿gQ?g{Gz?gp=
ףgq=
ףp?g      gQ?)
r   g?gɿg333333?皙gܿgzGgQοg(\gRQrh  rb  ro  g֔  ?r8   rH   r  gYJ?r  g*  ?r   r  r  r   g(?g(?gS?)r+   r   r   r=  r   r\   )rO   r]   r  r  ra  r1   r1   r2   test_exact_pvalJ  s    zTestWilcoxon.test_exact_pvalr]   )r   r   rx   )r   r   r  rJ  r{   )r   r   rx   rJ  r  ir;   r<   c                 C   s@   t |\}}t|}||dk  }t|| t|d d S Nr   r   )r   r=  r+   r   r  r	   )rO   r]   r   ra  wtruer1   r1   r2   test_exact_p_1f  s
   

zTestWilcoxon.test_exact_p_1c                 C   s   t ddd }t ddd}tt||tj||dd t dd}t }|jtdd	 t|\}}W d    n1 s@w   Y  ttj|d
d||f t dd}tt|tj|d
d d S )Nr   r  r  r   rb  r;  r=   zExact p-value calculationr  r@  r   r  )r+   r\   r	   r   r=  r
   r  r3  )rO   r]   r  r  r  r   ra  r1   r1   r2   	test_autop  s   zTestWilcoxon.test_autoN)rb   rc   rd   r>  rC  rG  rI  r`  ra  re  rh  rp  rs  ru  r*  rH  rI  rx  ry  r1   r1   r1   r2   r8    s$    	
	&
r8  c                   @   re   )
	TestKstatc                    s   t jd t jd  fdddD }g d}t||dd tj d	d
}tj dd
}tj dd
}t|||f|d d ddd d S )Ni}  r	  c                    s   g | ]}t  |qS r1   )r   kstat)rt   r   rP   r1   r2   ru     s    z>TestKstat.test_moments_normal_distribution.<locals>.<listcomp>r  )gHQ,?glZ)rI?g88}?gcJN?r  rA   r   )momentr   rx   r   g{Gz?{Gz?r  )r+   r>  rV  r  r   r   r}  )rO   momentsr  m1m2m3r1   r|  r2    test_moments_normal_distribution  s   "z*TestKstat.test_moments_normal_distributionc                 C   rU   r*   )rV   rW   r   r{  rX   r1   r1   r2   rY     rZ   zTestKstat.test_empty_inputc                 C   *   t d}t j|d< tt|t j d S NrC   r8   r+   r\   r   r	   r   r{  rm   r1   r1   r2   r        

zTestKstat.test_nan_inputc                 C   s*   t d}dD ]}tttj||d qd S )Nr   )r   g/$@ry  )r+   r\   rV   rW   r   r{  rO   rP   r   r1   r1   r2   test_kstat_bad_arg  s   
zTestKstat.test_kstat_bad_argN)rb   rc   rd   r  rY   r   r  r1   r1   r1   r2   rz    s
    rz  c                   @   r5   )TestKstatVarc                 C   rU   r*   rV   rW   r   kstatvarrX   r1   r1   r2   rY     rZ   zTestKstatVar.test_empty_inputc                 C   r  r  r  rm   r1   r1   r2   r     r  zTestKstatVar.test_nan_inputc                 C   s    dg}d}t ttj||d d S )Nr   r   ry  r  r  r1   r1   r2   rn     s   zTestKstatVar.test_bad_argN)rb   rc   rd   rY   r   rn   r1   r1   r1   r2   r    s    r  c                   @   sN   e Zd Zdd Zdd Zdd Zejje	 ddd	d
 Z
dd Zdd ZdS )TestPpccPlotc                 C      t ddddd | _d S Nr{     r"  r=  r3   r]   rX   r1   r1   r2   setup_method  r4   zTestPpccPlot.setup_methodc                 C   H   d}t j| jdd|d\}}g d}t|tjdd|d t|| d S )Nr{   r   N)g9%	?gxM _?g";9?g®XZ?g?num)r   	ppcc_plotr]   r   r+   r  )rO   r  svalsppccppcc_expectedr1   r1   r2   rT     
   zTestPpccPlot.test_basicc                 C   s   t j| jdddd\}}t j| jddt jd\}}t||dd t||dd t | jdd\}}t||dd t||dd d S )Nr  r   tukeylambdar   g#B;rA   )r   r  r]   r  r   )rO   svals1ppcc1svals2ppcc2svals3ppcc3r1   r1   r2   	test_dist  s   
zTestPpccPlot.test_distr!  r  c                 C   Z   t  }|d}tj| jddt d || |d}tj| jdd|d t   d S Nr  ir  r#  )r%  r&  r'  r   r  r]   delaxesr(  rO   r)  r0  r1   r1   r2   r1       


zTestPpccPlot.test_plot_kwargc                 C   s2   t ttj| jdd t ttjg ddddd d S )Nr   r   rw   r   r   )rV   rW   r   r  r]   rX   r1   r1   r2   test_invalid_inputs  s   
z TestPpccPlot.test_invalid_inputsc                 C   s@   t g dd\}}t|tjdddd t|tjdtd d S )Nr   r   r  r  r   )r   r  r   r+   r  r]  float)rO   r  r  r1   r1   r2   r4    s   zTestPpccPlot.test_emptyN)rb   rc   rd   r  rT   r  r*  rH  r  r7  r1  r  r4  r1   r1   r1   r2   r    s    
r  c                   @   re   )
TestPpccMaxc                 C   s   dg}t ttj|dd d S r   )rV   rW   r   ppcc_maxrm   r1   r1   r2   test_ppcc_max_bad_arg  s   z!TestPpccMax.test_ppcc_max_bad_argc                 C   s2   t jjddddddd }tt |dd	d
 d S )Nffffffr   r  '   r        @{r;   rH   )r   r  r.   r   r  r   r1   r1   r2   test_ppcc_max_basic  s   zTestPpccMax.test_ppcc_max_basicc                 C   sp   t jjddddddd }t j|dd	}t j|t jd	}t|d
dd t|d
dd t |}t|d
dd d S )Nr  r   r  r  r  r   r  r  r   r  r{   rH   )r   r  r.   r  r   )rO   r]   max1max2max3r1   r1   r2   r    s   
zTestPpccMax.test_distc                 C   s`   t jjddddddd }ttt j|dd	 tt j|d
d	ddd tt j|dd	ddd d S )Nr  r   r  r  r  r   r  )r   r   r  )brack)r   r   r  r;   rH   )r   r   )r   r  r.   rV   rW   r  r   r   r1   r1   r2   
test_brack  s   
zTestPpccMax.test_brackN)rb   rc   rd   r  r  r  r  r1   r1   r1   r2   r    s
    r  c                   @   r#  )TestBoxcox_llfc                 C   sT   t jjdddd}d}t ||}|j d tt| d  }t	|| d S )Nr  r   1  r   r   r   r   rL  r   )
r   r   r.   
boxcox_llfr   r+   r,   r  rS   r   )rO   r]   lmbdallfllf_expectedr1   r1   r2   rT     s
   $zTestBoxcox_llf.test_basicc                 C   D   t jjdddd}d}t ||}t |t|}t||dd d S Nr   r   r  r  r   r  rA   )r   r   r.   r  listr   rO   r]   r  r  llf2r1   r1   r2   test_array_like  
   zTestBoxcox_llf.test_array_likec                 C   P   t jjdddd}d}t ||}t |t||gj}t||g|dd d S r  )r   r   r.   r  r+   vstackr1  r   r  r1   r1   r2   test_2d_input  s
   zTestBoxcox_llf.test_2d_inputc                 C      t ttdg  d S rl   )r   r+   isnanr   r  rX   r1   r1   r2   r4  (     zTestBoxcox_llf.test_emptyc                 C   s"   g d}t d|}t|d d S )N)g     h@      m@r  g     x@ig#:x1)r   r  r   )rO   rP   r  r1   r1   r2   test_gh_6873+  s   zTestBoxcox_llf.test_gh_6873N)rb   rc   rd   rT   r  r  r4  r  r1   r1   r1   r2   r    s    
r  (  iU>  iϵ i i_
 i7 i is iv  i+ i* i-q ir: i  i i i i_ i4  i iD  i i, iE i\  i i| ig  iM iӢ i[  i! i if6 iP i i iI i/ iF& i i i`# i i  i& i i iQ ia i
 iXZ iw  i_ im  iD i- iB  i:C i	 i' i  i % ix i  iz_ iQ i!# i= i  ix i
 iU7 i_ i1k iW9 i& ix iI i  iT i
 i i  irv im i͟  i iO iv i) iP i4 i= ii  iU i i ig> ix	 i` i| i0 i` iM	 i ' i[ i8 i  i3[ i  i) i`N i
	 i  i[ i:+ ip i	% i
  iV i[$ iX  i_ ib ieZ i i`?	 it i  i iE i iT iy i imv i< iZ i! i i(  ir
 i i iɪ i| i i> ii  i i* i# ix i i& i i# i1g iP i) iۆ i{% i? i۽ i| i i& i- iư) i+c i7  iQ i	o i& i|_ iy
 i6 i  i i i:J i#j i8 i  i%  i i%W i( i` i: ina  i  i i% i$ i4 iA( iq  i  i" iy iXZ i izr i i$ i; is igY i i7 ip	 ip iL6" iDH i i!3 i iI( i|C iш i im  iS i  i{8 i  i i iP i2L i( i@ is iu i2 i_ ij i)m igZ i_ iu	 i9 ik  i'; i ii iK i܃ i i= iv i= i ia  ir  i i i}W i*  i  iJ i i  i2  iUF iNj	 iV," ip iA i% i iv9 iI5 i  iT- i  i8~ i?e` i1 iV  i/X i= i4 i iح i^ i i/ ig i( i i|B; ip i ic iE' iE i$
 i i{ ik i: iI iA_ i iP ioP  i6_ i iܲ iO	 iA7 iL  iA i iV  iwm i iZh i|V iԇC iw i& i~ i i'  i!  i i
 i} ie i3 i i i_ i  i@) i  iV ii iF i? i&
 i"Y iw i  iH iS ird ib i@ iK iA; iz| id+ i	 i
8 i"  i	 iO i_ i$l i!9 ii iT iuL iQ'% iqp iC i]K@ i! i i  iIp| i i?	B i i io0 i i

 i  i5% i i i* i i iv0 i* i* iw i  ive i$ i i- iD~ i- i i,4 ib i  i iP0 i) im i	 iR iQ i
 i i c                   @   s   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Ze	j
dg ddd Zdd Ze	j
dg ddd Zdd ZdS )
TestBoxcoxc                 C   s   t ddddd }tj|dd}t||d  tj|dd}t|dd|   tj|dd}t|t| tjt|dd}t|t| ttd	d
}t|t	d	 d S )Nr{   r   r	  r=  r   r  r   r   r   r   )
r3   r   r   r   r+   r,   r  r   r	   r]  )rO   r]   xtr1   r1   r2   test_fixed_lmbdak  s   zTestBoxcox.test_fixed_lmbdac                 C   sL   d}t jjdddd}|| d |  }t |\}}t|d| dd	 d S )
Ng      @r   iP  i  )r   r   r   r   r   r   rH   )r   r   r.   r   r   )rO   r  r]   x_invr  r   r1   r1   r2   test_lmbda_None}  s
   zTestBoxcox.test_lmbda_Nonec                 C   s   t jd}tdd|dd }tj|dd\}}}t|ddg tj|d	d\}}}t|d
dg tdd|dd }tj|dd\}}}t|ddg tj|dd\}}}t|ddg d S )NrT  r{   r   r=  rr  )alphagx@gv[%@r  gxˉEk?g8Bhk @r;   r  r[   r  g\?gx̡&@r    g7(U@gp:sY@)r+   r>  r   r3   r   r   r   )rO   rB  r]   r  ri   r1   r1   r2   
test_alpha  s   zTestBoxcox.test_alphac                 C   sR   t ddg}tttj| tttjt dg tttjt dgdgg d S )Nr   r   r   )r+   r   rV   rW   r   r   r   r1   r1   r2   test_boxcox_bad_arg  s    zTestBoxcox.test_boxcox_bad_argc                 C      t tg jdk d S N)r   )r   r   r   r  rX   r1   r1   r2   r4    r4   zTestBoxcox.test_emptyc                 C   s    t t\}}t|ddd d S )NgsHjdrr   rA   r   r   _boxcox_datar   )rO   r  lamr1   r1   r2   r    s   zTestBoxcox.test_gh_6873bounds)r   r   )r   r   )r   rt  c                    sH    fdd}t jtd |d\}} d |  k r d k s"J  J d S )Nc                       t j|  ddS Nbounded)r  r  r   minimize_scalarfunr  r1   r2   	optimizer     zBTestBoxcox.test_bounded_optimizer_within_bounds.<locals>.optimizerr  r  r   r   )r   r   r  )rO   r  r  r  r  r1   r  r2   $test_bounded_optimizer_within_bounds  s   (z/TestBoxcox.test_bounded_optimizer_within_boundsc                    sj   t jtd d\}}|d |d f ddi fdd}t jtd |d\}}||ks,J t| d	  d S )
Nr  r_  r   xatolr  c                    s   t j|  ddS )Nr  )r  r  optionsr  r  r  r  r1   r2   r    s   zPTestBoxcox.test_bounded_optimizer_against_unbounded_optimizer.<locals>.optimizerr  r   r  )rO   r  r  r  lmbda_boundedr1   r  r2   2test_bounded_optimizer_against_unbounded_optimizer  s   
z=TestBoxcox.test_bounded_optimizer_against_unbounded_optimizerr  )strr   r   r_  c                 C   sD   t jtdd tjtd |d W d    d S 1 sw   Y  d S )Nz`optimizer` must be a callabler(  r  r*  r   rW   r   r   r  rO   r  r1   r1   r2   $test_bad_optimizer_type_raises_error  s   "z/TestBoxcox.test_bad_optimizer_type_raises_errorc                 C   sP   dd }d}t jt|d tjtd |d W d    d S 1 s!w   Y  d S )Nc                 S   s   dS rl   r1   r  r1   r1   r2   r    s   zCTestBoxcox.test_bad_optimizer_value_raises_error.<locals>.optimizerz;`optimizer` must return an object containing the optimal...r(  r  r  )rO   r  r+  r1   r1   r2   %test_bad_optimizer_value_raises_error  s
   "z0TestBoxcox.test_bad_optimizer_value_raises_errorN)rb   rc   rd   r  r  r  r  r4  r  r*  rH  rI  r  r  r  r  r1   r1   r1   r2   r  i  s    
	
	
r  c                   @   sh   e Zd Zdd Zdd Zdd Zdd Zej	d	g d
ej	dg ddd Z
dd Zdd ZdS )TestBoxcoxNormmaxc                 C   r  Nr{   r   r	  r=  r  rX   r1   r1   r2   r    r4   zTestBoxcoxNormmax.setup_methodc                 C      t | j}t|ddd d S )Nm?r   rA   r   boxcox_normmaxr]   r   rO   r   r1   r1   r2   test_pearsonr     zTestBoxcoxNormmax.test_pearsonrc                 C   s<   t j| jdd}t|ddd t | j\}}t|| d S )Nmler  	.!?r   rA   )r   r  r]   r   r   )rO   r   r  maxlog_boxcoxr1   r1   r2   test_mle  s   zTestBoxcoxNormmax.test_mlec                 C   s&   t j| jdd}t|ddgdd d S )Nr   r  r  r  r   rA   r  )rO   
maxlog_allr1   r1   r2   test_all  s   zTestBoxcoxNormmax.test_allr  )r  pearsonrr   r  r  c                    sN    fdd}t j| j||d}t d |k sJ t| d k s%J d S )Nc                    r  r  r  r  r  r1   r2   r    r  zITestBoxcoxNormmax.test_bounded_optimizer_within_bounds.<locals>.optimizer)r  r  r   r   )r   r  r]   r+   r   )rO   r  r  r  r   r1   r  r2   r    s   
z6TestBoxcoxNormmax.test_bounded_optimizer_within_boundsc                    sv   t | j}t|d}t|d |d dG dd d  fdd}t j| j|d}||ks3J t||d	 d S )
Nr{   r~  i  c                   @   s   e Zd ZdS )z?TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.MyResultN)rb   rc   rd   r1   r1   r1   r2   MyResult  s    r  c                    s6   g }D ]	}| | | q  }t| |_|S r*   )appendr+   argminr]   )r  objsr  r_   r  lmbda_ranger1   r2   r    s   z@TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.optimizer)r  r   )r   r  r]   r+   roundr  r   )rO   r  lmbda_roundedr  lmbda2r1   r  r2   test_user_defined_optimizer  s   	z-TestBoxcoxNormmax.test_user_defined_optimizerc                 C   s^   t j}tj| jd |d tjtdd tj| jd|d W d    d S 1 s(w   Y  d S )N)r  r  z,`brack` must be None if `optimizer` is givenr(  )g       rL  )r   r  r   r  r]   r*  r   rW   r  r1   r1   r2   2test_user_defined_optimizer_and_brack_raises_error  s   
"zDTestBoxcoxNormmax.test_user_defined_optimizer_and_brack_raises_errorN)rb   rc   rd   r  r  r  r  r*  rH  rI  r  r  r  r1   r1   r1   r2   r    s    r  c                   @   sF   e Zd Zdd Zdd Zejje dddd Z	d	d
 Z
dd ZdS )TestBoxcoxNormplotc                 C   r  r  r  rX   r1   r1   r2   r  '  r4   zTestBoxcoxNormplot.setup_methodc                 C   r  )Nr{   r  r   r  )g!X4}?g΁i?g-	15?gW$?g\ܱ{?r  )r   boxcox_normplotr]   r   r+   r  )rO   r  lmbdasr  r  r1   r1   r2   rT   *  r  zTestBoxcoxNormplot.test_basicr!  r  c                 C   r  r  )r%  r&  r'  r   r  r]   r  r(  r  r1   r1   r2   r1  2  r  z"TestBoxcoxNormplot.test_plot_kwargc                 C   s.   t ttj| jdd t ttjddgdd d S )Nr   r   r   )rV   rW   r   r  r]   rX   r1   r1   r2   r  ?  s   z&TestBoxcoxNormplot.test_invalid_inputsc                 C   s   t tg ddjdk d S rv  )r   r   r  r   rX   r1   r1   r2   r4  E  r  zTestBoxcoxNormplot.test_emptyN)rb   rc   rd   r  rT   r*  rH  r  r7  r1  r  r4  r1   r1   r1   r2   r  &  s    
r  c                   @   r5   )TestYeojohnson_llfc                 C   r  )Nr   r   r  r  r   r  rA   )r   r   r.   yeojohnson_llfr  r   r  r1   r1   r2   r  K  r  z"TestYeojohnson_llf.test_array_likec                 C   r  r  )r   r   r.   r  r+   r  r1  r   r  r1   r1   r2   r  R  s
   z TestYeojohnson_llf.test_2d_inputc                 C   r  rl   )r   r+   r  r   r  rX   r1   r1   r2   r4  Y  r  zTestYeojohnson_llf.test_emptyN)rb   rc   rd   r  r  r4  r1   r1   r1   r2   r  I  s    r  c                   @   s   e Zd Zdd Zejdg ddd Zdd Zd	d
 Z	ejde
je
jgdd Zejde
je
je
je
jgdd Zdd ZdS )TestYeojohnsonc                 C   s  t jd}tdd|dd }t |dksJ tj|dd}t|| tj|dd}t|dd|d    tj|dd}t|t |d  tj|dd}t|| tdd|dd }t |dk seJ tj|d	d}t|t | d   tj|dd}t|| tj|d
d}t|d| d  d  tdd|dd	 }t |dk rJ t |dkrJ |dk}tj|dd}t|| ||  tj|dd}t|| dd|| d    tj|dd}t|| t || d  tj|dd}t|| ||  | }tj|d	d}t|| t ||  d   tj|dd}t|| ||  tj|d
d}t|| d||  d  d  d S )Nr	  r{   r   r=  r   r   r  r   r   rx   )	r+   r>  r   r3   r   r   
yeojohnsonr   r,   )rO   rB  r]   r  posnegr1   r1   r2   r  _  sL   


 $zTestYeojohnson.test_fixed_lmbdar  )r   r_  r  r   c                 C   s   dd }d}t jd t jjdd|d}|||}t|\}}t||dd	 tdt j	|| | d
d td|
 dd td| dd d S )Nc                 S   s   t j| j| jd}| dk}t|t dk r"t | | d ||< nt | | | d d| d ||< t|d t dkrZdt d|  | |   d dd|   || < |S dt | |    || < |S )Nr   r   r   r   r   )r+   r]  r  r   absspacingexppower)r]   r  r  r  r1   r1   r2   _inverse_transform  s   $
z:TestYeojohnson.test_lmbda_None.<locals>._inverse_transformi N  r  r   r   )r   r   r   r~  r  r   rH   )r+   r>  rV  normalr   r  r   r   linalgr   rQ   rS   )rO   r  r   	n_samplesr]   r  r  r   r1   r1   r2   r    s   
zTestYeojohnson.test_lmbda_Nonec                 C   r  r  )r   r   r  r  rX   r1   r1   r2   r4    r4   zTestYeojohnson.test_emptyc                 C   sD   t jjdddd}t |\}}t t|\}}t||dd d S )Nr   r   r  r  r  rA   )r   r   r.   r  r  r   )rO   r]   xt1r  xt2r1   r1   r2   r    s   zTestYeojohnson.test_array_liker   c                 C   sP   t jd|d}d}tjt|d t| W d    d S 1 s!w   Y  d S )Nr8   r   z>Yeo-Johnson transformation is not defined for complex numbers.r(  )r+   r\   r*  r   rW   r   r  )rO   r   r]   err_msgr1   r1   r2   test_input_dtype_complex  s
   "z'TestYeojohnson.test_input_dtype_complexc                 C   sZ   t jd|d}t jdt jd}t|\}}t|\}}t||dd t||dd d S )Nr<   r   r  rA   )r+   r\   float64r   r  r   )rO   r   x_intx_floatxt_int	lmbda_intxt_floatlmbda_floatr1   r1   r2   test_input_dtype_integer  s   z'TestYeojohnson.test_input_dtype_integerc                 C   sN   t g d}t|\}}t|d \}}t||dd t||dd d S )N)(HAG."A    5fAg(@Ag   VAgHzuoAg)\¿+Ar2  g(\^@Ag    RAr1  g    #Ag\µ=Ag
ףNAr0  gp=oHAg    6bAg{csAg)\tmhAg/BAg   '|Ag(\wipAg    ݃hAgQlBAgףp @Ar   r   rA   )r+   r   r   r  r   r   )rO   r]   xt_yeolam_yeoxt_boxlam_boxr1   r1   r2   test_input_high_variance  s
   z'TestYeojohnson.test_input_high_varianceN)rb   rc   rd   r  r*  rH  rI  r  r4  r  r+   	complex64
complex128r'  int8uint8int16int32r/  r7  r1   r1   r1   r2   r  ]  s    /
$

r  c                   @   r5   )TestYeojohnsonNormmaxc                 C   r  r  r  rX   r1   r1   r2   r    r4   z"TestYeojohnsonNormmax.setup_methodc                 C   r  )Ng?r   rA   )r   yeojohnson_normmaxr]   r   r  r1   r1   r2   r    r   zTestYeojohnsonNormmax.test_mlec                 C   s*   g d}t |}tj|dddsJ d S )N)gffffff@g r   rL  ffffff?g333333@rS  gffffff@g?g@rQ  rM  g"@r  g      gzG?r  r  )r   r?  r+   allclose)rO   r]   r  r1   r1   r2   test_darwin_example  s   
z)TestYeojohnsonNormmax.test_darwin_exampleN)rb   rc   rd   r  r  rB  r1   r1   r1   r2   r>    s    r>  c                   @   sB  e Zd Zejdejdfejdfej	dfgdd Z
dd Zejd	ejejfejejfej	ejfgd
d Zdd Zdd Zdd Zejdejdfejdfej	dfgdd Zejdejejej	gdd Zejdejejej	gdd Zejdejejdddfejejdddfej	ejdddfgd d! Zejdejd"eg d#ed$d%gdfejd&eg d'ejd( d)  edd*gdfej	d+eg d,ed-d.gdfgd/d0 Zejdejdfejdfej	dfgd1d2 Zejdejejej	gd3d4 Zejdejejej	gd5d6 Zejd7d8d9d)d:d;d<ejged8d9d)d:d;d<ejgd=d>d?d@ejdAdBgggejdejejej	gdCdD Zejd7d8d9d)d:d;d<ejged8d9d)d:d;d<ejgd=d>d?d@ejdAdBgggejdejejej	gdEdF Z dGdH Z!dIdJ Z"dKdL Z#dMS )NTestCircFuncsztest_func,expectedgA:wv?g @pz?gE|2@c                 C   s(   t g d}t||dd|dd d S Nc  r{   r   r  r   ^  h  r  r  rA   r+   r   r   rO   	test_funcr  r]   r1   r1   r2   test_circfuncs  s   zTestCircFuncs.test_circfuncsc                 C   s   t g d}| }tj|dd}t||dd |t j d  }|d }tj|dd}t||dd |	 }tj
|dd}t||dd d S )	N)r  r"  ry  rv     g     4@g3333333@rH  rI  r   rA   rM  rL  r  )r+   r   rQ   r   circmeanr   pirR   circvarrS   circstd)rO   r]   M1M2V1V2S1S2r1   r1   r2   test_circfuncs_small  s   z"TestCircFuncs.test_circfuncs_smallztest_func, numpy_funcc                 C   s<   t dgd dgd  }||}||}t||dd d S )Ng˛ɦv9?r   guv9?r   g:0yU>r  rJ  )rO   rL  
numpy_funcr]   circstatr!  r1   r1   r2   test_circfuncs_close	  s   z"TestCircFuncs.test_circfuncs_closec                       t g dg dg dg tj dd}tj  dd}t||dd tj ddd	} fd
dt jd D }t||dd tj ddd	} fddt jd D }t||dd d S )NrE  _  r;   rz   `  r9   ]  ie  r9   r<   if  rz   id  rH  rI  r@   rA   r   r  r  c                       g | ]}t j | d dqS rH  rI  r   rO  rt   r`   r]   r1   r2   ru   !	      z4TestCircFuncs.test_circmean_axis.<locals>.<listcomp>r   c                    &   g | ]}t j d d |f ddqS NrH  rI  rf  rg  rh  r1   r2   ru   %	     & )r+   r   r   rO  ravelr   r  r  )rO   rS  rT  r1   rh  r2   test_circmean_axis	  s   
z TestCircFuncs.test_circmean_axisc                    r]  )NrE  r^  rb  rH  rI  dy=rA   r   rc  c                    rd  re  r   rQ  rg  rh  r1   r2   ru   2	  ri  z3TestCircFuncs.test_circvar_axis.<locals>.<listcomp>r   c                    rj  rk  rp  rg  rh  r1   r2   ru   6	  rl  )r+   r   r   rQ  rm  r   r  r  )rO   rU  rV  r1   rh  r2   test_circvar_axis(	     
zTestCircFuncs.test_circvar_axisc                    r]  )NrE  r^  rb  rH  rI  ro  rA   r   rc  c                    rd  re  r   rR  rg  rh  r1   r2   ru   C	  ri  z3TestCircFuncs.test_circstd_axis.<locals>.<listcomp>r   c                    rj  rk  rs  rg  rh  r1   r2   ru   G	  rl  )r+   r   r   rR  rm  r   r  r  )rO   rW  rX  r1   rh  r2   test_circstd_axis9	  rr  zTestCircFuncs.test_circstd_axisc                 C   s"   g d}t ||dd|dd d S rD  )r   rK  r1   r1   r2   test_circfuncs_array_likeJ	  s   z'TestCircFuncs.test_circfuncs_array_likerL  c                 C   s   t t|g  d S r*   )r   r+   r  )rO   rL  r1   r1   r2   r4  R	  r   zTestCircFuncs.test_emptyc                 C   s.   ddddddt jg}tt ||dd d S )	NrF  r{   r   r  r   rG  rH  rI  r+   r   r   r  rO   rL  r]   r1   r1   r2   test_nan_propagateW	  s   z TestCircFuncs.test_nan_propagateg5:v@gfPb?rv  gee?gOQOv?g#q@gauy+@c                 C   s   t g dddddddt jgdt jt jt jt jt jt jgg}| D ]-}||d	|d
}|d u r8tt | q#t|d || dd tt |dd    q#d S )N)rF  r{   r   r  r   rG  r   r_  r;   rz   r`  r9   ra  r   rH  rc  r   r  rA   )r+   r   r   keysr   r  r   r   rO   rL  r  r]   r  outr1   r1   r2   test_nan_propagate_array]	  s   

z&TestCircFuncs.test_nan_propagate_arraygrQvv@)g     v@rP  rM  g     8v@r  g     u@g5@v?g2>Bjv@gG2?)gaF+?z?r}  gA ~@碲?r~  rM  r   gX҄?gJPm$@)gj  @竃 ?r  g)h;@{ ?r  gun2@g+fa @c                 C   s   t ddddddt jgddd	d
ddt jgt jt jt jt jt jt jt jgg}| D ].}||dd|d}|d u rBt||| dd q*t|d d || dd tt |d  q*d S )NrF  r{   r   r  r   rG  r_  r;   rz   r`  r9   ra  rH  omit)r  
nan_policyr  r  rA   r   )r+   r   r   ry  r   r   r  rz  r1   r1   r2   test_nan_omit_arrayr	  s   z!TestCircFuncs.test_nan_omit_arrayc                 C   s0   ddddddt jg}t||ddd	|d
d d S )NrF  r{   r   r  r   rG  rH  r  r  r  r  rA   )r+   r   r   rK  r1   r1   r2   test_nan_omit	  s   
zTestCircFuncs.test_nan_omitc                 C   s2   t jt jt jt jt jg}tt ||dd d S )Nr  r  rv  rw  r1   r1   r2   test_nan_omit_all	  s   zTestCircFuncs.test_nan_omit_allc                 C   sj   t t jt jt jt jt jgt jt jt jt jt jgg}||ddd}tt |  tt|dk d S )Nr  r   )r  r  r   )r+   r   r   r   r  r   r   )rO   rL  r]   r{  r1   r1   r2   test_nan_omit_all_axis	  s   z$TestCircFuncs.test_nan_omit_all_axisr]   rF  r{   r  r   rG  r_  r;   rz   r`  r9   ra  c                 C      t t||ddd d S )NrH  raiser  rV   rW   rw  r1   r1   r2   test_nan_raise	     zTestCircFuncs.test_nan_raisec                 C   r  )NrH  foobarr  r  rw  r1   r1   r2   test_bad_nan_policy	  r  z!TestCircFuncs.test_bad_nan_policyc                 C   s$   d}|}t |}t||dd d S )Nr   r   rA   )r   rO  r   )rO   r]   rS  rT  r1   r1   r2   test_circmean_scalar	  s   
z"TestCircFuncs.test_circmean_scalarc                 C   s@   t tdddtjtj }t|tjk  t|tj k d S )Nr   r   r_  )r   rO  r+   r\   rP  r   rF  r1   r1   r2   test_circmean_range	  s   z!TestCircFuncs.test_circmean_rangec                 C   sZ   t jddgdd}ttj|ddd ttj|dddd	d
 ttj|dddd	d
 d S )N   r   r;  r   rM  rI  g     @e@gLwqA?r  rA   g)#p_4@)r+   r   r	   r   rO  r   rQ  rR  r   r1   r1   r2   test_circfuncs_uint8	  s   z"TestCircFuncs.test_circfuncs_uint8N)$rb   rc   rd   r*  rH  rI  r   rO  rQ  rR  rM  rY  r+   rQ   rR   rS   r\  rn  rq  rt  ru  r4  rx  r   r|  r   rP  r  r  r  r  r  r  r  r  r  r1   r1   r1   r2   rC    s    








	



rC  c                   @   s~   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zdd Zejdddgdd ZdS )TestMedianTestc                 C   s   t ttjg d d S ru  rV   rW   r   median_testrX   r1   r1   r2   test_bad_n_samples	  rx  z!TestMedianTest.test_bad_n_samplesc                 C   s   t ttjg g d d S ru  r  rX   r1   r1   r2   r  	  r3  z TestMedianTest.test_empty_samplec                 C   s&   t ttjg dg dddgdd d S )N)r   r   r   r   )r   r   r   r   r   r   tiesr  rX   r1   r1   r2   test_empty_when_ties_ignored	  s   
z+TestMedianTest.test_empty_when_ties_ignoredc                 C   s8   t ttjg dg d t ttjg dg ddd d S )N)r   r   r   abover  r  rX   r1   r1   r2   test_empty_contingency_row	  s   
z)TestMedianTest.test_empty_contingency_rowc                 C       t ttjg dddgdd d S )Nrw   rz   r{   r  r  r  rX   r1   r1   r2   test_bad_ties	     
zTestMedianTest.test_bad_tiesc                 C   r  )Nrw   rz   r{   r  r  r  rX   r1   r1   r2   r  	  r  z"TestMedianTest.test_bad_nan_policyc                 C   r  )Nrw   rz   r{   r  )r  )rV   ra  r   r  rX   r1   r1   r2   rb  	  r  zTestMedianTest.test_bad_keywordc                 C   s\   g d}g d}t ||\}}}}t|d t|ddgddgg t|d t|d d S )Nrw   rL  r   r   r   )r   r  r	   r   )rO   r]   r  rf  ra  medtblr1   r1   r2   test_simple	  s   

zTestMedianTest.test_simplec                 C   s   g d}ddg}g d}t |||\}}}}t|d t|g dg dg t j|||dd\}}}}t|d t|g dg d	g t j|||d
d\}}}}t|d t|g dg d	g d S )Nr  r{   r8   )r;   r<   r9   )r   r   rx   )rz   r   r   r   r  )rz   r   r   r  )r   r   rx   )r   r  r	   )rO   r]   r  zrf  ra  r  r  r1   r1   r2   test_ties_options	
  s   


z TestMedianTest.test_ties_optionsc                 C   s   ddt jg}g d}tj||dd}tj||dd\}}}}t|t jt jt jd f t|d t|d t|d	 t|t d
dgddgg tttj||dd d S )Nr   r   )rz   r{   r8   	propagater  r  g     ?gإvao?rN  r   r  )	r+   r   r   r  r	   r   r   rV   rW   )rO   r]   r  mt1r   ra  r  r  r1   r1   r2   test_nan_policy_options
  s   


z&TestMedianTest.test_nan_policy_optionsc                 C   s2  g d}g d}t ||\}}}}t|d t|ddgddgg t |\}}}	}
t|| t|| t j||dd\}}}}t|d t|ddgddgg t j|dd\}}}	}
t|| t|| t j||dd	\}}}}t|d t|ddgddgg t j|dd	\}}}	}
t|| t|| d S )
Nry   )r   rz   r8   r<   rz   r   r   r   )lambda_Fr[  )r   r  r	   chi2_contingencyr   )rO   r]   r  rf  ra  r  r  exp_statexp_pdofer1   r1   r2   rT   *
  s(   







zTestMedianTest.test_basicr[  FTc                 C   s>   g d}g d}t j|||d}t|j|j|j|jf| d S )Nrw   r  )r   r  r	   rL   r   r[  table)rO   r[  r]   r  r_   r1   r1   r2   r  I
  s   zTestMedianTest.test_resultN)rb   rc   rd   r  r  r  r  r  r  rb  r  r  r  rT   r*  rH  rI  r  r1   r1   r1   r2   r  	  s    r  c                   @   sz   e Zd Zdd Zejdej d ejd gdfddej gdfgdd	 Z	d
d Z
dd Zdd Zdd Zdd ZdS )TestDirectionalStatsc                 C   s   t t g d }t t g d }t jt |t | t |t | t |fdd}t|}|j}t 	|d}t g d}t
|| d S )N)	g33333su@g      O@g33333sB@g      ;@g     pv@g@g333333I@gYv@rQ  )	gfffffP@g,Q@gfffffQ@gfffffT@g     S@g     @R@g33333SQ@gffffffM@g33333I@r   r  rz   )gHP?g&S:g[<)r+   deg2radr   stackcossinr   directional_statsmean_directionr  r   )rO   declinclrP   dirstatsdirectional_meanmean_roundedreference_meanr1   r1   r2   "test_directional_stats_correctnessT
  s   
z7TestDirectionalStats.test_directional_stats_correctnesszangles, refr   r   r   r   c                 C   sL   t |r||}tjt|t|gdd}dt|j }t|| d S )Nr   r  )	callabler+   r  r  r  r   r  mean_resultant_lengthr   )rO   anglesr   rP   r_   r1   r1   r2   'test_directional_stats_2d_special_casesh
  s
   z<TestDirectionalStats.test_directional_stats_2d_special_casesc           
      C   s   t jd}dt j |d }t jt |t |fdd}t|}|j	}t 
|d |d }|dt j  }t|}t|| d|j }t|}	t||	 d S )Nl	   }7Qh0P4u&4M r   )r<  r   r  r   )r+   r>  r?  rP  r  r  r  r   r  r  arctan2rO  r   r  rQ  )
rO   rB  testdatatestdata_vectorr  r  directional_mean_anglerO  directional_varcircular_varr1   r1   r2   test_directional_stats_2ds
  s$   





z.TestDirectionalStats.test_directional_stats_2dc                 C   sf   t g dg dg}t |d}t g dg dgg dg dgg}tj|dd}t||j d S )N+Pz?r  r   )r        r   )r   r   r   r   )r   r   r   r   r  )r+   r   tiler   r  r   r  )rO   rP   
full_arrayr  r  r1   r1   r2    test_directional_mean_higher_dim
  s   

z5TestDirectionalStats.test_directional_mean_higher_dimc                 C   sN   g dg dg}t |}t|}t|}t|j|j t|j|j d S )Nr  )r  r  r   )r+   r  r   r  r   r  r  )rO   rP   
data_arrayr_   r   r1   r1   r2   )test_directional_stats_list_ndarray_input
  s   


z>TestDirectionalStats.test_directional_stats_list_ndarray_inputc                 C   sR   t d}d}tjtt|d t| W d    d S 1 s"w   Y  d S )N)r{   zIsamples must at least be two-dimensional. Instead samples has shape: (5,)r(  )	r+   r   r*  r   rW   reescaper   r  )rO   rP   r+  r1   r1   r2   test_directional_stats_1d_error
  s
   
"z4TestDirectionalStats.test_directional_stats_1d_errorc                 C   sh   t g dg dg}tj|dd}|t jj|ddd }tj|dd}t|j|j t|j|j d S )Nr  )g+Pz?g      r   T)	normalizer   )r  keepdimsF)	r+   r   r   r  r"  r   r   r  r  )rO   rP   r_   normalized_datar   r1   r1   r2    test_directional_stats_normalize
  s   
z5TestDirectionalStats.test_directional_stats_normalizeN)rb   rc   rd   r  r*  rH  rI  r+   rP  r  r  r  r  r  r  r1   r1   r1   r2   r  R
  s    

r  c                   @   sp   e Zd Zdd Zdd Zejdg ddfg dd	fgd
d Zdd Z	ejdg ddd Z
dd ZdS )TestFDRControlc                 C   s  d}t jt|d tg d W d    n1 sw   Y  t jt|d tg d W d    n1 s9w   Y  t jt|d tddtjg W d    n1 sYw   Y  d}t jt|d tjg dd	d
 W d    n1 s{w   Y  d}t jt|d tjg ddd W d    n1 sw   Y  t jt|d tjg ddd W d    d S 1 sw   Y  d S )Nz.`ps` must include only numbers between 0 and 1r(  )r   r  r@  )r  r@  r   r  r@  zUnrecognized `method` 'YAK')r  r@  rg   YAKr  z#`axis` must be an integer or `None`r!  r  r  )r*  r   rW   r   false_discovery_controlr+   r   )rO   r+  r1   r1   r2   test_input_validation
  s*   "z$TestFDRControl.test_input_validationc                 C   s,   g d}t |}g d}t||dd d S )N)r  g;On?g~jt?gI+?g&1?gB`"۹?gv?S㥛?);On?r  I+?r  gˡE?gHzG?r  r  r  r  )r   r  r   )rO   psr_   r   r1   r1   r2   test_against_TileStats
  s   
z%TestFDRControl.test_against_TileStatscase)
go?gQ]
Y?ghD@?u?D?gn@t?g=T	[ے?g=!s?r  r  g	L?bh)
glߢ?gq?g.0?R]D?g)<Gh?g$͂?g2-p?r  r  g]:?byc                 C   sT   |\}}t jd}tjjddd|d}|d |d< tj||d}t||d	d
 d S )N   <U) r  r  r   r=  r;   rx   r  r   r  )r+   r>  r?  r   
loguniformr.   r  r   )rO   r  r   r  rB  r  r_   r1   r1   r2   test_against_R
  s   zTestFDRControl.test_against_Rc                 C   sJ   t jd}tjjddd|d}tj|d d}t| }t|| d S Nr  r  r  )rx   rz   r{   r=  r  )	r+   r>  r?  r   r  r.   r  rm  r	   )rO   rB  r  r_   r   r1   r1   r2   test_axis_None
  s
   zTestFDRControl.test_axis_Noner  )r   r   r   c                 C   sL   t jd}tjjddd|d}tj||d}t tj||}t|| d S r  )	r+   r>  r?  r   r  r.   r  apply_along_axisr	   )rO   r  rB  r  r_   r   r1   r1   r2   	test_axis
  s
   zTestFDRControl.test_axisc                 C   s8   t tdgdg t tdd t tg g  d S )Nr  )r   r   r  rX   r1   r1   r2   r  
  s   zTestFDRControl.test_edge_casesN)rb   rc   rd   r  r  r*  rH  rI  r  r  r  r  r1   r1   r1   r2   r  
  s"    

r  )Ur|   r  numpyr+   numpy.randomr   numpy.testingr   r   r   r   r   r   r	   r
   r*  r   rV   r  scipyr   r   scipy.stats._morestatsr   r   r   common_testsr   
_hypotestsr   r   scipy.stats._binomtestr   scipy.stats._distr_paramsr   dict
matplotlibrcParamsmatplotlib.pyplotpyplotr%  r7  	Exceptionr&  r'  r(  r)  r*  r+  r,  r-  r.  r/  r3   r6   rf   r   r   rK  r  r$  r?  rg  r  r  r  r  r  r8  rz  r  r  r  r  r  r  r  r  r  r  r>  rC  r  r  r  r1   r1   r1   r2   <module>   s   (
!{  I E= WG .^ U6%(3uH# e d