+
    ,i.                       ^ RI t ^ RIt^ RIt^ RIt^ RIt^ RIHt ^ RIt^ RI	H
t
 ^ RIHtHtHtHtHtHtHtHt ^ RIt^ RIHt ^ RIHtHtHt ^ RIHtHtHt ^RIH t  ^R	I!H"t"H#t# ^ R
I$H%t% ^ RI&H't' ^ RI(H)t)H*t*H+t+H,t, ^ RI-H.t. ^ RI/H0t0H1t1H2t2 ]Pf                  Ph                  t4]5! ]'4      t' ^ RI6t6R]6Pn                  R&   ^ RI8H9t: Rt;. EROt=. EROt>. EROt?. EROt@. EROtA. EROtB. EROtC. EROtD. EROtE. EROtFR tG ! R R4      tH ! R R4      tI ! R R4      tJ ! R R4      tK ! R R4      tL ! R R4      tM ! R  R!4      tN ! R" R#4      tO ! R$ R%4      tP ! R& R'4      tQR( tR ! R) R*4      tS ! R+ R,4      tT ! R- R.4      tU. EROtV ! R/ R04      tW ! R1 R24      tX ! R3 R44      tY ! R5 R64      tZ ! R7 R84      t[. R9NR:NR;NR<NR=NR>NR?NR@NRANRBNRCNRDNRENRFNRGNRHNRINRJNRKNRLNRMNRNNRONRPNRQNRRNRSNRTNRUNRVNRWNRXNRYNRZNR[NR\NR]NR^NR_NR`NRaNRbNRcNRdNReNRfNRgNRhNRiNRjNRkNRlNRmNRnNRoNRpNRqNRrNRsNRtNRuNRvNRwNRxNRyNRzNR{NR|NR}NR~NRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNER NERNERNERNERNERNERNERNERNER	NER
NERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNER NER!NER"NER#NER$NER%NER&NER'NER(NER)NER*NER+NER,NER-NER.NER/NER0NER1NER2NER3NER4NER5NER6NER7NER8NER9NER:NER;NER<NER=NER>NER?NER@NERANERBNERCNERDNERENERFNERGNERHNERINERJNERKNERLNERMNERNNERONERPNERQNERRNERSNERTNERUNERVNERWNERXNERYNERZNER[NER\NER]NER^NER_NER`NERaNERbNERcNERdNEReNERfNERgNERhNERiNERjNERkNERlNERmNERnNERoNERpNERqNERrNERsNERtNERuNERvNERwNERxNERyNERzNER{NER|NER}NER~NERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNERNt\ ! ER ER4      t] ! ER ER4      t^ ! ER ER4      t_ ! ER ER4      t` ! ER ER4      ta ! ER ER4      tb ! ER ER4      tc ! ER ER4      td ! ER ER4      te ! ER ER4      tf ! ER ER4      tg ! ER ER4      thR#   ]< d    Rt; ELi ; i(      N)partial)RandomState)assert_array_equalassert_almost_equalassert_array_lessassert_array_almost_equalassert_assert_allcloseassert_equalsuppress_warnings)raises)optimizestatsspecial)
_abw_state_get_As_weibull_Avals_weibull)check_named_results)_get_wilcoxon_distr_get_wilcoxon_distr2)_binary_search_for_binom_tst)distcont)SmallSampleWarningtoo_small_nd_omittoo_small_1d_omittoo_small_1d_not_omit)is_numpy)xp_assert_closexp_assert_equalxp_assert_lessAggbackendTFc                  j    \         P                  ! \        P                  P                  ! V / VB 4      # N)nplogr   gammarvs)argskwargss   *,B/usr/lib/python3/dist-packages/scipy/stats/tests/test_morestats.py_old_loggamma_rvsr,   E   s#    66%++//426233    c                   2   a  ] tR t^It o R tR tR tRtV tR# )TestBayes_mvsc                X   . R
Op\         P                  ! V4      w  r#p\        VP                  R4       \	        VP
                  RRR7       \        VP                  R4       \	        VP
                  RRR7       \        VP                  R^R7       \	        VP
                  RRR7       R	# )         "@ư>rtol      $@&.>_紫@decimal+=Nr1   	            r@      gURH#j@gV[%@gai	@gӇu8@gYs?g1O@)r   	bayes_mvsr   	statisticr
   minmaxselfdatameanvarstds   &    r+   
test_basicTestBayes_mvs.test_basicJ   s     '.3DNNC0%L!	# 	CMM40

$K"	$ 	CMM+=rJ

$L"	$r-   c                D    \        \        \        P                  . 4       R # r$   )assert_raises
ValueErrorr   rE   rI   s   &r+   test_empty_inputTestBayes_mvs.test_empty_input\   s    j%//26r-   c                    \         P                  ! ^4      pRp\        P                  ! V4      pV F  p\	        WB4       K  	  R# )   N)rF   rG   )r%   aranger   rE   r   )rI   x
attributesresis   &    r+   test_result_attributes$TestBayes_mvs.test_result_attributes_   s4    IIbM,
ooa A. r-    N)	__name__
__module____qualname____firstlineno__rN   rT   r]   __static_attributes____classdictcell____classdict__s   @r+   r/   r/   I   s     $$7/ /r-   r/   c                   8   a  ] tR t^ht o R tR tR tR tRtV t	R# )TestMvsdistc                   . R
Op\         P                  ! V4      w  r#p\        VP                  4       R4       \	        VP                  R4      RRR7       \        VP                  4       R4       \	        VP                  R4      RRR7       \        VP                  4       R^R7       \	        VP                  R4      RRR7       R	# )r1   r2   ?r;   r4   r6   r7   r8   r9   Nr<   rB   rC   rD   )r   mvsdistr   rK   r
   intervalrH   s   &    r+   rN   TestMvsdist.test_basici   s    &t,3DIIK-c* -AGL	N 	CHHJ-S) ,?EJ	L 	CHHJ(:BGS) ,@FK	Mr-   c                D    \        \        \        P                  . 4       R # r$   rQ   rR   r   rl   rS   s   &r+   rT   TestMvsdist.test_empty_inputx   s    j%--4r-   c                J    ^.p\        \        \        P                  V4       R#    Nrp   rI   rJ   s   & r+   test_bad_argTestMvsdist.test_bad_arg{   s    sj%--6r-   c                   \         P                  ! 4       ;_uu_ 4        \         P                  ! R \        4       \        P
                  ! . RO4       Uu. uF  qP                  4       NK  	   p\        P
                  ! . RO4       Uu. uF  qP                  4       NK  	   pRRR4       R# u upi u upi   + '       g   i     R# ; i)errorNrt         rt   r{   r|         )warningscatch_warningssimplefilterRuntimeWarningr   rl   rK   rI   rY   s   & r+   
test_warnsTestMvsdist.test_warns   s     $$&&!!'>:$}}Y787!VVX78$}}_=>=!VVX=> '&8> '&&s)   6B6B,,B6B1 B6,
B66C	r_   N)
r`   ra   rb   rc   rN   rT   rv   r   rd   re   rf   s   @r+   ri   ri   h   s     M57
? ?r-   ri   c                      a  ] tR t^t o R tR t]P                  P                  R. ^.^^.34      R 4       t	R t
R tR tRtV tR	# )
TestShapiroc                8   . ROp\         P                  ! V4      w  r#\         P                  ! V4      p\        VR^R7       \        VP                  R^R7       \        VR^R7       \        VP                  R^R7       . ROp\         P                  ! V4      w  r#\         P                  ! V4      p\        VR^R7       \        VP                  R^R7       \        VR^R7       \        VP                  R^R7       \         P
                  P                  ^^^dRR7      p\         P                  ! V4      w  r#\         P                  ! V4      p\        VR^R7       \        VP                  R^R7       \        VR	^R7       \        VP                  R	^R7       . ROpR
pRp	\         P                  ! V4      w  r#\         P                  ! V4      p\        W(^R7       \        VP                  V^R7       \        W9^R7       \        VP                  V	^R7       R# ))\(?   ?r9      Ɍ?IddY?Dio?iNa locscalesizerandom_stateg   E?g   Oٴ?ghyܝ?g0M?N)r   {Gz@q=
ףp@HzG$@@Q	@q=
ףp?Q@zG?      @Q?RQ?Q	@(\?zG?zG@(\#@333333@Gz?p=
ף@)(\?=
ףp=?\(\@ffffff@\(\?(\?Gz@(\ףp=
@皙?)\(?RQr   p=
ף?@
ףp=
?{Gz?\(\@{Gz@Q@)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 @)r   shapiror   rF   pvaluenormr(   )
rI   x1wpwshapiro_testx2x3x4
W_expected
p_expecteds
   &         r+   rN   TestShapiro.test_basic   s   8 b!}}R(A2A>L224GQRSB 4a@L//1EqQ& b!}}R(Ay!4L22IqIB3L//!D ZZ^^8^Lb!}}R(A11=L224FPQRB 3Q?L//1DaP? 

b!}}R(A15L22JJBA6L//QGr-   c                   . RO. RO.p\         P                  ! V4      w  r#\         P                  ! V4      p\        VR^R7       \        VP                  R^R7       \        VR^R7       \        VP                  R^R7       . R	O. R
O.p\         P                  ! V4      w  r#\         P                  ! V4      p\        VR^R7       \        VP                  R^R7       \        VR^R7       \        VP                  R^R7       R# )r   r   r9   r   r   r   N)
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   r   r   r   r   r   r   r   r   )r   r   r   rF   r   )rI   r   r   r   r   r   s   &     r+   test_2dTestShapiro.test_2d   s     "89 b!}}R(A2A>L224GQRSB 4a@L//1EqQ&' b!}}R(Ay!4L22IqIB3L//!Dr-   rY   c                L   \         P                  ! \        \        R 7      ;_uu_ 4        \        P
                  ! V4      p\        VP                  \        P                  4       \        VP                  \        P                  4       RRR4       R#   + '       g   i     R# ; imatchN)pytestwarnsr   r   r   r   r   rF   r%   nanr   rI   rY   r[   s   && r+   test_not_enough_values"TestShapiro.test_not_enough_values   sW    \\,4IJJ--"C/RVV, KJJJs   ABB#	c                   \         P                  ! R 4      p\         P                  V^	&   \        P                  ! V4      w  r#\        P                  ! V4      p\        V\         P                  4       \        VP                  \         P                  4       \        V\         P                  4       \        VP                  \         P                  4       R# r6   N)	r%   rX   r   r   r   r   rF   r   r   )rI   rY   r   r   r   s   &    r+   test_nan_inputTestShapiro.test_nan_input   s    IIcNvv!a }}Q'Q\++RVV4 	B'L//8r-   c                    \         P                  ! . RO4      w  r\         P                  ! V4      pRp\        W4RR7       R# ) h㈵>r4   N)r   i = i@ )gPq~?gF?)r   boxcoxr   r
   )rI   	trans_valmaxlogr[   refs   &    r+   test_gh14462TestShapiro.test_gh14462   s5     "LL)AB	mmI& 2t,r-   c                    \         P                  ! . RO4      pVP                  ^ 8  g   Q h. ROp\         P                  ! V4      p\        VP                  R4       \        VP                  RRR7       R# )9B.?gr??gel?r3   r4   N)r           r   )g1gIۿg0!?)r   r   r   r
   rF   )rI   r[   rY   s   &  r+   test_length_3_gh18322!TestShapiro.test_length_3_gh18322   sX     mm:;zzQ KmmA'78

O$?r-   r_   N)r`   ra   rb   rc   rN   r   r   markparametrizer   r   r   r   rd   re   rf   s   @r+   r   r      sZ     )HVE* [[S2sQF"34- 5-9-"@ @r-   r   c                      a  ] tR tRt o R tR tR tR tR tR t	R t
R	 tR
 t]P                  P                  R 4       t]P                  P#                  R. RO4      R 4       tR tRtV tR# )TestAndersoni  c                   \        R 4      pVP                  ^2R7      pVP                  ^2R7      p\        P                  ! V4      w  rEp\        VRR V4       \        P                  ! V4      w  rEp\        WERR 4       \        P                  ! ^
4      p^ V^ &   \        P                  ! V4      w  rEp\        VR4       R# )Ir   Ng	@)	r   standard_exponentialstandard_normalr   andersonr   r%   onesr
   )rI   rsr   r   Acritsigvs   &       r+   test_normalTestAnderson.test_normal  s    $$$"$-R(~~b)$s)Q'~~b)!"#Y'GGBK!~~a( 	8$r-   c                   \        R 4      pVP                  ^2R7      pVP                  ^2R7      p\        P                  ! VR4      w  rEp\        WERR 4       \        P                  ! RR7      ;_uu_ 4        \        P                  ! VR4      w  rEpRRR4       \        WER,          8  4       R#   + '       g   i     L&; i)r   r   exponNignoreallr   r   )	r   r   r   r   r   r   r%   errstater	   )rI   r   r   r   r   r   r   s   &      r+   
test_exponTestAnderson.test_expon  s    $$$"$-R(~~b'2!"#Y'[[X&& >>"g6LAS 'H '&s   7B00C 	c           	         \         P                  ! ^d4      pRV^ &   \        P                  ! VR4      w  r#p\	        V4      p\        P
                  P                  V4      w  rg\        P
                  P                  WV4      p\        P
                  P                  WV4      p	\         P                  ! ^V^,           4      p
V) \         P                  ! ^V
,          ^,
          WRRR1,          ,           ,          4      ,
          p\        W+4       R# )d   r   gumbelNr   )r%   r   r   r   lengumbel_lfitlogcdflogsfrX   rK   r
   )rI   r   a2r   r   nxbarsr  r	  r\   expected_a2s   &           r+   test_gumbelTestAnderson.test_gumbel#  s     GGCL!q(3#F..$$Q'&&q2$$Qa0IIa1b277AaC!Gtt0D#EFF(r-   c                J    \        \        \        P                  ^.RR7       R# rt   plate_of_shrimpdistN)rQ   rR   r   r   rS   s   &r+   rv   TestAnderson.test_bad_arg3  s    j%..1#<MNr-   c                    \        R 4      pVP                  ^2R7      p\        P                  ! V4      pRp\	        W44       R# )r   r   NrF   critical_valuessignificance_level)r   r   r   r   r   )rI   r   rY   r[   rZ   s   &    r+   r]   #TestAnderson.test_result_attributes6  s:    $###,nnQK
C,r-   c                    \        R 4      pVP                  ^dR7      p\        P                  ! VR4      w  r4p\        P                  ! VR4      w  rgp\	        Wc4       R# )r   r   r  r  N)r   r  r   r   r
   )	rI   r   rY   A1crit1sig1A2crit2sig2s	   &        r+   test_gumbel_lTestAnderson.test_gumbel_l=  sN     $II3I..H54..J74r-   c                "   \        R 4      pVP                  ^dR7      p\        P                  ! ^d4      pRV^ &   \        P
                  ! VR4      w  rEp\        P
                  ! VR4      w  rxp	\        WERR 4       \        WxR,          8  4       R# )r   r   Zd;?gumbel_rNr   r   )r   r  r%   r   r   r   r   r	   )
rI   r   r   r   r  r  r  r   r!  r"  s
   &         r+   test_gumbel_rTestAnderson.test_gumbel_rG  sz     $YYCY WWS\ 1..Z84..Z84"BCj)2Yr-   c                \   \         P                  ! . R	O4      p\        P                  ! VR4      pVP                  P
                  w  r4p\        W4V3R
RR7       \        VP                  RRR7       VP                  VP                  ^ ,          8  g   Q h^V,          p\        VRRR7       \        R,          p\        R,          pWvR,
          R,          W,
          ,          ,           p	\         P                  ! VP                  V	8  4      '       g   Q h\        VP                  V	RR7       R# )   weibull_minMb`?r4   gp=
ף?MbP?皙?atolN)
r+           r4        r5  u   r6  )g
ףp=
@gzGX@gQS@g?r   g?)r%   arrayr   r   
fit_resultparamsr
   rF   r  r   r   )
rI   rY   r[   mr   r   cAs40As45As_refs
   &         r+   test_weibull_min_case_A$TestAnderson.test_weibull_min_case_AV  s    HHGHnnQ.--)=DIu48}}s2215555E6-b!b!S:.$+>>vvc))F23333++V$?r-   c                    \         P                  ! . RO4      pRp\        P                  ! \        VR7      ;_uu_ 4        \
        P                  ! VR4       RRR4       R#   + '       g   i     R# ; i)J   z/Maximum likelihood estimation has converged to r   r,  N)rD  9   0        r>   F      rG    ;            iF  )r%   r9  r   r   rR   r   r   )rI   rY   messages   &  r+   test_weibull_min_case_B$TestAnderson.test_weibull_min_case_Bi  sH    HH 5 6C]]:W55NN1m, 6555s   A!!A2	c           	     x   \         P                  ! . RO4      ) pRpRp\        P                  ! \        VR7      p\        P
                  ! \        VR7      pT;_uu_ 4        V;_uu_ 4        \        P                  ! VR4       RRR4       RRR4       R#   + '       g   i     L; i  + '       g   i     R# ; i)r+  z:Critical values of the test statistic are given for the...z;An error occurred while fitting the Weibull distribution...r   r,  N)	r+  K   rE     k   r>   =   +   rG  )	r%   r9  r   r   UserWarningr   rR   r   r   )rI   rY   wmessageemessagewcontextecontexts   &     r+   test_weibull_warning_error'TestAnderson.test_weibull_warning_errorq  sl     XX=>>OP<<8<==8<XxxNN1m,  XXxxXXXs$    B(+B	B(B% B((B9	distnamec                $   \         P                  P                  R 4      pVR9   d   RMTp\        \        V4      p\
        V,          pVP                  ! VRRRV/ p\        P                  ! Wa4      pVP                  P                  '       g   Q hR# )l   'Dah$r  r     r   N>   r  extreme1)
r%   randomdefault_rnggetattrr   r   r(   r   r:  success)rI   r`  rngreal_distnamer  r;  rY   r[   s   &&      r+   test_anderson_fit_params%TestAnderson.test_anderson_fit_params}  s~    
 ii##$67'/3I'I& 	um,-(HHf:4:c:nnQ)~~%%%%%r-   c                    ^p\        \        ^V,          4      \        R,          4       \        P                  p\        \        ^V,          4      \        ^ ,          4       R# rt   Nr   )r   r   r   r%   inf)rI   r<  s   & r+   test_anderson_weibull_As%TestAnderson.test_anderson_weibull_As  s@    _QqS)>"+=>FF_QqS)>!+<=r-   r_   N)r   r   r  rc  r  r'  logisticr,  )r`   ra   rb   rc   r   r   r  rv   r]   r#  r(  rA  rQ  r   r   thread_unsafer^  r   rj  ro  rd   re   rf   s   @r+   r   r     s     %.) O-  @&- [[	- 	- [[ZOP	&P	&> >r-   r   c                      a  ] tR tRt o R tR t]P                  P                  R 4       t	R t
]P                  P                  R 4       tR tR tR	 tR
 tRtV tR# )TestAndersonKSampi  c                :   . R	Op\         P                  ! . R
O4      p\         P                  ! . RO4      p\         P                  ! . RO4      p\        P                  ! WW43RR7      w  rVp\	        VR^4       \        . ROVR,          ^4       \        VRRR7       R# )皙YC@Fmidrankg@r   r   NgJ4a?Mb0?r0  Nrv  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        F@g     F@r~  ffffffA@r  g33333A@gB@gfffffB@gD@gffffffE@gCl?gvq-?g,Ԛ?gl@gPk	@r%   r9  r   anderson_ksampr   r   r
   rI   t1t2t3t4Tktmps   &       r+   test_example1a TestAndersonKSamp.test_example1a  s}    
 >XXFGXXFGXXFG(("")95I	Bq)!"J"$S'1	.60r-   c                b   \         P                  ! . R	O4      p\         P                  ! . R
O4      p\         P                  ! . RO4      p\         P                  ! . RO4      p\        P                  ! WW43RR7      w  rVp\	        VR^4       \        . ROVR,          ^4       \        VRRR7       R# )rv  Trw  gQ@ry  r-  rz  r0  Nr{  r|  r}  r  r  r  r  s   &       r+   test_example1b TestAndersonKSamp.test_example1b  s    
 XXFGXXFGXXFGXXFG(("")94H	Bq)!"J"$S'1	.60r-   c                J   . ROp. ROp. ROp. ROp. ROp. ROp. ROp. ROp. ROp	. ROp
. ROp. ROp. ROp. ROpWW4WVWxWWW3p\         P                  ! VRR7      w  ppp\        VR^4       \        . ROVR,          ^4       \	        VRRR7       \
        P                  P                  R4      p\         P                  ! R	VR
