o
    tBh                    @   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mZ d dlmZ d dlmZ dd	lmZ d
dlmZ d dlmZ zd dlZdejd< d dlm Z! dZ"W n e#yr   dZ"Y nw g dZ$g dZ%g dZ&g dZ'g dZ(g dZ)g dZ*g dZ+g dZ,g dZ-G dd dZ.G dd dZ/G dd  d Z0G d!d" d"Z1G d#d$ d$Z2G d%d& d&Z3G d'd( d(Z4G d)d* d*Z5G d+d, d,Z6G d-d. d.e6Z7G d/d0 d0Z8G d1d2 d2Z9G d3d4 d4Z:G d5d6 d6Z;G d7d8 d8Z<G d9d: d:Z=G d;d< d<Z>G d=d> d>Z?G d?d@ d@Z@G dAdB dBZAg dCZBG dDdE dEZCG dFdG dGZDG dHdI dIZEG dJdK dKZFG dLdM dMZGG dNdO dOZHG dPdQ dQZIG dRdS dSZJdS )T    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   )check_named_results   )_get_wilcoxon_distr)_binary_search_for_binom_tst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                   @   $   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std rE   w/var/www/html/riverr-enterprise-integrations-main/venv/lib/python3.10/site-packages/scipy/stats/tests/test_morestats.py
test_basic0   s   
zTestBayes_mvs.test_basicc                 C      t ttjg  d S N)assert_raises
ValueErrorr   r<   r@   rE   rE   rF   test_empty_inputB      zTestBayes_mvs.test_empty_inputc                 C   s0   t d}d}t|}|D ]}t|| qd S )N   )r=   r>   )nparanger   r<   r   )r@   x
attributesresirE   rE   rF   test_result_attributesE   s   

z$TestBayes_mvs.test_result_attributesN)__name__
__module____qualname__rG   rM   rV   rE   rE   rE   rF   r'   /   s    r'   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 )Nr(   r/   g?r0   r1   r2   r4   r5   r6   r7   r8   r9   r;   )r   mvsdistr   rB   r   intervalr?   rE   rE   rF   rG   O   s   
zTestMvsdist.test_basicc                 C   rH   rI   rJ   rK   r   r\   rL   rE   rE   rF   rM   ^   rN   zTestMvsdist.test_empty_inputc                 C      dg}t ttj| d S Nr   r^   r@   rA   rE   rE   rF   test_bad_arga      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 rE   rB   .0rR   rE   rE   rF   
<listcomp>k       z*TestMvsdist.test_warns.<locals>.<listcomp>r   r      c                 S   re   rE   rf   rg   rE   rE   rF   ri   l   rj   r   r   rl         )warningscatch_warningssimplefilterRuntimeWarningr   r\   rL   rE   rE   rF   
test_warnsf   s
   
"zTestMvsdist.test_warnsN)rW   rX   rY   rG   rM   rb   rt   rE   rE   rE   rF   r[   N   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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=
ף@   ?r)   r9      Ɍ?)(\?=
ףp=?\(\@ffffff@\(\?(\?Gz@(\ףp=
@皙?)\(?RQr   p=
ף?@
ףp=
?{Gz?\(\@{Gz@Q@IddY?Dio?rl   ro   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?rn   )r   shapiror   r=   pvaluenormrvs)
r@   x1wpwshapiro_testx2x3x4
W_expected
p_expectedrE   rE   rF   rG   p   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)
rv   rw   rx   ry   rz   r{   r|   r}   r~   r   )
r   r   r   r   r   r   r   r   r   r   r   r)   r9   r   )
r   r   r   r   r   r   r   r   r   r   )
r   r   r   r   r   r   r   r   r   r   r   r   rl   )r   r   r   r=   r   )r@   r   r   r   r   r   rE   rE   rF   test_2d   s$   

zTestShapiro.test_2dc                 C   s&   t ttjg  t ttjg g g g d S rI   rJ   rK   r   r   rL   rE   rE   rF   rM         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)rJ   rK   r   r   rP   arrayobjectrL   rE   rE   rF   test_not_enough_values   s   "z"TestShapiro.test_not_enough_valuesc                 C   r_   r`   r   r@   rR   rE   rE   rF   rb      rc   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 )Nr4   r*   r   )	rP   rQ   nanr   r   r	   r=   r   r   )r@   rR   r   r   r   rE   rE   rF   test_nan_input   s   



zTestShapiro.test_nan_inputN)	rW   rX   rY   rG   r   rM   r   rb   r   rE   rE   rE   rF   ru   o   s    +ru   c                   @   D   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 )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   rP   onesr   )r@   rsr   r   AcritsigvrE   rE   rF   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   rP   errstater   )r@   r   r   r   r   r   r   rE   rE   rF   
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   gumbelr   r   r   )rP   r   r   r   lengumbel_lfitlogcdflogsfrQ   rB   r   )r@   r   a2r   r   nxbarsr   r   rU   expected_a2rE   rE   rF   test_gumbel   s   
*zTestAnderson.test_gumbelc                 C      t ttjdgdd d S Nr   plate_of_shrimpdist)rJ   rK   r   r   rL   rE   rE   rF   rb         zTestAnderson.test_bad_argc                 C   s0   t d}|jdd}t|}d}t|| d S )Nr   r   r   r=   critical_valuessignificance_level)r   r   r   r   r   )r@   r   rR   rT   rS   rE   rE   rF   rV      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   )	r@   r   rR   A1crit1sig1A2crit2sig2rE   rE   rF   test_gumbel_l  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   rP   r   r   r   r   r   )
r@   r   r   r   r   r   r   r   r   r   rE   rE   rF   test_gumbel_r  s   
zTestAnderson.test_gumbel_rN)
rW   rX   rY   r   r   r   rb   rV   r   r   rE   rE   rE   rF   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@r   g      F@g     F@r   ffffffA@r   g33333A@gB@gfffffB@gD@gffffffE@Fmidrankg@rl   gCl?gvq-?g,Ԛ?gl@gPk	@r   ro   rn   gJ4a?Mb0?atolrP   r   r   anderson_ksampr   r   r   r@   t1t2t3t4TktmprE   rE   rF   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 )Nr   r   r   r   Tr   gQ@rl   r   r   ro   rn   gMb`?r  r  r  r  rE   rE   rF   test_example1b1  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   rO   )      !      i  r8   :   %   r   A   r*      i     $      v   
"         r#  C   9   >   r,      r!  Z   r   <      =   1   r8      8      O   T   ,   ;   r  r        i6  L      r3     r&        F   e   r<  J   r%  0   r  i  r+   r=     r  i  r4     7   i@  r/  h         /            r  r9  i  W   r,   x   r8   r&  rI     G   rJ  rB  *   r0  ro   r+   rO     rl   r8   rQ  rS  r8   rS     r)  r   rT  4   _   a   3   rS  rn      r#     D   M   P   r   rT  j      R   6   r"     .   o   '   ?   r#     r#     r.  r   r3  f   H   r'  rf  rl   rO         r1  X   rd  ro   ro   r  r'        rX     r9  r.   r8   	g  r*   r+   i  i[  rl   rF  r   i  r      ro   i  #   r+   i  r#  r   r,   b   ro   U   [   +      rl   r;  rk     r8   r%  rb      r$  r4        rC  r8   r}  B   r,  r!  Fr   gM
@rl   g^I+?g~:?g+e?gH.!@gcZB>h@r   ro   rn   	^)p?r  r  r   r  r   r   r   r@   r  r  r	  r
  t5t6t7t8t9t10t11t12t13t14r  r  r  rE   rE   rF   test_example2aA  s0   z TestAndersonKSamp.test_example2ac                 C   r  )Nr  r  r   r(  r:  r?  rD  rM  rW  rj  rs  ru  rx  r~  Tr   g'1Z
@rl   r  r   ro   rn   r  r  r  r  r  rE   rE   rF   test_example2ba  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messageg     @D@Fr   g=
ףpD@rl   MbP?gGzD@zp-value capped      ?g?ܵrn         ?g?ܵ      @g!u?gףp=
׳?g{Gzt?r   r  r3   r)   gBiq?g鷯?      '@g:M@g46<Rv?gMb@?g      +@g6>W[1@g1ZGUS?-C6?)
rP   linspacer
   filterUserWarningr   r  r   r	   r   )r@   r   supr   _r  rE   rE   rF   test_R_kSamples  sJ   !



