
    (ph                         S SK rS SKJr  S SKJs  Jr  S SKJ	r	J
r
JrJr  S rS rS rS rS rS r " S	 S
5      rS rg)    N)assert_equalassert_almost_equalassert_assert_allclosec                      [         R                  " S5      n U S-   n[        [        R                  " X5      S5        [         R
                  " SSSS9n[        [        R                  " X5      S5        g )N   
   r      )numgX?)nparanger   mscompare_medians_mslinspace)xyy2s      W/var/www/html/venv/lib/python3.13/site-packages/scipy/stats/tests/test_mstats_extras.pytest_compare_medians_msr   	   sT    
		!A	BA--a3Q7	Qr	"B--a4nE    c                     [         R                  " S5      n [        [        R                  " U 5      SSS9  [         R
                  " U 5      U l        SU R                  S S& [        [        R                  " U 5      SSS9  [        [        R                  " U SS	9R                  S
:H  5        [         R                  " S5      R                  S5      n[        [        R                  " USS9SS/5        [         R
                  " U5      Ul        SUR                  S S2S S 24'   [        [        R                  " USS9SS/5        g )N      g+=)rtolFr      Tvar      )r   r   r   axisr	      )
mar   r   r   hdmedian	make_maskmaskr   sizereshape)r   x2s     r   test_hdmedianr*      s    
		"ABKKNAE2\\!_AFAFF2AJBKKNAE2 BKKt$))Q./ 
2		w	'BBKK+b"X6ll2BGBGGBQBENBKK+aV4r   c                  N   [         R                  R                  S5        [         R                  R                  S5      n [        R
                  " U 5      n[        UR                  U R                  :H  5        [        R
                  " U SS/S9n[        UR                  S:H  5        g )Ni d   r         ?)pointsr   )	r   randomseedrandnr   rshr   shaper'   )r   ress     r   test_rshr5   &   sn    IINN6
		A
&&)C CII ! &&Ar7
#CCHHMr   c                  z    [         R                  " / SQ5      n [        [        R                  " U 5      / SQS5        g )N)M   W   X   r                  i  i(  i+  i2  ix  i  i    i  i3
  )gTK@gϽKF@gQh@r   )r#   arrayr   r   mjcidatas    r   	test_mjcirE   3   s*    88 < =D&CAFr   c                      [         R                  " / SQ5      n [        [        R                  " U S5      SS5        [        [        R                  " [        R                  " U S5      S5      SS/5        g )N)i!  i+  i.  i<  i?  i@  iB  iD  iR  i]  i{  i  i  i  r@   皙?g@r
   )rG   rG   gfffff@g̴@)	r#   rA   r   r   trimmed_meanr   r   roundtrimmed_mean_cirC   s    r   test_trimmed_mean_cirK   :   sU    88 2 3DS15!<",,T)<Q?!r   c                  6   [         R                  " S5      n [        [         R                  " [        R
                  " U 5      5      SS/S5        U R                  S5      R                  SS5      n[        [        R
                  " USS9/ S	Q/ S
Q/S5        [        [        R
                  " USS9U R                  S5      R                  SS5      5        SS/n [        R
                  " U 5      n[        [         R                  " U5      R                  5       5        g )Nr,   RBj8@]kSUR@r"   r   r   r    )rM   rM   rM   )rN   rN   rN   r
   r   )r   r   r   asarrayr   idealfourthsrepeatr(   r   isnanall)testtest_2D_results      r   test_idealfourthsrX   C   s    99S>D