7      p\         P                  ! VRVR7      p\        VP                  V4       \        VP                  V4       \	        VP                  VRR7       R# )   Frw  gM
@ry  	^)p?rz  r0  l   <Ic i'  n_resamplesrh  )rx  methodga2U0*C?Nr  rW   )   rG  !      i     :   %   r  A   r=      i     $      v   
"         r  C   rE  >   r?      r  Z   
   <      rW  1   r     8      O   T   ,   rL  rG  r        i6  L   rO  r     r        rI  e   r  rD  rE  rF  rG  rH  r>   rI  rJ  rG  rK  rL  rM  7   i@  r  h         /            r  r  i  W   r?   x   r  r  r  r+  G   r  rJ  *   r  r   r>   r     r|   r  r  r  r  r     r  rt   r  4   _   a   3   r  r~      r     D   M   P   rt   r  j      R   6   r     .   o   '   ?   r     r     r  2   r  f   H   r  r  r|   rW         r  X   r  r   r   r  r        r     r  rA   r  	g  r=   r>   i  i[  r|   r  r{   i  r     r   i  #   r>   i  r  r  r?   b   r   U   [   rX     r|   r  r     r  rE  r      r  rL        rM  r  r  B   rW  r  g^I+?g~:?g+e?gH.!@gcZB>h@)r   r  r   r   r
   r%   rd  re  PermutationMethodr   rF   r  r   )rI   r  r  r  r  t5t6t7t8t9t10t11t12t13t14samplesr  r  r  rh  r  r[   s   &                     r+   test_example2a TestAndersonKSamp.test_example2a  s   
 (G45%?@IB14'? 2222CcO((%@	BBq)!"J"$S'1	.60ii##$78((TsC""7E&I3=="-3..3

AD1r-   c                   . R	Op. R
Op. ROp. ROp. ROp. ROp. ROp. ROp. ROp	. ROp
. ROp. ROp. ROp. ROp\         P                  ! WW4WVWxWWW3RR7      w  ppp\        VR^4       \        . ROVR,          ^4       \	        VRRR7       R# )r  Trw  g'1Z
@ry  r  rz  r0  Nr  r  r  r  r  r  r  r  r  r  r  r  r  r  r  )r   r  r   r   r
   )rI   r  r  r  r  r  r  r  r  r  r	  r
  r  r  r  r  r  r  s   &                 r+   test_example2b TestAndersonKSamp.test_example2b  s     (G45%?@IB14'? ((""""*,3S*G157	B 	Bq)!"J"$S'1	.60r-   c                2   \         P                  ! ^^d^d4      p\        4       ;_uu_ 4       pVP                  \        RR7       \
        P                  ! WR,           .RR7      w  r4pRRR4       \        XR^4       \        XR4       \        4       ;_uu_ 4       pVP                  \        RR7       \
        P                  ! WR,           .4      w  r4pRRR4       \        VR	^4       \        VR4       \        4       ;_uu_ 4       pVP                  \        R
R7       \
        P                  ! WR,           .RR7      w  r4pRRR4       \        VR^4       \        VR4       \        4       ;_uu_ 4       pVP                  \        R
R7       \
        P                  ! WR,           .4      w  r4pRRR4       \        VR^4       \        VR4       \
        P                  ! WR,           .RR7      w  r4p\        VR^4       \        VRR^ R7       \
        P                  ! W^,           .4      w  r4p\        VR^4       \        VRR^ R7       \
        P                  ! WR,           .RR7      w  r4p\        VR^4       \        VRR^ R7       \
        P                  ! WR,           .RR7      w  r4p\        VR^4       \        VRR^ R7       R#   + '       g   i     ELl; i  + '       g   i     EL; i  + '       g   i     EL; i  + '       g   i     ELb; i)rt   zp-value floored)rP  g     @D@Frw  Ng=
ףpD@r.  gGzD@zp-value capped      ?      ?      @g!u?gףp=
׳?{Gzt?r1  r5   gBiq?g鷯?      '@g:M@g46<Rv?gMb@?g      +@g6>W[1@g1ZGUS?-C6?g?ܵg?ܵ)
r%   linspacer   filterrY  r   r  r   r   r
   )rI   r   supr  _r  s   &     r+   test_R_kSamples!TestAndersonKSamp.test_R_kSamples  sx   D [[C%   CJJ{,=J>**BT	?EJGA! ! 	Avq)Q  CJJ{,=J>**BT	?;GA! ! 	Avq)Q   CJJ{,<J=**BR=%HGA! ! 	Aw*Q  CJJ{,<J=**BR=9GA! ! 	Aw*Q &&H~uEaAvq)6A6 &&F|4aAvq)6A6 &&IFaAvq)7a8&&IFaAvq)7a8U !   !   !   !  s/   ;K	9K;K139L	K	K.	1L	L	c                l    \        \        \        P                  \        P
                  ! ^4      4       R# r   NrQ   rR   r   r  r%   r   rS   s   &r+   test_not_enough_samples)TestAndersonKSamp.test_not_enough_samplesL  s    j%"6"6
Cr-   c                    \        \        \        P                  \        P
                  ! ^4      \        P
                  ! ^4      34       R# r#  r$  rS   s   &r+   test_no_distinct_observations/TestAndersonKSamp.test_no_distinct_observationsO  s,    j%"6"6wwqz2771:.	0r-   c                p    \        \        \        P                  \        P
                  ! ^4      . 34       R# r#  r$  rS   s   &r+   test_empty_sample#TestAndersonKSamp.test_empty_sampleS  s!    j%"6"6R8HIr-   c                    . ROp\         P                  ! . RO4      p\        P                  ! W3RR7      pRp\	        W44       \        VP                  VP                  4       R# )rv  Frw  Nr{  r|  r  )r%   r9  r   r  r   r   r  r   )rI   r  r  r[   rZ   s   &    r+   r]   (TestAndersonKSamp.test_result_attributesV  sK    =XXFG""B8U;K
C,S++SZZ8r-   r_   N)r`   ra   rb   rc   r  r  r   r   xslowr  r  rr  r   r%  r(  r+  r]   rd   re   rf   s   @r+   rt  rt    sp     1"1  [[$2 $2L1@ [[N9 N9`D0J	9 	9r-   rt  c                   N  a  ] tR tRt o R tR tR t]P                  P                  R. ^.3^.. 3.4      R 4       t
R tR tR	 t]P                  P                  R
. RO. RORR3. RO. RORR3. RO. RORR3. RO. RORR3. RO. RORR3. RO. RORR3.4      R 4       tR tRtV tR# )
TestAnsariib  c                   . ROp. ROp\        4       ;_uu_ 4       pVP                  \        R4       \        P                  ! W4      w  rERRR4       \        XR^4       \        XR^4       R#   + '       g   i     L,; i)rt   %Ties preclude use of exact statistic.Ng     7@g8#oG?rt   r{   r|   r|   r~   r|   r{   r1   rt   r1   rt   r~   rt   )r   r  rY  r   ansarir   )rI   rY   yr  Wpvals   &     r+   
test_smallTestAnsari.test_smalld  s_    $  CJJ{$KLll1(GA ! 	AtR(D"5r:	 ! s   /A..A>	c                J   \         P                  ! R4      p\         P                  ! R4      p\        4       ;_uu_ 4       pVP                  \        R4       \
        P                  ! W4      w  rERRR4       \        XR^4       \        XR^4       R#   + '       g   i     L,; i)r  r3  Ng     0g@g|:?r  rV  r  c   r  r  m   l   r  r>  r  `   r  r  rV  q   t   rB  n   r  rV  r@  r  r  i   g   rD  rF  r  r  rA  r@  rG  r  r   rH  rB  r@  r  r>  )r%   r9  r   r  rY  r   r6  r   )rI   ramsayparekhr  r8  r9  s   &     r+   test_approxTestAnsari.test_approxm  s     J K $ %   CJJ{$KLll62GA ! 	Aub)D"5r: ! s   /BB"	c                t    \         P                  ! . RO. RO4      w  r\        VR^4       \        VR^4       R# )rt   r6   g?Nrt   r{   r|   r~   )rW   r   r  r@   r  r>   )r   r6  r   )rI   r8  r9  s   &  r+   
test_exactTestAnsari.test_exact{  s.    ,,|-CDAtR(D"6:r-   r)   c                H   \         P                  ! \        \        R 7      ;_uu_ 4        \        P
                  ! V!  p\        VP                  \        P                  4       \        VP                  \        P                  4       RRR4       R#   + '       g   i     R# ; ir   )r   r   r   r   r   r6  r   rF   r%   r   r   )rI   r)   r[   s   && r+   rv   TestAnsari.test_bad_arg  sW    \\,4IJJ,,%C/RVV, KJJJs   ABB!	c                    . ROp. ROp\        4       ;_uu_ 4       pVP                  \        R4       \        P                  ! W4      pRRR4       Rp\        XV4       R#   + '       g   i     L ; i)rt   r3  Nr4  r5  rF   r   )r   r  rY  r   r6  r   )rI   rY   r7  r  r[   rZ   s   &     r+   r]   !TestAnsari.test_result_attributes  sS    $  CJJ{$KL,,q$C ! -
C,	 ! s   -A  A0	c                    . ROp. ROpRp\        \        VR7      ;_uu_ 4        \        P                  ! WRR7       RRR4       R#   + '       g   i     R# ; i)rt   z!'alternative' must be 'two-sided'r   fooalternativeNrN  r   r1   r?   r@   )rQ   rR   r   r6  )rI   r   r   r   s   &   r+   test_bad_alternativeTestAnsari.test_bad_alternative  s;    3:U33LLU3 4333s   AA	c                   . R	Op. R
Op\         P                  ! W4      w  r4\         P                  ! WRR7      P                  p\         P                  ! WRR7      P                  pVR8  g   Q hVR8  g   Q h\        P                  P                  V\        V4      \        V4      4      p\        We,           ^V,           RR7       \        Wd^,          RR7       \        V^V,           V^,          ,
          RR7       \         P                  ! W!RR7      P                  p\         P                  ! W!RR7      P                  p	VR8  g   Q hV	R8  g   Q hR# )r   lessrX  greaterffffff?皙?-q=r0  N)rt   r   r  rW   r  r  )r  g      !@      #@g      %@r  g      )@)r   r6  r   r   apmfr  r
   )
rI   r   r   rF   r9  pval_lpval_gprobpval_l_reversepval_g_reverses
   &         r+   test_alternative_exact!TestAnsari.test_alternative_exact  s   '.  ,,r.	b&9@@b)<CC}}}} ||	3r7CG<T> 	QU3$tAvE: b&AHHb)DKK$$$$$$r-   zx, y, alternative, expectedr^  g@A?r_  gG$I?g%I$I?c                d    \         P                  ! WVR 7      P                  p\        WTRR7       R# )rX  rb  r0  N)r   r6  r   r
   )rI   rY   r7  rY  expectedr9  s   &&&&& r+   test_alternative_exact_with_R(TestAnsari.test_alternative_exact_with_R  s$    8 ||Ak:AAU3r-   c                   \         P                  P                  ^ ^^d^{R7      p\         P                  P                  ^ ^^d^{R7      p\         P                  ! WRR7      P                  p\         P                  ! WRR7      P                  p\        VRRR7       \        VRRR7       \         P                  P                  ^ ^^<^{R7      p\         P                  P                  ^ R	^<^{R7      p\         P                  ! W4      P                  p\         P                  ! WRR7      P                  p\         P                  ! WRR7      P                  p\        WE^,          RR7       \        V^V^,          ,
          RR7       R
# )r   r   r   r^  rX  r_        ?rb  r0  r         ?N)r   r   r(   r6  r   r
   )rI   r   r   rg  rh  r9  s   &     r+   test_alternative_approx"TestAnsari.test_alternative_approx  s   ZZ^^Aqs^=ZZ^^Aqs^= b&9@@b)<CC%0%0 ZZ^^Aqr^<ZZ^^As#^>||B#**b&9@@b)<CCQU3$q&u5r-   r_   NrN  rZ  rz   )r~   r   r1   r?   r@   r}   )r1   r?   r@   )r`   ra   rb   rc   r:  rK  rO  r   r   r   rv   r]   r[  rl  rp  rv  rd   re   rf   s   @r+   r1  r1  b  s     ;;;
 [[VrA3i1#r%;<- =--4%6 [[% fo	>
i	A
_fo	>
_i	A
9fo	>
9i	A	C44&6 6r-   r1  c                   z   a  ] tR tRt o R tR tR t]P                  P                  R4      R 4       t
R tRtV tR	# )
TestBartletti  c           
     J   \         \        \        \        \        \
        \        \        \        \        .
pV Uu. uF  q1P                  V4      NK  	  pp\        P                  ! V!  w  rE\        WAP                  R 4      4       \        WQP                  R4      4       R# u upi )gzU@ /4@gG?N)g1g2g3g4g5g6g7g8g9g10asarrayr   bartlettr   )rI   xpr)   argTr9  s   &&    r+   	test_dataTestBartlett.test_data  sn    BBBBC8+/04C

340..$'::&789jj9: 1s   B c                    R p\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! VP                  R.4      4       RRR4       R#   + '       g   i     R# ; i)z-Must enter at least two input sample vectors.r   rt  N)r   r   rR   r   r  r  )rI   r  rP  s   && r+   test_too_few_argsTestBartlett.test_too_few_args  s>    A]]:W55NN2::rd+, 6555s   'AA)	c           
         \         \        \        \        \        \
        \        \        \        \        .
pV Uu. uF  q1P                  V4      NK  	  pp\        P                  ! V!  pRp\        WEVR7       R# u upi )rF   )r  NrT  )r{  r|  r}  r~  r  r  r  r  r  r  r  r   r  r   )rI   r  r)   r  r[   rZ   s   &&    r+   r]   #TestBartlett.test_result_attributes  sX    BBBBC8+/04C

340nnd#,
C3 1s   A9z*ignore:invalid value encountered in dividec                J   \         \        \        \        \        \
        \        \        \        \        . 3pV Uu. uF  q1P                  V4      NK  	  pp\        V4      '       dG   \        P                  ! \        \        R 7      ;_uu_ 4        \         P"                  ! V!  pRRR4       M\$        P&                  P)                  4       ;_uu_ 4       pVP+                  \,        R4       VP+                  \.        R4       VP+                  \,        R4       \         P"                  ! V!  pRRR4       VP                  VP0                  4      p\3        XP4                  V4       \3        VP6                  V4       R# u upi   + '       g   i     L^; i  + '       g   i     Lq; i)r   Ninvalid value encounteredz$var\(\): degrees of freedom is <= 0.z!Degrees of freedom <= 0 for slice)r{  r|  r}  r~  r  r  r  r  r  r  r  r   r   r   r   r   r   r  r%   testingr   r  r   rY  r   r   rF   r   )rI   r  r)   r  r[   r  NaNs   &&     r+   test_empty_argTestBartlett.test_empty_arg  s   BBBBC<+/04C

340B<<08MNNnnd+ ON --//3

>+FG

;(OP

>+NOnnd+ 0 jj s+

C( 1NN 0/s   E:E?AF?F	F"	c                :   VP                  . ROVP                  R7      pVP                  . ROVP                  R7      pVP                  . ROVP                  R7      p\        P                  ! W#V4      pVP	                  VP
                  ^ 8  4      '       g   Q hR# )333333$@dtypeN)r  gffffff$@g$@g$@)L$@     $@33333$@gfffff$@)g$@r  r  r  )r  float32r   r  r   rF   )rI   r  re  br=  r[   s   &&    r+   test_negative_pvalue_gh21152)TestBartlett.test_negative_pvalue_gh21152  su    JJ/rzzJBJJ32::JFJJ32::JFnnQ1%vvcmmq())))r-   r_   N)r`   ra   rb   rc   r  r  r]   r   r   filterwarningsr  r  rd   re   rf   s   @r+   ry  ry    sC     ;-
4 [[ LM) N)"* *r-   ry  c                   V   a  ] tR tRt o R tR tR tR tR tR t	R t
R	 tR
 tRtV tR# )
TestLevenei  c           
         \         \        \        \        \        \
        \        \        \        \        .
p\        P                  ! V!  w  r#\        VR ^4       \        VR^4       R# )gYpK?gI{]?N)r{  r|  r}  r~  r  r  r  r  r  r  r   levener   )rI   r)   r8  r9  s   &   r+   r  TestLevene.test_data  sE    BBBBC8,,%A115D/15r-   c                    \         P                  ! \        \        \        R R7      w  r\         P                  ! \        \        \        RRR7      w  r4\        W4       \        W$4       R# )rK   centertrimmedr   r  proportiontocutN)r   r  r{  r|  r}  r   )rI   W1pval1W2pval2s   &    r+   test_trimmed1TestLevene.test_trimmed1   sE     LLRF;	LLRI146	B#E)r-   c                   . ROp. R	Op\         P                  P                  R4       \         P                  P                  V4      p\        P
                  ! WRRR7      w  rE\        P
                  ! W2RRR7      w  rg\        P
                  ! V^R
 V^R
 RR7      w  r\        WH4       \        Wh4       \        Wy4       R# )333333?  r        ?r  rK   r  Nr         @      @      @      @      @      @g      Y@r   r        @r  g      @r  g      @      i@r   )r%   rd  seedpermutationr   r  r   )
rI   rY   r7  r   W0pval0r  r  r  r  s
   &         r+   test_trimmed2TestLevene.test_trimmed2)  s    66
		tYY""1% LLi168	LLy168	 LL1R!Ab'&A	B#B#E)r-   c                X   \         P                  ! R^^4      p\         P                  P                  R4       \         P                  P	                  V4      pV^,          p\
        P                  ! WRR7      w  rE\
        P                  ! W#RR7      w  rg\        WF4       \        WW4       R# )rt   r  rK   r  medianNr   )r%   r  rd  r  r  r   r  r   )rI   rY   r   r7  r  r  r  r  s   &       r+   test_equal_mean_median!TestLevene.test_equal_mean_median;  su    KKAr"
		tYY""1%qDLLf5	LLx8	B#E)r-   c                x    \         P                  ! R^^4      p\        \        \        P
                  WRR7       R# rt   皙?)portiontocutNr   )r%   r  rQ   	TypeErrorr   r  r   s   & r+   test_bad_keywordTestLevene.test_bad_keywordE  s&    KKAr"iq#Fr-   c                x    \         P                  ! R^^4      p\        \        \        P
                  WRR7       R# rt   trimr  Nr   )r%   r  rQ   rR   r   r  r   s   & r+   test_bad_center_value TestLevene.test_bad_center_valueI  s&    KKAr"j%,,VDr-   c                F    \        \        \        P                  ^.4       R# rs   )rQ   rR   r   r  rS   s   &r+   r  TestLevene.test_too_few_argsM  s    j%,,4r-   c           
         \         \        \        \        \        \
        \        \        \        \        .
p\        P                  ! V!  pRp\        W#4       R# )rF   NrT  )r{  r|  r}  r~  r  r  r  r  r  r  r   r  r   )rI   r)   r[   rZ   s   &   r+   r]   !TestLevene.test_result_attributesP  s7    BBBBC8llD!,
C,r-   c                    \         P                  ! ^^.^^..4      p\        \        \        P
                  \        V4       R# rs   )r%   r9  rQ   rR   r   r  r{  r   s   & r+   test_1d_inputTestLevene.test_1d_inputW  s/    HHq!fq!f%&j%,,A6r-   r_   N)r`   ra   rb   rc   r  r  r  r  r  r  r  r]   r  rd   re   rf   s   @r+   r  r    s9     6**$*GE5-7 7r-   r  c                     a  ] tR tRt o RtR tR tR tR t]	P                  P                  R. RO4      R 4       t]	P                  P                  R. R O4      R	 4       t]	P                  P                  R
. R!O4      R 4       t]	P                  P                  R. R"O4      R 4       t]	P                  P                  R. R#O. R$O. R%O. R&O. R'O. R(O. R)O. R*O. R+O. R,O. R-O. R.O. R/O. R0O. R1O. R2O. R3O. R4O. R5O. R6O. R7O. R8O. R9O. R:O.4      R 4       tR t]	P                  P                  RR;R<.4      R 4       tR tR tR tR tR tR tR t]	P                  P5                  ]P8                  R=8*  RR7      R 4       tRtV tR# )>TestBinomTesti\  zTests for stats.binomtest.c                   R p\         P                  ! RRR4      p\        VP                  RVR7       \         P                  ! RRR4      p\        VP                  RVR7       \         P                  ! RRR4      p\        VP                  R	VR7       \         P                  ! R
RR4      p\        VP                  RRR7       R# )绽|=iΙ i@o@Q?rt  r4   iΙ gX7?i	ϙ g=?iϙ gybw?r7   Nr   	binomtestr
   r   rI   r5   r[   s   &  r+   test_two_sided_pvalues1%TestBinomTest.test_two_sided_pvalues1c  s     ooh$7

Cd3ooh$7

O$?ooh$7

O$?ooh$7

O$?r-   c                   R p\         P                  ! ^	^RR7      p\        VP                  RVR7       \         P                  ! ^^R4      p\        VP                  RVR7       \         P                  ! ^^R4      p\        VP                  RVR7       \         P                  ! ^^R4      p\        VP                  RVR7       \         P                  ! ^^
R4      p\        VP                  R	VR7       \         P                  ! ^^R
4      p\        VP                  RVR7       \         P                  ! ^^R4      p\        VP                  RVR7       R# )r  r  r  r  gX-h?r4   g㤶y?g	y}?g+,?r  g      ?r/  g{Gz?g333333?g<,?Nr  r  s   &  r+   test_two_sided_pvalues2%TestBinomTest.test_two_sided_pvalues2q  s    ooa2.

O$?ooaT*

$5DAoob"d+

O$?ooaT*

O$?ooaR(

G$7ooaB'

Dt4ooaB'

F6r-   c                J   R p\         P                  ! RRR4      p\        VP                  ^VR7       \         P                  ! ^^/R4      p\        VP                  ^VR7       \         P                  ! ^^.R4      p\        VP                  ^VR7       \         P                  ! ^^,R4      p\        VP                  ^VR7       \         P                  ! ^^R4      p\        VP                  ^VR7       \         P                  ! ^^R4      p\        VP                  ^VR7       R# )	r  i  i  r  r4   NgAL &W?g!YB?gE]t?r  r  s   &  r+   test_edge_casesTestBinomTest.test_edge_cases  s    ooc3,

AD1ooaT*

AD1oob"e,

AD1oob"e,

AD1ooaS)

AD1ooaS)

AD1r-   c           	       aa ^
oRo^p\         P                  ! \         P                  ! SS,          4      S^,           4      p\        P                  P                  VSS4      p\         P                  ! \        P                  P                  VSS4      V8*  ^ R7      p\        VV3R lV) \         P                  ! SS,          4      S4      pSV,
          \        V\        P                  P                  VSS4      8H  4      ,           p\        WFRR7       ^p\         P                  ! \         P                  ! SS,          4      ^,           4      p\        P                  P                  VSS4      p\         P                  ! \        P                  P                  VSS4      V8*  ^ R7      p\        VV3R lV^ \         P                  ! SS,          4      4      pV^,           p\        WFRR7       R# )r  r  axisc                 H   < \         P                  P                  V SS4      ) # r$   r   binomrf  r   r  r  s   &r+   <lambda>>TestBinomTest.test_binary_srch_for_binom_tst.<locals>.<lambda>  s    +0;;??2q!+D*Dr-   r7   r4   c                 F   < \         P                  P                  V SS4      # r$   r  r  s   &r+   r  r    s    */++//"a*Cr-   N)r%   rX   ceilr   r  rf  sumr   intr
   floor)	rI   kr\   dy1ixy2r  r  s	   &      @@r+   test_binary_srch_for_binom_tst,TestBinomTest.test_binary_srch_for_binom_tst  sd    IIbgga!enac*KKOOAq!$VVEKKOOAq!,1:) +E+,"bgga!enaA Vc!u{{r1a8899T*IIbhhq1uo)*KKOOAq!$VVEKKOOAq!,1:) +D*+QQA !VT*r-   z"alternative, pval, ci_low, ci_highc                   \         P                  ! ^^dRVR7      p\        VP                  VRR7       \	        VP
                  R4       VP                  RR7      p\        VP                  VP                  3W43RR7       R# )	r  r  r  r  rY  rb  r4   皙?r`  confidence_levelN	r   r  r
   r   r   rF   proportion_cilowhighrI   rY  r9  ci_lowci_highr[   cis   &&&&&  r+   test_confidence_intervals1(TestBinomTest.test_confidence_intervals1  se     oobC4[I

Du5S]]C(5)F+<5Ir-   c                   \         P                  ! ^^2RVR7      p\        VP                  VRR7       \	        VP
                  R4       VP                  RR7      p\        VP                  VP                  3W43RR7       R# )	r|   r  r  r3   r4   gQ?Gz?r  Nr
  r  s   &&&&&  r+   test_confidence_intervals2(TestBinomTest.test_confidence_intervals2  se     ooa2+F

Dt4S]]D)5)F+<4Hr-   zalternative, pval, ci_highc                    \         P                  ! ^ ^
RVR7      p\        VP                  VRR7       VP	                  RR7      p\        VP                  R4       \        VP                  VRR7       R# )	r   r  r  rY  r3   r4   r`  r  r   N)r   r  r
   r   r  r   r  r  )rI   rY  r9  r  r[   r  s   &&&&  r+   !test_confidence_interval_exact_k0/TestBinomTest.test_confidence_interval_exact_k0  sY     ooatE

Dt45RVVS!t4r-   zalternative, pval, ci_lowc                    \         P                  ! ^
^
RVR7      p\        VP                  VRR7       VP	                  RR7      p\        VP                  R4       \        VP                  VRR7       R# )	r  r  r  r3   r4   r`  r  rt  N)r   r  r
   r   r  r   r  r  )rI   rY  r9  r  r[   r  s   &&&&  r+   %test_confidence_interval_exact_k_is_n3TestBinomTest.test_confidence_interval_exact_k_is_n  sY     oob"+F

Dt45RWWc"T2r-   z+k, alternative, corr, conf, ci_low, ci_highc                    \         P                  ! V^
RVR7      pV'       d   RpMRpVP                  WHR7      p	\        V	P                  V	P
                  3WV3RR7       R# )	r  r  r  wilsonccwilson)r	  r  r3   r4   N)r   r  r  r
   r  r  )
rI   r  rY  corrconfr  r  r[   r  r  s
   &&&&&&&   r+   test_ci_wilson_method#TestBinomTest.test_ci_wilson_method  sT    > ooa2+FFFD)F+<4Hr-   c                    \         P                  ! ^^R4      p\        VP                  R4       \        VP                  R4       R# )r~   r  rt  N)r   r  r   rF   r   rI   r[   s   & r+   &test_estimate_equals_hypothesized_prop4TestBinomTest.test_estimate_equals_hypothesized_prop  s2     ooaT*S]]D)SZZ%r-   zk, nc                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! W4       RRR4       R#   + '       g   i     R# ; i)z must be an integer not less thanr   Nr   r   rR   r   r  )rI   r  r  s   &&&r+   test_invalid_k_nTestBinomTest.test_invalid_k_n  s?    ]]:!CE EOOA!E E E Es   AA	c                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! ^^
R4       RRR4       R#   + '       g   i     R# ; i)z+k \(11\) must not be greater than n \(10\).r   r  Nr+  rS   s   &r+   test_invalid_k_too_big$TestBinomTest.test_invalid_k_too_big  sC    ]]:!OQ QOOBD)Q Q Q Q   AA	c                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! ^
^.^R4       RRR4       R#   + '       g   i     R# ; i)zk must be an integer.r   r  N)r   r   r  r   r  rS   s   &r+   test_invalid_k_wrong_type'TestBinomTest.test_invalid_k_wrong_type#  sA    ]]9!8: :OORHb$/: : : :   A

A	c                t   R p\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! ^2^RR7       RRR4       Rp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! ^2^RR7       RRR4       R#   + '       g   i     LZ; i  + '       g   i     R# ; i)zp \(-0.5\) must be in range...r   r  Nzp \(1.5\) must be in range...ru        r+  rI   rP  s   & r+   test_invalid_p_range"TestBinomTest.test_invalid_p_range(  sp    3]]:W55OOBt, 62]]:W55OOBs+ 65 65 655s   B/B&B#	&B7	c                    \         P                  ! ^^
RR7      pRp\        P                  ! \        VR7      ;_uu_ 4        VP                  RR7       RRR4       R#   + '       g   i     R# ; i)r|   r  r  z/confidence_level \(-1\) must be in the intervalr   r  Nr   r   r  r   r   rR   r  )rI   r[   rP  s   &  r+   test_invalid_confidence_level+TestBinomTest.test_invalid_confidence_level0  sJ    ooa2-D]]:W55r2 6555s    AA.	c                    \         P                  ! ^^
RR7      p\        P                  ! \        RR7      ;_uu_ 4        VP                  RR7       RRR4       R#   + '       g   i     R# ; i)r|   r  r  z$method \('plate of shrimp'\) must ber   zplate of shrimpr  Nr=  r'  s   & r+   test_invalid_ci_method$TestBinomTest.test_invalid_ci_method6  sH    ooa2-]]:-TUU%67 VUUUs   AA,	c           	         \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! ^^
RRR7       RRR4       R#   + '       g   i     R# ; i)zalternative \('ekki'\) not...r   r  ekkir  Nr+  rS   s   &r+   test_invalid_alternative&TestBinomTest.test_invalid_alternative;  s6    ]]:-MNNOOAs? ONNNr5  c                x    \         P                  ! ^^
RR7      p\        VP                  VP                  4       R# )r|   r  r  N)r   r  r   proportion_estimaterF   r'  s   & r+   
test_aliasTestBinomTest.test_alias?  s(    ooa2-S,,cmm<r-   z32-bit does not overflowreasonc                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! ^^\
        P                  P                  R7       RRR4       R#   + '       g   i     R# ; i)zError in function...r   r7  N)r   r   OverflowErrorr   r  sys
float_infominrS   s   &r+   test_boost_overflow_raises(TestBinomTest.test_boost_overflow_raisesC  s@     ]]=0FGGOOAqCNN$6$67 HGGGs   2A!!A2	r_   N))r^  g$W?r   g?)r_  g}?gGI~?rt  )	two-sidedg}G?g>?6?g4M?))r^  ga&+w?r   gny?)r_  g5!Rx?g5%?rt  )rU  g~th?gSuR&|?g|KnG?))r^  gլ?g̓A?)r_  rt  rt  )rU  gEow?g;Kڞj?))r^  rt  r   )r_  6>g@>}_?)rU  rV  gcڒ ?)r|   rU  Tr`  g`?g?)r|   rU  Tr  grt!?g/?)r|   rU  Fr`  gVK5?g5dM?)r|   rU  Fr  gXCVu^?g*r[?)r|   r^  Tr`  r   gP-V?)r|   r^  Tr  r   g^'?)r|   r^  Fr`  r   gS5a?)r|   r^  Fr  r   g?)r|   r_  Tr`  g!Շ`ⰸ?rt  )r|   r_  Tr  g$x"?rt  )r|   r_  Fr`  gW~=?rt  )r|   r_  Fr  g	L?rt  )r   rU  Tr`  r   gX?)r   rU  Fr`  r   gN?)r   r^  Tr`  r   g/<7#9?)r   r^  Fr`  r   gAA?)r   r_  Tr`  r   rt  )r   r_  Fr`  r   rt  )r  rU  Tr`  g7?rt  )r  rU  Fr`  g؍s?rt  )r  r^  Tr`  r   rt  )r  r^  Fr`  r   rt  )r  r_  Tr`  ghadn?rt  )r  r_  Fr`  gAB/?rt  )r   r   )r   r{   l        ) r`   ra   rb   rc   __doc__r  r  r  r  r   r   r   r  r  r  r  r$  r(  r,  r/  r3  r:  r>  rB  rF  rJ  skipifrP  maxsizerS  rd   re   rf   s   @r+   r  r  \  s    $@7"2+@ [[AHIJIJ [[AFGIGI [[9DE5	E5 [[8FG3	G3 [[5	=	;	<	=	0	0	1	1	4	4	4	5	5	6	0	0	-	.	6	7	+	,	4	43	6:I;:I& [[Vfg%67" 8"
*
0
,38
@= [[u,5OP8 Q8r-   r  c                   J   a  ] tR tRt o R tR tR tR tR tR t	R t
R	tV tR
# )TestFligneriJ  c                    \         P                  ! ^4      p\        \        P                  ! W^,          4      R^4       R# )r   N)gnf	@gDq?)r%   rX   r   r   fligner)rI   r   s   & r+   r  TestFligner.test_dataL  s,    YYq\!%--E":"L"$	&r-   c                4  a	 \         P                  P                  ^{4      o	V	3R lpV! \        4      pV! \        4      pV! \
        4      p\        P                  ! W#VRR7      w  rV\        P                  ! W#VRRR7      w  rx\        WW4       \        Wh4       R# ){   c                    < \         P                  ! V 4      R SP                  \        V 4      4      ,          ,           P	                  4       # )r  )r%   r  randnr  tolist)gr   s   &r+   _perturb+TestFligner.test_trimmed1.<locals>._perturbX  s1    JJqMEBHHSV,<$<<DDFFr-   rK   r  r  r   r  N)	r%   rd  r   r{  r|  r}  r   r]  r   )
rI   re  g1_g2_g3_Xsq1r  Xsq2r  r   s
   &        @r+   r  TestFligner.test_trimmed1S  s{     YY""3'	G rlrlrl mmCc&AmmCc)479D'E)r-   c                    . ROp. ROp\         P                  ! WRRR7      w  r4\         P                  ! V^R	 V^R	 RR7      w  rV\        W54       \        WF4       R# )
r  r  r  r  rK   r  Nr  r  r   )r   r]  r   )rI   rY   r7  rj  r  rk  r  s   &      r+   r  TestFligner.test_trimmed2f  sU    66mmA49; mmAaGQqWVDD'E)r-   c                x    \         P                  ! R^^4      p\        \        \        P
                  WRR7       R# r  )r%   r  rQ   r  r   r]  r   s   & r+   r  TestFligner.test_bad_keyword  s&    KKAr"i3Gr-   c                x    \         P                  ! R^^4      p\        \        \        P
                  WRR7       R# r  )r%   r  rQ   rR   r   r]  r   s   & r+   r  !TestFligner.test_bad_center_value  s&    KKAr"j%--fEr-   c                F    \        \        \        P                  ^.4       R# rs   )rQ   rR   r   r]  rS   s   &r+   test_bad_num_argsTestFligner.test_bad_num_args  s    j%--!5r-   c                   \         P                  ! ^4      p\        P                  ! \        \
        R7      ;_uu_ 4        \        P                  ! W^,          . 4      p\        VP                  \         P                  4       \        VP                  \         P                  4       RRR4       R#   + '       g   i     R# ; i)r   r   N)r%   rX   r   r   r   r   r   r]  r   rF   r   r   r   s   &  r+   r  TestFligner.test_empty_arg  sj    IIaL\\,4IJJ--a4,C/RVV, KJJJs   A'B00C	r_   N)r`   ra   rb   rc   r  r  r  r  r  rt  r  rd   re   rf   s   @r+   r[  r[  J  s/     &*&