z!TestAndersonKSamp.test_R_kSamplesc                 C   s   t ttjtd d S Nro   rJ   rK   r   r  rP   r   rL   rE   rE   rF   test_not_enough_samples  r   z)TestAndersonKSamp.test_not_enough_samplesc                 C   s"   t ttjtdtdf d S r  r  rL   rE   rE   rF   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  rL   rE   rE   rF   test_empty_sample     z#TestAndersonKSamp.test_empty_samplec                 C   s:   g d}t g d}tj||fdd}d}t|| d S )Nr   r   Fr   r   )rP   r   r   r  r   )r@   r  r  rT   rS   rE   rE   rF   rV     s
   z(TestAndersonKSamp.test_result_attributesN)rW   rX   rY   r  r  r  r  r  r  r  r  rV   rE   rE   rE   rF   r     s      P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	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   rl   rl   rn   rl   r   r)   r   r)   r   rn   r   %Ties preclude use of exact statistic.g     7@rS  g8#oG?)r
   r  r  r   ansarir   )r@   rR   yr  WpvalrE   rE   rF   
test_small  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re  k   r   c   rk  r_  m   l   rF  r  r>  `   rX  rk  r  q   t   r  n   ry  r  r  r_  ry  i   g   r  r  rF  r   r  r  r  rF  r   r  r  r  r_  r  r  g     0g@rS  g|:?)rP   r   r
   r  r  r   r  r   )r@   ramsayparekhr  r  r  rE   rE   rF   test_approx  s   

zTestAnsari.test_approxc                 C   s4   t g dg d\}}t|dd t|dd d S )Nr   r   rl   rn   )rO   ro   r0  r-   r   r+   r4   rS  g?r,   )r   r  r   )r@   r  r  rE   rE   rF   
test_exact  s   zTestAnsari.test_exactc                 C   s(   t ttjg dg t ttjdgg  d S r`   rJ   rK   r   r  rL   rE   rE   rF   rb     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  r=   r   )r
   r  r  r   r  r   )r@   rR   r  r  rT   rS   rE   rE   rF   rV     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  ro   r)   r,   r-   z!'alternative' must be 'two-sided'matchfooalternativer  )r@   r   r   r  rE   rE   rF   test_bad_alternative  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   ro   r   rO   r0  r5  )r  g      !@      #@g      %@r  g      )@lessr  greaterffffff?皙?r   -q=r  r   )r   r  r   r   pmfr   r   )
r@   r   r   r=   r  pval_lpval_gprobpval_l_reversepval_g_reverserE   rE   rF   test_alternative_exact  s   z!TestAnsari.test_alternative_exactzx, y, alternative, expectedr  r  r  g@A?r  rk   )rn   ro   r)   r,   r-   gG$I?g%I$I?rm   )r)   r,   r-   c                 C   s$   t j|||dj}t||dd d S )Nr  r  r  )r   r  r   r   )r@   rR   r  r  expectedr  rE   rE   rF   test_alternative_exact_with_R3  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   ro   r   {   r   r   r   r  r  r  r   r  r  r   r*  g      ?r   )r   r   r   r  r   r   )r@   r   r   r  r  r  rE   rE   rF   test_alternative_approxR  s   z"TestAnsari.test_alternative_approxN)rW   rX   rY   r  r  r  rb   rV   r  r  pytestmarkparametrizer  r  rE   rE   rE   rF   r    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   )r@   argsTr  rE   rE   rF   	test_dataj     zTestBartlett.test_datac                 C      t ttjdg d S r`   )rJ   rK   r   r  rL   rE   rE   rF   rb   q     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   r  r   r@   r  rT   rS   rE   rE   rF   rV   u     
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 rI   )r  r  r  r  r  r  r  r  r  r  r	   rP   r   r   r  )r@   r  rE   rE   rF   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   rl   rn   )rP   r   rJ   rK   r   r  r  r   rE   rE   rF   test_1d_input     zTestBartlett.test_1d_inputN)rW   rX   rY   r  rb   rV   r  r  rE   rE   rE   rF   r  h  s    r  c                   @   r   )
TestLevenec              
   C   r  )NgYpK?r,   gI{]?)r  r  r  r  r  r  r  r  r  r  r   levener   )r@   r  r  r  rE   rE   rF   r    r  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 )NrB   centertrimmedr   r  proportiontocut)r   r  r  r  r  r   )r@   W1pval1W2pval2rE   rE   rF   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   r        @r  g      @r  g      @g      i@  r        ?r  r   r   rB   r  )rP   randomseedpermutationr   r  r   )
