
    doi D              "          S SK rS SKrS SKJr  S SKJr  S SKJrJ	r	J
r
  S SKJrJr  S SKJr  S SKJrJr  S SKJr  \R&                  " S	5      r\R*                  R-                  \(       + S
S9/r\R*                  R0                  rS\SSS.S jr\" SSS9\" SSS9\" SSS9\R*                  R7                  S\" \R8                  0 5      \" \R:                  0 5      \" \R<                  SS05      /5      \R*                  R7                  SS S/5      S 5       5       5       5       5       r\" SSS9\" SSS9\" SSS9\R*                  R7                  S/ SQ5      \R*                  R7                  SSS /5      S! 5       5       5       5       5       r \" SSS9\" SSS9\" SS"S9\R*                  R7                  S\" \RB                  S#S05      \" \RD                  0 5      \" \RD                  S$S05      \" \RF                  0 5      \" \RF                  S$S05      \" \RH                  0 5      \" \RJ                  S%S05      \" \RJ                  S%S05      \" \RJ                  S%S&05      \" \RJ                  S%S'05      \" \RL                  S%S05      \" \RL                  S%S05      \" \RN                  0 5      \" \RP                  0 5      \" \RR                  0 5      \" \RT                  0 5      \" \RV                  0 5      \0 4\" \RX                  S(S)05      \" \RZ                  S(S)05      \" \R\                  S*S+05      \" \R^                  S,S+05      \" \R`                  S(S)05      \" \Rb                  S(S)05      /5      \R*                  R7                  S/ SQ5      S- 5       5       5       5       5       r2\" \Rf                  5      \" SSS9\" SSS9\" SS"S9\R*                  R7                  SS S/5      \R*                  R7                  S.0 /5      S/ 5       5       5       5       5       5       r4\" SSS9\" SSS9\" SS"S9\R*                  R7                  S0\" \Rj                  5      \" \Rl                  5      \" \Rn                  5      /5      \R*                  R7                  S/ SQ5      S1 5       5       5       5       5       r8\" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R7                  S4\" \Rr                  5      \" \Rt                  5      \" \Rv                  5      /5      \R*                  R7                  S/ SQ5      S5 5       5       5       5       5       5       r<\" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R{                  S65      \R*                  R7                  S4\" \R|                  5      \" \R~                  5      \" \R                  5      \" \R                  5      /5      \R*                  R7                  S/ SQ5      S7 5       5       5       5       5       5       5       rBS8 rC\" S9 S \CS:9SOS S;.S< jj5       rD\" \R                  \R                  5      \" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R7                  S=/ S>Q5      \R*                  R7                  S?S S/5      \R*                  R7                  S/ SQ5      S@ 5       5       5       5       5       5       5       5       rG\" \R                  5      \" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R7                  SA/ SBQ5      \R*                  R7                  S/ SQ5      SC 5       5       5       5       5       5       5       rI\" \R                  5      \" SSS9\" SSS9\" SS"S9\" S2S3S9SD 5       5       5       5       5       rK\R*                  R-                  \R                  " \R                  5      \R                  " SE5      :  SFS9SG 5       rN\" \R                  5      \" SSS9\" SSS9\" SS"S9\R*                  R{                  SH5      SI 5       5       5       5       5       rP\" \R                  5      \" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R7                  S\R                  0 4\R                  SJS 04\R                  SJS04/5      \R*                  R7                  S/ SQ5      SK 5       5       5       5       5       5       5       rS\" SSS9\" SSS9\" SS"S9\" S2S3S9\R*                  R7                  S4\" \R                  5      \" \R                  5      /5      SL 5       5       5       5       5       rV\" \R                  5      \" SSS9\" SSS9\R*                  R7                  SMSS /5      \R*                  R7                  S/ SQ5      SN 5       5       5       5       5       rXg)P    N)stats)version)xp_assert_closexp_assert_equal_length_nonmasked)make_xp_pytest_parammake_xp_test_case)SCIPY_ARRAY_API)_xp_mean_xp_var)_axis_nan_policy_factorymarrayzcspecial function dispatch to marray required for these tests is hidden behind SCIPY_ARRAY_API flag.)reasonfloat64)      l   !wN4 )dtypexpshapeseedc                j   [         R                  U5      n[        R                  R	                  U5      n/ / p[        U 5       H  n	UR                  US9n
UR                  S5      (       a  SU
-  S-  UR                  US9-  n
U
R                  U5      n
UR                  U
5        UR                  US9S:  nUR                  U5        M     / n/ n[        Xx5       HR  u  pUR                  UR                  XS95        UR                  5       nUR                  X'   UR                  U5        MT     X\U4$ )N)sizecomplex
   y              $@g      ?mask)r   _get_namespacenprandomdefault_rngrange
startswithstandard_normalastypeappendzipasarraycopynan)n_arraysr   r   r   r   mxprngdatasmasksidatar   marrays
nan_arraysarray	nan_arrays                   Y/var/www/html/land-ocr/venv/lib/python3.13/site-packages/scipy/stats/tests/test_marray.py
get_arraysr6      s   



#C
))