*4HF6- -r-   r[  c               #     "   . ROp . ROp\        V4       F  w  r#\        P                  P                  V4      pVP                  ^d4      pVP	                  ^ ^c^R7      pVP	                  ^^^R7      p\        Wg4       F/  w  r\        V^,           W,           4       F  p
WX,          WZ&   K  	  K1  	  VP                  V4       \        P                  ! V^4      w  rWR.W,          O5x  K  	  R# 5i)g]YC?)r  r  r   r^  N))g]YCgQɣ?)gme5gϵ5?)g;jgH?)iFeiHcMi%:)		enumerater%   rd  re  integersziprangeshufflesplit)expected_resultsseedssir  rh  xytie_indnum_ties_per_indr\   r  jrY   r7  s                r+   mood_cases_with_tiesr    s     @ .Ee$ii##D)ZZ_,,12A,6<<AAA<>2DA1q5!%( ) 3 	BxxAF1-111 %s   C.C0c                   "  a  ] tR tRt o ]P
                  P                  R]! 4       4      R 4       t]P
                  P                  R. RO4      R 4       t	R t
R tR tR	 tR
 tR t]P
                  P                  R. RO4      R 4       tRtV tR# )TestMoodi  z$x,y,alternative,stat_expect,p_expectc                n    \         P                  ! WVR7      w  rg\        WFRR7       \        WWRR7       R# )a  
Example code used to generate SAS output:
DATA myData;
INPUT X Y;
CARDS;
1 0
1 1
1 2
1 3
1 4
2 0
2 1
2 4
2 9
2 16
ods graphics on;
proc npar1way mood data=myData ;
   class X;
    ods output  MoodTest=mt;
proc contents data=mt;
proc print data=mt;
  format     Prob1 17.16 Prob2 17.16 Statistic 17.16 Z 17.16 ;
    title "Mood Two-Sample Test";
proc print data=myData;
    title "Data for above results";
  run;
rX  gؗҜ<r0  Nr   moodr
   )rI   rY   r7  rY  stat_expectp_expectrF   r   s   &&&&&&  r+   test_against_SASTestMood.test_against_SAS  s+    < "JJqE	U;u5r-   zalternative, expectedc                \    . ROp. ROp\         P                  ! W4VR7      p\        WR4       R# )r  rX  Nr=  rE  r  )rI   rY  ro  rY   r7  r[   s   &&&   r+   test_against_SAS_2TestMood.test_against_SAS_2  s(    <=jj;7&r-   c                D   \         P                  P                  R 4       \         P                  P                  ^
^4      p\         P                  P                  ^^4      p\        P
                  ! W4      w  r4\        P
                  ! W!4      w  rV\        W4.V) V.4       R# )r  N)r%   rd  r  rb  r   r  r   )rI   r   r   z1p1z2p2s   &      r+   test_mood_order_of_args TestMood.test_mood_order_of_args  sn     			tYY__R#YY__R#B#B#!2(bS"I6r-   c                    . ROp. ROp\         P                  ! V4      p\         P                  ! V4      pRVn        RVn        \        \        P
                  ! WRR7      RR.4       R# )	g?Nr  g^}E	?)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{   )rW   r{   g4)r%   r9  shaper   r   r  )rI   r   r   s   &  r+   test_mood_with_axis_none!TestMood.test_mood_with_axis_none  sY    4	H XXb\XXb\!%**R$"?#1="A	Cr-   c           
        ^p\         P                  P                  R4       \         P                  P                  ^
V4      p\         P                  P                  ^V4      p\        P
                  ! W#4      w  rE\        V4       FD  p\        WF,          WV,          .\        P
                  ! VRV3,          VRV3,          4      4       KF  	  VP                  4       pVP                  4       p\        P
                  ! W#^R7      w  rE\        V4       FB  p\        WG,          WW,          .\        P
                  ! W'R3,          W7R3,          4      4       KD  	  R# )r   r  NNNr  N)	r%   rd  r  rb  r   r  r|  r   	transpose)rI   nyr   r   	z_vectestpval_vectestr  r\   s   &       r+   test_mood_2dTestMood.test_mood_2d  s     
		tYY__R$YY__R$"'**R"4	rA%y|\_&E&+jjAqD2ad8&DF 
 \\^\\^"'**R!"<	rA%y|\_&E&+jjqD2d8&DF r-   c           
        Rp\         P                  P                  R4       \         P                  P                  ! V!  p\         P                  P                  ! V!  p\	        ^4       F  p\
        P                  ! W#VR7      w  rV^^.^ ^.^ ^.3p\	        WV,          ^ ,          ,          4       F  p\	        WV,          ^,          ,          4       F  p	V^ 8X  d   VRW3,          p
VRW3,          pM4V^8X  d   W(RV	3,          p
W8RV	3,          pMW(V	R3,          p
W8V	R3,          p\        WXV	3,          WhV	3,          .\
        P                  ! W4      4       K  	  K  	  K  	  R# )r  r  r  r  N)r  r   r1   )r%   rd  r  rb  r|  r   r  r   )rI   r  r   r   r  r  r  axes_idxr\   r  slice1slice2s   &           r+   test_mood_3dTestMood.test_mood_3d  s5   
		tYY__e$YY__e$!HD&+jjd&C#I AAA/H5$!234ud^A%678Aqy!#AqG!#AqG!#q!G!#q!G!#q!G!#q!G-yA/;qD/A/C.3jj.HJ 9 5 r-   c                P   \         P                  ! \        \        R 7      ;_uu_ 4        \        P
                  ! ^.. 4      p\        VP                  \        P                  4       \        VP                  \        P                  4       RRR4       R#   + '       g   i     R# ; ir   )r   r   r   r   r   r  r   rF   r%   r   r   r'  s   & r+   test_mood_bad_argTestMood.test_mood_bad_arg8  s[    \\,4IJJ**aS"%C/RVV, KJJJs   A!BB%	c                   \         P                  P                  ^ 4      p\        P                  P                  R^dVR7      p\        P                  P                  R^dVR7      p\        P                  ! W#RR7      w  rE\        P                  ! W#RR7      w  rg\        P                  ! W#RR7      w  rYFu;8X  d	   V8X  g   Q h Q h\        V^ RR	7       \        Wu^,          4       \        V	^V^,          ,
          4       \        P                  ! \        R
R7      ;_uu_ 4        \        P                  ! W#RR7       RRR4       R#   + '       g   i     R# ; i)r         ?)r   r   r   g      ?rU  rX  r^  r_  Hz>r0  z`alternative` must be...r   z	ekki-ekkiN)r%   rd  r   r   r   r(   r  r
   r   r   rR   )
rI   rh  rY   r7  stat1r  stat2r  stat3p3s
   &         r+   test_mood_alternativeTestMood.test_mood_alternative?  s    ii##A&JJNNCcNBJJNNCcNBJJq=	JJq8	JJq;	&&&&&&AD)qD!A1H%]]:-GHHJJq5 IHHHs   #EE	rY  c                    \         P                  P                  R 4      pVP                  R4      pVP                  R4      p\        P
                  ! W4VR7      p\        VP                  VP                  3V4       R# )l	   hTt fU6H~ rX  N)r  rt   )rW   rt   )	r%   rd  re  r   r   r  r   rF   r   )rI   rY  rh  r   r   r[   s   &&    r+   test_resultTestMood.test_resultQ  s^    ii##$KL  )  )jj[9cmmSZZ0#6r-   r_   N))rU  )3oXQ?gpL?)r^  )r  gԣl?)r_  )r  gpL?)rU  r^  r_  )r`   ra   rb   rc   r   r   r   r  r  r  r  r  r  r  r  r  r  rd   re   rf   s   @r+   r  r    s     [[C135656@ [[4@A'A'7C8F.J8-6$ [[],LM7 N7r-   r  c                      a  ] tR tRt o R tR tR t]P                  P                  ]P                  P                  ]'       * RR7      R 4       4       tR tR	 tR
 tRtV tR# )TestProbploti[  c                ,   \         P                  P                  ^RR7      p\         P                  ! VRR7      w  r#. ROp\	        V\
        P                  ! V4      4       \	        W$4       \         P                  ! VRR7      w  rV. ROp\	        Wg4       R# )	r  90  rs  Fr  TN)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+/?)gƝV?g?gà?)r   r   r(   probplotr
   r%   sort)rI   rY   osmosrosm_expectedr[   res_fitres_fit_expecteds   &       r+   rN   TestProbplot.test_basic]  sn    JJNNN7>>!/G
 	RWWQZ(*~~aT2?2r-   c                x   \         P                  P                  ^dRR7      p\         P                  ! VRRR7      w  r#\         P                  ! V^ RR7      w  rE\         P                  ! VRRR7      w  rg\	        W$4       \	        W&4       \	        W54       \	        W74       \         P                  ! VRRR7      w  rR# )r  i@ rs  NFsparamsr  r_   )r   r   r(   r  r
   )