r@   rR   r  r   W0pval0r  r  r	  r
  rE   rE   rF   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   rB  r  rl   rB   r  median)rP   r  r  r  r  r   r  r   )r@   rR   r   r  r  r  r	  r
  rE   rE   rF   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   rB  皙?)portiontocut)rP   r  rJ   	TypeErrorr   r  r   rE   rE   rF   test_bad_keyword  r   zTestLevene.test_bad_keywordc                 C   r  Nr   r   rB  trimr  )rP   r  rJ   rK   r   r  r   rE   rE   rF   test_bad_center_value  r   z TestLevene.test_bad_center_valuec                 C   r  r`   )rJ   rK   r   r  rL   rE   rE   rF   test_too_few_args  s   zTestLevene.test_too_few_argsc              
   C   r  r  )r  r  r  r  r  r  r  r  r  r  r   r  r   r  rE   rE   rF   rV     r  z!TestLevene.test_result_attributesc                 C   r  r  )rP   r   rJ   rK   r   r  r  r   rE   rE   rF   r    r  zTestLevene.test_1d_inputN)rW   rX   rY   r  r  r  r  r$  r'  r(  rV   r  rE   rE   rE   rF   r     s    	
r   c                   @   s`   e Zd ZdZeej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 )
TestBinomPzTests for stats.binom_test.c                 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-"^?rS  r  i  g`w?i           ?r  gz?binom_test_funcr   )r@   r  rE   rE   rF   r    s   zTestBinomP.test_datac                 C   s   t t| jg d d S Nrk   rJ   rK   r/  rL   rE   rE   rF   test_bad_len_x     zTestBinomP.test_bad_len_xc                 C   s(   t t| jdg t t| jdgdd d S )Nr   r   r   r1  rL   rE   rE   rF   
test_bad_n  s   zTestBinomP.test_bad_nc                 C   s   t t| jddgdd d S )Nr   r  r-  r1  rL   rE   rE   rF   
test_bad_p  s   
zTestBinomP.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 )NrY     gUUUUUU?r  r  r  g-$Ǻl?r  gs0,j.?	two-sidedgJfJ"f?r.  r@   rT   rE   rE   rF   test_alternatives  s   

zTestBinomP.test_alternativesl        z32-bit does not overflowreasonc                 C   s   t t| jddtjjd d S )Nr  r)   r-  )rJ   OverflowErrorr/  sys
float_infominrL   rE   rE   rF   test_boost_overflow_raises  s   z%TestBinomP.test_boost_overflow_raisesN)rW   rX   rY   __doc__staticmethodr   
binom_testr/  r  r2  r5  r6  r;  r  r  skipifr?  maxsizerB  rE   rE   rE   rF   r)    s    

r)  c                   @   s   e Zd ZdZedddZdS )TestBinomTestPzJ
    Tests for stats.binomtest as a replacement for stats.binom_test.
    Nr  r9  c                 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.r8  )	rP   
atleast_1dastypeint_r   rK   r   	binomtestr   )rR   r   r  r  resultrE   rE   rF   r/    s   
zTestBinomTestP.binom_test_func)Nr  r9  )rW   rX   rY   rC  rD  r/  rE   rE   rE   rF   rH    s    rH  c                   @   sx  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dAS )B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   r2   iΙ gX7?i	ϙ g=?iϙ gybw?r6   r   rL  r   r   r@   r3   rT   rE   rE   rF   test_two_sided_pvalues1  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 )NrO  r*   rB  rP  r   r  gX-h?r2   rn   g㤶y?rS  g	y}?r,   g+,?rl   r   r  g      ?r   g?g{Gz?g333333?g<,?rQ  rR  rE   rE   rF   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 )NrO  i  i  r  r   r2   rl   rI  gAL &W?r.   rd  g!YB?rO   r3  gE]t?r,   r)   rS  rQ  rR  rE   rE   rF   test_edge_cases3  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  rl   r   r   axisc                    s   t j|   S rI   r   binomr  r   rT  rE   rF   <lambda>O  rN   z>TestBinomTest.test_binary_srch_for_binom_tst.<locals>.<lambda>r6   r2   r,   c                    s   t j|  S rI   rY  r[  rT  rE   rF   r\  [     )rP   rQ   ceilr   rZ  r  sumr   intr   floor)r@   krU   dy1ixy2rE   rT  rF   test_binary_srch_for_binom_tstB  s(    z,TestBinomTest.test_binary_srch_for_binom_tstz"alternative, pval, ci_low, ci_high))r  g$W?r   g?)r  g}?gGI~?r   )r9  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 )
Nr0  r   r  r   r  r  r  r2   皙?r  confidence_level	r   rL  r   r   r	   proportion_estimateproportion_cilowhighr@   r  r  ci_lowci_highrT   cirE   rE   rF   test_confidence_intervals1b  
   z(TestBinomTest.test_confidence_intervals1))r  ga&+w?r   gny?)r  g5!Rx?g5%?r   )r9  g~th?gSuR&|?g|KnG?c                 C   rh  )
Nrl   r   rj  ri  ư>r2   gQ?r"   rk  rm  rr  rE   rE   rF   test_confidence_intervals2q  rw  z(TestBinomTest.test_confidence_intervals2zalternative, pval, ci_high))r  gլ?g̓A?)r  r   r   )r9  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  r8  rx  r2   r  rk  r   )r   rL  r   r   ro  r	   rp  rq  )r@   r  r  rt  rT   ru  rE   rE   rF   !test_confidence_interval_exact_k0  
   z/TestBinomTest.test_confidence_interval_exact_k0zalternative, pval, ci_low))r  r   r   )r  6>g@>}_?)r9  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  r8  rx  r2   r  rk  r   )r   rL  r   r   ro  r	   rq  rp  )r@   r  r  rs  rT   ru  rE   rE   rF   %test_confidence_interval_exact_k_is_n  r{  z3TestBinomTest.test_confidence_interval_exact_k_is_nz+k, alternative, corr, conf, ci_low, ci_high)rl   r9  Tr  g`?g?)rl   r9  Tr"   grt!?g/?)rl   r9  Fr  gVK5?g5dM?)rl   r9  Fr"   gXCVu^?g*r[?)rl   r  Tr  r   gP-V?)rl   r  Tr"   r   g^'?)rl   r  Fr  r   gS5a?)rl   r  Fr"   r   g?)rl   r  Tr  g!Շ`ⰸ?r   )rl   r  Tr"   g$x"?r   )rl   r  Fr  gW~=?r   )rl   r  Fr"   g	L?r   )r   r9  Tr  r   gX?)r   r9  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   r9  Tr  g7?r   )r   r9  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!  ri  wilsonccwilson)rl  methodrx  r2   )r   rL  ro  r   rp  rq  )
