o
    GeS                      @   sV  d dl Zd dlmZmZ d dlmZmZ d dlZG dd dZ	G dd dZ
g g dg fg d	g fg d
g fg dg fg dg fdgddgfdgd	dgfdgd
dgfdgddgfdgddgfg ddg dfg dd	g dfg dd
g dfg ddg dfg ddg dfg ddg dfg dd	g dfg dd
g dfg ddg dfg ddg dfg ddg dfg dd	g dfg dd
g dfg ddg dfg ddg dfg ddg dfg dd	g dfg dd
g dfg ddg d fg ddg d!fd"gd# dedd$fR Zd%d& ZdS )'    N)assert_equalassert_array_equal)rankdata
tiecorrectc                   @   s4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )TestTieCorrectc                 C   s&   t jg t jd}t|}t|d dS )z9An empty array requires no correction, should return 1.0.dtype      ?Nnparrayfloat64r   r   selfranksc r   e/var/www/html/SPIRAL.web.tool/spiral_venv/lib/python3.10/site-packages/scipy/stats/tests/test_rank.py
test_empty
   s   zTestTieCorrect.test_emptyc                 C   s(   t jdgt jd}t|}t|d dS )z;A single element requires no correction, should return 1.0.r	   r   Nr
   r   r   r   r   test_one   s   zTestTieCorrect.test_onec                 C   s<   t d}t|}t|d t d}t|}t|d dS )z*Arrays with no ties require no correction.       @r	         @N)r   aranger   r   r   r   r   r   test_no_correction   s   


z!TestTieCorrect.test_no_correctionc                 C   s,  t g d}t|}d}|j}d|d | |d |   }t|| t g d}t|}d}|j}d|d | |d |   }t|| t g d}t|}d}|j}d|d | |d |   }t|| t g d}t|}d}d}|j}d|d | |d |  |d |   }t|| d	S )
z8Check a few basic examples of the tie correction factor.)r	         @r   r   r	      )      ?r   r   )r	   r   r   r   r   )r   r         @r   r   N)r   r   r   sizer   )r   r   r   TNexpectedT1T2r   r   r   
test_basic   s2   


(zTestTieCorrect.test_basicc                 C   sZ   d\}}t t ||}|j}tt|}t|d||d |  t|d |    d S )N)i     r	   r   )r   repeatr   r   r   r   r   float)r   ntiekanoutr   r   r   test_overflowB   s
   .zTestTieCorrect.test_overflowN)__name__
__module____qualname__r   r   r   r$   r-   r   r   r   r   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	g dZ
ejgejgd  Zejdddgejdee
edd ZdS )TestRankDatac                 C   sN   t jg td}t|}t|t jg t jd tg }t|t jg t jd dS )z0stats.rankdata([]) should return an empty array.r   Nr   r   intr   r   r   )r   r*   rr   r   r   r   L   s
   zTestRankData.test_emptyc                 C   sX   dg}t j|td}t|}t|t jdgt jd t|}t|t jdgt jd dS )z/Check stats.rankdata with an array of length 1.d   r   r	   Nr2   )r   datar*   r4   r   r   r   r   T   s   zTestRankData.test_onec                 C   s  g d}t jg dt jd}t j|td}t|}t|| t|}t|| g d}t jg dt jd}t j|td}t|}t|| t|}t|| g d}t jg dt jd}t j|td}t|}t|| t|}t|| |dd	}t|}t|| d
S )zBasic tests of stats.rankdata.)r5   
   2   )r   r	   r   r   )(   r7      r7   r8   )r   r   r   r         @)   r<   r<   r7   r7   r7   )r;   r;   r;   r   r   r      r   N)r   r   r   r3   r   r   reshape)r   r6   r!   r*   r4   a2dr   r   r   r$   ]   s0   