rI   rY   osm1osr1osm2osr2osm3osr3r  r  s
   &         r+   test_sparams_keyword!TestProbplot.test_sparams_keywordl  s    JJNN&N9 ^^At?
^^Aqe<
^^Aru=
####>>!RU;Sr-   c                8   \         P                  P                  ^RR7      p\         P                  ! VRRRR7      w  r#\         P                  ! VR\         P                  RR7      w  rE\        W$4       \        W54       \        \        \         P                  VRR7       \        \        \         P                  V. R7        ! R R	4      p\         P                  ! VRRR
7      w  r#\         P                  ! W! 4       RR7      w  rE\        W$4       \        W54       R# )r  r  rs  Ft)r  r  r  zwrong-dist-namer  c                   *   a  ] tR tRt o RtR tRtV tR# )3TestProbplot.test_dist_keyword.<locals>.custom_disti  z6Some class that looks just enough like a distribution.c                D    \         P                  P                  V^R7      # )r{   )r   )r   r   ppf)rI   qs   &&r+   r  7TestProbplot.test_dist_keyword.<locals>.custom_dist.ppf  s    zz~~aQ~//r-   r_   N)r`   ra   rb   rc   rW  r  rd   re   rf   s   @r+   custom_distr    s     H0 0r-   r  r  )r  r  N)r|   )r{   )	r   r   r(   r  r  r
   rQ   rR   AttributeError)rI   rY   r  r  r  r  r  s   &      r+   test_dist_keywordTestProbplot.test_dist_keywordz  s    JJNNN7^^A5sDI
^^A5uwwM
##j%..!:KLnennabA	0 	0
 ^^At?
^^AKMuE
##r-   no matplotlibrL  c                   \         P                  ! 4       pVP                  ^o4       \        P                  P                  ^^dRR7      p\        P                  ! V\         R7      w  r4\         P                  ! 4        \        P                  ! VRR7      w  rV\        P                  ! VR\         R7      p\         P                  ! 4        \        P                  ! VRRR7      p\        \        V4      \        V4      u;8H  ;'       d4    \        V4      u;8H  ;'       d    \        V4      u;8H  ;'       d    ^8H  Mu 4       \        W54       \        W74       \        W84       \        WF4       \         P                  ! 4       pVP                  ^o4      p	\        P                  ! VRV	R7       \         P                  ! 4        R# )r  t rs  plotNF)r  r  )pltfigureadd_subplotr   r  r(   r  closer	   r  r
   )
rI   figrY   res1fitres1res2fitres2res3res4axs
   &         r+   test_plot_kwargTestProbplot.test_plot_kwarg  s    jjlGGKK'K:qs3		qt4~~aU5		~~aU6 	D	SYEE#d)EEs4yEEAEF###) jjl__S!qe"-		r-   c                J    \        \        \        P                  ^.RR7       R# r  )rQ   rR   r   r  rS   s   &r+   test_probplot_bad_args#TestProbplot.test_probplot_bad_args  s    j%..1#<MNr-   c                x   \        \        P                  ! . R R7      \        P                  ! . 4      \        P                  ! . 4      34       \        \        P                  ! . RR7      \        P                  ! . 4      \        P                  ! . 4      3\        P
                  \        P
                  R334       R# )Fr  Tr   N)r   r   r  r%   r9  r   rS   s   &r+   
test_emptyTestProbplot.test_empty  ss    U^^BE2hhrlBHHRL1	3U^^BD1xx|RXXb\2vvrvvs+-	.r-   c                   R p\         P                  ! RR7      ;_uu_ 4        \        P                  ! \        VR7      ;_uu_ 4        \        \        P                  ! ^.RR7      \         P                  ! R.4      \         P                  ! ^.4      3\         P                  \         P                  \         P                  334       RRR4       RRR4       R#   + '       g   i     L; i  + '       g   i     R# ; i),One or more sample arguments is too small...r   invalidr   Tr  r   N)
r%   r   r   r   r   r   r   r  r9  r   r9  s   & r+   test_array_of_size_one#TestProbplot.test_array_of_size_one  s    @kk(++ll-W==688RD>288QC=966266266245 > ,+== ,++s$   %C+A?C	C+C(#C++C<	r_   N)r`   ra   rb   rc   rN   r  r  r   r   rr  rX  have_matplotlibr  r  r  r  rd   re   rf   s   @r+   r  r  [  sg     3<$( [[[[O+OD E 0O.5 5r-   r  c                     a  ] tR tRt o R tR tR tR tR tR t	R t
R	 tR
 tR tR tR t]P"                  P%                  R. RO. RO. RO.4      R 4       tR t]P"                  P*                  R 4       t]P"                  P%                  R. R O4      R 4       tR t]P"                  P%                  RRR.4      R 4       t]P"                  P%                  RR]P6                  ! 4       34      R 4       tR tRtV tR# )!TestWilcoxoni  c                    \        \        \        P                  ^^.^^.R4       \        \        \        P                  ^^.^^.RR7       \        \        \        P                  ^.^
,          RR7       R# )rt   dummyrX  xyzrA  N)rQ   rR   r   wilcoxonrS   s   &r+   test_wilcoxon_bad_arg"TestWilcoxon.test_wilcoxon_bad_arg  sW     	j%..1a&1a&'Jj%..1a&1a&")	+j%..1#b&Gr-   c           	        \         P                  ! ^4      p\         P                  ! RR7      ;_uu_ 4        \        P                  ! WRRR7      w  r#\        W#3R\         P                  34       \        P                  ! WRRR7      w  r#\        W#3R\         P                  34       RRR4       \        \        P                  ! WR	RR7      R
4       R#   + '       g   i     L5; i)r  r   r  wilcox
asymptoticrA  r   prattNzsplit)g     @Z@rt  )r%   rX   r   r   r  r   r   )rI   rY   r   r  s   &   r+   test_zero_diffTestWilcoxon.test_zero_diff  s    IIbM [[**>>!FDA!#rvv/>>!EDA!#rvv/	 + 	U^^A(<H#	% +*s   A1CC$	c                b    . ROp. ROp\         P                  ! WRRRR7      p\        VR4       R# )	rt   r  r
  F)zero_methodr  
correctionNrN  )rt   r{   r|   r   )r   g|N?)r   r  r
   )rI   rY   r7  r[   s   &   r+   
test_prattTestWilcoxon.test_pratt  s0     nnQw|(-/78r-   c                    . ROp\         P                  ! VRRR7      p\         P                  ! VRRR7      p\         P                  ! VRRR7      pR# )rt   r  r
  )r  r  r  r	  N)rt   r{   r|   r   r   r|   rt   r{   rt   rt   r{   )r   r  )rI   arrr  s   &  r+   test_wilcoxon_arg_type#TestWilcoxon.test_wilcoxon_arg_type  s?     1NN3GLINN3H\JNN3H\Jr-   c                T   . ROp\        R^4      p\        P                  ! \        W!4       UUu. uF  w  r4V.V,          NK  	  upp4      p\        P                  ! VP
                  4      p\        P                  ! WVRRRR7      w  rx\        VR4       \        VR4       \        P                  ! WVRRRR7      w  rx\        VR4       \        VR	4       \        P                  ! WVR
RRR7      w  rx\        VR4       \        VR4       \        P                  ! . RO4      p\        P                  ! . RO4      p\        P                  ! WVRRR7      w  rx\        V^"4       \        VRRR7       \        P                  ! WVRRR7      w  rx\        V^"4       \        VRRR7       R# u uppi )rt   r  r
  Fr  r  i  g.i?r  i  gFYv]Uj?r	  iG  g%^ Ez?r  r  g^҂<?r3   r4   Tg:b+?N)	rt   r~   r  rW   r@   r~   r   rt   r{   r  rH  r  r        y   r     rB     p         w   r4  r"     r3  rB  s   r  r     )r|  r%   concatenater{  zerosr   r   r  r
   r9  r   )	rI   freqnumsur   rY   r7  r  r  s	   &        r+   test_accuracy_wilcoxon#TestWilcoxon.test_accuracy_wilcoxon  sS   ,R|NND@QC!GG@AHHQVV~~aGL).0301~~aH\).0301~~aH\).03+, HHQRHHQR~~au\JQ940~~atLIQ9407 As   F$
c                d   \         P                  ! . RO4      p\         P                  ! . RO4      p\        P                  ! WRRR7      p\        P                  ! WRRR7      p\        P                  ! WRRR7      p\        P                  ! WRRR7      pY4u;8X  d   Tu;8X  d	   V8X  g   Q h Q hR# )	r|   r	  approxrA  r
  )modeN)r|   r   r  r?      r  r  r{   r@   r8  r=   r  )r{   r   rt   r  r   r   r>   r  r>  r>      rM  )r%   r9  r   r  )rI   rY   r7  r  r  r  r  s   &      r+   test_approx_modeTestWilcoxon.test_approx_mode  s     HH@AHHAB~~aHX>~~aH\B~~aH8<~~aH<@+t+t+++++r-   c                    \         P                  ! . RO4      p\         P                  ! . RO4      p\        P                  ! WRRR7      pRp\	        W44       R# )r  Fr
  r  Nr  r#  rT  )r%   r9  r   r  r   )rI   rY   r7  r[   rZ   s   &    r+   test_wilcoxon_result_attributes,TestWilcoxon.test_wilcoxon_result_attributes  s?    HHQRHHQRnnQeLI,
C,r-   c                   \         P                  P                  R 4      pVP                  ^4      VP                  ^4      r2\        P                  ! W#RR7      p\        P
                  P                  VP                  ^,          4      p\        VP                  V4       \        P                  ! W#RR7      p\        VR4      '       d   Q h\        P                  ! W#4      p\        VR4      '       d   Q hR# )l   	m$S r
  rA  exact
zstatisticN)r%   rd  re  r   r  r   r  r   r
   r=  hasattr)rI   rh  rY   r7  r[   r   s   &     r+   test_wilcoxon_has_zstatistic)TestWilcoxon.test_wilcoxon_has_zstatistic  s    ii##K0zz"~szz"~1nnQ,7jjnnSZZ\*,nnQ'23----nnQ"3-----r-   c                    \         P                  ! R .^
,          RRR7      w  rRp\        V^ 4       \        W#RR7       \         P                  ! R .^
,          RRR7      w  rR	p\        V^ 4       \        W#RR7       R
# )r  r
  Fr  g4IťY?r3   r4   Tr  gc2_?N)r   r  r   r
   )rI   statr  
expected_ps   &   r+   test_wilcoxon_tieTestWilcoxon.test_wilcoxon_tie$  sp     ..#L,13 
T1D1..#(46 
T1D1r-   c                   . ROp. ROp\         P                  ! WRRRR7      w  r4\        V^4       \        VR^R7       \         P                  ! WRRRR7      w  r4\        V^4       \        VR	^R7       \         P                  ! WR
RRR7      w  r4\        V^4       \        VR^R7       \         P                  ! WR
RRR7      w  r4\        V^4       \        VR^R7       R# )}   r^  r
  F)rY  r  r  g13}?r9   T)rY  r  r  g7a%?r_  g?gڏoJT?N)
rG  r(  r     rH  r(  rH  rG  rH  r5  )
rD  z   rG  r  rH  |   r`     r5     )r   r  r   r   )rI   rY   r7  r   r  s   &    r+   test_onesidedTestWilcoxon.test_onesided9  s     ?>~~a|).0QAy!4~~a4%13QAy!4~~a	%1eEQAy!4~~a	)-lDQAy!4r-   c                    \        ^^34       Fi  p\        V4      p\        V4      p\        W^,           ,          ^,          ^,           \	        V4      4       \        \        V4      ^4       \        W#4       Kk  	  R# rs   )r|  r   r   r   r  r  r   )rI   r  pmf1pmf2s   &   r+   test_exact_basicTestWilcoxon.test_exact_basicY  sX    q"A&q)D'*DaCQD	2TA&%d1 r-   c                   \         P                  ! . RO4      p\         P                  ! . RO4      p\        P                  ! WRRR7      w  r4\	        VR^R7       \        P                  ! WRRR7      w  r4\	        VR^R7       \        P                  ! WRRR7      w  r4\	        VR	^R7       \         P
                  ! ^ ^4      R
,           p\         P
                  ! ^^ R4      p\        P                  ! WRRR7      w  r4\	        VR^R7       \        P                  ! WRRR7      w  r4\	        VR^R7       \        P                  ! WRRR7      w  r4\	        VR^R7       R# )(\?rU  r<  rY  r  g֔  ?r9   r^  gYJ?r_  g*  ?r  g(?g(?gS?N)
rU  g=
ףp=?g(\?gQ޿gQ?g{Gz?gp=
ףgq=
ףp?g      gQ?)
r   g?gɿg333333?皙gܿgzGgQοg(\gRQr   )r%   r9  r   r  r   rX   )rI   rY   r7  r  r  s   &    r+   test_exact_pvalTestWilcoxon.test_exact_pvala  s   HH # $HH $ %~~aGLAy!4~~awGAy!4~~a	'JAz15IIas"IIb!R ~~aGLAy!4~~awGAy!4~~a	'JAy!4r-   rY   c                    \         P                  ! V4      w  r#\        P                  ! V4      pW^ 8  ,          P	                  4       p\        W$4       \        V^4       R# r   N)r   r  r%   r9  r  r   )rI   rY   r   r  wtrues   &&   r+   test_exact_p_1TestWilcoxon.test_exact_p_1}  sD     ~~a HHQKa%QQr-   c           	     l   \         P                  ! ^ ^24      R,           p\         P                  ! ^2^ R4      p\        \        P                  ! W4      \        P                  ! WRR7      4       \        P
                  ! 4       p\         P                  ! R^4      p\        P                  ! V4      w  rV\        WV3\        P                  ! WCR7      4       \         P                  ! R^	4      p\        P                  ! V4      w  rV\        WV3\        P                  ! VRR7      4       ^WD^ 8H  &   \        P                  ! V4      w  rV\        WV3\        P                  ! VRR7      4       \         P                  ! ^^44      p\        \        P                  ! V4      \        P                  ! VRR7      4       R# )	r   r  r<  rA  r
  Nr   r   rc  )r%   rX   r   r   r  r  )rI   rY   r7  pmr  r   r  s   &      r+   	test_autoTestWilcoxon.test_auto  s,   IIas"IIb!R U^^A)^^A9	; $$&IIb!~~a 
 	aVU^^A9: IIb!~~a aVU^^AlCDq&	~~a aVU^^AlCD IIaU^^A&q(NOr-   c                    \         P                  ! R^4      p\        P                  ! V4      pRp\	        W#4       ^W^ 8H  &   \        P                  ! V4      pRp\	        W#4       R# )r   Nrc  )g     ;@g     P?)r  g     0?)r%   rX   r   r  r   )rI   r  r[   r   s   &   r+   test_auto_permutation_edge_case,TestWilcoxon.test_auto_permutation_edge_case  sU     IIb!nnQ"Sq&	nnQ Sr-   r   c                   \         P                  P                  R 4      pVP                  VR7      p\        P                  ! V\        P
                  ! 4       R7      p\        P                  ! VRR7      p\        VP                  VP                  4       \        VP                  VP                  4       VP                  V^
,          R7      p\         P                  P                  R4      p\        P
                  ! ^cVR7      p\        P                  ! W6R7      p\         P                  P                  R4      p\        P
                  ! ^cVR7      p\        P                  ! W6R7      p\        \         P                  ! VP                  ^4      VP                  4       \        VP                  VP                  4       R# )l   aQG1X@
r   rA  r<  l   VsWFr  )r  r   N)
r%   rd  re  r   r  r  r   rF   r   round)rI   r   rh  rY   r[   r   r`  s   &&     r+   test_permutation_method$TestWilcoxon.test_permutation_method  s!   ii##$56JJDJ!nnQu'>'>'@AnnQw/S]]CMM2SZZ,JJDGJ$ii##$56$$=nnQ*ii##$56$$#FnnQ*RXXcjj!,cjj9SZZ,r-   c                P   \         P                  P                  R 4      pVP                  RR7      p\         P                  VR&   \
        P                  ! V4      p\
        P                  ! VRR7      p\        W44       \        VR4      '       g   Q h\        VR4      '       d   Q hR# )l   HF )Z r   r
  rA  r=  N)r  r{   )r   rt   )	r%   rd  re  normalr   r   r  r
   r>  )rI   rh  r   r[   r   s   &    r+   6test_method_auto_nan_propagate_ND_length_gt_50_gh20591CTestWilcoxon.test_method_auto_nan_propagate_ND_length_gt_50_gh20591  s     ii##O4JJGJ$&&$nnQnnQ|4!sL))))3-----r-   r  r<  r
  c                   . ROp. ROp\         P                  ! W#RVR7      p\         P                  ! W2RVR7      p\        V4      \        V4      ^,           ,          ^,          p\        VP                  4      VP                  8w  g   Q h\        WeP                  ,
          VP                  RR7       \        VP                  VP                  RR7       R# )	r  r^  rV  r_  V瞯<r4   N)
r  r  rW  r  rD  r  r  r  r  rL  )
r  r  E   rW  rT  r  r  r  r  r  )r   r  r  r  rF   r
   r   )rI   r  var1var2r   r[   max_statistics   &&     r+   test_symmetry_gh19872_gh20752*TestWilcoxon.test_symmetry_gh19872_gh20752  s     87nnTVFKnnTYvND	SY]3a73==!S]]2225s}}5Q