r@   rb  r  corrconfrs  rt  rT   r  ru  rE   rE   rF   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 )Nrn   rT  r  r   )r   rL  r	   rn  r   r:  rE   rE   rF   &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   rK   r   rL  )r@   rb  r   rE   rE   rF   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  rS  r   r  r  rL   rE   rE   rF   test_invalid_k_too_big  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 )	Nrl   r   r!  rT  zmust be in the intervalr  r   rk  r   rL  r  r   rK   ro  r:  rE   rE   rF   test_invalid_confidence_level     "z+TestBinomTest.test_invalid_confidence_levelc                 C   r  )	Nrl   r   r!  rT  zmethod must ber  zplate of shrimpr  r  r:  rE   rE   rF   test_invalid_ci_method  r  z$TestBinomTest.test_invalid_ci_methodN)rW   rX   rY   rC  rS  rU  rV  rg  r  r  r  rv  ry  rz  r}  r  r  r  r  r  r  rE   rE   rE   rF   rN    sr     


	
	


rN  c                   @   r   )TestFlignerc                 C   &   t d}tt||d dd d S )Nro   r   )gnf	@gDq?rS  )rP   rQ   r   r   flignerr@   r   rE   rE   rF   r    
   
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 )NrO  )rP   asarrayrandnr   tolist)gr   rE   rF   r\    s     z+TestFligner.test_trimmed1.<locals>.<lambda>rB   r  r  r   r  )	rP   r  r   r  r  r  r   r  r   )	r@   _perturbg1_g2_g3_Xsq1r  Xsq2r
  rE   r  rF   r    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 )
Nr  r  r  r  r  r   r   rB   r  )r   r  r   )r@   rR   r  r  r  r  r
  rE   rE   rF   r    s   

$
zTestFligner.test_trimmed2c                 C   r  r   )rP   r  rJ   r#  r   r  r   rE   rE   rF   r$    r   zTestFligner.test_bad_keywordc                 C   r  r%  )rP   r  rJ   rK   r   r  r   rE   rE   rF   r'    r   z!TestFligner.test_bad_center_valuec                 C   r  r`   )rJ   rK   r   r  rL   rE   rE   rF   test_bad_num_args  r  zTestFligner.test_bad_num_argsc                 C   s.   t d}tt jt jft||d g  d S )Nro   r   )rP   rQ   r	   r   r   r  r   rE   rE   rF   r    s   
$zTestFligner.test_empty_argN)
rW   rX   rY   r  r  r  r$  r'  r  r  rE   rE   rE   rF   r    s    r  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	e
jjdde
jdg ddd ZdS )TestMoodc                 C   r  )Nro   r   )g8!giT?rS  )rP   rQ   r   r   moodr  rE   rE   rF   	test_mood$  r  zTestMood.test_moodc                 C   s`   t jd t jdd}t jdd}t||\}}t||\}}t||g| |g d S )Nr  r   r   rO   )rP   r  r  r  r   r  r   )r@   r   r   z1p1z2p2rE   rE   rF   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   )rO   r   rW  g4g^}E	?)rP   r   shaper   r   r  )r@   r   r   rE   rE   rF   test_mood_with_axis_none5  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 )Nro   r  r   rO   r   rW  )	rP   r  r  r  r   r  ranger   	transpose)r@   nyr   r   	z_vectestpval_vectestjrU   rE   rE   rF   test_mood_2dQ  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   ro   r)   r  rl   rW  r   r   r   )rP   r  r  r  r  r   r  r   )r@   r  r   r   rX  r  r  axes_idxrU   r  slice1slice2rE   rE   rF   test_mood_3dh  s2   

zTestMood.test_mood_3dc                 C   s   t ttjdgg  d S r`   )rJ   rK   r   r  rL   rE   rE   rF   test_mood_bad_arg     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   r,  r   )r   r   g      ?r9  r  r  r  Hz>r  r   r   zalternative must be...r  z	ekki-ekki)rP   r  r  r   r   r   r  r   r  r   rK   )	r@   rR   r  stat1r  stat2r  stat3p3rE   rE   rF   test_mood_alternative  s   "zTestMood.test_mood_alternativez,SciPy needs tie correction like R (gh-13730)r<  zalternative, expected))r9  )!}?g-?)r  )g~?g4?)r  )r  g-?c                 C   s.   g d}g d}t j|||d}t|| d S )Nr  r  r  )r   r  r   )r@   r  r  rR   r  rT   rE   rE   rF   test_mood_alternative_against_R  s   z(TestMood.test_mood_alternative_against_RN)rW   rX   rY   r  r  r  r  r  r  r  r  r  xfailr  r  rE   rE   rE   rF   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 )	Nr0  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   rP   sort)r@   rR   osmosrosm_expectedrT   res_fitres_fit_expectedrE   rE   rF   rG     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   rE   )r   r   r   r  r   )
r@   rR   osm1osr1osm2osr2osm3osr3r  r  rE   rE   rF   test_sparams_keyword  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 )Nr0  r  r  Ft)rl   )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)r@   qrE   rE   rF   r    r]  z7TestProbplot.test_dist_keyword.<locals>.custom_dist.ppfN)rW   rX   rY   rC  r  rE   rE   rE   rF   custom_dist  s    r  )r   r  )r   r   )	r   r   r   r  r  r   rJ   rK   AttributeError)r@   rR   r  r  r  r  r  rE   rE   rF   test_dist_keyword  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 )
Nre  rl   r   t r  plotF)r   r  r   )pltfigureadd_subplotr   r  r   r  closer   r   r   )
r@   figrR   res1fitres1res2fitres2res3res4axrE   rE   rF   test_plot_kwarg  s$   
@