2??4#89"9-a1kk!n$$R*Ga888::;= a8A..r!46q6Dood#GBHHW!!#$r   c                   .    \ rS rSr/ SQrS rS rS rSrg)TestQuantilesS   )dgxc%?gF4vE?g݌Y?gn?gJ]X?g'h?g|f2?g!}_z`?g3!?g[Q,?gW1y?g+?gZy7В?g?g|GA?gPw2?g0?g3;??g,7tZ?gxދ?gޠt$O?giL?g|rmq?gsl?g݋Y?g)چ?g
zS?g'	A?g9io?gOB?g5G~?g9 ?g*sD?g\?gN?g6??gƟ	'?g$	?gLw5?g2J(?g2^?g?]?gͬ?g8~n7?gxY?gE!?ga@'?gʇ3?gF~\?gB誽?gX/u?g/:L@{?gb9?gvY?g5XEs?gd"2?g-=X.?giE"?g/E?g4J3H?ǒ?g+[{?go?gbm?g\-?g'W?@?g0?tz?g=z5?g}v?gXPT?g<*R?gX?gO{F?g::2?gǺ?g/i?g6?gRfU?gu9?g!G[?gT=~?g?gAM07?g,jKV?g"ЄI?g΅?g\?gRSI?g	o?g2QJ)d~?g4n/?g+--C?gX?g|wiƘ?gttњc,?g{˻;Ib?gّ?gv+:?gR <e?c           	         U R                   n[        [        R                  " USS/5      SS/5        [        R                  " U/ SQ5      n[        U/ SQ5        [        R
                  " U5      R                  SS5      n[        R                  " U/ SQSS	9n[        US S 2S4   [        R                  " US S 2S4   / SQ5      5        [        US S 2S
4   [        R                  " US S 2S
4   / SQ5      5        [        R                  " U/ SQSSS9n[        US   [        R                  " US S 2S4   / SQSS95        [        US   [        R                  " US S 2S
4   / SQSS95        g )Ng        r-   r]   r\   )g      ?      ?g      ?)g]Ar4?g5$?i?gJJ5d?r	   r   r    rO   T)r!   r   ).r   r   ).rO   )rD   r   r   hdquantilesr   rA   r(   )selfrD   hdqs      r   test_hdquantilesTestQuantiles.test_hdquantilesi   s    yyBNN4R9(+6	8nnT"34C!IJxx~%%b,nnT/q9C!HbnnT!A#Y&OPC"Ir~~d1R4j'QRnnT/qTBCJNN4!9_N	PCKNN4":o4P	Rr   c                 P   [         R                  " U R                  5      n[        U R                  5      n[        R
                  " U R                  X"45      n[        R                  " [        R                  " U5      5      nX4   R                  X"S-
  5      n[         R                  " USS9n[        R                  " USS9nUS-
  U-  [        R                  " XV-
  S-  SS9-  S-  n[        X5        [        U/ SQ5        [         R                  " SS/5      n[        U/ SQ5        g )Nr
   r    r   r   r_   )giYk?gK4U}?gN5i	t?)r_   r_   r_   )r   hdquantiles_sdrD   lenr   broadcast_tological_noteyer(   r`   meansumr   )	ra   hd_std_errsnjdata	jselectorjdist
jdist_meanjstdtwo_data_pointss	            r   test_hdquantiles_sd!TestQuantiles.test_hdquantiles_sdz   s    
 ''		2 		N		A62NN266!9-	 ((aC0u1-WWU+
1a"&&%"4q!8qAABFK.K)JK++QF3O_=r   c                     [         R                  " U R                  5      u  p[        UR                  UR                  s=:H  =(       a
    S:H  5        g s  5        g )Nr   )r   mquantiles_cimjrD   r   r'   )ra   ci_lowerci_uppers      r   test_mquantiles_cimj"TestQuantiles.test_mquantiles_cimj   s<    //		:33!3434r    N)	__name__
__module____qualname____firstlineno__rD   rc   ru   r{   __static_attributes__r}   r   r   rZ   rZ   S   s    ID*R">.5r   rZ   c                      [         R                  R                  S5      n U R                  SS9n[        [        R
                  " U5      S5        [        [        R
                  " US5      S5        g )Nl   WE~-S    )r'   )gl?gz0JL?g?)g/?gL0?)r   r/   default_rngr   r   median_cihs)rngr   s     r   test_median_cihsr      sU     ))

 3
4C


ABNN1%'?@ BNN1c*,JKr   )numpyr   numpy.mar#   scipy.stats.mstatsstatsmstatsr   numpy.testingr   r   r   r   r   r*   r5   rE   rK   rX   rZ   r   r}   r   r   <module>r      sN       , ,F5(
G!% A5 A5HLr   