CJJU;r-   c                z    \         P                  ! \        P                  ! ^4      VR7      p\	        V^ ^.4       R# )r   rA  N)r   r  r%   r+  r
   )rI   r  r[   s   && r+   test_all_zeros_exact!TestWilcoxon.test_all_zeros_exact  s*    
 nnRXXa[8aV$r-   c           	     R   R p\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RO^^.4       RRR4       Rp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RO^^.RR7       RRR4       \        \        RR4      ;'       g    \        P                  P                  pRp\         P                  ! W!R7      ;_uu_ 4        \        P                  ! . RO. RO^R7       RRR4       R	p\         P                  ! W!R7      ;_uu_ 4        \        P                  ! . RO. RO^RR
7       RRR4       R#   + '       g   i     EL!; i  + '       g   i     L; i  + '       g   i     L; i  + '       g   i     R# ; i)z/Array shapes are incompatible for broadcasting.r   Nz0operands could not be broadcast together with...T)_no_deco	AxisErrorz8source: axis 3 is out of bounds for array of dimension 1r  z+`axis` must be compatible with the shape...)r  rz  rz   r~   r   r1   )	r   r   rR   r   r  rf  r%   
exceptionsr{  )rI   rP  r{  s   &  r+   .test_wilcoxon_axis_broadcasting_errors_gh22051;TestWilcoxon.test_wilcoxon_axis_broadcasting_errors_gh22051  s
    D]]:W55NN9q!f- 6 E]]:W55NN9q!ft< 6 BT2MMbmm6M6M	L]]944NN9ia8 5 @]]944NN9ia$G 54 655 65
 54 544s/   E1E/+F2FE,	/E?	F	F&	r_   N)r   r   r|   )r   r{   r8  r  r   )r   r   r|   r  rc  ir?   r@   )r|   r   r  ) r`   ra   rb   rc   r  r  r  r  r/  r6  r9  r?  rD  rM  rR  rX  r   r   r   r]  ra  r/  rd  rh  rl  rt  r   r  rw  r~  rd   re   rf   s   @r+   r  r    s3    H%9K1@	,-.2*5@258 [[S;#5#@#B CCP< [[  [[VZ0- 1-(.  [[X'>?
< @
< [[X1H1H1J'KL% M%H Hr-   r  c                      a  ] tR tRt o R t]P                  P                  R4      R 4       tR t	]P                  P                  R^ R.4      R 4       t]P                  P                  R	. RO4      R
 4       tRtV tR# )	TestKstati  c           	     >   \         P                  P                  R 4      pVP                  VP	                  R4      VP
                  R7      pTP                  R
 Uu. uF  p\        P                  ! W44      NK  	  up4      pVP                  . ROVP                  R7      p\        WVRR7       \        P                  ! V^R7      p\        P                  ! V^R7      p\        P                  ! V^R7      p	\        VP                  WxV	34      VRR RRR	7       R# u upi )i}  r  r  r  r4   )orderNg{Gz?{Gz?r  rN  )gHQ,?glZ)rI?g88}?gcJN?r   )r%   rd  r   r  rb  float64stackr   kstatr  r   moment)
rI   r  rh  rJ   r  momentsro  m1m2m3s
   &&        r+    test_moments_normal_distribution*TestKstat.test_moments_normal_distribution  s    ii##E*zz#))E*"**z=((,G,QEKK0,GH::I$(JJ  05 \\$a(\\$a(\\$a(""."DtT Hs   D1ignore:invalid value encountered in scalar dividec                   \        V4      '       dX   \        P                  ! \        \        R 7      ;_uu_ 4        \
        P                  ! VP                  . 4      4      pRRR4       MM\        P                  ! RR7      ;_uu_ 4        \
        P                  ! VP                  . 4      4      pRRR4       \        XVP                  VP                  4      4       R#   + '       g   i     L7; i  + '       g   i     LJ; ir   Nr   r  )r   r   r   r   r   r   r  r  r%   r   r   r   )rI   r  r[   s   && r+   rT   TestKstat.test_empty_input  s    B<<08MNNkk"**R.1 ON X..kk"**R.1 /RZZ/0 ON /.s   &C&C0C-	0D 	c                    VP                  R 4      pVP                  V^8H  VP                  V4      p\        \        P
                  ! V4      VP                  VP                  4      4       R# r   rX   wherer   r   r   r  r  rI   r  rJ   s   && r+   r   TestKstat.test_nan_input%  I    yy~xx	26640D)2::bff+=>r-   r  g/$@c                    VP                  ^
4      pRp\        P                  ! \        VR7      ;_uu_ 4        \        P
                  ! W1R7       RRR4       R#   + '       g   i     R# ; i)r  z'k-statistics only supported for 1<=n<=4r   r  N)rX   r   r   rR   r   r  )rI   r  r  rJ   rP  s   &&&  r+   test_kstat_bad_argTestKstat.test_kstat_bad_arg+  sB     yy};]]:W55KK" 6555s   AA+	casec                    Vw  r4\         P                  ! VP                  \        4      V4      p\	        WRP                  V4      4       R # r$   )r   r  r  x_kstatr   )rI   r  r  r  r   r[   s   &&&   r+   test_against_RTestKstat.test_against_R3  s4     kk"**W-q1ZZ_-r-   r_   N))rt   gM~,@)r{   2.L)@)r|   gP'')r~    -"ba)r`   ra   rb   rc   r  r   r   r  rT   r   r   r  r  rd   re   rf   s   @r+   r  r    s     U [[ ST1 U1? [[S1e*-# .# [[V &? @.	@.r-   r  c                      a  ] tR tRt o ]P
                  P                  R4      R 4       tR t]	! RRR7      R 4       t
R	 tR
tV tR# )TestKstatVariE  r  c                   VP                  . 4      p\        V4      '       dI   \        P                  ! \        \
        R 7      ;_uu_ 4        \        P                  ! V4      pRRR4       M>\        P                  ! RR7      ;_uu_ 4        \        P                  ! V4      pRRR4       \        XVP                  VP                  4      4       R#   + '       g   i     L7; i  + '       g   i     LJ; ir  )r  r   r   r   r   r   r   kstatvarr%   r   r   r   )rI   r  rY   r[   s   &&  r+   rT   TestKstatVar.test_empty_inputF  s    JJrNB<<08MNNnnQ' ON X..nnQ' /RZZ/0 ON /.s   C
C#C 	#C3	c                    VP                  R 4      pVP                  V^8H  VP                  V4      p\        \        P
                  ! V4      VP                  VP                  4      4       R# r   r  r  s   && r+   r   TestKstatVar.test_nan_inputQ  r  r-   Tz2input validation of `n` does not depend on backendnp_onlyrM  c                    ^.p^
pRp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! W#R7       RRR4       R#   + '       g   i     R# ; i)rt   zOnly n=1 or n=2 supported.r   r  N)r   r   rR   r   r  )rI   r  rJ   r  rP  s   &&   r+   rv   TestKstatVar.test_bad_argW  sA     s.]]:W55NN4% 6555s   AA	c                   \        \        4      pR pRp\        P                  ! VP	                  \        4      ^4      pW2,          p\        WQP	                  V4      4       \        P                  ! VP	                  \        4      ^4      p^V^,          ,          V,          V^,
          V,          ,           W"^,           ,          ,          p\        WQP	                  V4      4       R# )r  Nr  )r  r  r   r  r  r   )rI   r  r  k2k4r[   r   s   &&     r+   test_against_R_mathworld%TestKstatVar.test_against_R_mathworlda  s    
 LnnRZZ0!4fZZ_-nnRZZ0!4QwqyAaC8#qS	2ZZ_-r-   r_   N)r`   ra   rb   rc   r   r   r  rT   r   skip_xp_backendsrv   r  rd   re   rf   s   @r+   r  r  E  sW     [[ ST1 U1? dQS&S&. .r-   r  c                      a  ] tR tRt o R tR tR t]P                  P                  ]
'       * RR7      R 4       tR tR	 tR
tV tR# )TestPpccPlotit  c                :    \        ^RRR7      ^,           V n        R# r     r  rs  Nr,   rY   rS   s   &r+   setup_methodTestPpccPlot.setup_methodu      "13WEIr-   c           	         ^p\         P                  ! V P                  R^
VR7      w  r#. ROp\        V\        P
                  ! R^
VR7      4       \        W44       R# )r   NnumN)g9%	?gxM _?g";9?g®XZ?g?)r   	ppcc_plotrY   r
   r%   r  )rI   r  svalsppccppcc_expecteds   &    r+   rN   TestPpccPlot.test_basicx  sG    oodffc2;%r{{3:;,r-   c                j   \         P                  ! V P                  R^
RR7      w  r\         P                  ! V P                  R^
\         P                  R7      w  r4\	        WRR7       \	        W$RR7       \         P                  ! V P                  R^
4      w  rV\	        WRR7       \	        W&RR7       R# )r  tukeylambdar  g#B;r4   Nr  )r   r  rY   r  r
   )rI   svals1ppcc1svals2ppcc2svals3ppcc3s   &      r+   	test_distTestPpccPlot.test_dist  s    RmLR-2->->@U351R8U351r-   r  rL  c                X   \         P                  ! 4       pVP                  ^o4      p\        P                  ! V P
                  R^\         R7       VP                  V4       VP                  ^o4      p\        P                  ! V P
                  R^VR7       \         P                  ! 4        R# r  r  Ni)r  r  r  r   r  rY   delaxesr  rI   r  r  s   &  r+   r  TestPpccPlot.test_plot_kwarg  sm     jjl__S!Rc2B __S!Rb1		r-   c           	         \        \        \        P                  V P                  ^^ 4       \        \        \        P                  . RO^ ^RR7       R# )rt   r  r  Nrz   )rQ   rR   r   r  rY   rS   s   &r+   test_invalid_inputs TestPpccPlot.test_invalid_inputs  s5    j%//4661a@ 	j%//9a,	.r-   c           	         \         P                  ! . ^ ^4      w  r\        V\        P                  ! ^ ^^PR7      4       \        V\        P
                  ! ^P\        R7      4       R# )r   r  r  N)r   r  r
   r%   r  r+  float)rI   r  r  s   &  r+   r  TestPpccPlot.test_empty  sD     oob!Q/r{{1aR89bhhr78r-   rY   N)r`   ra   rb   rc   r  rN   r  r   r   rX  r  r  r  r  rd   re   rf   s   @r+   r  r  t  sP     J-
2 [[O+OD
 E
.9 9r-   r  c                   8   a  ] tR tRt o R tR tR tR tRtV t	R# )TestPpccMaxi  c                N    ^.p\        \        \        P                  VRR7       R# r  )rQ   rR   r   ppcc_maxru   s   & r+   test_ppcc_max_bad_arg!TestPpccMax.test_ppcc_max_bad_arg  s    sj%..$=NOr-   c                    \         P                  P                  R^RRRR7      R,           p\        \         P                  ! V4      R	^R7       R# )
ffffff?r  '   r        @r9   Nffffff{)r   r  r(   r   r  r   s   & r+   test_ppcc_max_basicTestPpccMax.test_ppcc_max_basic  sF    !!$ASu/6 " 8:=>ENN1-/CQOr-   c                X   \         P                  P                  R
^RRRR7      R,           p\         P                  ! VRR7      p\         P                  ! V\         P                  R7      p\	        VR^R7       \	        VR^R7       \         P                  ! V4      p\	        VR^R7       R	# )r  r  r  r  r   r  r  r  r9   Nr  r  )r   r  r(   r  r   )rI   rY   max1max2max3s   &    r+   r  TestPpccMax.test_dist  s    !!$ASu/6 " 8:=> ~~am4~~ae&7&78D"6BD"6B ~~a D"6Br-   c                .   \         P                  P                  R	^RRRR7      R,           p\        \        \         P
                  VR
R7       \        \         P
                  ! VRR7      R^R7       \        \         P
                  ! VRR7      R^R7       R# )r  r  r  r  r   r  brackr9   Nr  )r   rt  r  )r   rt   r  r   r{   )r   r  r(   rQ   rR   r  r   r   s   & r+   
test_brackTestPpccMax.test_brack  sx    !!$ASu/6 " 8:=>j%..!?KENN1F;0!	= 	ENN1G<0!	=r-   r_   N)
r`   ra   rb   rc   r  r  r  r  rd   re   rf   s   @r+   r  r    s!     P
P
C	= 	=r-   r  c                      a  ] tR tRt o ]P
                  P                  RRR.4      R 4       t]! RRR7      R	 4       t	R
 t
]P
                  P                  R 4       tR tR tR tRtV tR# )TestBoxcox_llfi  r  r  r  c                   \        W!4      p\        P                  P                  R ^
RR7      p^p\        P                  ! WRP                  WCR7      4      pVP                  ) R,          \        P                  ! \        P                  ! VP                  4       ^,          4      4      ,          p\        WbP                  WsR7      4       R# )r  1  r   r   r   r  r  N)rf  r   r   r(   
boxcox_llfr  r   r%   r&   r  rM   r   )rI   r  r  dtrY   lmbdallfllf_expecteds   &&&     r+   rN   TestBoxcox_llf.test_basic  s    RJJNN2ENBujjj&=>w|bffRVVAEEGQJ-?&@@ZZZ?@r-   Tz,array-likes only accepted for NumPy backend.r  c                    \         P                  P                  ^d^
RR7      p^p\         P                  ! W24      p\         P                  ! V\	        V4      4      p\        WERR7       R# r  r  r  rb  r4   N)r   r   r(   r  listr   rI   r  rY   r   r  llf2s   &&    r+   test_array_likeTestBoxcox_llf.test_array_like  sQ     JJNN%N@u(tAw/.r-   c                @   \         P                  P                  ^d^
RR7      p^p\         P                  ! W24      p\         P                  ! V\        P
                  ! W".4      P                  4      p\        VP                  WD.4      VP                  V4      RR7       R# r  )	r   r   r(   r  r%   vstackr  r   r  r  s   &&    r+   test_2d_inputTestBoxcox_llf.test_2d_input  st     JJNN%N@u(ryy!'8':':;

C:.

40@uMr-   c                v   R p\        V4      '       d   \        P                  ! \        VR7      M\        P
                  ! 4       pT;_uu_ 4        VP                  VP                  \        P                  ! ^VP                  . 4      4      4      4      '       g   Q h RRR4       R#   + '       g   i     R# ; i)r  r   N)
r   r   r   r   
contextlibnullcontextisnanr  r   r  )rI   r  rP  contexts   &&  r+   r  TestBoxcox_llf.test_empty  ss    @FNrll6<< 2'B"..0 	W88BJJu'7'72::b>'JKLLLL WWWs   AB''B8	c                    VP                  . RO4      p\        P                  ! RV4      p\        W1P                  R4      4       R# )     h@N)r        m@r  g     x@ig#:x1)r  r   r  r   rI   r  rJ   r  s   &&  r+   test_gh_6873TestBoxcox_llf.test_gh_6873  s7     zz67r4(ZZ(:;<r-   c                    VP                  . ROVP                  R7      p\        P                  ! RV4      p\	        W1P                  RVP                  R7      RR7       R# )  r  g:0yE>gƠ>r4   N)r  i  i  i  i  gl.)r  r  r   r  r   r  s   &&  r+   test_instability_gh20021'TestBoxcox_llf.test_instability_gh20021  sL    zz8

zKtT* 	ZZ(=RZZZP!	#r-   c                   VP                  ^d^.RR..4      p\        P                  ! ^V^ R7      pVP                  \        P                  ! ^VR,          4      \        P                  ! ^VR,          4      .4      p\	        W44       \        P                  ! ^V^R7      pVP                  \        P                  ! ^VR,          4      \        P                  ! ^VR,          4      .4      p\	        WV4       R# )	r  i,  i  r  N)r  r   )r  rt   )r   r  )rt   r  )r  r   r  r  r   )rI   r  rJ   
llf_axis_0llf_0
llf_axis_1llf_1s   &&     r+   	test_axisTestBoxcox_llf.test_axis  s    zzC:Sz23%%aA6
QT
+QT
+
  	
*%%aA6
QT
+QT
+
  	
*r-   r_   N)r`   ra   rb   rc   r   r   r   rN   r  r	  r  rr  r  r  r  r$  rd   re   rf   s   @r+   r  r    s     [[Wy)&<=A >A dKM/M/N [[M M=#+ +r-   r  iU>  iϵ i i_
 i7 i is iv  i+ i* i-q ir: i  i i i i_ i4  i iD  i i, iE i\  i i| ig  iM iӢ i[  i! i if6 iP i i iI i/ iF& i i i`# i i  i& i i iQ ia i
 iXZ iw  i_ im  iD i- iB  i:C i	 i' i  i % ix i  iz_ iQ i!# i= i  ix i
 iU7 i_ i1k iW9 i& ix iI i  iT i
 i i  irv im i͟  i iO iv i) iP i4 i= ii  iU i i ig> ix	 i` i| i0 i` iM	 i ' i[ i8 i  i3[ i  i) i`N i
	 i  i[ i:+ ip i	% i
  iV i[$ iX  i_ ib ieZ i i`?	 it i  i iE i iT iy i imv i< iZ i! i i(  ir
 i i iɪ i| i i> ii  i i* i# ix i i& i i# i1g iP i) iۆ i{% i? i۽ i| i i& i- iư) i+c i7  iQ i	o i& i|_ iy
 i6 i  i i i:J i#j i8 i  i%  i i%W i( i` i: ina  i  i i% i$ i4 iA( iq  i  i" iy iXZ i izr i i$ i; is igY i i7 ip	 ip iL6" iDH i i!3 i iI( i|C iш i im  iS i  i{8 i  i i iP i2L i( i@ is iu i2 i_ ij i)m igZ i_ iu	 i9 ik  i'; i ii iK i܃ i i= iv i= i ia  ir  i i i}W i*  i  iJ i i  i2  iUF iNj	 iV," ip iA i% i iv9 iI5 i  iT- i  i8~ i?e` i1 iV  i/X i= i4 i iح i^ i i/ ig i( i i|B; ip i ic iE' iE i$
 i i{ ik i: iI iA_ i iP ioP  i6_ i iܲ iO	 iA7 iL  iA i iV  iwm i iZh i|V iԇC iw i& i~ i i'  i!  i i
 i} ie i3 i i i_ i  i@) i  iV ii iF i? i&
 i"Y iw i  iH iS ird ib i@ iK iA; iz| id+ i	 i
8 i"  i	 iO i_ i$l i!9 ii iT iuL iQ'% iqp iC i]K@ i! i i  iIp| i i?	B i i io0 i i

 i  i5% i i i* i i iv0 i* i* iw i  ive i$ i i- iD~ i- i i,4 ib i  i iP0 i) im i	 iR iQ i
 i i c            	       <  a  ] tR tRt o R tR tR tR tR tR t	]
P                  P                  R. RO4      R	 4       tR
 t]
P                  P                  R. RO4      R 4       tR t]
P                  P                  R]P$                  ! . RO4      ]P$                  ! ]P&                  ^*^.4      .4      R 4       t]
P                  P                  R]P$                  ! . RO4      ]P$                  ! . RO4      ]P$                  ! . RO4      .4      R 4       tRtV tR# )
TestBoxcoxiJ  c                p   \        ^^2RR7      ^,           p\        P                  ! V^R7      p\        W!^,
          4       \        P                  ! VRR7      p\        V^^V,          ,
          4       \        P                  ! V^ R7      p\        V\        P
                  ! V4      4       \        P                  ! \        V4      ^ R7      p\        V\        P
                  ! V4      4       \        P                  ! \        P                  ! ^
4      ^4      p\        V\        P                  ! ^
4      4       R# )r   r  rs  r   Nr   )
r,   r   r   r
   r%   r&   r  r   r   r+  )rI   rY   xts   &  r+   test_fixed_lmbdaTestBoxcox.test_fixed_lmbdaL  s    abu=A\\!1%E"\\!2&A!G$\\!1%BFF1I& \\$q'+BFF1I& \\"''"+q)R"&r-   c                    R p\         P                  P                  ^
RRR7      pW!,          ^,           V) ,          p\         P                  ! V4      w  rE\	        VRV,          ^R7       R# )g      @iP  i  )r   r   r   r9   Nr   )r   r   r(   r   r   )rI   r   rY   x_invr*  r   s   &     r+   test_lmbda_NoneTestBoxcox.test_lmbda_None^  sS     JJNNrDNAQ5&)\\%(
FBJ:r-   c                   \         P                  P                  R 4      p\        ^^2VR7      ^,           p\        P
                  ! VRR7      w   r4\        VRR.4       \        P
                  ! VRR7      w   r4\        VRR.4       \        ^R	VR7      ^,           p\        P
                  ! VR
R7      w   r4\        VRR.4       \        P
                  ! VRR7      w   r4\        VRR.4       R# )r  rs  r  )alphagx@gv[%@ra  gxˉEk?g8Bhk @r  r.  g\?gx̡&@+?g7(U@gp:sY@N)r%   rd  r   r,   r   r   r
   )rI   rh  rY   r  rm   s   &    r+   
test_alphaTestBoxcox.test_alphah  s    ii##D)abs;a? at41#46G"HIat41#57H"IJ ac<rAau519k":;au51:z":;r-   c                J   \         P                  ! R^.4      p\        \        \        P
                  V4       \        \        \        P
                  \         P                  ! ^.4      4       \        \        \        P
                  \         P                  ! ^.^..4      4       R# rm  )r%   r9  rQ   rR   r   r   r   s   & r+   test_boxcox_bad_argTestBoxcox.test_boxcox_bad_argy  s_    HHb!Wj%,,2j%,,!>j%,,1#s0DEr-   c                ^    \        \        P                  ! . 4      P                  R8H  4       R# r   N)r   )r	   r   r   r  rS   s   &r+   r  TestBoxcox.test_empty  s    R &&$./r-   c                Z    \         P                  ! \        4      w  r\        VRRR7       R# )gsHjdr?r   r4   NgsHjdrr   r   _boxcox_datar
   )rI   r7  lams   &  r+   r  TestBoxcox.test_gh_6873  s     l+ 	YT2r-   boundsc                   a V3R  lp\         P                  ! \        RVR7      w  r4S^ ,          Tu;8  d   S^,          8  g   Q h Q hR# )c                 6   < \         P                  ! V SR R7      # bounded)rA  r  r   minimize_scalarfunrA  s   &r+   	optimizerBTestBoxcox.test_bounded_optimizer_within_bounds.<locals>.optimizer      ++C3<> >r-   Nr   rJ  )r   r   r>  )rI   rA  rJ  r  r   s   &f   r+   $test_bounded_optimizer_within_bounds/TestBoxcox.test_bounded_optimizer_within_bounds  s?    	> <<DINay5,6!9,,,,,r-   c                   aa \         P                  ! \        R R7      w  rVR,           V^,           3oRR/oVV3R lp\         P                  ! \        R VR7      w  rWB8w  g   Q h\        VS^ ,          4       R # )Nr)  r  xatolrb  c                 8   < \         P                  ! V SR SR7      # )rE  )rA  r  optionsrF  )rI  rA  rS  s   &r+   rJ  PTestBoxcox.test_bounded_optimizer_against_unbounded_optimizer.<locals>.optimizer  s!    ++C3<gO Or-   rM  r=  )rI   r  r   rJ  lmbda_boundedrA  rS  s   &    @@r+   2test_bounded_optimizer_against_unbounded_optimizer=TestBoxcox.test_bounded_optimizer_against_unbounded_optimizer  so    
 <<D9 #+uqy)E"	O
 !<<D2;=%%%vay1r-   rJ  c                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! \
        RVR7       RRR4       R#   + '       g   i     R# ; i)z`optimizer` must be a callabler   NrM  r   r   rR   r   r   r>  rI   rJ  s   &&r+   $test_bad_optimizer_type_raises_error/TestBoxcox.test_bad_optimizer_type_raises_error  s6     ]]:-MNNLLTYG ONNNs   AA	c                    R  pRp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! \
        RVR7       RRR4       R#   + '       g   i     R# ; i)c                     ^# )rt   r_   )rI  s   &r+   rJ  CTestBoxcox.test_bad_optimizer_value_raises_error.<locals>.optimizer  s    r-   z/return an object containing the optimal `lmbda`r   NrM  rY  )rI   rJ  rP  s   &  r+   %test_bad_optimizer_value_raises_error0TestBoxcox.test_bad_optimizer_value_raises_error  s?    
	 D]]:W55LLTYG 6555s   AA#	bad_xc                    Rp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! V4       RRR4       R#   + '       g   i     R# ; i)zHTest boxcox_normmax raises ValueError if x contains non-positive values.z#only positive, finite, real numbersr   N)r   r   rR   r   boxcox_normmax)rI   rb  rP  s   && r+   "test_negative_x_value_raises_error-TestBoxcox.test_negative_x_value_raises_error  s7    
 8]]:W55  ' 6555s   AA	rY   c                   \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! V4      w  r#\
        P                  ! \
        P                  ! V4      4      '       g   Q h RRR4       R#   + '       g   i     R# ; i)The optimal lambda isr   N)r   r   rY  r   r   r%   r   isfinite)rI   rY   xt_bclam_bcs   &&  r+   test_overflowTestBoxcox.test_overflow  sS     \\+-DEE!LLOME66"++e,---- FEEEs   AA;;B	r_   N)r   rt   )r   r{   )r   rW  )strrt   r{   r  )rt   ig@
     L@     x@     4@     @@     d@rv  g     @g     <@g     \@g     @rr  rs  rt  ru  rv  gFq$jgS	#jgGewa$jgVBi$jg[(>2$j)r`   ra   rb   rc   r+  r/  r4  r7  r  r  r   r   r   rN  rV  r[  r`  r%   r9  r   re  rl  rd   re   rf   s   @r+   r'  r'  J  s     '$;<"F03 [[X'FG- H-2* [[[*>?H @H