zTestProbplot.test_plot_kwargc                 C   r   r   )rJ   rK   r   r  rL   rE   rE   rF   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  rP   r   r   rL   rE   rE   rF   
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   )rP   r   r	   r   r  r   r   rL   rE   rE   rF   test_array_of_size_one  s   "z#TestProbplot.test_array_of_size_oneN)rW   rX   rY   rG   r  r  r  r  rF  have_matplotlibr  r  r  r  rE   rE   rE   rF   r    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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)rJ   rK   r   wilcoxonrL   rE   rE   rF   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 )Nr0  wilcoxapproxr  prattzsplit)g     @Z@r   )rP   rQ   rJ   rK   r   r  r	   r   rE   rE   rF   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   rl   ro   Sample size too smallr  r  r  zero_methodr  )r   g|N?)r
   r  r  r   r  r   )r@   rR   r  r  rT   rE   rE   rF   
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   rl   r   r   rl   r   r   r   r   r   r  r  r  r  r  )r   r  )r@   arrr  rE   rE   rF   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   rn   rT  rO   r-   rn   ro   r   r   ro   c                 S   s   g | ]	\}}|g| qS rE   rE   )rh   ur   rE   rE   rF   ri   :  s    z7TestWilcoxon.test_accuracy_wilcoxon.<locals>.<listcomp>r  r  r  i  g.i?r  i  gFYv]Uj?r  iG  g%^ Ez?rO  r  r  rn        y   rh     r     p         w      r	        r  s   r  r[     F
correctionr  r!  g^҂<?rx  r2   Tg:b+?)r  rP   concatenatezipzerosr   r   r  r   r   r	   )r@   freqnumsrR   r  r  r  rE   rE   rF   test_accuracy_wilcoxon7  s*   








z#TestWilcoxon.test_accuracy_wilcoxonc                 C   s@   t g d}t g d}tj||ddd}d}t|| d S )Nr  r
  Fr  r  r  )rP   r   r   r  r   )r@   rR   r  rT   rS   rE   rE   rF   test_wilcoxon_result_attributesT  s
   z,TestWilcoxon.test_wilcoxon_result_attributesc                 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   rx  r2   Tr  gc2_?)r   r  r	   r   )r@   statr  
expected_prE   rE   rF   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)
}   r  r;     r   r  r   r  r      )
r  z   r  rO  r   |   r     r!     r  r  r  r  r  r  rC  g13}?r)   r9   T)r  r  r  g7a%?r  g?gڏoJT?)r
   r  r  r   r  r	   r   )r@   rR   r  r  r   r  rE   rE   rF   test_onesidedn  s<   	



zTestWilcoxon.test_onesidedc                 C   sL   t ddD ]}t|}t||d  d d t| tt|d|  qd S )Nr   r8  r   )r  r   r	   r   r_  )r@   r   cntrE   rE   rF   test_exact_basic  s
   zTestWilcoxon.test_exact_basicc                 C   s6  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 t dd }tttj|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Qr9  exactr&  g֔  ?r)   r9   r  gYJ?r  g*  ?r   r0  r  r   g(?g(?gS?r8  r   r  )rP   r   r   r  r   rQ   rJ   rK   )r@   rR   r  r  r  rc  rE   rE   rF   test_exact_pval  s$   zTestWilcoxon.test_exact_pvalrR   )r   r   rl   )r   r   r  ro   )r   r   rl   r  r  ir,   r-   c                 C   s@   t |\}}t|}||dk  }t|| t|d d S Nr   r   )r   r  rP   r   r_  r	   )r@   rR   r   r  wtruerE   rE   rF   test_exact_p_1  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   r5  r  r   r+  r  r.   zExact p-value calculationr  r  r   rC  )rP   rQ   r	   r   r  r
   r  r  )r@   rR   r  rc  r  r   r  rE   rE   rF   	test_auto  s   zTestWilcoxon.test_autoN)rW   rX   rY   r  r  r  r   r  r  r  r'  r)  r,  r  r  r  r0  r1  rE   rE   rE   rF   r    s"    	
	&
r  c                   @   rZ   )
	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 rE   )r   kstat)rh   r   rA   rE   rF   ri     s    z>TestKstat.test_moments_normal_distribution.<locals>.<listcomp>r  )gHQ,?glZ)rI?g88}?gcJN?r  r2   r   )momentr   rl   r   g{Gz?{Gz?r  )rP   r  r  r  r   r   r5  )r@   momentsr  m1m2m3rE   r4  rF    test_moments_normal_distribution  s   "z*TestKstat.test_moments_normal_distributionc                 C   rH   rI   )rJ   rK   r   r3  rL   rE   rE   rF   rM     rN   zTestKstat.test_empty_inputc                 C   *   t d}t j|d< tt|t j d S Nr4   r)   rP   rQ   r   r	   r   r3  ra   rE   rE   rF   r        

zTestKstat.test_nan_inputc                 C   s*   t d}dD ]}tttj||d qd S )Nr   )r   g/$@r4  )rP   rQ   rJ   rK   r   r3  r@   rA   r   rE   rE   rF   test_kstat_bad_arg  s   
zTestKstat.test_kstat_bad_argN)rW   rX   rY   r;  rM   r   rA  rE   rE   rE   rF   r2    s
    r2  c                   @   r&   )TestKstatVarc                 C   rH   rI   rJ   rK   r   kstatvarrL   rE   rE   rF   rM     rN   zTestKstatVar.test_empty_inputc                 C   r<  r=  r>  ra   rE   rE   rF   r     r?  zTestKstatVar.test_nan_inputc                 C   s    dg}d}t ttj||d d S )Nr   r   r4  rC  r@  rE   rE   rF   rb     s   zTestKstatVar.test_bad_argN)rW   rX   rY   rM   r   rb   rE   rE   rE   rF   rB    s    rB  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 jjddddd | _d S Nro     r  r  r   loggammar   rR   rL   rE   rE   rF   setup_method  r  zTestPpccPlot.setup_methodc                 C   H   d}t j| jdd|d\}}g d}t|tjdd|d t|| d S )Nro   r   N)g9%	?gxM _?g";9?g®XZ?g?num)r   	ppcc_plotrR   r   rP   r  )r@   rO  svalsppccppcc_expectedrE   rE   rF   rG   	  
   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 )NrM  r   tukeylambdar   g#B;r2   )r   rR  rR   rW  r   )r@   svals1ppcc1svals2ppcc2svals3ppcc3rE   rE   rF   	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 Nre  ir0  r  )r  r  r  r   rR  rR   delaxesr  r@   r  r  rE   rE   rF   r       


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   rk   r   r   )rJ   rK   r   rR  rR   rL   rE   rE   rF   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^  rP  r   )r   rR  r   rP   r  r  float)r@   rS  rT  rE   rE   rF   r  2  s   zTestPpccPlot.test_emptyN)rW   rX   rY   rK  rG   r^  r  r  rF  r  r  rd  r  rE   rE   rE   rF   rE    s    
rE  c                   @   rZ   )
TestPpccMaxc                 C   s   dg}t ttj|dd d S r   )rJ   rK   r   ppcc_maxra   rE   rE   rF   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,   r9   )r   rW  r   r   rg  r   rE   rE   rF   test_ppcc_max_basicA  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 )Nri  r   r  rj  rk  r   rl  rW  r   rm  ro   r9   )r   rW  r   rg  r   )r@   rR   max1max2max3rE   rE   rF   r^  F  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 )Nri  r   r  rj  rk  r   rl  )r   r   r  )brack)r   r   rm  r,   r9   )r   r   )r   rW  r   rJ   rK   rg  r   r   rE   rE   rF   
test_brackT  s   
zTestPpccMax.test_brackN)rW   rX   rY   rh  rn  r^  rs  rE   rE   rE   rF   rf  ;  s
    rf  c                   @   r  )TestBoxcox_llfc                 C   sT   t jjdddd}d}t ||}|j d tt| d  }t	|| d S )Nrj  r   1  r   r   r   r   r  r   )