%Cr58_zzuz%I&&d7S=!4!4%!4!@@D{{5!Tzzuz%,T  GJ5(s{{5{45JJL	&&	)$	 ) ##    z
dask.arrayz/Arrays need `device` attribute: dask/dask#11711z	jax.numpyz"JAX doesn't allow item assignment.torchz	marray#99zfun, kwargsp   axis   c                     [        SUS9u  pEnU " US   4US   US.UD6nU " US   4US   SUS.UD6n[        UR                  UR                  U5      5        g )Nr:   r   r   r<   )weightsr;   omit)r?   
nan_policyr;   )r6   r   r0   r'   	funkwargsr;   r   r+   r1   narraysresrefs	            r5   
test_xmeanrH   3   si     'qR0C'
gaj
B'!*4
B6
BC
gaj
U'!*d
Uf
UCCHHbjjo.r7   )r   r<   NkeepdimsFTc                     [        SUS9u  p4n[        XS9n[        US   4SUS   0UD6n[        US   4US   SS.UD6n[        UR                  UR                  U5      5        g )	Nr:   r>   r;   rI   r   r?   r<   r@   r?   rA   )r6   dictr   r   r0   r'   )	r;   rI   r   r+   r1   rE   rD   rF   rG   s	            r5   test_xp_meanrN   A   so     'qR0C't/F
71:
<wqz
<V
<C
71:
Owqzf
O
OCCHHbjjo.r7   zarray-api-compat#242orderbiasn      limits)g?g?
lowerlimitg      ?
upperlimitc                     [        SUS9u  pEn[        US9U-  nU " US   40 UD6nU " US   4SS0UD6n[        UR                  UR	                  U5      5        g Nr<   r>   r;   r   rA   r@   )r6   rM   r   r0   r'   rB   s	            r5   test_severalrZ   N   sd    > 'qR0C't_v%F
gaj
#F
#C
gaj
6V
6v
6CCHHbjjo.r7   rD   c                    [        SUS9u  p4n[        U S9U-  n[        R                  " US   40 UD6n[        R                  " US   4SS0UD6n[	        UR
                  R                  UR                  UR
                  5      5        [	        UR                  S   R                  UR                  UR                  S   R                  5      5        [	        UR                  S   R                  UR                  UR                  S   R                  5      5        [	        UR                  R                  UR                  UR                  R                  5      5        [	        UR                  R                  UR                  UR                  R                  5      5        [	        UR                  R                  UR                  UR                  R                  5      5        g rX   )r6   rM   r   describer   nobsr0   r'   minmaxvarianceskewnesskurtosis)r;   rD   r   r+   r1   rE   rF   rG   s           r5   test_describerb   t   sH    'qR0C't_v%F
..
.v
.C
..
A
A&
ACCHHMM2::chh#78CJJqM&&

3::a=3E3E(FGCJJqM&&

3::a=3E3E(FGCLL%%rzz#,,2C2C'DECLL%%rzz#,,2C2C'DECLL%%rzz#,,2C2C'DEr7   rC   c                 L   U [         R                  :X  a
  [        SUS9O	[        SUS9u  p4nU " USU06nUR                  U " USUS.65      n[	        UR
                  UR                  )    XrR                  U5      )    5        [        UR                  US   R                  5        g )Nr:   r>   r<   r;   r@   rA   r;   r   )	r   zmapr6   r'   r   r0   r   isnanr   )rC   r;   r   r+   r1   rE   rF   rG   s           r5   test_zscorerg      s     69EJJ5FZb1",Q2"6 C'
w
"T
"C
**S'f4@
ACCHHchhY'hhsm^)<=CHHgajoo.r7   cupyzspecial functions won't workfc                    U R                   n[        SUS9u  pEnUS:X  a5  UR                  US   XS LS9US'   [        R                  " US   XS LS9US'   U " USU06nU " USUS.6n[        UR                  R                  UR                  UR                  5      5        [        UR                  R                  UR                  UR                  5      5        UR                  5       n	UR                  5       n
[        U	R                  R                  UR                  U
R                  5      5        [        U	R                  R                  UR                  U
R                  5      5        g )	Nr:   r>   ttest_1sampr<   rK   r;   r@   rd   )__name__r6   meanr   nanmeanr   	statisticr0   r'   pvalueconfidence_intervallowhigh)ri   r;   r   f_namer+   r1   rE   rF   rG   res_ciref_cis              r5   
test_ttestrw      s    ZZF&qR0C'XXgajt$>NXO
ZZ
D@PQ

W
 4
 C
Wd
3CCMM&&

3==(ABCJJOORZZ

%;<$$&F$$&FFJJOORZZ

%;<FKK$$bjj&=>r7   z7ignore::scipy.stats._axis_nan_policy.SmallSampleWarningc                 &   [        SUSS9u  p4nU " USU06nU " USUS.6n[        UR                  R                  UR	                  UR                  5      5        [        UR
                  R                  UR	                  UR
                  5      5        g )Nr<   r      r   r   r;   r@   rd   r6   r   ro   r0   r'   rp   )ri   r;   r   r+   r1   rE   rF   rG   s           r5   test_normality_testsr}      st     'qRx@C'
W
 4
 C
Wd
3CCMM&&

3==(ABCJJOORZZ

%;<r7   c                 0    U R                  SS 5      b  S$ S$ )Nf_expr:   r<   )get)rD   s    r5   pd_nsamplesr      s    

7D)51<1<r7   c                      [        U 5      $ N)tuple)argss    r5   <lambda>r      s    dr7   )paired	n_samplesrY   c                .    [         R                  " XXBUS9$ )N)r;   ddoflambda_r   power_divergence)f_obsr   r   r   r;   s        r5   power_divergence_refr      s    !!%TgVVr7   r   )pearsonzlog-likelihoodzfreeman-tukeyzmod-log-likelihoodneymanzcressie-read	chisquarer   c                 ~  ^  [        SUSS9u  pEn[        X!S9nT S:X  a  Sm S nOU 4S jnU" US	   40 UD6n	[        US	   4S
T S.UD6n
[        U	R                  R
                  UR                  U
S	   5      5        [        U	R                  R
                  UR                  U
S   5      5        [        R                  " US	   5      [        R                  " US   5      -  n[        R                  " US   U) -  USS9[        R                  " US	   U) -  USS9-  nUS	==   UR                  U5      -  ss'   US	==   U-  ss'   U" U0 UD6n	[        US
T S.UD6n
[        U	R                  R
                  UR                  U
S	   5      5        [        U	R                  R
                  UR                  U
S   5      5        g )Nr:   )      r{   )r;   r   r   r   c                  .    [         R                  " U 0 UD6$ r   )r   r   )r   rD   s     r5   ri   *test_power_divergence_chisquare.<locals>.f   s    ??D3F33r7   c                  4   > [         R                  " U ST0UD6$ )Nr   r   )r   rD   r   s     r5   ri   r      s    ))4KKFKKr7   r   r@   )rA   r   r<   TrK   )r6   rM   r   r   ro   r0   r'   rp   r   rf   nansum)r   r   r;   r   r+   r1   rE   rD   ri   rF   rG   common_mask	normalizes   `            r5   test_power_divergence_chisquarer      s    'qRv>C't'F+	4	L GAJ
!&
!C
wqz
Xfg
XQW
XCCMM&&

3q6(:;CJJOORZZA%78 ((71:&'!*)==K71:44$O99WQZ;,6TDQRIAJ"**Y''JAJ)J
W

C
FG
Vv
VCCMM&&

3q6(:;CJJOORZZA%78r7   method)fisherr   mudholkar_georgetippettstoufferc                    [        SUSS9u  p4n[        XS9n[        R                  " US   40 UD6n[        R                  " US   4SS0UD6n[	        UR
                  R                  UR                  UR
                  5      5        [	        UR                  R                  UR                  UR                  5      5        U S:w  a  g [        R                  " US   4S	US
   0UD6n[        R                  " US   4US
   SS.UD6n[	        UR
                  R                  UR                  UR
                  5      5        [	        UR                  R                  UR                  UR                  5      5        g )Nr:   ry   r{   )r   r;   r   rA   r@   r   r?   r<   rL   )	r6   rM   r   combine_pvaluesr   ro   r0   r'   rp   )	r   r;   r   r+   r1   rE   rD   rF   rG   s	            r5   test_combine_pvaluesr      s3    'qRx@C'+F




5f
5C




Hv
H
HCCMM&&

3==(ABCJJOORZZ

%;<




IGAJ
I&
IC



 =GAJ+1=5;=C CMM&&

3==(ABCJJOORZZ

%;<r7   c           
      .   Sn[        SXS9u  p#n[        R                  " [        R                  " U Vs/ s H  n[        R                  " U5      PM     sn5      SS9R                  [        5      nU Vs/ s H  oUU)    PM
     nnUS   S-  US   S-  sUS'   US'   US   S-  US   S-  sUS'   US'   [        R                  " U6 n[        R                  " U6 nU R                  U5      nU R                  U5      (       a  U R                  U) 5      (       d   e[        UR                  R                  U)    U R                  UR                  5      5        [        UR                  R                  U)    U R                  UR                  5      5        [        UR                  R                  U5        [        UR                  R                  U5        UR                  R                   U:X  d   eUR                  R                   U:X  d   eg s  snf s  snf )	Nry   r   r{   r   rY   r:   d   r   )r6   r   sumstackrf   r$   boolr   ttest_ind_from_statsr'   anyr   ro   r0   rp   r   r   )	r   r   r+   r1   rE   argr   rF   rG   s	            r5   test_ttest_ind_from_statsr     s    E&qR=C'66"((W=WcBHHSMW=>QGNNtTD%,-WcD5zWG-$QZ#-wqzC/?GAJ
$QZ#-wqzC/?GAJ


$
$g
.C

$
$g
.C::dD66$<<BFFD5MM))CMM&&u-rzz#--/HICJJOOTE*BJJszz,BCCMM&&-CJJOOT*==%'''::u$$$ >-s    H5H2z/Call to _getnamespace fails with AttributeErrorc                      [         R                  [        5      n SS/SS//nSS/SS//nU R                  XS9n[        R
                  " [        SS	9   [        US
U S9  S S S 5        g ! , (       d  f       g = f)Ng      ?g       @g      @g      @FTr   z7`axis` must be an integer or None for use with `MArray`)match)r   r<   )r;   r   )r   r   r   r'   pytestraisesNotImplementedErrorr   )r   r0   r   marrs       r5   1test_length_nonmasked_marray_iterable_axis_raisesr   .  s|     
		r	"B#Jc
#DENT5M*D::d:&D 
*G
I$V3
I 
I 
Is   A''
A5zignore::RuntimeWarningc                 l   [        SSU S9u  pn[        R                  " U6 nUS   [        R                  " [        R
                  " US   5      SS9)    US'   [        R                  " U6 n[        UR                  R                  U R                  UR                  5      5        [        UR                  R                  U R                  UR                  5      5        U R	                  UR                  R                  5      (       a   eU R	                  UR                  R                  5      (       a   eg )Nr<   )r   rR   r   r   r   rY   )r6   r   directional_statsr   r   rf   r   mean_directionr0   r'   mean_resultant_lengthr   )r   r+   r1   rE   rF   rG   s         r5   test_directional_statsr   >  s     'qR@C'

!
!7
+CRVVBHHWQZ$8qAABGAJ

!
!7
+CC&&++RZZ8J8J-KLC--22JJs889;vvc((--....vvc//4455555r7   	equal_varc                 ,   [        SUS9u  pEnU " USU0UD6nU " USUS.UD6n[        UR                  R                  UR	                  UR                  5      5        [        UR
                  R                  UR	                  UR
                  5      5        g )NrR   r>   r;   r@   rd   r|   rB   s	            r5   test_k_sample_testsr   O  s|     'qR0C'
w
,T
,V
,C
w6
?
?CCMM&&

3==(ABCJJOORZZ

%;<r7   c                    [        SSUS9u  p#nU " U6 nUu  pg[        R                  " U5      [        R                  " U5      -  nU " Xh)    Xx)    5      n	[        UR                  R
                  UR                  U	R                  5      5        [        UR                  R
                  UR                  U	R                  5      5        U [        R                  :X  ao  UR                  5       u  pU	R                  5       u  p[        U
R
                  UR                  U5      5        [        UR
                  UR                  U5      5        g g )Nr:   )   r   )r6   r   rf   r   ro   r0   r'   rp   r   pearsonrrq   )ri   r   r+   r1   rE   rF   xyr   rG   
res_ci_lowres_ci_high
ref_ci_lowref_ci_highs                 r5   test_pearsonrr   b  s     'q"=C'
W+CDA88A;!$D
AeHah
CCMM&&

3==(ABCJJOORZZ

%;<ENN"%"9"9";
"%"9"9";

J)?@(("**[*AB	 r7   qkc                     [        U (       a  SOSUS9u  p4n[        R                  " USU06n[        R                  " USUS.6n[        UR                  UR                  U5      5        g )Nr:   r<   r>   r;   r@   rd   )r6   r   entropyr   r0   r'   )r   r;   r   r+   r1   rE   rF   rG   s           r5   test_entropyr   y  sT     'BqA"=C'
--
,t
,C
--V$
?CCHHbjjo.r7   r   )Ynumpyr   r   scipyr   	packagingr   scipy._lib._array_apir   r   r   r   r	   r
   scipy.stats._stats_pyr   r   scipy.stats._axis_nan_policyr   importorskipr   markskipif
pytestmarkskip_xp_backendsskip_backendr6   parametrizegmeanhmeanpmeanrH   rN   momentskewra   semkstatkstatvarcircmeancircvarcircstdgstd	variationtmeantvartmintmaxtstdtsemrZ   r\   rb   zscoregzscorere   rg   rk   	ttest_rel	ttest_indrw   filterwarningsskewtestkurtosistest
normaltestjarque_berar}   r   r   r   r   r   r   r   r   r   parse__version__r   r   r   bartlettf_onewayr   r   pointbiserialrr   r   r    r7   r5   <module>r      s
       U U I 1 3 A 
		X	&
KK6  
 {{++"+&~ $2 l#TUk"FGgk*)=ekk2)N)=ekk2)N)=ekkCQR8)T)V W !Q(/ )W + H V/ l#TUk"FGgk*.eT]3/ 4 / + H V
/ l#TUk"FGg45%,,!5%**b)%**vuo6%.."-%..65/:%))R(%++Qx0%++Qx0%++Qx0%++Qx0%..3(3%..3(3%.."-%--,%--,%**b)%//2.r]%++*'=>%**x&<=%**|S&9:%**|S&9:%**x&<=%**x&<=/4 ./ /5 6 H V</ 5>>"l#TUk"FGg45!Q(B4(
F ) ) 6 H V #
F l#TUk"FGg45!5ell!C!5emm!D!5ejj!A!C D ./ /D 6 H V/ l#TUk"FGg45f;<3E4E4EF3EOOD3EOODF G .? /G = 6 H V?  l#TUk"FGg45f;<UV3ENNC3E4F4FG3E4D4DE3E4E4EFH I .= /	I W = 6 H V== 3DKXWA W YW 5??E$:$:;l#TUk"FGg45f;< %2 3 !Q(.9 / )3 = 6 H V <9B 5(()l#TUk"FGg45f;< $; <.= /< = 6 H V *=* 5--.l#TUk"FGg45f;<% = 6 H V /
%* GMM"..1GMM#4FFL  N4N4 5**+l#TUk"FGg4545	6 6 6 H V ,
	6 5>>"l#TUk"FGg45f;<
^^R
^^k4()
^^k5)*) 
 .= / = 6 H V #= l#TUk"FGg45f;<3ENNC3E4H4HIK LCL = 6 H VC" 5==!l#TUk"FGt}-./ / . H V "
/r7   