H [[bhh01288RVVRO3LM
(
( [[S
 : 	; 	9:
MN# ..r-   r'  c                     a  ] tR tRt o R tR tR tR t]P                  P                  R. RO4      ]P                  P                  R. RO4      R	 4       4       t]P                  P                  R
 4       tR t]P                  P                  R. RO. RO34      R 4       tR t]P                  P                  R]P$                  ! . RO]P&                  R7      ]P$                  ! . RO]P&                  R7      ]P$                  ! . RO]P(                  R7      ]P$                  ! . RO]P(                  R7      .4      ]P                  P                  R. RO4      ]P                  P                  RR.4      R 4       4       4       t]P                  P                  R. RO. RO.4      ]P                  P                  RR.4      R 4       4       tRtV tR# )TestBoxcoxNormmaxi  c                :    \        ^^2RR7      ^,           V n        R# r   r  rs  Nr  rS   s   &r+   r  TestBoxcoxNormmax.setup_method      "12EBQFr-   c                b    \         P                  ! V P                  4      p\        VR RR7       R# )m?r3   r4   Nr   rd  rY   r
   rI   r   s   & r+   test_pearsonrTestBoxcoxNormmax.test_pearsonr  s"    %%dff-t4r-   c                    \         P                  ! V P                  R R7      p\        VRRR7       \         P                  ! V P                  4      w  r#\        W14       R# )mlerA  	.!?r3   r4   N)r   rd  rY   r
   r   )rI   r   r  maxlog_boxcoxs   &   r+   test_mleTestBoxcoxNormmax.test_mle  sC    %%dffU;t4 !<</.r-   c                j    \         P                  ! V P                  R R7      p\        VRR.RR7       R# )r   rA  r  r  r3   r4   Nr  )rI   
maxlog_alls   & r+   test_allTestBoxcoxNormmax.test_all  s)    ))$&&?

Xx$8tDr-   r  r  rA  c                   a V3R  lp\         P                  ! V P                  VVR7      p\        P                  ! S^ ,          V8  4      '       g   Q h\        P                  ! VS^,          8  4      '       g   Q hR# )c                 6   < \         P                  ! V SR R7      # rD  rF  rH  s   &r+   rJ  ITestBoxcoxNormmax.test_bounded_optimizer_within_bounds.<locals>.optimizer  rL  r-   )r  rJ  N)r   rd  rY   r%   r   )rI   r  rA  rJ  r   s   &&f  r+   rN  6TestBoxcoxNormmax.test_bounded_optimizer_within_bounds  sa    	> %%dffV09;vvfQi&())))vvfvay())))r-   c                R  aa \         P                  ! V P                  4      p\        P                  ! V^4      p\        P
                  ! VR,
          VR,           R4      o ! R R4      oVV3R lp\         P                  ! V P                  VR7      pWA8w  g   Q h\        WAR4       R# )	r   r  i  c                       ] tR tRtRtR# )?TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.MyResulti  r_   N)r`   ra   rb   rc   rd   r_   r-   r+   MyResultr    s    r-   r  c                    < . pS F  pVP                  V ! V4      4       K  	  S! 4       pS\        P                  ! V4      ,          Vn        V# r$   )appendr%   argminrY   )rI  objsr   r[   r  lmbda_ranges   &   r+   rJ  @TestBoxcoxNormmax.test_user_defined_optimizer.<locals>.optimizer	  sB    D$CJ' %*C		$0CEJr-   rJ  r   N)r   rd  rY   r%   rg  r  r
   )rI   r   lmbda_roundedrJ  lmbda2r  r  s   &    @@r+   test_user_defined_optimizer-TestBoxcoxNormmax.test_user_defined_optimizer  s     $$TVV,*kk-"4mD6H$O	 		 %%dff	Bt,r-   c                4   \         P                  p\        P                  ! V P                  R VR7       \
        P                  ! \        RR7      ;_uu_ 4        \        P                  ! V P                  RVR7       R R R 4       R #   + '       g   i     R # ; i)N)r  rJ  z,`brack` must be None if `optimizer` is givenr   )g       r  )r   rG  r   rd  rY   r   r   rR   rZ  s   & r+   2test_user_defined_optimizer_and_brack_raises_errorDTestBoxcoxNormmax.test_user_defined_optimizer_and_brack_raises_error	  sv    ,,	 	TVV49E
 ]]: .D E E   {+46E E E Es   $BB	rY   c                   R p\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! VRR7      pRRR4       \
        P                  ! \        P                  ! VX4      4      P                  4       '       g   Q h\
        P                  ! \
        P                  4      P                  R,          pV^ 8  d   \
        P                  ! V4      M\
        P                  ! V4      p\        P                  ! WS4      p\        Wd\
        P                  ! V4      ,          4       R#   + '       g   i     L; i)zThe optimal lambda is...r   r  rA  Nr  )r   r   rY  r   rd  r%   ri  r   r   r   finfor  maxrR  r
   sign)rI   rY   rP  r   ymaxx_treme	y_extremes   &&     r+   rl  TestBoxcoxNormmax.test_overflow	  s     -\\+W55((59E 6{{7>>!U3488::::xx

#''%/$qy"&&)bffQiNN72		"''%.#89 65s   D//D?	c                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! V P
                  RR7       RRR4       R#   + '       g   i     R# ; i)z `ymax` must be strictly positiver   )r  Nr   )r   r   rR   r   rd  rY   rS   s   &r+   test_negative_ymax$TestBoxcoxNormmax.test_negative_ymax,	  s8    ]]:-OPP  b1 QPPPs   #AA#	r  r  Nc           
        \         P                  ! \        R R7      ;_uu_ 4        Ve   RV/M/ p\        P                  ! V3RV/VB p\
        P                  ! V4      \
        P                  ! V4      .p\        \        \        P                  ! We4      4      4      pVf2   \
        P                  ! VP                  4      P                  R,          p\        W'RR7       RRR4       R#   + '       g   i     R# ; i)rh  r   Nr  r  r  r   r4   )r   r   rY  r   rd  r%   rR  r  absr   r  r  r
   )rI   rY   r  r  kwarglmbr  ymax_ress   &&&&    r+   'test_user_defined_ymax_input_float64_329TestBoxcoxNormmax.test_user_defined_ymax_input_float64_320	  s    & \\+-DEE&*&6VTNBE&&qAA5ACvvay"&&),G3u||G9:;H|xx(,,u4D6 FEEEs   B:C))C:	c                0   \         P                  ! V\         P                  R 7      p\         P                  ! V\         P                  R 7      p\        P
                  ! \        RR7      ;_uu_ 4        \        P                  ! W2R7       RRR4       \        P                  ! WBR7       \        P                  ! V\         P                  VR7      p\        P                  ! V\         P                  VR7      p\        WVRR7       R#   + '       g   i     L; i)r  rh  r   rA  N)r  r  r  r4   )r%   r  r  r  r   r   rY  r   rd  rn  r
   )rI   rY   r  x_32x_64lmb_32lmb_64s   &&&    r+   test_user_defined_ymax_inf,TestBoxcoxNormmax.test_user_defined_ymax_infM	  s     zz!2::.zz!2::. \\+-DEE  5 FT1 %%dG%%dGT2 FEs   1DD	r  )r  pearsonrr   rn  rw  )gx	 ?gN~jh ?g{ ?gXkI ?g}9" ?)g	i@g     `h@gfffffh@r  gi@)gKH9KH9r  r  r  r  )g    _Bgꌠ9Y>)FN)r`   ra   rb   rc   r  r  r  r  r   r   r   rN  slowr  r  rl  r  r%   r9  r  r  r  r  rd   re   rf   s   @r+   rz  rz    s    G5/E [[X'AB[[X'FG	* H C	* [[- -,6  [[6JLM	:M	:2 [[S
9zz	# 	Mzz	# 	4zz	# 	;zz	##  [[V%78[[Xw/
7 0 9"
7 [[S+2	#  [[Xw/3 03r-   rz  c                      a  ] tR tRt o R tR t]P                  P                  ]	'       * RR7      R 4       t
R tR tR	tV tR
# )TestBoxcoxNormplotid	  c                :    \        ^RRR7      ^,           V n        R# r  r  rS   s   &r+   r  TestBoxcoxNormplot.setup_methode	  r  r-   c           	         ^p\         P                  ! V P                  R^
VR7      w  r#. ROp\        V\        P
                  ! R^
VR7      4       \        W44       R# )r   r  r  Nr  )g!X4}?g΁i?g-	15?gW$?g\ܱ{?)r   boxcox_normplotrY   r
   r%   r  )rI   r  lmbdasr  r  s   &    r+   rN   TestBoxcoxNormplot.test_basich	  sI    ,,TVVS"B%C ;<,r-   r  rL  c                X   \         P                  ! 4       pVP                  ^o4      p\        P                  ! V P
                  R^\         R7       VP                  V4       VP                  ^o4      p\        P                  ! V P
                  R^VR7       \         P                  ! 4        R# r  )r  r  r  r   r  rY   r  r  r  s   &  r+   r  "TestBoxcoxNormplot.test_plot_kwargp	  sq     jjl__S!dffc2C8B __S!dffc2B7		r-   c                    \        \        \        P                  V P                  ^^ 4       \        \        \        P                  R^.^ ^4       R# rm  )rQ   rR   r   r  rY   rS   s   &r+   r  &TestBoxcoxNormplot.test_invalid_inputs}	  s6    j%"7"7AFj%"7"7"a!QGr-   c                b    \        \        P                  ! . ^ ^4      P                  ^ 8H  4       R# r[  )r	   r   r  r   rS   s   &r+   r  TestBoxcoxNormplot.test_empty	  s$    %%b!Q/449:r-   r  N)r`   ra   rb   rc   r  rN   r   r   rX  r  r  r  r  rd   re   rf   s   @r+   r  r  d	  sL     J- [[O+OD
 E
H; ;r-   r  c                   2   a  ] tR tRt o R tR tR tRtV tR# )TestYeojohnson_llfi	  c                    \         P                  P                  ^d^ RR7      p^p\         P                  ! W!4      p\         P                  ! V\	        V4      4      p\        W4RR7       R# r  )r   r   r(   yeojohnson_llfr  r
   rI   rY   r   r  r  s   &    r+   r	  "TestYeojohnson_llf.test_array_like	  sO    JJNNN?""5,##E473.r-   c                   \         P                  P                  ^d^
RR7      p^p\         P                  ! W!4      p\         P                  ! V\        P
                  ! W.4      P                  4      p\        W3.VRR7       R# r  )r   r   r(   r  r%   r  r  r
   r  s   &    r+   r   TestYeojohnson_llf.test_2d_input	  s`    JJNN%N@""5,##E299aV+<+>+>?
Du5r-   c                n    \        \        P                  ! \        P                  ! ^. 4      4      4       R# rs   )r	   r%   r  r   r  rS   s   &r+   r  TestYeojohnson_llf.test_empty	  s     --a456r-   r_   N)	r`   ra   rb   rc   r	  r  r  rd   re   rf   s   @r+   r  r  	  s     /67 7r-   r  c                     a  ] tR tRt o R t]P                  P                  R. RO4      R 4       tR t	R t
]P                  P                  R]P                  ]P                  .4      R 4       t]P                  P                  R]P                  ]P                   ]P"                  ]P$                  .4      R	 4       tR
 t]P                  P                  R]P*                  ! R]! R4      R.4      ]P*                  ! R]! R4      R.4      ]P*                  ! R]! R4      ) R.4      ]P*                  ! R]! R4      ]! R4      ]! R4      ) R.4      .4      R 4       t]P                  P                  R]P*                  ! . RO4      ]P*                  ! . RO4      ]P*                  ! . RO4      .4      R 4       t]P                  P                  R]P*                  ! . RO4      ]P*                  ! . RO4      .4      ]P                  P                  R. RO4      ]P                  P                  R^R.4      R 4       4       4       t]P                  P                  R]P*                  ! . R O4      ]P*                  ! . R!O4      ]P*                  ! . R"O4      .4      ]P                  P                  R^R.4      ]P                  P                  RRR#.4      R 4       4       4       tRtV tR# )$TestYeojohnsoni	  c                   \         P                  P                  R 4      p\        ^^2VR7      ^,           p\         P                  ! V^ 8  4      '       g   Q h\
        P                  ! V^R7      p\        W24       \
        P                  ! VRR7      p\        V^^V^,           ,          ,
          4       \
        P                  ! V^ R7      p\        V\         P                  ! V^,           4      4       \
        P                  ! V^R7      p\        W24       \        ^^2VR7      ^,
          p\         P                  ! V^ 8  4      '       g   Q h\
        P                  ! V^R7      p\        V\         P                  ! V) ^,           4      ) 4       \
        P                  ! V^R7      p\        W24       \
        P                  ! V^R7      p\        V^V) ^,           ,          ^,
          4       \        ^^2VR7      ^,
          p\         P                  ! V^ 8  4      '       d   Q h\         P                  ! V^ 8  4      '       d   Q hV^ 8  p\
        P                  ! V^R7      p\        W4,          W$,          4       \
        P                  ! VRR7      p\        W4,          ^^W$,          ^,           ,          ,
          4       \
        P                  ! V^ R7      p\        W4,          \         P                  ! W$,          ^,           4      4       \
        P                  ! V^R7      p\        W4,          W$,          4       V( p\
        P                  ! V^R7      p\        W5,          \         P                  ! W%,          ) ^,           4      ) 4       \
        P                  ! V^R7      p\        W5,          W%,          4       \
        P                  ! V^R7      p\        W5,          ^W%,          ) ^,           ,          ^,
          4       R# )r  rs  r)  Nr   )	r%   rd  r   r,   r   r   
yeojohnsonr
   r&   )rI   rh  rY   r*  posnegs   &     r+   r+  TestYeojohnson.test_fixed_lmbda	  s   ii##E* abs;a?vva!e}}}aq)ar*AQUO,aq)BFF1q5M*aq) abs;a?vva!e}}}aq)RVVQBF^O,aq)aq)A!aL1,- abs;a?66!a%==  66!q&>>!!1faq)(ar*Q!&1*%5!56aq)
!34aq)(daq)"&&!&1"5!56aq)(aq)qvgk!2Q!67r-   r   c                   R  pRp\         P                  P                  R4      pVP                  ^ ^VR7      pV! WQ4      p\        P
                  ! V4      w  rx\        WRR7       \        ^ \         P                  P                  WW,
          4      V,          ^R7       \        ^ VP                  4       ^R7       \        ^VP                  4       ^R7       R# )c                    \         P                  ! V P                  V P                  R 7      pV ^ 8  p\	        V4      \         P
                  ! R4      8  d'   \         P                  ! W,          4      ^,
          W#&   M;\         P                  ! W,          V,          ^,           ^V,          4      ^,
          W#&   \	        V^,
          4      \         P
                  ! R4      8  dO   ^\         P                  ! ^V,
          ) W( ,          ,          ^,           ^^V,
          ,          4      ,
          W#( &   V# ^\         P                  ! W( ,          ) 4      ,
          W#( &   V# )r  rt  )r%   r+  r  r  r  spacingexppower)rY   r   r.  r  s   &&  r+   _inverse_transform:TestYeojohnson.test_lmbda_None.<locals>._inverse_transform	  s    HHQWWAGG4Eq&C 5zBJJrN*VVAF^a/
XXafunq&8!e)DqH
 519~

2."((QY<!D'+AA+E+,E	?#< <d
 L  "&&!D'"22dLr-   i N  r  )r   r   r   r  r0  r9   N)r%   rd  r   rk  r   r  r
   r   linalgr   rK   rM   )	rI   r   r  	n_samplesrh  rY   r.  r*  r   s	   &&       r+   r/  TestYeojohnson.test_lmbda_None	  s    
	& 	ii##G,JJ1AYJ8"1,%%e,