r   r   r   
boxcox_llfr   rP   logr_  rD   r   )r@   rR   lmbdallfllf_expectedrE   rE   rF   rG   b  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   ru  rv  r   r  r2   )r   r   r   rw  listr   r@   rR   ry  rz  llf2rE   rE   rF   test_array_likei  
   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   rw  rP   vstackr  r   r  rE   rE   rF   test_2d_inputp  s
   zTestBoxcox_llf.test_2d_inputc                 C      t ttdg  d S r`   )r   rP   isnanr   rw  rL   rE   rE   rF   r  z     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   rw  r   )r@   rA   rz  rE   rE   rF   test_gh_6873}  s   zTestBoxcox_llf.test_gh_6873N)rW   rX   rY   rG   r  r  r  r  rE   rE   rE   rF   rt  `  s    
rt  (  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 jj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| d S )	Nro   r   r  r  r   ry  r   r   )r   rJ  r   boxcoxr   rP   rx  r~  )r@   rR   xtrE   rE   rF   test_fixed_lmbda  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   r9   )r   r   r   r  r   )r@   ry  rR   x_invr  maxlogrE   rE   rF   test_lmbda_None  s
   zTestBoxcox.test_lmbda_Nonec                 C   s   t jd}tjjdd|dd }tj|dd\}}}t|ddg tj|d	d\}}}t|d