zTestRankData.test_basicc                    s   dd dd  fdd} fdd}dd }t  |||d	fd
d}g d}|tj|d |tj|dd tjg ddd}|tj|dd d S )Nc                        fdd D S )Nc                    &   g | ] d t  fddD  qS )   c                 3       | ]}| k V  qd S Nr   .0ijr   r   	<genexpr>|       XTestRankData.test_rankdata_object_string.<locals>.<lambda>.<locals>.<listcomp>.<genexpr>sumrF   r*   rH   r   
<listcomp>|      & NTestRankData.test_rankdata_object_string.<locals>.<lambda>.<locals>.<listcomp>r   rP   r   rP   r   <lambda>|       z:TestRankData.test_rankdata_object_string.<locals>.<lambda>c                    r@   )Nc                    s"   g | ] t  fd dD qS )c                 3   s    | ]}| kV  qd S rD   r   rE   rH   r   r   rJ   }   rK   rL   rM   rO   rP   rH   r   rQ   }   s   " rS   r   rP   r   rP   r   rT   }   rU   c                    s    dd t | D S )Nc                 S   s   g | ]\}}||fqS r   r   )rF   rG   xr   r   r   rQ   ~   s    rS   )	enumeraterP   )min_rankr   r   rT   ~   s    c                    s   dd t |  | D S )Nc                 S   s   g | ]
\}}|| d  qS )r   r   )rF   rG   rI   r   r   r   rQ      s    zRTestRankData.test_rankdata_object_string.<locals>.average_rank.<locals>.<listcomp>)ziprP   )max_rankrX   r   r   average_rank   s   z>TestRankData.test_rankdata_object_string.<locals>.average_rankc                    s   t |   fdd| D S )Nc                    rA   )rB   c                 3   rC   rD   r   rE   rH   r   r   rJ      rK   zZTestRankData.test_rankdata_object_string.<locals>.dense_rank.<locals>.<listcomp>.<genexpr>rM   rO   brH   r   rQ      rR   zPTestRankData.test_rankdata_object_string.<locals>.dense_rank.<locals>.<listcomp>)r   uniquerP   r   r\   r   
dense_rank   s   
z<TestRankData.test_rankdata_object_string.<locals>.dense_rank)minmaxordinalaveragedensec                    s,   dD ]}t | |d}t| | |  qd S )N)r`   ra   rd   rb   rc   methodr   r   )r*   rf   r,   )rankfr   r   check_ranks   s   z=TestRankData.test_rankdata_object_string.<locals>.check_ranks)	foobarquxxyzabcefgaceqweqaz   object)r   rB   r=   gX9v@r   gT㥛 	@r   )dictr   randomchoiceastyper   )r   ordinal_rankr[   r_   ri   valr   )rZ   rX   rh   r   test_rankdata_object_string{   s   z(TestRankData.test_rankdata_object_stringc                 C   s   t jddgt jd}t|}t|ddg t jddgt jd}t|}t|ddg t jddgt jd}t|}t|ddg d S )Nl            l           r   r	   r   l)r   r   uint64r   r   int64)r   r6   r4   r   r   r   test_large_int   s   zTestRankData.test_large_intc                 C   sD   dD ]}t j|td}t|}d|d  }t||| d|  qd S )N)i'  i i@B r   g      ?rB   ztest failed with n=%d)r   onesr3   r   r   )r   r+   r6   r4   expected_rankr   r   r   test_big_tie   s   
zTestRankData.test_big_tiec                 C   s`   g dg dg}g dg dg}t |dd}t|| g dg dg}t |d	d}t|| d S )
N)r   r=   rB   )   r=   r=   )r	   r   r	   )r   r   r   r   )axisr	   r   r   )r   r   r   rB   rg   )r   r6   	expected0r0	expected1r1r   r   r   	test_axis   s   
zTestRankData.test_axis)rc   r`   ra   rd   rb   r   r   r   rB   zmethod, dtypec                 C   s8   d}t |}t|||d}t|j| t|j| d S )N)r   r   )rf   r   )r   zerosr   r   shaper   )r   r   rf   r   r   r6   r4   r   r   r   test_size_0_axis   s
   
zTestRankData.test_size_0_axisN)r.   r/   r0   r   r   r$   r{   r~   r   r   methodsr   r   int_dtypespytestmarkparametrizerY   r   r   r   r   r   r1   J   s    	r1   rc   r`   ra   rd   rb   r5   r	   )r5   r5   r5   )r   r   r   )r	   r	   r	   )r   r   r   )r	   r   r   )r5   ,  rs   r   )r5   rs   r   rs   )r	   r   r   r   )r	   r   r   r   )r	   r   r   r   )r	   r   r   r   )r	   r   r   r   )r5   rs   r   rs   r5   )r         @r;   r   r   )r	   r   r;   r   r	   )r   r   r;   r   r   )r	   r   r   r   r	   )r	   r   r;   r   r   r7   r:   g      ?@c                  C   s*   t D ]\} }}t| |d}t|| qd S )Nre   )_casesr   r   )valuesrf   r!   r4   r   r   r   
test_cases   s   r   )numpyr   numpy.testingr   r   scipy.statsr   r   r   r   r1   r   r   r   r   r   r   r   <module>   s    Bz	
 !"#$&*