D1Aryy~~af5	A1MArwwy!4Arvvx3r-   c                ^    \        \        P                  ! . 4      P                  R8H  4       R# r:  )r	   r   r  r  rS   s   &r+   r  TestYeojohnson.test_empty	  s       $**d23r-   c                    \         P                  P                  ^d^ RR7      p\         P                  ! V4      w  r#\         P                  ! \	        V4      4      w  rC\        W$RR7       R# r  )r   r   r(   r  r  r
   )rI   rY   xt1r  xt2s   &    r+   r	  TestYeojohnson.test_array_like	  sL    JJNNN?!!!$!!$q'*u-r-   r  c                    \         P                  ! ^VR7      pRp\        P                  ! \        VR7      ;_uu_ 4        \
        P                  ! V4       RRR4       R#   + '       g   i     R# ; i)r1   r  z>Yeo-Johnson transformation is not defined for complex numbers.r   N)r%   rX   r   r   rR   r   r  )rI   r  rY   err_msgs   &&  r+   test_input_dtype_complex'TestYeojohnson.test_input_dtype_complex	  sG    IIau%]]:W55Q 6555s   A  A1	c                   \         P                  ! ^VR7      p\         P                  ! ^\         P                  R7      p\        P                  ! V4      w  rE\        P                  ! V4      w  rg\        WFRR7       \        WWRR7       R# )r@   r  r  r4   N)r%   rX   r  r   r  r
   )rI   r  x_intx_floatxt_int	lmbda_intxt_floatlmbda_floats   &&      r+   test_input_dtype_integer'TestYeojohnson.test_input_dtype_integer
  s`    		!5)))ARZZ0!,,U3 % 0 0 9t4	T:r-   c                    \         P                  ! . RO4      p\        P                  ! V4      w  r#\        P                  ! V^,           4      w  rE\        W$RR7       \        W5RR7       R# )(HAr3   r4   N)r  G."A    5fAg(@Ag   VAgHzuoAg)\¿+Ar   g(\^@Ag    RAr  g    #Ag\µ=Ag
ףNAr  gp=oHAg    6bAg{csAg)\tmhAg/BAg   '|Ag(\wipAg    ݃hAgQlBAgףp @A)r%   r9  r   r  r   r
   )rI   rY   xt_yeolam_yeoxt_boxlam_boxs   &     r+   test_input_high_variance'TestYeojohnson.test_input_high_variance
  sQ    HH " #  **1-,,q1u-T2t4r-   rY   rt  r   r  rn  c                    \         P                  ! \        R R7      ;_uu_ 4        \        P                  ! V4      w  r#RRR4       R#   + '       g   i     R# ; i)z Yeo-Johnson input must be finiter   N)r   r   rR   r   r  )rI   rY   r  r  s   &&  r+   test_nonfinite_input#TestYeojohnson.test_nonfinite_input
  s8     ]]:-OPP#..q1OF QPPPr1  c           
        R  p\         P                  ! RR7      ;_uu_ 4        \        P                  ! V4      w  r4\        P                  ! V^,           \        W$R7      R7      w  rV\         P                  ! \         P                  ! V4      4      '       g   Q h\         P                  ! \         P                  ! V4      4      '       g   Q h\        WFRR7       \        W5RR7       RRR4       R#   + '       g   i     R# ; i)	c                 p    \         P                  ! W) VR R7      p\         P                  ! 4       pW#n        V# )g`sbO>)xtol)r   	fminboundOptimizeResultrY   )rI  r  outresults   &&  r+   rJ  /TestYeojohnson.test_overflow.<locals>.optimizer/
  s0    $$S(G(KC,,.FHMr-   raiser   )r  r  r3   r4   r  N)	r%   r   r   r  r   r   ri  rL   r
   )rI   rY   rJ  r  r  r  r  s   &&     r+   rl  TestYeojohnson.test_overflow#
  s    	 [[W%%#..q1OF#llA!DFOF;;rvvf~....;;rvvf~....G48F6 &%%%s   B%C/	C//D 	r   r  c                   \         P                  ! R R7      ;_uu_ 4        \        P                  ! W1,          V,          4      w  rE\         P                  ! \         P
                  ! W1,          4      \         P
                  ! V4      8H  4      '       g   Q h\         P                  ! V4      '       g   Q h\         P                  ! \         P                  ! V4      4      '       g   Q h RRR4       R#   + '       g   i     R# ; i)r  r   N)r%   r   r   r  r   r  ri  rL   )rI   rY   r   r  r  r  s   &&&&  r+   #test_overflow_underflow_signed_data2TestYeojohnson.test_overflow_underflow_signed_data>
  s     [[W%%#..tx%/?@OF66"''$(+rwwv>????;;w'''';;rvvf~....	 &%%%s   BC011C0#C00D	r  Nc                   \         P                  ! R R7      ;_uu_ 4        W!,          pVP                  \         P                  4      p\        P
                  ! WCR7      p\        P                  ! WFR7      p\        P
                  ! WSR7      p\        P                  ! WXR7      p	\         P                  ! \         P                  ! V4      \         P                  ! V4      8H  4      '       g   Q h\         P                  ! V4      '       g   Q h\         P                  ! \         P                  ! V4      4      '       g   Q hWh8X  g   Q h\         P                  ! Wy8H  4      '       g   Q h RRR4       R#   + '       g   i     R# ; i)r  r   r  r)  N)r%   r   astyper  r   yeojohnson_normmaxr  r   r  ri  rL   )
rI   rY   r  r  r  r  lam_yeo_int
xt_yeo_intlam_yeo_floatxt_yeo_floats
   &&&&      r+   test_integer_signed_data'TestYeojohnson.test_integer_signed_dataM
  s     [[W%%HEll2::.G225FK))%CJ!44WJM ++GIL66"''%.BGGJ,??@@@@;;{++++;;rvvj12222///66*45555 &%%%s   C(E0	1E0;*E00F	r_   )r   r  r  r{         rq  rw  rx  )rt   rb  g3#I9gu?j/ gnFgZbtir   )r   rt   r{   r|   )r   r   r{   r8  )r   r   r   r  )r`   ra   rb   rc   r+  r   r   r   r/  r  r	  r%   	complex64
complex128r  int8uint8int16int32r  r  r9  r  r  rl  r  r  rd   re   rf   s   @r+   r  r  	  s    -8^ [[Wn5"4 6"4H4. [[Wr||R]]&CD  E  [[Wrww"((BHH&MN; O;5 [[S
#uU|S)*
#uU|S)*
#e}c*+
$eeElU5\M3GH	# 22 [[S
 : 	; 	9:
MN# 77$ [[S
 : 	;
9:# 
 [[W&LM[[VaW-/ . N/ [[S

 
# 
 [[VaW-[[WtWo66 7 .6r-   r  c                   2   a  ] tR tRt o R tR tR tRtV tR# )TestYeojohnsonNormmaxic
  c                :    \        ^^2RR7      ^,           V n        R# r|  r  rS   s   &r+   r  "TestYeojohnsonNormmax.setup_methodd
  r~  r-   c                b    \         P                  ! V P                  4      p\        VR RR7       R# )g?r3   r4   N)r   r  rY   r
   r  s   & r+   r  TestYeojohnsonNormmax.test_mleg
  s"    ))$&&1t4r-   c                |    . ROp\         P                  ! V4      p\        P                  ! VRRR7      '       g   Q hR# )ffffff@gzG?r.  r0  N)r.  g rt  r  r  g333333@r  gffffff@g?g@r  r  g"@r  g      )r   r  r%   allclose)rI   rY   r   s   &  r+   test_darwin_example)TestYeojohnsonNormmax.test_darwin_examplek
  s2    ((+{{5%d3333r-   r  N)	r`   ra   rb   rc   r  r  r0  rd   re   rf   s   @r+   r(  r(  c
  s     G54 4r-   r(  c                     a  ] tR tRt o ]P
                  P                  R]P                  R3]P                  R3]P                  R3.4      R 4       tR t]P
                  P                  R]P                  ]P                  3]P                  ]P                  3]P                  ]P                   3.4      R	 4       t]P
                  P                  R
]P                  ]P                  ]P                  .4      R 4       t]P
                  P                  R]P                  R3]P                  R3]P                  R3.4      R 4       t]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       t]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       t]P
                  P                  R]P                  R]P,                  ^ R^R/3]P                  R]P,                  ^ R^R/3]P                  R]P,                  ^ R^R/3.4      R 4       tR tR tR tR tR tR tRtV tR# ) TestCircFuncsit
  test_func,expectedA:wv?gpz?E|2@c                r    VP                  . RO4      p\        V! VRR7      VP                  V4      4       R# )     0v@h  r  N)r8  r  r       pv@r6        u@r  r   rI   	test_funcro  r  rY   s   &&&& r+   test_circfuncsTestCircFuncs.test_circfuncs{
  s-    
 JJ67	!#.

80DEr-   c                   VP                  . R
OVP                  R7      pVP                  V4      p\        P                  ! VRR7      p\        WCRR7       VP                  W!P                  ,          ^,          ^ R7      pVR,          p\        P                  ! VRR7      p\        WeRR7       VP                  V^ R7      p\        P                  ! VRR7      p\        WRR7       R	# )r  r  r9  r:  r   r4   r  r  r  N)r  rJ  r  r     g     4@g3333333@)r  r  rK   r   circmeanr   rL   picircvarrM   circstd)	rI   r  rY   M1M2V1V2S1S2s	   &&       r+   test_circfuncs_small"TestCircFuncs.test_circfuncs_small
  s     JJ7rzzJJWWQZ^^AC(T*VVAeeGCKAV. "W]]13'T*VVA!V$]]13'T*r-   ztest_func, numpy_funcc                    \         P                  ! R .^
,          R.^d,          ,           4      pV! VP                  V4      4      pVP                  V! V4      4      p\        WVRR7       R# )g˛ɦv9?guv9?g:0yU>r0  N)r%   r  r   )rI   r?  
numpy_funcr  rY   circstatrk  s   &&&&   r+   test_circfuncs_close"TestCircFuncs.test_circfuncs_close
  sW     JJ+,r15H4IC4OOPRZZ]+JqM*t4r-   circfuncc                H   VP                  . RO. RO. RO.4      pV! VRR7      pV! VP                  VR	4      RR7      p\        WAP                  V4      4       V! VR^R7      p\        VP                  ^ ,          4       Uu. uF  qb! W6R3,          RR7      NK  	  pp\        WAP                  V4      4       V! VR^ R7      p\        VP                  ^,          4       Uu. uF  qb! VRV3,          RR7      NK  	  pp\        WAP                  V4      4       R# u upi u upi )
c  r9  r:  r  r  r  N)rX  r   r{   r  r  ^  )_  r?   r~   `  r=   ]  )ie  r=   r@   if  r~   g     @v@r   )r  reshaper   r|  r  r  )rI   r  rV  rY   r[   r   r\   s   &&&    r+   test_circmean_axis TestCircFuncs.test_circmean_axis
  s     JJ1013 4 qs#rzz!U+#6ZZ_-qs+49!''!*4EF4EqxQ$c*4EFXXc]+qs+49!''!*4EF4Eqx!Q$c*4EFXXc]+ G Gs   D!D @pz?c                r    VP                  . RO4      p\        V! VRR7      VP                  V4      4       R# )rX  r9  r:  N)rX  r   r{   r  r  r<  r=  r>  s   &&&& r+   test_circfuncs_array_like'TestCircFuncs.test_circfuncs_array_like
  s-    
 JJ12	!#.

80DEr-   r?  c                :   VP                   pVP                  . VR 7      p\        V4      '       d;   \        P                  ! \
        \        R7      ;_uu_ 4        V! V4      pRRR4       Mc\        P                  P                  4       ;_uu_ 4       pVP                  \        R4       VP                  \        R4       V! V4      pRRR4       \        XVP                  VP                  VR 7      4       R#   + '       g   i     L9; i  + '       g   i     LL; i)r  r   NzMean of empty slicer  )r  r  r   r   r   r   r   r%   r  r   r  r   r   r   )rI   r?  r  r  rY   r[   r  s   &&&    r+   r  TestCircFuncs.test_empty
  s     

JJrJ'B<<08MNNl ON --//3

>+@A

>+FGl	 0
 	RZZeZ<= ON 0/s   	C75D
7D	
D	c           	         VP                  R ^^R^
R\        P                  .4      p\        V! VRR7      VP                  VP                  4      4       R# )rX  r  rZ  r9  r:  N)r  r%   r   r   )rI   r?  r  rY   s   &&& r+   test_nan_propagate TestCircFuncs.test_nan_propagate
  sC     JJQ3C89	!#.

2660BCr-   Ng5:v@gfPb?gee?gOQOv?g#q@gauy+@c                b   VP                  . ROR^^R^	R\        P                  .^\        P                  \        P                  \        P                  \        P                  \        P                  \        P                  ..4      pVP                  4        F  pV! VRVR7      pVf'   \	        WcP                  VP                  4      4       K8  \        V^ ,          VP                  W%,          4      4       \	        VR,          VP                  VR,          VP                  4      4       K  	  R# )	rX  r[  r\  r]  r9  rY  N:rt   NN)rX  r   r{   r  r  rZ  rt   )r  r%   r   keysr   r   	full_like)rI   r?  ro  r  rY   r  r  s   &&&&   r+   test_nan_propagate_array&TestCircFuncs.test_nan_propagate_array
  s     JJ4aCC8BFFBFFBFFBFFBFFBFFKM N MMODACd3C|ZZ%78A

8>(BCBc"grvv)FG $r-   c                    VP                  R 4      R,          pTp\        P                  ! V4      p\        WCRR7       R# )rt  r   r4   Nr_   )r  r   rE  r   )rI   r  rY   rI  rJ  s   &&   r+   test_circmean_scalar"TestCircFuncs.test_circmean_scalar
  s1    JJrN2^^AT*r-   c                   \         P                  ! VP                  ^ ^R4      VP                  VP                  ) 4      p\	        W!P                  VP                  4      4       \	        V) VP                  VP                  4      4       R# )r   r  N)r   rE  rX   rF  r    r  )rI   r  r<  s   && r+   test_circmean_range!TestCircFuncs.test_circmean_range
  s[     NN299Q3/?q**RUU+,r2::bee,-r-   c                j   VP                  ^^
.VP                  R7      p\        \        P                  ! V^R7      VP                  R4      4       \        \        P
                  ! V^R7      VP                  R4      4       \        \        P                  ! V^R7      VP                  R4      4       R# )   r  r:  g     @e@gLwqA?g)#p_4@N)r  r$  r   r   rE  rG  rH  )rI   r  rY   s   && r+   test_circfuncs_uint8"TestCircFuncs.test_circfuncs_uint8
  sx     JJRyJ1qs3RZZ5FGac2BJJ4OPac2BJJ{4KLr-   c                    \         P                  ! VP                  ^ .4      4      p\        P                  ! RV4      R8X  g   Q hR# )r   rt  N)r   rH  r  mathcopysign)rI   r  r7  s   && r+   test_circstd_zeroTestCircFuncs.test_circstd_zero
  s3    MM"**aS/*}}S!$+++r-   c                   VP                  R R^24      pVP                  VP                  V4      V8H  4      '       d*   VP                  VP                  V4      R8H  4      '       g   Q hV^VP                  ,          ,          ^VP                  ,          ,          V8g  pVP                  V4      '       g   Q hW#,          p\        P                  ! VR,          ^R7      pVP                  WB8H  4      '       g   Q hR# )r7   gA:9>rt  Nr  )r  N)r  r   sincosrF  anyr   rE  )rI   r  rY   r<  r7  s   &&   r+   !test_circmean_accuracy_tiny_input/TestCircFuncs.test_circmean_accuracy_tiny_input  s     KKdB'vvbffQi1n%%"&&c1A*B*BBB!bee)_BEE	*q0vvayyyDNN1W:A.vvaf~~~r-   c                X   VP                  R VP                  R7      p\        P                  ! VP	                  V4      VP                  V4      4      pVP                  W1P                  R7      p\        P                  ! W!P                  VP                  ) R7      p\        WTRRR7       R# )g ؅W4vCr  )r  r  ro  r   )r5   r1  N)
r  r  r{  atan2r  r  r   rE  rF  r   )rI   r  rY   r7  ro  actuals   &&    r+   !test_circmean_accuracy_huge_input/TestCircFuncs.test_circmean_accuracy_huge_input  st     JJt2::J.JJrvvay"&&),::azz:2BEE6:u3?r-   r_   ) r`   ra   rb   rc   r   r   r   r   rE  rG  rH  r@  rO  r%   rK   rL   rM   rT  r`  rd  r  ri  r   rn  rq  rt  rx  r}  r  r  rd   re   rf   s   @r+   r3  r3  t
  s     [[1$~~{;$}}.BC$}}k:<=F	=F+( [[4$~~rww7$}}bff5$}}bff5785	85 [[Z%..*/--*/--*9 :,:,  [[1$~~{;$}}.BC$}}k:<=F	=F [[[5>>5==+0==+: ;>;> [[[5>>5==+0==+: ;D;D [[1$~~ $bffaq*MO$}} $bff !#7 !#7 9:  %}} $bffaQ
KMNO
HO
H+.M,
*@ @r-   r3  c                   .  a  ] tR tRt o ]P
                  P                  R]P                  RR^ ]	P                  ! . RO4      ^]	P                  ! RR.4      /3]P                  RR^ ]	P                  ! . RO4      ]	P                  ^,          ^,          ,          ^]	P                  ! RR	.4      /3]P                  RR
^ ]	P                  ! . R O4      ^]	P                  ! RR.4      /3.4      R 4       t]P
                  P                  R]P                  R3]P                  R3]P                  R3.4      R 4       t]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       t]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       t]P
                  P                  RR^^R^
R]	P$                  .]	P                  ! R^^R^
R]	P$                  .R^^R]	P$                  ^	R..4      .4      ]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       4       t]P
                  P                  RR^^R^
R]	P$                  .]	P                  ! R^^R^
R]	P$                  .R^^R]	P$                  ^	R..4      .4      ]P
                  P                  R]P                  ]P                  ]P                  .4      R 4       4       tRtV tR# )!TestCircFuncsNanPolicyi"  r4  NgrQvv@g5@v?g2>Bjv@gG2?rb  gX҄?gJPm$@gun2@g+fa @c                   \         P                  ! R ^^R^
R\         P                  .R^^R^	R\         P                  .\         P                  \         P                  \         P                  \         P                  \         P                  \         P                  \         P                  ..4      pVP                  4        F  pVf#   V! VRRVR	7      p\	        WRV,          R
R7       K)  \
        P                  ! \        \        R7      ;_uu_ 4        V! VRRVR	7      p\	        VRR W$,          R