dg tjjdd|dd }tj|dd\}}}t|ddg tj|dd\}}}t|ddg d S )Nr  ro   r   r  r,  )alphagx@gv[%@r  gxˉEk?g8Bhk @r,   rH  rO   r  g\?gx̡&@r   g7(U@gp:sY@)rP   r  r   r   rJ  r   r  r   )r@   rngrR   r  r]   rE   rE   rF   
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   )rP   r   rJ   rK   r   r  r   rE   rE   rF   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  rL   rE   rE   rF   r    r   zTestBoxcox.test_emptyc                 C   s    t t\}}t|ddd d S )NgsHjdrh㈵>r2   r   r  _boxcox_datar   )r@   r  lamrE   rE   rF   r    s   zTestBoxcox.test_gh_6873bounds)r   r   )r   r   )r   r*  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  rE   rF   	optimizer     zBTestBoxcox.test_bounded_optimizer_within_bounds.<locals>.optimizerry  r  r   r   )r   r  r  )r@   r  r  r  ry  rE   r  rF   $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  rE   rF   r    s   zPTestBoxcox.test_bounded_optimizer_against_unbounded_optimizer.<locals>.optimizerr  r   r  )r@   r  ry  r  lmbda_boundedrE   r  rF   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   rK   r   r  r  r@   r  rE   rE   rF   $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 r`   rE   r  rE   rE   rF   r  #  s   zCTestBoxcox.test_bad_optimizer_value_raises_error.<locals>.optimizerz;`optimizer` must return an object containing the optimal...r  r  r  )r@   r  r  rE   rE   rF   %test_bad_optimizer_value_raises_error  s
   "z0TestBoxcox.test_bad_optimizer_value_raises_errorN)rW   rX   rY   r  r  r  r  r  r  r  r  r  r  r  r  r  rE   rE   rE   rF   r    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   rF  Nro   r   r  r  rI  rL   rE   rE   rF   rK  ,  r  zTestBoxcoxNormmax.setup_methodc                 C      t | j}t|ddd d S )Nm?rx  r2   r   boxcox_normmaxrR   r   r@   r  rE   rE   rF   test_pearsonr/     zTestBoxcoxNormmax.test_pearsonrc                 C   s<   t j| jdd}t|ddd t | j\}}t|| d S )Nmler  	.!?rx  r2   )r   r  rR   r   r  )r@   r  r  maxlog_boxcoxrE   rE   rF   test_mle3  s   zTestBoxcoxNormmax.test_mlec                 C   s&   t j| jdd}t|ddgdd d S )Nr   r  r  r  rx  r2   r  )r@   
maxlog_allrE   rE   rF   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  rE   rF   r  C  r  zITestBoxcoxNormmax.test_bounded_optimizer_within_bounds.<locals>.optimizer)r  r  r   r   )r   r  rR   rP   r   )r@   r  r  r  r  rE   r  rF   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 )
Nro   r6  i  c                   @   s   e Zd ZdS )z?TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.MyResultN)rW   rX   rY   rE   rE   rE   rF   MyResultR  s    r  c                    s6   g }D ]	}| | | q  }t| |_|S rI   )appendrP   argminrR   )r  objsry  rT   r  lmbda_rangerE   rF   r  U  s   z@TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.optimizer)r  r  )r   r  rR   rP   roundr  r   )r@   ry  lmbda_roundedr  lmbda2rE   r  rF   test_user_defined_optimizerL  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)rr  r  z,`brack` must be None if `optimizer` is givenr  )g       r  )r   r  r   r  rR   r  r   rK   r  rE   rE   rF   2test_user_defined_optimizer_and_brack_raises_errorb  s   
"zDTestBoxcoxNormmax.test_user_defined_optimizer_and_brack_raises_errorN)rW   rX   rY   rK  r  r  r  r  r  r  r  r  r  rE   rE   rE   rF   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   rF  rG  rI  rL   rE   rE   rF   rK  t  r  zTestBoxcoxNormplot.setup_methodc                 C   rL  )Nro   rM  r   rN  )g!X4}?g΁i?g-	15?gW$?g\ܱ{?rP  )r   boxcox_normplotrR   r   rP   r  )r@   rO  lmbdasrT  rU  rE   rE   rF   rG   w  rV  zTestBoxcoxNormplot.test_basicr  r<  c                 C   r_  r`  )r  r  r  r   r  rR   ra  r  rb  rE   rE   rF   r    rc  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   )rJ   rK   r   r  rR   rL   rE   rE   rF   rd    s   z&TestBoxcoxNormplot.test_invalid_inputsc                 C   s   t tg ddjdk d S r.  )r   r   r  r   rL   rE   rE   rF   r    r  zTestBoxcoxNormplot.test_emptyN)rW   rX   rY   rK  rG   r  r  rF  r  r  rd  r  rE   rE   rE   rF   r  s  s    
r  c                   @   r&   )TestYeojohnson_llfc                 C   r|  )Nr   r   ru  rv  r   r  r2   )r   r   r   yeojohnson_llfr~  r   r  rE   rE   rF   r    r  z"TestYeojohnson_llf.test_array_likec                 C   r  r}  )r   r   r   r  rP   r  r  r   r  rE   rE   rF   r    s
   z TestYeojohnson_llf.test_2d_inputc                 C   r  r`   )r   rP   r  r   r  rL   rE   rE   rF   r    r  zTestYeojohnson_llf.test_emptyN)rW   rX   rY   r  r  r  rE   rE   rE   rF   r    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S )TestYeojohnsonc                 C   s  t jd}tjj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jjdd|dd }t |dk siJ tj|d	d}t|t 	| d   tj|dd}t|| tj|d
d}t|d| d  d  tjj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  ro   r   r  r   r   r  r   r   rl   )
rP   r  r   r   rJ  r   r   
yeojohnsonr   rx  )r@   r  rR   r  posnegrE   rE   rF   r    sL   


 $zTestYeojohnson.test_fixed_lmbdary  )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   )rP   r  r  r   absspacingexppower)rR   ry  r  r  rE   rE   rF   _inverse_transform  s   $
z:TestYeojohnson.test_lmbda_None.<locals>._inverse_transformi N  rk  r   r   )r   r   r   r6  r  r   r9   )rP   r  r  normalr   r  r   r   linalgr   rB   rD   )r@   ry  r  	n_samplesrR   r  r  r  rE   rE   rF   r    s   
zTestYeojohnson.test_lmbda_Nonec                 C   r  r  )r   r   r  r  rL   rE   rE   rF   r     r   zTestYeojohnson.test_emptyc                 C   sD   t jjdddd}t |\}}t t|\}}t||dd d S )Nr   r   ru  rv  r  r2   )r   r   r   r  r~  r   )r@   rR   xt1r  xt2rE   rE   rF   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 )Nr)   r   z>Yeo-Johnson transformation is not defined for complex numbers.r  )rP   rQ   r  r   rK   r   r  )r@   r   rR   err_msgrE   rE   rF   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  r2   )rP   rQ   float64r   r  r   )r@   r   x_intx_floatxt_int	lmbda_intxt_floatlmbda_floatrE   rE   rF   test_input_dtype_integer  s   z'TestYeojohnson.test_input_dtype_integerN)rW   rX   rY   r  r  r  r  r  r  r  rP   	complex64
complex128r  int8uint8int16int32r  rE   rE   rE   rF   r    s    /
$
r  c                   @   r&   )TestYeojohnsonNormmaxc                 C   rF  r  rI  rL   rE   rE   rF   rK    r  z"TestYeojohnsonNormmax.setup_methodc                 C   r  )Ng?rx  r2   )r   yeojohnson_normmaxrR   r   r  rE   rE   rF   r    r  zTestYeojohnsonNormmax.test_mlec                 C   s*   g d}t |}tj|dddsJ d S )N)gffffff@g r   r  gffffff?g333333@r  gffffff@g?g@r  r  g"@r  g      gzG?r  r  )r   r  rP   allclose)r@   rR   ry  rE   rE   rF   test_darwin_example#  s   
z)TestYeojohnsonNormmax.test_darwin_exampleN)rW   rX   rY   rK  r  r  rE   rE   rE   rF   r    s    r  c                   @   s4  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dd(gdfej	d)eg d*ed+d,gdfgd-d. Zejdejdfejdfej	dfgd/d0 Zejdejejej	gd1d2 Zejdejejej	gd3d4 Zejd5d6d7d8d9d:d;ejged6d7d8d9d:d;ejgd<d=d>d?ejd@dAgggejdejejej	gdBdC Zejd5d6d7d8d9d:d;ejged6d7d8d9d:d;ejgd<d=d>d?ejd@dAgggejdejejej	gdDdE ZdFdG Z dHdI Z!dJdK Z"dLS )MTestCircFuncsztest_func,expectedgA:wv?gk~_ЀBE@gE|2@c                 C   s(   t g d}t||dd|dd d S Nc  ro   r   rt  r   ^  h  rq  r  r2   rP   r   r   r@   	test_funcr  rR   rE   rE   rF   test_circfuncs-  s   zTestCircFuncs.test_circfuncsc                 C   s~   t g d}| }tj|dd}t||dd | }tj|dd}t||dd | }tj	|dd}t||dd d S )N)r0  rB  r'  r#     g     4@g3333333@r  r  r  r2   r  )
rP   r   rB   r   circmeanr   rC   circvarrD   circstd)r@   rR   M1M2V1V2S1S2rE   rE   rF   test_circfuncs_small5  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  r  )r@   r   
numpy_funcrR   circstatr  rE   rE   rF   test_circfuncs_closeC  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 )Nr  _  r,   rn   `  r*   ]  ie  r*   r-   if  rn   id  r  r  r1   r2   r   rq  rX  c                       g | ]}t j | d dqS r  r  r   r  rh   rU   rR   rE   rF   ri   W      z4TestCircFuncs.test_circmean_axis.<locals>.<listcomp>r   c                    &   g | ]}t j d d |f ddqS Nr  r  r  r  r  rE   rF   ri   [     & )rP   r   r   r  ravelr   r  r  )r@   r  r  rE   r  rF   test_circmean_axisN  s   
z TestCircFuncs.test_circmean_axisc                    r  )Nr  r  r  r  r  dy=r2   r   r  c                    r  r  r   r  r  r  rE   rF   ri   h  r  z3TestCircFuncs.test_circvar_axis.<locals>.<listcomp>r   c                    r  r  r#  r  r  rE   rF   ri   l  r  )rP   r   r   r  r   r   r  r  )r@   r  r	  rE   r  rF   test_circvar_axis^     
zTestCircFuncs.test_circvar_axisc                    r  )Nr  r  r  r  r  r"  r2   r   r  c                    r  r  r   r  r  r  rE   rF   ri   y  r  z3TestCircFuncs.test_circstd_axis.<locals>.<listcomp>r   c                    r  r  r&  r  r  rE   rF   ri   }  r  )rP   r   r   r  r   r   r  r  )r@   r
  r  rE   r  rF   test_circstd_axiso  r%  zTestCircFuncs.test_circstd_axisc                 C   s"   g d}t ||dd|dd d S r  )r   r  rE   rE   rF   test_circfuncs_array_like  s   z'TestCircFuncs.test_circfuncs_array_liker   c                 C   s   t t|g  d S rI   )r   rP   r  )r@   r   rE   rE   rF   r    r  zTestCircFuncs.test_emptyc                 C   s.   ddddddt jg}tt ||dd d S )	Nr  ro   r   rt  r   r  r  r  rP   r   r   r  r@   r   rR   rE   rE   rF   test_nan_propagate  s   z TestCircFuncs.test_nan_propagateg5:v@gfPb?r.  gڴV0@gLԿAB@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)r  ro   r   rt  r   r  r   r  r,   rn   r  r*   r  r   r  r  r   r  r2   )rP   r   r   keysr   r  r   r   r@   r   r  rR   rX  outrE   rE   rF   test_nan_propagate_array  s   
z&TestCircFuncs.test_nan_propagate_arraygrQvv@)g     v@r  r  g     8v@r  g     u@g5@v?g2>Bjv@gtUuYK@)gR @zw<5 ?r0  g~/ϔ(@K ?r1  g sQP@gJPm$@)gj  @竃 ?r2  g)h;@{ ?r3  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 )Nr  ro   r   rt  r   r  r  r,   rn   r  r*   r  r  omit)rq  
nan_policyrX  r  r2   r   )rP   r   r   r,  r   r   r  r-  rE   rE   rF   test_nan_omit_array  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 )Nr  ro   r   rt  r   r  r  r4  rq  r5  r  r2   )rP   r   r   r  rE   rE   rF   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 )Nr4  r5  r)  r*  rE   rE   rF   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 )Nr4  r   )r5  rX  r   )rP   r   r   r   r  r   r   )r@   r   rR   r.  rE   rE   rF   test_nan_omit_all_axis  s   z$TestCircFuncs.test_nan_omit_all_axisrR   r  ro   r   rt  r   r  r  r,   rn   r  r*   r  c                 C      t t||ddd d S )Nr  raiser7  rJ   rK   r*  rE   rE   rF   test_nan_raise     zTestCircFuncs.test_nan_raisec                 C   r<  )Nr  foobarr7  r>  r*  rE   rE   rF   test_bad_nan_policy  r@  z!TestCircFuncs.test_bad_nan_policyc                 C   s$   d}|}t |}t||dd d S )Nr   r  r2   )r   r  r   )r@   rR   r  r  rE   rE   rF   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   r  rP   rQ   pir   )r@   mrE   rE   rF   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   r  r  g     @e@g.}}VW{@r  r2   g)#p_4@)rP   r   r	   r   r  r   r  r  r   rE   rE   rF   test_circfuncs_unit8  s   z"TestCircFuncs.test_circfuncs_unit8N)#rW   rX   rY   r  r  r  r   r  r  r  r  r  rP   rB   rC   rD   r  r!  r$  r'  r(  r  r+  r   r/  r   r6  r8  r:  r;  r?  rB  rC  rF  rH  rE   rE   rE   rF   r  ,  s    












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d Z
dd Zdd Zdd ZdS )TestMedianTestc                 C   s   t ttjg d d S r0  rJ   rK   r   median_testrL   rE   rE   rF   test_bad_n_samples	  r3  z!TestMedianTest.test_bad_n_samplesc                 C   s   t ttjg g d d S r0  rJ  rL   rE   rE   rF   r  
	  r  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   tiesrJ  rL   rE   rE   rF   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   aboverM  rJ  rL   rE   rE   rF   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 )Nrk   rn   ro   r  rM  rJ  rL   rE   rE   rF   test_bad_ties"	     
zTestMedianTest.test_bad_tiesc                 C   rR  )Nrk   rn   ro   rA  r9  rJ  rL   rE   rE   rF   rB  &	  s    z"TestMedianTest.test_bad_nan_policyc                 C   rR  )Nrk   rn   ro   r  )r  )rJ   r#  r   rK  rL   rE   rE   rF   r$  )	  rT  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 )Nrk   r  r   r   r   )r   rK  r	   r   )r@   rR   r  r  r  medtblrE   rE   rF   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  ro   r)   )r,   r-   r*   )r   r   rl   )rn   r   r   r   rM  )rn   r   r   rP  )r   r   rl   )r   rK  r	   )r@   rR   r  zr  r  rE  rV  rE   rE   rF   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   )rn   ro   r)   	propagater9  r4  g     ?gإvao?r  r   r=  )	rP   r   r   rK  r	   r   r   rJ   rK   )r@   rR   r  mt1r   r  rE  r  rE   rE   rF   test_nan_policy_optionsP	  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 )
Nrm   )r   rn   r)   r-   rn   r   r   r   )lambda_F)r  )r   rK  r	   chi2_contingencyr   )r@   rR   r  r  r  rE  rV  exp_statexp_pdoferE   rE   rF   rG   ]	  s(   







zTestMedianTest.test_basicN)rW   rX   rY   rL  r  rO  rQ  rS  rB  r$  rW  rY  r\  rG   rE   rE   rE   rF   rI  	  s    rI  )Krp   r?  numpyrP   numpy.randomr   numpy.testingr   r   r   r   r   r   r	   r
   r  r   rJ   scipyr   r   scipy.stats._morestatsr   common_testsr   
_hypotestsr   scipy.stats._binomtestr   
matplotlibrcParamsmatplotlib.pyplotpyplotr  r  	Exceptionr  r  r  r  r  r  r  r  r  r  r'   r[   ru   r   r   r  r  r   r)  rH  rN  r  r  r  r  r2  rB  rE  rf  rt  r  r  r  r  r  r  r  r  rI  rE   rE   rE   rF   <module>   s   (
!ZV G E, SD ^ G6%(2qH#q Y