R7       \        \         P                  ! VR,          4      4       RRR4       K  	  R#   + '       g   i     K  ; i)rX  r  rZ  r[  r\  r]  Nr9  omit)r  
nan_policyr  r  r4   r   r   )r%   r9  r   rl  r
   r   r   r   r   r	   r  )rI   r?  ro  rY   r  r  s   &&&   r+   test_nan_omit_array*TestCircFuncsNanPolicy.test_nan_omit_array(  s	   ( HHsAq#r37Aq#q#rvv6vvrvvrvvrvvrvvrvvrvvNP Q MMOD|TJd^$?\\"4<MNN#ACFNC#CHhn4HBHHSW-. ON $
 ONNs   A
E&&E7r5  r6  c                b    R ^^R^
R\         P                  .p\        V! VRRR7      VRR7       R# )	rX  r  rZ  r9  r  r  r  r  r4   N)r%   r   r
   )rI   r?  ro  rY   s   &&& r+   test_nan_omit$TestCircFuncsNanPolicy.test_nan_omitI  s4    
 !QRbff-	!#&A t	-r-   r?  c                x   \         P                  \         P                  \         P                  \         P                  \         P                  .p\        P                  ! \        \
        R 7      ;_uu_ 4        \        \         P                  ! V! VRR7      4      4       RRR4       R#   + '       g   i     R# ; i)r   r  r  N)r%   r   r   r   r   r   r	   r  rI   r?  rY   s   && r+   test_nan_omit_all(TestCircFuncsNanPolicy.test_nan_omit_allR  sc     VVRVVRVVRVVRVV4\\,4EFFBHHYqV<=> GFFFs   6(B((B9	c           
        \         P                  ! \        \        R 7      ;_uu_ 4        \        P
                  ! \        P                  \        P                  \        P                  \        P                  \        P                  .\        P                  \        P                  \        P                  \        P                  \        P                  ..4      pV! VR^R7      p\        \        P                  ! V4      P                  4       4       \        \        V4      ^8H  4       RRR4       R#   + '       g   i     R# ; i)r   r  )r  r  N)r   r   r   r   r%   r9  r   r	   r  r   r  )rI   r?  rY   r  s   &&  r+   test_nan_omit_all_axis-TestCircFuncsNanPolicy.test_nan_omit_all_axisY  s     \\,4EFF266266266266266B66266266266266BD EAA&q9CBHHSM%%'(CHM" GFFFs   C>D11E	rY   rX  r  rZ  r[  r\  r]  c                ,    \        \        WR RR7       R# )r9  r  r  NrQ   rR   r  s   &&&r+   test_nan_raise%TestCircFuncsNanPolicy.test_nan_raisec  s     	j)SWMr-   c                ,    \        \        WR RR7       R# )r9  foobarr  Nr  r  s   &&&r+   test_bad_nan_policy*TestCircFuncsNanPolicy.test_bad_nan_policyl  s     	j)SXNr-   r_   )g     v@r  r  g     8v@rd  g     u@)gaF+?z?r  gA ~@碲?r  )gj  @竃 ?r  g)h;@{ ?r  )r`   ra   rb   rc   r   r   r   r   rE  r%   r9  rG  rF  rH  r  r  r  r  r   r  r  rd   re   rf   s   @r+   r  r  "  sK     [[1$~~ $&7 !288 -4 $5 !288Z,F#G IJ
  %}} $&: !288 -N $./1uuSy1n$= !"288-A-@-B $C	 DE  %}} $&7 !288 -9 $: !"288Z,D#E	 GHIJ&/'J&/ [[1$~~{;$}}.BC$}}k:<=-	=-
 [[[5>>5==+0==+: ;?;?
 [[[5>>5==+0==+: ;#;# [[S"Aq#r3? hhaCS"&&(I),aCC(H(J KLM [[[5>>5==+0==+: ;N;	MN [[S"Aq#r3? hhaCS"&&(I),aCC(H(J KLM [[[5>>5==+0==+: ;O;	MOr-   r  c                      a  ] tR tRt o R tR tR tR tR tR t	R t
R	 tR
 tR tR t]P                   P#                  RRR.4      R 4       tRtV tR# )TestMedianTestiv  c                H    \        \        \        P                  . RO4       R# rt   Nrz   rQ   rR   r   median_testrS   s   &r+   test_bad_n_samples!TestMedianTest.test_bad_n_samplesx  s    j%"3"3Y?r-   c                J    \        \        \        P                  . . RO4       R# r  r  rS   s   &r+   r+   TestMedianTest.test_empty_sample|  s    j%"3"3RCr-   c           	     X    \        \        \        P                  . RO. RO^^ .RR7       R# )rt   r   tiesN)rt   rt   rt   rt   )r{   r   rt   r  rS   s   &r+   test_empty_when_ties_ignored+TestMedianTest.test_empty_when_ties_ignored  s$    
 	j%"3"3"I1vH	Fr-   c                    \        \        \        P                  . RO. RO4       \        \        \        P                  . RO. RORR7       R# )rt   abover  N)rt   rt   rt   r  rS   s   &r+   test_empty_contingency_row)TestMedianTest.test_empty_contingency_row  s3    
 	j%"3"3Y	J 	j%"3"3Y	"	$r-   c                R    \        \        \        P                  . RO^^.RR7       R# )rt   rW  r  Nrz   r  rS   s   &r+   test_bad_tiesTestMedianTest.test_bad_ties  s    j%"3"3YA 	"r-   c                R    \        \        \        P                  . RO^^.RR7       R# )rt   r  r  Nrz   r  rS   s   &r+   r  "TestMedianTest.test_bad_nan_policy  s    j%"3"3YA!)	+r-   c                R    \        \        \        P                  . RO^^.RR7       R# )rt   rW  )rW  Nrz   )rQ   r  r   r  rS   s   &r+   r  TestMedianTest.test_bad_keyword  s    i!2!2I1v	!r-   c                    . ROp. ROp\         P                  ! W4      w  r4rV\        VR4       \        V^^.^^..4       \        V^ 4       \        V^4       R# )rt   r  Nrz   )r   r  r   r   )rI   rY   r7  rB  r  medtbls   &      r+   test_simpleTestMedianTest.test_simple  sZ    !--a3 	S#3!Q!Q 01 	T1Qr-   c                v   . ROp^^.p. ROp\         P                  ! WV4      w  rErg\        V^4       \        V. RO. RO.4       \         P                  ! WVRR7      w  rErg\        V^4       \        V. RO. R	O.4       \         P                  ! WVRR7      w  rErg\        V^4       \        V. R
O. R	O.4       R# )rt   r   r  r  NrN  )r?   r@   r=   )r   rt   r|   )r~   rt   r   )r~   r   r   )r   r{   r|   )r   r  r   )rI   rY   r7  zrB  r  r<  r  s   &       r+   test_ties_options TestMedianTest.test_ties_options  s    F  ++A!4QS9i01++A!(CQS9i01++A!'BQS9i01r-   c                   ^^\         P                  .p. R	Op\        P                  ! WRR7      p\        P                  ! WRR7      w  rErg\	        V\         P                  \         P                  \         P                  R34       \        VR4       \        VR4       \	        VR4       \	        V\         P                  ! ^ ^.^^..4      4       \        \        \        P                  WRR7       R# )
rt   	propagater  r  Ng     ?gإvao?r  r  r|  )	r%   r   r   r  r   r
   r9  rQ   rR   )rI   rY   r7  mt1r  r  r<  r  s   &       r+   test_nan_policy_options&TestMedianTest.test_nan_policy_options  s    266N=&&q?
aS266266266489././QQ1a&1a&!123j%"3"3QgNr-   c                   . ROp. ROp\         P                  ! W4      w  r4rV\        V^4       \        V^^.^^..4       \         P                  ! V4      w  rxr\	        W74       \	        WH4       \         P                  ! W^ R7      w  r4rV\        V^4       \        V^^.^^..4       \         P                  ! V^ R7      w  rxr\	        W74       \	        WH4       \         P                  ! WRR7      w  r4rV\        V^4       \        V^^.^^..4       \         P                  ! VRR7      w  rxr\	        W74       \	        WH4       R# )rt   )lambda_FrC  Nr}   )r{   r~   r1   r@   )r   r  r   chi2_contingencyr
   )rI   rY   r7  rB  r  r<  r  exp_statexp_pdofes   &          r+   rN   TestMedianTest.test_basic  s$    ++A1QSAq6Aq6*+"'"8"8"='!++A!<QSAq6Aq6*+"'"8"8a"H'!++AUCQSAq6Aq6*+"'"8"8"O'!r-   r  FTc                    . ROp. ROp\         P                  ! W#VR7      p\        VP                  VP                  VP
                  VP                  3V4       R# )rt   rC  Nrz   )r   r  r   rF   r   r  table)rI   r  rY   r7  r[   s   &&   r+   r  TestMedianTest.test_result  sA    <cmmSZZSYYGMr-   r_   N)r`   ra   rb   rc   r  r+  r  r  r  r  r  r  r  r  rN   r   r   r   r  rd   re   rf   s   @r+   r  r  v  sk     @DF
$"+!2(O"> [[\E4=9N :Nr-   r  c                   `  a  ] tR tRt o R t]P                  P                  R]P                  ) ^,          ]P                  ^,          .R3^ ^]P                  ,          .R3.4      R 4       t
R tR t]! R	R
R7      R 4       tR t]P                  P                  RRR.4      R 4       tRtV tR# )TestDirectionalStatsi  c                   \         P                  ! \         P                  ! . RO4      4      ) p\         P                  ! \         P                  ! . RO4      4      ) p\         P                  ! \         P                  ! V4      \         P                  ! V4      ,          \         P                  ! V4      \         P
                  ! V4      ,          \         P
                  ! V4      3^R7      pVP                  VP                  4       4      pVP                  VP                  4       4      pVP                  VP                  4       4      p\        P                  ! V4      pVP                  pVP                  . RO4      p\        WgRR7       R# )33333su@r  r  r0  N)	r  g      O@g33333sB@g      ;@r;  g@g333333I@gYv@r  )	gfffffP@g,Q@gfffffQ@gfffffT@g     S@g     @R@g33333SQ@gffffffM@g33333I@)gHP?g&S:g[<)r%   deg2radr9  r  r  r  r  rc  r   directional_statsmean_directionr   )rI   r  declinclrJ   dirstatsdirectional_meanreference_means   &&      r+   "test_directional_stats_correctness7TestDirectionalStats.test_directional_stats_correctness  s    

288 %< = > >

288 %< = > >xxt4t4'  
 zz$++-(zz$++-(zz$++-(**40#22$>?(tDr-   zangles, refrt  r   c                   VP                  V4      pVP                  V4      pVP                  VP                  V4      VP                  V4      .^R7      p^\        P
                  ! V4      P                  ,
          p\        WR4       R# )rt   r  N)r  r  r  r  r   r  mean_resultant_lengthr   )rI   anglesr   r  rJ   r[   s   &&&&  r+   'test_directional_stats_2d_special_cases<TestDirectionalStats.test_directional_stats_2d_special_cases  si    
 F#jjoxx8qxA%))$/EEE!r-   c                |   \         P                  P                  R 4      pVP                  ^VP                  ,          VP                  R4      ,          4      pVP                  VP                  V4      VP                  V4      3^R7      p\        P                  ! V4      pVP                  pVP                  V^,          V^ ,          4      pV^VP                  ,          ,          p\        P                  ! V4      p\        Wx4       RVP                  ,
          p	\        P                  ! V4      p
\        W4       R# )l	   }7Qh0P4u&4M r  rt  N)rb  )r%   rd  re  r  rF  r  r  r  r   r  r  r  rE  r   r  rG  )rI   r  rh  testdatatestdata_vectorr  r  directional_mean_anglerE  directional_varcircular_vars   &&         r+   test_directional_stats_2d.TestDirectionalStats.test_directional_stats_2d  s     ii##$FG::a"%%i#**X*>>?((BFF8$4$&FF8$4$6() # + **?;#22!#*:1*=?OPQ?R!S!71ruu9!E>>(+.9x===}}X.6r-   c                   VP                  . RO. RO.4      pVP                  VP                  VR4      4      pVP                  . RO. RO.. RO. RO..4      p\        P                  ! V^R7      p\	        VP
                  V4       R# )+Pz?r  Nr  r  r   )r  r8  r   )r{   r{   r{   rt   )rt  r   r   )r  tiler   r  r   r  )rI   r  rJ   
full_arrayro  r  s   &&    r+    test_directional_mean_higher_dim5TestDirectionalStats.test_directional_mean_higher_dim-  s     zz/02 3ZZl ;<
:: , . , , ./ 0 **:A>//:r-   Tzchecking array-like inputr  c                b   . RO. RO.pVP                  W!P                  R7      p\        P                  ! V4      p\        P                  ! V4      p\	        VP
                  VP                  VP
                  4      4       \	        VP                  VP                  VP                  4      4       R# )r  r  Nr  )r  r8  r   )r  r  r   r  r   r  r  )rI   r  rJ   
data_arrayr   r[   s   &&    r+   )test_directional_stats_list_ndarray_input>TestDirectionalStats.test_directional_stats_list_ndarray_input:  s     %&:;ZZJJZ7
%%d+%%j1**

3#5#56	811

3#<#<=	?r-   c                   VP                  R4      pRp\        P                  ! \        \        P
                  ! V4      R7      ;_uu_ 4        \        P                  ! V4       RRR4       R#   + '       g   i     R# ; i)r   zIsamples must at least be two-dimensional. Instead samples has shape: (5,)r   N)r   )r   r   r   rR   reescaper   r  )rI   r  rJ   rP  s   &&  r+   test_directional_stats_1d_error4TestDirectionalStats.test_directional_stats_1d_errorF  sN    wwu~6]]:RYYw-?@@##D) A@@@s   A--A>	r  r  r  c                   \         P                  ! . RO. RO.VR7      p\        P                  ! VP	                  V4      RR7      pV\         P
                  P                  VR	RR7      ,          p\        P                  ! VRR7      p\        VP                  VP	                  VP                  4      4       \        VP                  VP	                  VP                  4      4       R# )
r  r  T)	normalize)r  keepdimsFNr  )g+Pz?r   r   r   )
r%   r9  r   r  r  r  r   r   r  r  )rI   r  r  rJ   r[   normalized_datar   s   &&&    r+    test_directional_stats_normalize5TestDirectionalStats.test_directional_stats_normalizeN  s    
 xx--/6;=%%bjj&6$G29= "0 "? ?%%oG**

3#5#56	811

3#<#<=	?r-   r_   N)r`   ra   rb   rc   r  r   r   r   r%   rF  r  r  r  r  r   r  r  rd   re   rf   s   @r+   r  r    s     E, [[]55&(BEE!G	b!
QY- "	"7&; d+FG	? H	?* [[Wy)&<=? >?r-   r  c                      a  ] tR tRt o R tR t]P                  P                  R. ROR3. ROR3.4      R 4       t	R t
]P                  P                  R	. RO4      R
 4       tR tRtV tR# )TestFDRControli_  c                d   R p\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RO4       RRR4       \         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RO4       RRR4       \         P                  ! \        VR7      ;_uu_ 4        \        P                  ! RR\
        P                  .4       RRR4       Rp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RORR7       RRR4       Rp\         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . ROR	R
7       RRR4       \         P                  ! \        VR7      ;_uu_ 4        \        P                  ! . RORR
7       RRR4       R#   + '       g   i     EL; i  + '       g   i     ELT; i  + '       g   i     EL; i  + '       g   i     L; i  + '       g   i     L; i  + '       g   i     R# ; i)z.`ps` must include only numbers between 0 and 1r   r  r  NzUnrecognized `method` 'YAK'YAKrA  z#`axis` must be an integer or `None`ru  r  )r   r  r  )r  r  r{   )r  r  rk   rp  )r   r   rR   r   false_discovery_controlr%   r   r9  s   & r+   test_input_validation$TestFDRControl.test_input_validation`  s=   B]]:W55)).9 6]]:W55))-8 6]]:W55))3RVV*<= 6 0]]:W55))/%H 6 8]]:W55))/D 6]]:W55))/G 65 655555555 65 65555sG   F<,G1(G$G8HH<G	G!	$G5	8H	H	H/	c                \    . ROp\         P                  ! V4      p. ROp\        W#RR7       R# )r  r.  r0  N)r  g;On?g~jt?gI+?g&1?gB`"۹?gv?S㥛?);On?r  I+?r  gˡE?gHzG?r  r  )r   r  r
   )rI   psr[   r   s   &   r+   test_against_TileStats%TestFDRControl.test_against_TileStatss  s&    E++B/Ft,r-   r  bhbyc                    Vw  r#\         P                  P                  R 4      p\        P                  P                  RR^
VR7      pV^,          V^&   \        P                  ! WSR7      p\        WbRR7       R# )   <U) r.  r  rs  rA  r3   r0  N)r%   rd  re  r   
loguniformr(   r  r
   )rI   r  r   r  rh  r  r[   s   &&     r+   r  TestFDRControl.test_against_Rz  sf     ii##$78!!$"3!G11++B>t,r-   c                   \         P                  P                  R 4      p\        P                  P                  RRRVR7      p\        P                  ! VRR7      p\        P                  ! VP                  4       4      p\        W44       R# )r  r.  r  rs  Nr  r|   r~   r   )	r%   rd  re  r   r  r(   r  ravelr   )rI   rh  r  r[   r   s   &    r+   test_axis_NoneTestFDRControl.test_axis_None  sf    ii##$78!!$)#!N++BT:++BHHJ7Sr-   r  c                   \         P                  P                  R 4      p\        P                  P                  RRRVR7      p\        P                  ! W1R7      p\         P                  ! \        P                  W4      p\        WE4       R# )r  r.  r  rs  r  Nr"  )	r%   rd  re  r   r  r(   r  apply_along_axisr   )rI   r  rh  r  r[   r   s   &&    r+   r$  TestFDRControl.test_axis  sf    ii##$78!!$)#!N++B:!!%"?"?JSr-   c                    \        \        P                  ! R .4      R .4       \        \        P                  ! R 4      R 4       \        \        P                  ! . 4      . 4       R# )r  N)r   r   r  rS   s   &r+   r  TestFDRControl.test_edge_cases  sE    588$@4&I588>E588<bAr-   r_   N)
go?gQ]
Y?ghD@?u?D?gn@t?g=T	[ے?g=!s?r+  r+  g	L?)
glߢ?gq?g.0?R]D?g)<Gh?g$͂?g2-p?r,  r,  g]:?)r   rt   r   )r`   ra   rb   rc   r  r  r   r   r   r  r$  r$  r  rd   re   rf   s   @r+   r  r  _  s     H&- [[V 89=? 89=?@A	-A	- [[VZ0 1B Br-   r  c                      a  ] tR tRt o ]P
                  P                  R]P                  / 3]P                  R^/3]P                  R^/3]P                  / 3.4      R 4       tRtV tR# )TestCommonAxisi  r  r  c           
        Vw  r4\         P                  P                  R 4      pVP                  VP                  R4      4      pV! V3/ VBR^ /B pTP	                  \        VP                  ^,          4       Uu. uF  q! VRV3,          3/ VB NK  	  up4      p	\        Wy4       V! V3/ VBR^/B pTP	                  \        VP                  ^ ,          4       Uu. uF  q! WhR3,          3/ VB NK  	  up4      p	\        Wy4       V! V3/ VBRR/B pV! VP                  VR4      3/ VB p	\        Wy4       R# u upi u upi )l   @ms|Yr  r  N)r1   r?   r^  )	r%   rd  re  r  r  r|  r  r   r_  )
rI   r  r  rI  r*   rh  rY   r[   r\   r   s
   &&&       r+   r$  TestCommonAxis.test_axis  s*   
 ii##N3JJszz&)*!&v&A&hhaggaj8IJ8I1AadG.v.8IJK!!&v&A&hhaggaj8IJ8I1AdG.v.8IJK!!)v)D)"**Q&1&1! K Ks   :D:D?r_   N)r`   ra   rb   rc   r   r   r   r   semr  	variationr$  rd   re   rf   s   @r+   r.  r.    se      [[Vuyy"o',{{S!H&=',{{S!H&=',&;&= >"	>"r-   r.  )
jt?r&  V-?rt  gX9v?-?x&1?r3  +?rt  )
r4  r3  rt  r6  灕Cl?r4  r&  rt  r3  V-?)
ʡE?g/$?r8  r3  ףp=
?r7  rt  r3  r&  r&  )
gGz?r6  r7  rt  r;  r7  r4  r&  r6  r&  )
r4  r4  gCl?r  r6  rh|?r&  r5  g\(\?r&  )
羟/$?gh|?5?r=  r8  r9  r6  r;  r4  g1Zd?r&  )
r  Mb?r&  gjt?r4  rt  g}?5^I?g)\(?r&  r6  )
r4  rt  r3  rt  r6  r&  r4  r&  r6  r3  )
r6  r4  r&  r;  r&  r>  r>  r4  r3  r:  )
r:  r;  r<  r7  r8  r8  r:  r4  r>  r8  )gףp=
W0@gQ%@gGz'@g+@g33333/@g3333332@g
ףp=
@gq=
ףp$@g
ףp=
)@g\(\-@gGz0@gGzn3@g=
ףp= @gffffff&@gfffff)@g
ףp=-@Gz0@g3@g!@g)\('@g333333(@g
ףp=
.@r?  g{G0@gQ3@gq=
ף"@g\(\'@gR*@g33333.@g\(3@)ir  r{  r  rP  r   	functoolsr   numpyr%   numpy.randomr   numpy.testingr   r   r   r   r	   r
   r   r   r   r   rQ   scipyr   r   r   scipy.stats._morestatsr   r   r   common_testsr   
_hypotestsr   r   scipy.stats._binomtestr   scipy.stats._distr_paramsr   scipy.stats._axis_nan_policyr   r   r   r   scipy._lib._array_apir   scipy._lib._array_api_no_0dr   r   r    r   r  dict
matplotlibrcParamsmatplotlib.pyplotpyplotr  r  	Exceptionr{  r|  r}  r~  r  r  r  r  r  r  r,   r/   ri   r   r   rt  r1  ry  r  r  r[  r  r  r  r  r  r  r  r  r  r  r>  r'  rz  r  r  r  r(  r3  r  r  r  r  r.  r_   r-   r+   <module>rS     s  
   	 
    $. . .  * * * N N - B ? .T T +  ;;// >%*J	"#O LKKKKKKKKL
4/ />? ?Bu@ u@pN> N>bM9 M9`C6 C6L,* ,*^B7 B7Jk8 k8\G- G-T20l7 l7^^5 ^5BDH DHT
P
7. 7.t,. ,.^39 39l"= "=JC+ C+P0	000"0$*0,204:0<A0CI0
000!0#)0+10390;A0CH0 0 0 0 "0 $*0 ,10 3:0 <C0 EJ0 
	0 	0 	0 !	0 #)	0 +1	0 39	0 ;A	0 CI	0
 0
 0
 0
 %0
 '-0
 /60
 8>0
 @E0
 GM0 0 0 0 #0 %+0 -20 490 ;A0 CH0 0 0 0 "0 $*0 ,30 5:0 <B0 DK0 
0 0 0 $0 &,0 .50 7=0 ?D0 FL0 0 0 0 $0 &-0 /50 7<0 >D0 FL0 0 0 0 !0 #(0 *00 280 :@0 BH0 0 0 0 "0 $)0 +10 390 ;A0 CI0 0 0 0 "0 $*0 ,30 5:0 <A0 CI0 
0 0 0 "0 $*0 ,10 380 :A0 CI0 0 0 0 "0 $)0 +10 3:0 <C0 EK0 0 0 0 #0 %*0 ,20 4;0 =C0 EL0  
!0  !0  !0  #!0  %+!0  -2!0  4:!0  <B!0  DJ!0" #0" #0" #0" $#0" &+#0" -2#0" 4:#0" <C#0" EK#0$ %0$ %0$ %0$ #%0$ %+%0$ -3%0$ 5<%0$ >E%0$ GN%0& 
'0& '0& '0& #'0& %*'0& ,3'0& 5<'0& >C'0& EJ'0( )0( )0( )0( $)0( &,)0( .5)0( 7<)0( >E)0( GN)0* +0* +0* +0* #+0* %*+0* ,3+0* 5:+0* <C+0* EK+0, -0, -0, -0, #-0, %,-0, .5-0, 7=-0, ?F-0, HM-0. 
/0. /0. /0. !/0. #)/0. +1/0. 39/0. ;A/0. CI/00 100 100 100 #100 %+100 -4100 6=100 ?F100 HN102 302 302 302 %302 '-302 /5302 7<302 >D302 FK304 504 504 504 $504 &,504 .3504 5<504 >C504 EK506 706 706 706 $706 &+706 -4706 6<706 >D706 FL708 
908 908 908 !908 #)908 +2908 4:908 <A908 CI90: ;0: ;0: ;0: ";0: $*;0: ,2;0: 4:;0: <A;0: CH;0< 
=0< =0< =0<  =0< ")=0< +1=0< 3:=0< <B=0< DJ=0> ?0> ?0> ?0> %?0> ',?0> .5?0> 7=?0> ?D?0> FM?0@ A0@ A0@ A0@ !A0@ #)A0@ +2A0@ 4;A0@ =CA0@ EKA0B C0B C0B C0B $C0B &-C0B /5C0B 7=C0B ?FC0B HOC0D E0D E0D E0D %E0D ',E0D .4E0D 6=E0D ?FE0D HNE0F G0F G0F G0F #G0F %+G0F -3G0F 5<G0F >CG0F EJG0H I0H I0H I0H #I0H %+I0H -4I0H 6=I0H ?EI0H GMI0J K0J K0J K0J "K0J $)K0J +1K0J 38K0J :@K0J BHK0L M0L M0L M0L $M0L &-M0L /4M0L 6<M0L >DM0L FLM0N O0N O0N O0N %O0N ',O0N .5O0N 7=O0N ?EO0N GLO0P Q0P Q0P Q0P %Q0P '-Q0P /5Q0P 7>Q0P @EQ0P GMQ0R S0R S0R S0R #S0R %,S0R .4S0R 6<S0R >ES0R GLS0T U0T U0T U0T %U0T ',U0T .5U0T 7>U0T @GU0V W0V W0V W0V &W0V (.W0V 05W0V 7>W0V @FW0V HNW0X Y0X Y0X Y0X %Y0X '.Y0X 07Y0X 9>Y0X @EY0X GNY0Z [0Z [0Z [0Z #[0Z %,[0Z .5[0Z 7<[0Z >D[0Z FK[0\ ]0\ ]0\ ]0\ $]0\ &-]0\ /4]0\ 6=]0\ ?E]0\ GM]0^ _0fI. I.XK3 K3\ ;  ;F7 7(E6 E6P4 4"k@ k@\QO QOhN NBe? e?P?B ?BD" "ac  Os   ,Y 	Y,+Y,