(phA SSKrSSKrSSKJrJrJr SSKJr SSK J r SSK J r J r SSKJrJr SSKJrJr SSKJr SS KJrJrJrJr SSKJs Jr S r\R<R?S /S Q5\R<R?S /SQ5S55r \R<R?S /S Q5S5r!\R<R?S /S Q5\R<R?S /SQ5\R<R?SSS/5S555r"\R<RF\R<RIS5\R<R?S /S Q5S555r%SSSS .r&\R<R?S\&RO55S5r(SSSS .r)Sr*Sr+\R<RF\R<R?S\)RO55S 55r,S!S"S#.r-\R<RF\R<R?S\-RO55S$55r.\R<R?S S%S&/5\R<R?S SS'/5S(55r/\R<RIS)5\R<R?S /S Q5\R<R?S SS'/5S*555r0\R<R?S /S Q5S+5r1\R<R?S /S Q5S,5r2S-r3\R<R?S.SS//5S05r4\R<RIS15\R<R?S /S Q5S255r5S3r6\R<R?S4S5S6/5S75r7\R<R?S8/S9Q5\R<R?S /SQ5S:55r8S;r9\R<R?S /SQ5S<5r:\R<RF\R<R?S /S Q5S=55r;\R<RF\R<RyS>5S?55r=S@r>"SASB5r?"SCSD5r@"SESF5rASGrB\R<R?SH/SIQ5SJ5rC"SKSL5rDg)MN)assert_allclose assert_equalsuppress_warnings)array_api_compatible) rng_integers)array_namespaceis_numpy)xp_assert_closexp_assert_equal)statsspecial)root) bootstrapmonte_carlo_testpermutation_testpowerc Sn[R"[US9 [S[R 5 SSS5 Sn[R"[US9 [[ 5[R 5 SSS5 Sn[R"[US9 [/SQS/4[R 5 SSS5 Sn[R"[US9 [/SQ/SQ4[R S S 9 SSS5 S n[R"[US9 [S[R S S 9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R S S9 SSS5 Sn[R"[US9 [/SQ4[R SS9 SSS5 Sn[R"[US9 [/SQ4[R S S9 SSS5 S!n[R"[US9 [/SQ4[R S"S#9 SSS5 g!,(df  GN=f!,(df  GNe=f!,(df  GN4=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GN|=f!,(df  GNN=f!,(df  GN =f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GNh=f!,(df  GN:=f!,(df  g=f)$Nz%`data` must be a sequence of samples.matchz(`data` must contain at least one sample.z>each sample in `data` must contain two or more observations...rz=When `paired is True`, all samples must have the same length rrrT)paired0`vectorized` must be `True`, `False`, or `None`.ekki vectorized`axis` must be an integer.?axisz!could not convert string to floatni)confidence_levelz-`n_resamples` must be a non-negative integer. n_resamplesD@+`batch` must be a positive integer or None.batchz`method` must be inmethodz?`bootstrap_result` must have attribute `bootstrap_distribution' )bootstrap_resultz5Either `bootstrap_result.bootstrap_distribution.size`r,SeedSequence expects int or sequence of intsherringrng)pytestraises ValueErrorrnpmeantuple TypeError)messages T/var/www/html/venv/lib/python3.13/site-packages/scipy/stats/tests/test_resampling.pytest_bootstrap_ivr?s"5G z 1!RWW 29G z 1%'277# 2OG z 19qc"BGG, 2OG z 19l+RWWTB 2AG z 1!RWW0 2+G z 19,c2 22G z 19,$? 2>G z 19,U; 2>G z 19,V< 2test_bootstrap_batchrSPsIINN1 r2r"A aT277$v"#$C ADIINN1 aT277"VC 1D))--t/G/G/K/KL))..0H0H0M0MN$$d&9&9:c^^^[RRS5 Sn[RRU5m[RRU5mSSjmSUUU4Sjjn[R"[ T55n[ U4USS9n[ TT4TSSS9n[URUR5 [URUR5 g) NrrGc*X- S-RUS9$)Nrr#r:rPyr$s r> my_statistic+test_bootstrap_paired..my_statisticksqD))rTc,>TUnTUnT"X#5nU$N)ir$abresrZrPrYs r>my_paired_statistic2test_bootstrap_paired..my_paired_statisticns$ aD aD1  rTr4T)rr5) r9rIrJrKarangelenrrrLrO) r/nrcr_rQrRrZrPrYs @@@r>test_bootstrap_pairedrjcsIINN1 A qA qA* #a&A aT.A 6D aV\$A >DD,,d.F.FGD'')<)<=rTrTFc [RRS5 SSjnSnXAn[RRU5n[RRU5n[RRU5n[ XgU4X2USSSS9n U R R U RR S-:Xde/SQn XZU'[R"URU 5U5n[R"URU 5U5n[R"URU 5U5n[ XgU4X2USUSS9n [U RRU RR5 [U RRU RR5 [U RU R5 [U5n U RU5 [!U RRR U 5 [!U RRR U 5 [!U RR U 5 g) Nrc`URUS9URUS9-URUS9-$Nr#rW)rPrYzr$s r>rZ/test_bootstrap_vectorized..my_statistics3vv4v 166t6#44qvv4v7HHHrT)r0rErFrG)rr/r5r$r))rG)rrrre)r9rIrJrKrbootstrap_distributionshaperO broadcast_toreshaperrLrMrNlistpopr) r/r$rrZrq n_samplesrPrYrnrQrsrR result_shapes r>test_bootstrap_vectorizedrx}sIINN1I E I y!A y!A y!A aAY F 5D  ' ' - -""((61 23 2GDM  '*E2A  '*E2A  '*E2A aAY F3 8DD,,00,,002D,,11,,113D'')<)<=;LT))--33\B))..44lC$$**L9rTz#MemoryError with BCa observed in CIcb[RRS5n[RR SSSUS9nSn[R "[U5S- [R"U5[R"U5S9nURUS 9nUR5n[U4[RSXS 9n[S0UDS U0D6n[URUS S 9 [UR USS9 UR#[SUS95 [S0UDUSS.D6n[URR$UR'U5S S 9 UR#[SUS95 [S0UDUSS.D6n[URR(UR'SU- 5S S 9 g)Nl)ICbX r)locscalesizerHffffff?r)dfr|r}) confidence)data statisticr)r/r5r&gMb@?rtolga2U0*3?atolrr)r1lessr& alternativegreaterr^)r9rI default_rngr normrvstrhr:semintervalstddictrrrLrOupdaterNppfrM) r/r5ralphadistexpected_interval expected_seconfigrbs r>test_bootstrap_against_theoryrsa ))   3 4C ::>>aqt#> FD E 77c$ikrwwt}EIIdO LD  7((*K w"''t*F  5f 5u 5CC++->TJC&& $? MM$1s;<  If Iu& ICC++00$((5/M MM$1s;<  Lf Lu) LCC++//!E'1BNrT)gQ7@gHzS@)g\(<@g= ףp U@)gHz'@@gQV@zmethod, expectedc[R"/SQ5n[U4[RSUSS9n[ UR USS9 g)N)r0rF)@rg+@"2QYyi@Br)r)r/r5g{Gzt?r)r9arrayrr:rrL)r/expectedrPrbs r>test_bootstrap_against_RrsD" 5 6A QD"''wv CC++XEBrTiicT/SQn/SQnSn[RRS5n[R"X4USSUS9n[R"X4USSUS9n[R"X4US SUS9n[R "UR 5n[R "UR 5n[R "UR 5n S n Xz- [U 5- n X- [U 5- n X- [U 5- n U S :deU S :de[U 5S :deg)N)gbF?gp?g|Hܿgxy[e}u׿g*g_?ggi6g2?)gF?gK?g&F}/?g?g(iC?g}$ץ?ghu?gfj}!?cZ[R"XS9n[R"XS9nX4- $rm)r skew)rPrYr$s1s2s r>r1test_multisample_BCa_against_R..statistics% ZZ % ZZ %wrTltg y 4r@rG)r/r-r5rAbcag1*g333333ÿg333333?Q?)r9rIrr rr:rLabs)rPrYrr5 res_basic res_percentres_bca mid_basic mid_percentmid_bca mid_wboot diff_basic diff_percentdiff_bcas r>test_multisample_BCa_against_Rrs'  .A 9A  ))   2 3C '&)s4I//1&)L(+6Kooqfi$'S2G 556I''+99:Kggg112GI'Y7J+S^;L#S^3H    $   x=4  rTc[R"/SQ5n[R"/SQ5nSSjnX/n[R"X25nSnSnURnSn[ R "X2XVXx5u p[U S5 g) N) r0(.r4h)r&^ccT[R"XS9[R"XS9- $rmr9r:)rnrYr$s r>r:test_BCa_acceleration_against_reference..statistic,swwq$rwwq'<<'test_BCa_acceleration_against_referencer"s 78A /0A= 6D //$ *C D E,,K E++DT,7@KAqE/0rTc [RRS5 SnSnSn[R"SSS9nUR 5nSnUR Xr4S 9n[U4[R UUS US S 9n U Rn [R"U SU:XjS:-5n X:Xde[R"XU5Rn U S :deg)NrrG?rzrr|r}ir~rrfrr&r)r-r/r$皙? r9rIrJr rr:rrrLsum binomtestpvalue) r/rrir)r&r stat_truen_replicationsrrbcici_contains_truers r>#test_bootstrap_against_itself_1sampr;sIINN1 AK ::!1 %D IN 88.,8 -D TG gg%5 +!  C Bvvr!uy0YA5FGH  '' '__--//5v C<my_stat4test_bootstrap_against_itself_2samp..my_statws&))}rTrrrrrrrer)r/rn1n2r)r&rdist1dist2rrrrrbrrrs r>#test_bootstrap_against_itself_2samprgsIINN1 B BK JJ1A &E JJ3a (E uzz|+IN IIN/I 0E IIN/I 0E UN%%5 +!  C Bvvr!uy0YA5FGH  '' '__--//5v C<# UHoRT5v M g7fr]rW.0sampler$s r> Etest_bootstrap_vectorized_3samp..statistic..84;;t$$4 rr$rs` r>r2test_bootstrap_vectorized_3samp..statistic84888rTcJ>UHnURS:XaMe T"USS06$Nrr$rndimrrrs r> statistic_1d5test_bootstrap_vectorized_3samp..statistic_1d/F;;!# ##$'Q''rTrrzTrG)r r$r)r/r5Fr9rIrJrKrrrLrO) r/r$r rPrYrnrQrRrs @r>test_bootstrap_vectorized_3samprs 9( IINN1 q!A q!A q!A aAY dCA GD aAY CA GDD,,d.F.FGD'')<)<=rTz'Failure is not concerning; see gh-14107c ^^S SjmU4Sjn[RRS5 [RRSS5n[ U4TSUSSUSS9n[ U4US USS USS9n[ UR UR 5 [ URUR5 g) Nrc URUS9$rmrWrPr$s r>r2test_bootstrap_vectorized_1samp..statisticsvv4v  rTc6>URS:XdeT"USS9$Nrrr#r)rPrs r>r 5test_bootstrap_vectorized_1samp..statistic_1ds vv{{##rTrrzTrG)r r$r)r-r/r5Fr0rr)r/r$r rPrQrRrs @r>test_bootstrap_vectorized_1samprs!$ IINN1 q!A aT9D!$D D aTtest_bootstrap_degenerater"s =D  [[ *DCe99E"''&AS44rvvrvv6FGF+ * "''&9S,,.>?##Q' FE+ *s$&C?AC.C?. C< 8C?? D c J[RRS5n[R"SSS9nUR SUS9nU4n[ U[RUS[RRS5S9n[ U[RUS[RRS5S S 9n[URUR5 [URUR5 [UR[R5(deg) NlE<1 rrrrGr~rHi[%r/r)r5F)r/r)r5r ) r9rIrr rrrrrrLrO isinstancefloat64)r/r5rrrbrefs r>test_bootstrap_gh15678r)s ))   0C ::!1 %D 8838 /D 7D D%**V --d3 5C D%**V --d3 GCC++S-D-DEC&&(:(:; c(("** 5 55 5rTc [RRS5n[R"SSS9nUR SUS9n[R "U5nU4n[U[R SS[RRS5S 9nX4RR:Xde[[R"U5*[RSS[RRS5S 9n[URR*URR5 [URR*URR5 g) Nln!|'x`5rrrrGr$rBifr%)r9rIrr rrminrrLrMrmaxrrN)r5rrtrue_minrbrRs r>test_bootstrap_minr.s ))   0 1C ::!1 %D 8838 /Dvvd|H 7D D"&&C --d3 5C ..22 22 2 bhhtn_bffU..t4 6DS,,000,,113S,,111,,002rTadditional_resamplesrc T[RRS5nURSS9nSnUnX0-n[RRS5n[R"U4[R X1SSS9n[R"U4[R XAS S US 9n[RRS5n[R"U4[R XQS S S9n[ URURS S 9 [ URURS S 9 g)Nl ]4MrGrrlmdpi4 rrA)r)r5r&r/rrB)r)r5r&r/r1+=r) r9rIrr rr:rrOrL)r/r5rPrrn3rbr(s r>test_re_bootstrapr3s ))   0 1C  A B B "B ))   2 3C //1$R+/  FC //1$R+/+. 0C ))   2 3C //1$R+/ ?CC&&(:(:GC++S-D-D rTzSensible to machine precisionc[RRS5n[R"SSS9nUR SUS94n[ U[RUSS9n[R"S0UDS S 0D6nUR[ S US 95 [R"S0UDS SS.D6n[R"S0UDS SS.D6n[URRURRSS9 [URRURRSS9 [R"URR5(de[R"URR5(de[ R""[$SS9 [R"S0UDSS0D6 SSS5 g!,(df  g=f)Nl@QixtrrrrGr$rf)rrr5r$r&rrrrrrrr1rz`alternative` must be one ofrrz ekki-ekkir^)r9rIrr rrrrrrrrLrNrMisneginfisposinfr6r7r8)r/r5rrrrlgs r>test_bootstrap_alternativer9sk ))   3 4C ::!1 %D HH3cH 2 4D trvv3R @F 7&737A MM$1q9: L&L4VLA O&O4YOAA))..0E0E0J0J A))--q/D/D/H/H ;;q,,00 1 11 1 ;;q,,11 2 22 2 z)G H :&:k: I H Hs *G Gc Sn[RRS5 [RR"U6n[ [ R "U55n[US5H?nUSUSS24n[R"XSS9n[R"XE5(aM?e [R"[[ R "USS95SS9n[R"Xb5(deg) Nrrrzrrf.r#rr,) r9rIrJrKnextr_jackknife_resamplerangedelete array_equal concatenatert)rqrPrYr_slcry2s r>test_jackknife_resamplerF1s EIINN1 A [ , ,Q /0A 59 Q l99Q+~~c,,,,  [<>"  rTrng_name RandomStaterc[[RUS5nUc[R"US35 U"S5nU"S5nSnSn[RR S5 [RR "U6n[R"XdUS9n[U5HEnUSUSS24n [USUSUS5n USU 4n [R"X5(aMEe g)Nz not available.rr0r;r4.rf) getattrr9rIr6skiprJrKr_bootstrap_resampler@rrB) rGr5rng1rng2r)rqrPrYr_rDjsrs r>test_bootstrap_resamplerPDs "))Xt ,C { xj01 q6D q6DK EIINN1 A''DAA ; Q l $59eBi 8S"W:~~c,,,, rTscore)r?rcSn[RRS5 [RR"U6n[R "X0SS9nSnU"X0SS9S- n[ XFS5 g)N)r0rrfr#cN[R"[RX U5$r])r9apply_along_axisr percentileofscore)r`rQr$s r>vectorized_pos0test_percentile_of_score..vectorized_posfs""5#:#:DUKKrTrGV瞯<)r9rIrJrKr_percentile_of_scorer)rQr$rqrPprYp2s r>test_percentile_of_scorer_^s` EIINN1 A((;AL r *3 .BA5!rTcpSn[RRS5 [RR"U6n[RR"USS6S-n[R "X5n[ US5H-nX4n[R"XX$SS9n[XVS5 M/ g)N)r0rTrrfrGr#r[) r9rIrJrKr_percentile_along_axisr@rAr)rqrPqrYr_rbrs r>test_percentile_along_axisrcns EIINN1 A cr #c)A**10A 58_d==qt"5u-rTcz^SmU4Sjn[R"U5n[RR S5 [RR SSS5n[RR SSS5n[RR SSS5nT"X4XPS9nU"X4XPS9n[ Xg5 g) Nc.^[U4SjU55$)Nc3D># UHoRT5v M g7fr]rWrs r>r>test_vectorize_statistic..statistic..rrrrs` r>r+test_vectorize_statistic..statisticrrTcJ>UHnURS:XaMe T"USS06$rrr s r>r .test_vectorize_statistic..statistic_1dr rTrrrzr<rr#)r_vectorize_statisticr9rIrJrKr) r$r  statistic2rPrYrnrQrRrs @r>test_vectorize_statisticrms9(11,?JIINN1 q!QA q!QA q!QA Q1 (D aA )DDrTc [RRS5nSn[RR "USUS.6nSn[ U4X@SSSS 9n[R"URRRU:URRRU:-S S 9n[R"US :5(de[R"US :*5(deURRRS:XdeURRRS:XdeURRS:XdeURRS:Xdeg)Nls:)rrR)rGrGr$c [R"[R"X5[R"XSS9/5$)Nr)ddof)r9asarrayr:r)rr$s r>r/test_vector_valued_statistic..statistics1zz2774.66$1578 8rTrf')r/r$r)r-rr#ZrGrrG)rrGrs)r9rIrr rrrrrLrMTrNallrqrOrp)r/r5paramsrrrbcountss r>test_vector_valued_statisticr{sE ))   +C F ZZ^^V*3 GF8 VIyb $C 1CVVS,,0022V;..3355>@F 66&B,    66&C-   " " & & , , 88 8  " " ' ' - - 99 9    # #x // /  % % + +~ == =rTzignore::RuntimeWarningc^^[RRS5nSmU4SjmU4Sjn/SQ/SQ/SQ/SQ/S Q/n[R"U5Rn[ X!US S 9n[ UTUS S 9n[ URRS URRS S9 [ URRS URRS S9 g)Nl\QcURU5nURU5nXS- XS- -RU5nSU-URU5URU5-X4- S--- $)N).Nr)r:var)rPrYr$xmymcovs r> concordance9test_vector_valued_statistic_gh17715..concordancesn VVD\ VVD\yM!aY-&78>>tDCAEE$K!%%+5QFGGrTcN>X-nX-n[R"T"XVU55$r])r9 nan_to_num)tptnfpfnr$actualrrs r>r7test_vector_valued_statistic_gh17715..statistics'7}}[4@AArTc>>T"/UQUP76[RS4$)N.)r9newaxis)r$argsrs r>statistic_extradim@test_vector_valued_statistic_gh17715..statistic_extradims#%$%%bjj#o66rT)rrrr)rrrr)rr<rr)rr<rr)rrrT)r5rrr[r) r9rIrrrwrrrLrMrN)r5rrrbr(rrs @@r>$test_vector_valued_statistic_gh17715rs ))   2 3CH B 7       D 88D>  D D#d CC D)T :CC++//2++//e=C++003++00u>rTc4[RRS5nURS5nURS5nSn[R"X4U5 [R"UR UR 4USS9 Sn[ R"[US9 [R"XSS 2S 44U5 SSS5 [ R"[US9 [R"XSS 2S S244U5 SSS5 [ R"[US9 [R"UR UR S S2SS 244USS9 SSS5 g!,(df  N=f!,(df  Nw=f!,(df  g=f) NiF|)r0r)rErc@[R"XUS9R$rmr ttest_indrrXs r>r test_gh_20850..statisticsq$/999rTrr#z-Ignoring the dimension specified by `axis`...rr0r) r9rIrr rrwr6r FutureWarning)r5rPrYrr=s r> test_gh_20850rs) ))   +C 7A 7A: OOQFI& OOQSS!##J 2>G m7 3 crc1fI 2 4 m7 3 crc1Q3hK()4 4 m7 3 acc!A#ss(m,ia@ 4 3 4 3 3 3 3 3s$$ E'%"E8(6F ' E58 F Fc p\rSrSrSrS\4SjrSr\S5r \S5r \ RR\S55r\\ RRS \"S S 55S 55r\\ RRS S5S55r\ RR&\ RRS S5\ RRS\R("SSS55S555r\ RRS\R.\R045\ RRS S5\ RRS\R("SSS55S555r\ RRS\R4"SS 55S5r\ RR\ RRS\R("SSS55S55r\ RR&\ RRSS5\ RRS\"S55S 555rS!rS"rS#r \ RRCS5\ RRES$5S%55r#S&r$g)'TestMonteCarloHypothesisTesti皙?Nc^^^^UUUU4Sj$)Nc4>TRT"UST0UD6TS9$)NrHdtype)rq)rkwdsrrsrvs_inxps r>6TestMonteCarloHypothesisTest.get_rvs..s(RZZ0V20VQU0V6;&0&=rTr^)selfrrrrs ````r>get_rvs$TestMonteCarloHypothesisTest.get_rvss = =rTc^U4SjnU$)Nct>TRXS9nTRXSS9nURUnX#U- S-- $)Nr#r)r$ correctionrR)r:r~rq)rPr$mvrirs r>r=TestMonteCarloHypothesisTest.get_statistic..statisticsB%Aq2A A!cz> !rTr^)rrrs ` r> get_statistic*TestMonteCarloHypothesisTest.get_statistics " rTc ^TR/SQ5nS#U4SjjnSnTRS5TRS54n[RR[RR4n[ R "[US9 [XVSSS 9 SSS5 S n[ R "[US9 [U[RRUS S 9 SSS5 S n[ R "[US9 [U[RRUS S 9 SSS5 Sn[ R "[US9 [USU5 SSS5 [ R "[US9 TRSS/SS//5n[USS/U5 SSS5 Sn[ R "[US9 TR/SQ/5n[USS/U5 SSS5 Sn[ R "[US9 [U[RRS5 SSS5 Sn[ R "[US9 [U[RRUSS9 SSS5 Sn[ R "[US9 [U[RRUSS9 SSS5 Sn[ R "[US9 [U[RRUSS9 SSS5 Sn[ R "[US9 [U[RRUSS9 SSS5 Sn[ R "[US9 [U[RRUSS 9 SSS5 S!nU4S"jn[X&TR5 g!,(df  GN3=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GNq=f!,(df  GN:=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GNa=f!,(df  GN+=f!,(df  N=f![a+n[U5RU5(deSnAgSnAff=f)$N?@@c">TRXS9$rmrW)rPr$rs r>stat@TestMonteCarloHypothesisTest.test_input_validation..stat s7717( (rTz/Array shapes are incompatible for broadcasting.rrz)rrzrcgNrr^rXs r>rDTestMonteCarloHypothesisTest.test_input_validation..s1rTrfr#r!r"rr0`rvs` must be callable or sequence of callables.rrrg@cU$r]r^rPs r>rr#arTIf `rvs` is a sequence...cU$r]r^rs r>rr(rrTcU$r]r^rs r>rr(s1rTz`statistic` must be callable.)`n_resamples` must be a positive integer.r'r(r*r+r,`alternative` must be in...rrz!Signature inspection of statisticc\>TR[RRUS95$)Nr)rqr rr)r~rs r>r?TestMonteCarloHypothesisTest.test_input_validation..rvsDs!::ejjnn$n78 8rTr]) rqzerosr rrr6r7r8rr<r:str startswith)rrrrr=tempres ` r>test_input_validation2TestMonteCarloHypothesisTest.test_input_validationsNzz,' )D "((6"23zz~~uzz~~. ]]:W 5 T(<2 F6/ ]]:W 5 T5::>>4c B6E ]]:W 5 T5::>>4C H6E ]]9G 4 T4 .5 ]]9G 4::Bx"b23D TK#6 =5. ]]:W 5::|n-D TK#=t D62 ]]9G 4 T5::>>4 85> ]]:W 5 T5::>>4U K6> ]]:W 5 T5::>>4V L6@ ]]:W 5 T5::>>4u E6@ ]]:W 5 T5::>>4v F60 ]]:W 5 T5::>>4V L66 9 . T 0i6 56 56 55 4 4 4 6 5 5 46 56 56 56 56 5 .q6$$W-- -- .s O>%O%O+O==(P&P!&P3%Q"%Q*%Q)2%Q;:%R 0R O O(+ O:= P  P! P03 Q Q Q&) Q8; R  R S(!SSc^U4SjnSnTR/SQ5n[T5(a&[U[RR U5 g[ R"[US9 [U[RR U5 SSS5 [ R"[US9 [U[RR TRSS9 SSS5 g!,(df  N_=f!,(df  g=f)Nc&>TRU5$r]rW)rPrs r>non_vectorized_statisticWTestMonteCarloHypothesisTest.test_input_validation_xp..non_vectorized_statisticMs771: rTz!`statistic` must be vectorized...rrFr) rqr rr rrr6r7r8r:)rrrr=rs ` r>test_input_validation_xp5TestMonteCarloHypothesisTest.test_input_validation_xpKs 6L) B<< VUZZ^^5M N  ]]:W 5 VUZZ^^5M N6 ]]:W 5 VUZZ^^RWW O6 56 5 5 5s,&C+3/C<+ C9< D cx^^^ [RRS5nURUR SS95n[ U5m UUU 4SjmSTlSTlUTSSS.nTR[RR[RRS 5US 9US '[SS S 0UD6n[TR S5 [TRS 5 TR[RR[RRS 5US 9US 'STl[SS S0UD6n[TR S5 [TRS5 TR[RR[RRS 5US 9US 'STl[S0UD6n[TR S5 [TRS5 [URUR5 [URUR5 g)N)pr0rc>URS:XaSOURSn[UTR5TlT=RS- slTR T5"XS9$r)r rqr, batch_sizecounterr)rPr$rrrxp_tests r>r:TestMonteCarloHypothesisTest.test_batch..statisticdsXffkqwwqzJ#&z93G3G#HI    " %%g.q< test_batch'TestMonteCarloHypothesisTest.test_batch[s ii##K0 JJs***3 4!!$ =   )#49ll5::>>2993H3H3PUWlXU 0a040Y&&-Y))1-ll5::>>2993H3H3PUWlXU  1b1D1Y&&+Y))2.ll5::>>2993H3H3PUWlXU  '$'Y&&*Y))40 T[[1 T[[1rTr$rc [RRS5n/SQnSXA'URS5RnS[ U5;aSOSn[ [U5n[ X&5nURXGS9n [R"U S US 9n URXS 9n [U 5n URU 5n UR[RRX8US 9n [XU S SUS9nURU R US 9nURU R"US 9n[%UR U5 [%UR"UUR&S9 g)Nit$rrrrGr64float32r')r~r)popmeanr$r)rrTi N)r r)r$r)r9rIrrqrrrJrr ttest_1samprrrrrrrrr r)rr$rr5r~ dtype_default dtype_strdtype_nprrPrrrrrb ref_statistic ref_pvalues r>test_axis_dtype,TestMonteCarloHypothesisTest.test_axis_dtypes8 ii##G,  2,, !%]);!;I) 2y)&   T  :$$Q> JJqJ &!!$&&w/ ll5::>>3lCqyT+0t= 8#5#5U C ZZuZ=  }5 JTYY?rTr) two-sidedrrc[RRS5nURSS9n[R "USUS9nUR U5n[U5nURU5nUR[RRX2S9n[XHXqS9n [U RUR UR55 [U RUR UR5UR S9 g)N)rUrrrrr)r9rIrrr rrqrrrrrrr rrr) rrrr5rPr(rrrrbs r>test_alternative-TestMonteCarloHypothesisTest.test_alternativesii##H-   R  (2;? JJqM!!$&&w/ ll5::>>3l6qyJ rzz#--'@A BJJszz$:KrT)rrr`grRrzc ^^[RRS5n[RR USUS9n[R "U[RRTS9nU4SjnUR[RR U5n[XGUSSTS9n[URUR5 TS :Xa*[URURURS 9 gTS :Xa-[S UR- URURS 9 gg) NrrUr`r~rHrcv>[R"U[RRSTS9R$Nasymp)moder)r ks_1samprcdfr)rPrs r> statistic1dGTestMonteCarloHypothesisTest.test_against_ks_1samp..statistic1ds+>>!UZZ^^'.9;;D9 ErTrFr)r rrrrr)r9rIrr skewnormrrrrrrrrrr) rrr`r5rPrrnorm_rvsrbs ` r>test_against_ks_1samp2TestMonteCarloHypothesisTest.test_against_ks_1samps ii##H- NN  #  >>>!UZZ^^M E<< 4qK+/E+68  x'9'9: ) # CJJdii H F " AcjjL(// J#rThypotestrrrr=rc^[RRS5n[RR USUS9nT"XRS9nU4SjnUR [RR U5n[XXUSUS9n [U RUR5 [U RURURS9 g) Ni rrc$>T"XS9R$rm)r)rPr$rs r>rLTestMonteCarloHypothesisTest.test_against_normality_tests..statisticsA)33 3rTTr rr) r9rIrr rrrrrrrrr) rrrr`r5rPrrrrbs ` r>test_against_normality_tests9TestMonteCarloHypothesisTest.test_against_normality_testss ii##H- NN  3  ?A7 4<< 4qI$+68  x'9'9: HOO$))DrTc[RRS5n[RR USUS9n[R "U5nSnUR[RR U5n[X6USSS9n[URUR5 [URURURS9 g) Ni!Mr rc>[R"XS9R$rm)r normaltestrrs r>rGTestMonteCarloHypothesisTest.test_against_normaltest..statistics##A1;; ;rTTrr r)r9rIrr rrrrrrrrrr)rr`r5rPrrrrbs r>test_against_normaltest4TestMonteCarloHypothesisTest.test_against_normaltestsii##H- NN  3  ?##A& <<< 4qI$+46  x'9'9: HOO$))DrTc [RRS5n[RR USUS9n[R "U[RR5nSnUR[RR U5n[X6USSSS9n[URUR5 [URURURS 9 g) Ni rUrct[R"U[RR5R$r])r cramervonmisesrrrrs r>rMTestMonteCarloHypothesisTest.test_against_cramervonmises..statistic1ds#''5::>>:DD DrTrFrrr)r9rIrr rrrrrrrrrrr)rr`r5rPrrrrbs r>test_against_cramervonmises8TestMonteCarloHypothesisTest.test_against_cramervonmisessii##I. NN  #  >''5::>>: E<< 4qK+/E+46  x'9'9: HOO$))DrT dist_name)rlogisticr_c $^^UU4Sjn[5nUR[5 [USS9nSSS5 WR(deUR Sn[ RRS5n[RRUSUS9n[R"UT5n U Rn U RTS- n U4Sjn UR[![T5RU5n [5nUR[5 [#XU SS S S 9nSSS5 [%WRU 5 [%UR&U S UR(-S 9 g!,(df  GNB=f!,(df  NY=f)Nc>[RRS5n[RR USUS9n[R "UT5nURURT- $)NxrGr$) r9rIrr tukeylambdarandersonrcritical_values)r`r5rPrrr_s r>fun?TestMonteCarloHypothesisTest.test_against_anderson..funs`))'' 2C!!%%ac%DA~~a3H%%(@(@(CC CrTr)x0r rGr$cF>[R"UT5R$r])r r"r)rPrs r>rGTestMonteCarloHypothesisTest.test_against_anderson..statistic1d/s>>!Y/99 9rTrFrrrr)rfilterRuntimeWarningrsuccessrPr9rIrr r!rr"rsignificance_levelrrJrrrr)rrr_r$supsolr`r5rPr expected_stat expected_prdist_rvsrbs `` r>test_against_anderson2TestMonteCarloHypothesisTest.test_against_andersons; D  C JJ~ &sq/C!{{{ EE!Hii##I.    ! !!#C ! @>>!Y/ ** 003C7  :<<y 9 = =sC  C JJ~ &"1#.D.3LC!  }5 JQtyy[A5! &! s E/ #F/ E> Fc[RRS5n[R"S5n[ X!R[R SSS9nUR S:Xdeg)NlIL~QH(srGTrr g-C6?)r9rIrrrr:r)rr5rPrbs r>test_p_never_zero.TestMonteCarloHypothesisTest.test_p_never_zero<sTii##$78 HHSMq**bgg*.FDzzV###rTc[RRS5nURSS9URSS94nURUR4nSn[R "X#USS9n[R "USUS /SS9n[URUR5 [URURS S 9 g) NF:fmPrrc@[R"XUS9R$rmrrXs r>rFTestMonteCarloHypothesisTest.test_against_ttest_ind..statisticJs??1d3== =rTrfr#rr{Gz?r) r9rIrnormalr rrrrrrr5rrrrbr(s r>test_against_ttest_ind3TestMonteCarloHypothesisTest.test_against_ttest_indEsii##$67zzvz&  (::jj#**$ >$$T Cood1gQyr: s}}5 CJJT:rTc[RRS5nURSS9URSS9URSS9URSS94nURURURUR4nSn[R "X#USS S 9n[R "US S06n[URUR5 [URURS S 9 g)Nr8rvr)re)rf)rgcB[R"USU06R$)Nr$)r f_onewayr)r$rs r>rETestMonteCarloHypothesisTest.test_against_f_oneway..statisticYs>>43d3== =rTrfrr$rr${Gz?r) r9rIrr=r rrFrrrr>s r>test_against_f_oneway2TestMonteCarloHypothesisTest.test_against_f_onewayRsii##$67  )3::8:+D  )3::8:+DFjj#**cjj#**< >$$T 1:<nnd,, s}}5 CJJT:rTz2Statistic may not depend on sample order on 32-bitc^[RRS5mSnU4SjnU"S5n[R"X2[R SUS9n[R "URUR:*5n[R "URURS-:*5nXV:wdeURUS-US-- :Xdeg) NlU5S[R"SS9RUTS9-$)NrgZd;O?)r]r$)r bernoullirr~r5s r>rITestMonteCarloHypothesisTest.test_finite_precision_statistic..rvsks'%044$S4QQ QrTrGr)rr)g?r) r9rIrr rr~rnull_distributionrr)rr)rrPrbc0c1r5s @r>test_finite_precision_statisticVVC))S]]: ; VVC))S]]G-DD Exxzzb1f{Q7777rTr^)%__name__ __module__ __qualname____firstlineno__rr9rrrrrr6markxslowr parametrizer@rrslowlinspacerr skewtest kurtosistestrrgrrr2r5r?rJ fail_slowxfail_on_32bitrT__static_attributes__r^rTr>rrs D(,=@.@.D P P [[$2$2L [[VU2q\2@3@B [[],LMLNL( [[ [[],?@ [[S"++dC";<K=AK, [[Z%..%:L:L)MN [[],LM [[S"++b!Q"78E9NOE" [[S"))B"23E4E" [[ [[S"++dC";<E=E$ [[ [[[*>? [[S%(+(B,@(BT$ ;;" [[1 [[ TU8V8rTrc\rSrSrSr\R RS5r\R RS5r Sr Sr Sr g) TestPowerixc [RRS5n[RnUR UR 4nSnSn[ R"[US9 [X#USS9 SSS5 Sn[ R"[US9 [USU5 SSS5 [ R"[US9 [X!R S4U5 SSS5 S n[ R"[US9 [X!R 4U5 SSS5 [ R"[US9 [X#S 5 SSS5 S n[ R"[US9 [X#US S 9 SSS5 [ R"[US9 [X#U[R"SS5S 9 SSS5 Sn[ R"[US9 [X#USS9 SSS5 Sn[ R"[US9 [X#SS//SQ45 SSS5 [ R"[US9 [X#SS/SS/4S/SQ0S9 SSS5 Sn[ R"[US9 [SX45 SSS5 Sn[ R"[US9 [X#USS9 SSS5 [ R"[US9 [X#USS9 SSS5 Sn[ R"[US9 [X#USS 9 SSS5 [ R"[US9 [X#USS 9 SSS5 g!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GNi=f!,(df  GNM=f!,(df  GN0=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GN|=f!,(df  GN_=f!,(df  GND=f!,(df  GN'=f!,(df  g=f)!NlyB01)r0rFrrr"rrrr)r0z3`significance` must contain floats between 0 and 1.r significancerfrz`kwargs` must be a dictionaryr)kwargsz+shape mismatch: objects cannot be broadcastr0rE)rF rFrjrPz`test` must be callablez(`n_resamples` must be a positive integerir(g%@z"`batch` must be a positive integerr,) r9rIrr rr=r6r7r8rr<r^)rr5testrn_observationsr=s r>rTestPower.test_input_validationysii##$78zz3::&!D ]]:W 5 $^ <6E ]]9G 4 $n -5 ]]9G 4 $V,n =5. ]]:W 5 $ ~ 66 ]]:W 5 $U #6H ]]:W 5 $^! <6 ]]:W 5 $^"++b!:L M62 ]]9G 4 $^I >5@ ]]:W 5 $r2h 5 66 ]]:W 5 $r2hR13 :J K6, ]]9G 4 $ ,5= ]]:W 5 $^ =6 ]]:W 5 $^ >67 ]]:W 5 $^3 76 ]]:W 5 $^4 86 5Y6 55 4 4 46 5 5 56 5 5 55 46 5 5 55 46 5 5 56 5 5 5s% MMM'?M98 N ( N!N/ OO<O%5 O7% P  P P-. P? M M$' M69 N N N,/ N> O O"% O47 P P P*- P<? Q c|^[RRS5nU4SjmSTlSTl[ TSSS9n[RRS5n[ S 0UDURSS.D6n[TRS5 [TRS5 [RRS5nSTl[ S 0UDURS S.D6n[TRS 5 [TRS 5 [RRS5nSTl[ S 0UDURSS.D6n[TRS5 [TRS5 [UR UR 5 [UR UR 5 g) Nrc>URS:XaSO [U5n[UTR5TlT=RS- sl[ R "USUS9R$r)r rhr,rrr rr)rPr$rrls r>rl"TestPower.test_batch..testsSffks1vJ!*doo>DO LLA L$$Q5<< rTestPower.test_batchsFii##K0 =  bdCii##K05t515T\\4(T__a(ii##K0 6t626T\\2&T__b)ii##K0 8t848T\\1%T__d+TZZ,TZZ,rTc ^ ^ ^[RRS5n[R"SS/5mSU4Sjjn[R"/SQ5SS2[R 4nS/SQ0n[R"/SQ5n[ R"X!RX5US 9n/nUHenUH\n USHPm THGm SU U 4S jjn [ R"XRU US 9n URU R5 MI MR M^ Mg [R"XvRR5n[URUS S S9 g)Nlx7'\r皙?c >[R"T[[SURS-555n[ R "XUUS9$)Nr)rr$)r9 expand_dimsr;r@r r r)rPrr$popmeans_expandedpopmeanss r>rl*TestPower.test_vectorization..testsC "xuQQR ?S9T U $$Q{*.0 0rT)r0rrTrrGrr)rIr皙?rrhric2>[R"UTUTS9$)NrH)r r)rPr$ alternative_i popmean_is r>test2+TestPower.test_vectorization..test2s!#(#4#4Q AN$PPrTrgr<rI)rrre) r9rIrrrqrr rr=appendrsrqr)rr5rlnxrirhrbr(significance_inx_irtmpr}r~rxs @@@r>test_vectorizationTestPower.test_vectorizationsii##N388QH% 0ZZ- .q"**} =!ABzz":; kk$ B!')*N%+M%:M%- PP$kk%T7EG 399-&.&;+jjiioo.  3T=rTc[RRS5n[RnUR SSSS9nUR UR 4n[R"/SQ5n[R"X$X5S9n[R"USS2[R4URR5n[URUSS 9 g) N&(r0rG)rr0r)rIrzrrgrIr) r9rIrr rintegersr=rqrrrrrqr)rr5rlrmrrhrbs r>test_ttest_ind_nullTestPower.test_ttest_ind_nullsii##O4b#G<jj#**$zz"34 kk$^O|ArzzM'BCIIOOT  test_ttest_1samp_power TestPower.test_ttest_1samp_power sii##O45464467 _-$d| kk%++SZZ'3D 3T2rTr^N) rVrWrXrYrr6rZr]rrrrrcr^rTr>rerexsP69p [[!-!-F [[%>%>N <3rTrec \rSrSrSrSrSr\RRS\ RR\ RR/5\RRS/SQ5S55r\RRS\ RR\ RR/5\RRS S \R "S 5S -4S S\R""SS 54/5S55rSr\RR)5S5rSr\RRSS5S5r\RRSS5S5r\RRSS5S5rSr\RR75\RRSS5S55r\RR(\RRSS5S55r\RRSS5S5rSr\RRSS5S 5r \RR75\RRSS!5S"55r!/S#Q/S$QS%S&S'S(S)S*S+.r"/S,Q/S-QS.S/S0S1S2S3S+.r#\RR75\RRS4\"\#45S555r$\RR(\RRS6S75S855r%\RRS9S:5S;5r&\'"S<5S /4\'"S5S /S=Q/S>Q/4\'"S?5S /S=Q/S>QSS@//4/r(\RRSA\(5SB5r)\RRUS 5SC5r+SDr,gE)FTestPermutationTesti r1cL[RRS5Ulg)Nl F)r9rIrr5)rs r> setup_method TestPermutationTest.setup_method$s99(()<=rTcSnSn[R"[US9 [/SQS/4U5 SSS5 Sn[R"[US9 [SU5 SSS5 [R"[US9 [SU5 SSS5 Sn[R"[US9 [/SQ/SQ4US S 9 SSS5 S n[R"[US9 [/SQ/SQ4US S 9 SSS5 Sn[R"[US9 [/SQ/SQ4US S9 SSS5 Sn[R"[US9 [/SQ/SQ4USS9 SSS5 Sn[R"[US9 [/SQ/SQ4USS9 SSS5 Sn[R"[US9 [/SQ/SQ4USS9 SSS5 Sn[R"[US9 [/SQ/SQ4USS9 SSS5 Sn[R"[US9 [/SQ/SQ4US S9 SSS5 Sn[R"[US9 [/SQ/SQ4USS9 SSS5 g!,(df  GNK=f!,(df  GN-=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GN=f!,(df  GNb=f!,(df  GN?=f!,(df  GN=f!,(df  N=f!,(df  g=f)NcF[R"X4U5R$r]rrXs r>r:TestPermutationTest.test_permutation_test_iv..stat+s??A640:: :rTz2each sample in `data` must contain two or more ...rrrz6`data` must be a tuple containing at least two samples)rr!r"r#z `permutation_type` must be in...r)permutation_typerrrr'r(r*r+r,rrr2r3r4)r6r7r8rr<)rrr=s r>test_permutation_test_iv,TestPermutationTest.test_permutation_test_iv)sm ;G ]]:W 5 i!-t 46K ]]:W 5 T4 (6 ]]9G 4 Q %5/ ]]:W 5 i3T D65 ]]:W 5 i3T.4 66E ]]:W 5 i3Tc J6> ]]:W 5 i3Tu M6> ]]:W 5 i3Tv N6@ ]]:W 5 i3T G6@ ]]:W 5 i3T H60 ]]:W 5 i3Tv N6A ]]9G 4 i3Ty I5 4W6 56 5 4 46 56 5 6 56 56 56 56 56 55 4sI5 J J2J+'J=KK!K3;L0L%L)L:5 J J J(+ J:= K  K! K03 L L L&) L7: MrHr)pairingssamples independentc^ URRS5nURRS5nU 4Sjm ST lST lSUSS.n[R "X44T 4SU"S5S.UD6n[ T RS 5 [ T RS5 ST l[R "X44T 4S U"S5S.UD6n[ T RS 5 [ T RS 5 ST l[R "X44T 4SU"S5S.UD6n[ T RS 5 [ T RS5 [ URUR5 [ URUR5 g) Nr0c>URS:XaSO [U5n[UTR5TlT=RS- sl[ R "XS9[ R "XS9- $Nrr#)r rhr,rrr9r:)rPrYr$rrs r>r1TestPermutationTest.test_batch..statisticis[ffks1vJ#&z93G3G#HI    " 771(2771+@@ @rTrrT)r)rr r)r-rHrrrr)r5rIrrr rrr) rrrHrPrYrrQrRrrs @r>rTestPermutationTest.test_batch_sr HHOOB  HHOOB  A   #9I"$%%qfiLq3??LFJLY&&-Y))1- %%qfiLr3??LFJLY&&+Y))2. %%qfiLt3??LFJLY&&*Y))40T[[$++.T[[$++.rTzpermutation_type, exact_sizerrr)rrrr<cvURRS5nURRS5nSnUSS.n[R"XE4U4SU"S5S.UD6n[ UR R S5 [R"XE4U40UD6n[ UR R U5 g)NrcT[R"XS9[R"XS9- $rmrrXs r>r8TestPermutationTest.test_permutations..statistic771(2771+@@ @rTT)rr r)r)rH)r5rIr rrrQr~) rr exact_sizerHrPrYrrrbs r>test_permutations%TestPermutationTest.test_permutationss HHOOA  HHOOA  A%5"$$$aVYKA2>q/KEIKS**//3$$aVY?$?S**//rNTestPermutationTest.test_randomized_test_against_exact_both..statisticrrTTrrGr rr-rr5r)rIr) r binomr rrrr9infrrr rrr5rny permutationsrPrYrrrrbrRs r>'test_randomized_test_against_exact_both;TestPermutationTest.test_randomized_test_against_exact_boths% * }}RWb)L888 JJNNN # JJNNN #t A# [Ft'L'!%'LRVVLtL}}... DKKd;rTcSupSup4nSU-U:de[RRUS9n[RRUS9nXg4nSn SSSXS .n [X4S U0U D6n [X4S [R 0U D6n U R U R :Xde[U RU RS S 9 g) N)rN)rri}rrc0[R"X- US9$rmrrXs r>rQTestPermutationTest.test_randomized_test_against_exact_samples..statistics7715t, ,rTTrrGrr)rIr) r rrrr9rrrrrs r>*test_randomized_test_against_exact_samples>TestPermutationTest.test_randomized_test_against_exact_sampless + , "u|### JJNNN # JJNNN #t -# [Ft'L'!%'LRVVLtL}}... DKKd;rTc^ SURp!Sup4n[R"U5U:de[RR US9n[RR US9m U/nU 4Sjn[ R"U5n SSSXS.n [Xy4S U0U D6n [Xy4S [R0U D6n U RU R:Xde[U RU RS S 9 g) Nr)rri@rc8>[R"UT5S$Nrr pearsonrrPrYs r>rTTestPermutationTest.test_randomized_test_against_exact_pairings..statistic1ds>>!Q'* *rTTrrGrr)rIr)r5r factorialr rrrrkrr9rrrr)rrr5rrrrPrrrrrbrRrYs @r>+test_randomized_test_against_exact_pairings?TestPermutationTest.test_randomized_test_against_exact_pairingss'S*   $|333 JJNNN # JJNNN #s + 44[A " [Ft'L'!%'LRVVLtL}}... DKKd;rTrrc^URRSSS9nURRSSSS9n[R"X#TSS9nU4S jn[ X#4U[ R S URS 9n[URURURS 9 [URURURS 9 g) Nrrr~r}rzr)r~r|r}exact)rrcD>[R"XSTS9R$r)r ks_2samprrPrYrs r>r>TestPermutationTest.test_against_ks_2samp..statistic1ds!>>!W.9;;D9 ErTrr)rr5r) r5r=r rrr9rrrrr)rrrPrYrrrbs ` r>test_against_ks_2samp)TestPermutationTest.test_against_ks_2samps HHOO!O , HHOOO 3>>!KgN E v{+4$((D  x'9'9 J HOO$))DrTcURRSSS9nURRSSS9nSSSS .nXAn[R"X#US 9nS n[ X#4U[ R XRS 9n[URURURS 9 [URURURS 9 g)NrrrrzrrrrrrcB[R"X5R$r]r ansarirrs r>r.statistic1d<<%// /rTrr) r5r=r rrr9rrrrr) rrrPrYalternative_correspondencealternative_scipyrrrbs r>test_against_ansari'TestPermutationTest.test_against_ansari s HHOO!O , HHOO!O ,/8173>&@"7C<<2CD 0v{+6HHF  x'9'9 J HOO$))DrTc [RRSSURS9n[RRSSURS9n[R"X#SUS9nSn[ X#4US[ RUSURS 9n[URURURS 9 [URURURS 9 g) N)rrzrrr~r|rHrzrrHc@[R"XUS9R$rm)r mannwhitneyurrXs r>r@TestPermutationTest.test_against_mannwhitneyu..statistic*s%%a6@@ @rTT)r r)rr$r5r) r uniformrr5rrr9rrrrr)rrrPrYrrrbs r>test_against_mannwhitneyu-TestPermutationTest.test_against_mannwhitneyu"s MM  9!$((  K MM  9$TXX  N%%a L AvyT+-66{$%4885  x'9'9 J HOO$))DrTc[RRSSURS9n[RRSSSURS9n[R"XSS9nS n[ X4U[ RS URS 9n[URURURS 9 [URURURS 9 g) Nrr)r~r}rHrzr)r~r|r}rHrr.c@[R"XSS9R$N asymptoticr.)r cramervonmises_2samprrs r>r9TestPermutationTest.test_against_cvm..statistic1d;s!--a5ACCL9 MrTrrr) r rrr5rrr9rrrrr)rrPrYrrrbs r>test_against_cvm$TestPermutationTest.test_against_cvm4s JJNNN B JJNNqN I--a7C M v{+4$((D  x'9'9 J HOO$))DrTr$)rfrc[RRS5nURSS9nURSS9nURSS9nX4U4nSnSn[R"U5n [R"U5n [R "US5n [R "US 5n [R "US 5n U "XXS 9nU "XXS 9nS US SUR S.n[Xg4S[R0UD6n[Xg4SS0UD6n[URXRS9 [URURURS9 [URUSS9 [URURSS9 g)NlY}y!dhrr)rrr)rrrc<[R"U6R$r])r kruskalrrs r>rITestPermutationTest.test_vectorized_nsamp_ptype_both..statistic1dXs==$'11 1rTc<[R"U6R$r])r rrrs r>pvalue1dFTestPermutationTest.test_vectorized_nsamp_ptype_both..pvalue1d[s==$'.. .rT)rrr)rrrrr#Frrr r$rrr5r)rrgQ?rr) r9rIrrrkrrr5rrrrrr)rr$r5rPrYrnrrrrrx2rEz2expected_statisticexpected_pvaluerrbrRs r> test_vectorized_nsamp_ptype_both4TestPermutationTest.test_vectorized_nsamp_ptype_bothFs`ii##$78 JJQJ  JJIJ & JJIJ &ay 2 / 44[A 11(;__Q * __Q * __Q *&rr= 7$T)$1$((DtMbffMMLtLtL  '9 J t~~DIIF O$? DKKd;rTc^[RRSSURS9n[RRSSURS9nSnSnU4Sjn[R "U5nU"X#SS 9nUSn USn S STS UR[ RS .n [X#- 4U40U D6n [X#4U40U D6n [U RU RURS 9 TS:wa[U RXRS 9 [U RXRS 9 [U RU RURS 9 g)N)rr<rrrrzc@[R"USS9R$Nrrr wilcoxonr)rns r>statistic_1samp_1dETestPermutationTest.test_against_wilcoxon..statistic_1samp_1ds>>!8BB BrTc@[R"XSS9R$rrrs r>statistic_2samp_1dETestPermutationTest.test_against_wilcoxon..statistic_2samp_1ds>>!F;EE ErTc.>[R"XTS9$)Nr)r rrs r>test_1d:TestPermutationTest.test_against_wilcoxon..test_1ds>>!K@ @rTrr#Fr)r r$rrr5r)rr) r rrr5rrkr9rrrrrr)rrrPrYrrrrlrr/r0rrQrRs ` r>test_against_wilcoxon)TestPermutationTest.test_against_wilcoxonvs- MM  9!$((  K MM  9$TXX  N C F A//81%  a[ #Q{$-dhh!vv' (:CdC(:CdC TYYG + % DNNM J Zii@ T[[tyyArTc JURRSSSS9nSX"S:H'S SjnU"U5SSpen[R"XEXaS9n[R"U4US S [ R URUS 9n[URURURS 9 g)Nrrr0rrfc2[R"US:US9$)Nrr#r9rrs r>r=TestPermutationTest.test_against_binomtest..statistics66!a%d+ +rTrRrTr)r rr)r5rrr) r5rr rrr9rrrr) rrrPrkrir]rrbs r>test_against_binomtest*TestPermutationTest.test_against_binomtests HH  a  ,q&  ,A,Ca??1D$$aT96?13TXX1<>  HOO$))DrTc^URRSS9nXRRSS9-m[R"UTSS9nU4Sjn[ U4US[ R URS9n[URURURS9 [URURURS9 g) Nr<rrr.cD>[R"UTSS9R$r)r kendalltaurrs r>r@TestPermutationTest.test_against_kendalltau..statistic1ds##Aq>HH HrTr)rr)r5r) r5r=r rrr9rrrrr)rrPrrrbrYs @r>test_against_kendalltau+TestPermutationTest.test_against_kendalltaus HHOOO # Q' '##Aq9 It[:+-66txxA  x'9'9 J HOO$))DrTc ^U4Sjn[RRS5nURS5S:R[5nURS5SU--S:R[5m[ R RUT5Sn[U4US[RUUS9n[ R"XQS 9n[URUS5 g) NcB>[R"US:HTS:H-5$rr rs r>r@TestPermutationTest.test_against_fisher_exact..statistics6616a1f-. .rTl (yM4r9g333333?g?rr)rr)rr5r) r9rIrastypefloatr contingencycrosstabrr fisher_exactrr) rrrr5rPtabrbrRrYs @r>test_against_fisher_exact-TestPermutationTest.test_against_fisher_exacts /ii##$78 ZZ]S ( ( / ZZ]T!V #c ) 1 1% 8((A.q1tY+-66{#&(!!#? DG,rT)r=rcXURRSS9nURRSS9nURRSS9n[R"X!S9n[R"X1S9n[R"XAS9nUSnX#U4nSnSn[R "U5n[R "U5n U"[ R"U6S U06n U "[ R"U6S U06n S US S SS .n [XV4S[ R0U D6n [XV4SS0U D6n[U RXRS9 [U RURURS9 [U RXRS9 [U RURSS9 g)N)rrrr)rrr)rrrr#rcB[R"USSS9R$)NTrrankedr/)r page_trend_testrrs r>rLTestPermutationTest.test_vectorized_nsamp_ptype_samples..statistic1ds"((d0<>>Gi HrTcB[R"USSS9R$)NTrr!)r r#rrs r>rITestPermutationTest.test_vectorized_nsamp_ptype_samples..pvalue1ds"((d0799? @rTr$Frrrr)r{rrr)r5rIr rankdatarrkr9broadcast_arraysrrrrrr)rr$rPrYrnrrrrrrrrrbrRs r>#test_vectorized_nsamp_ptype_samples7TestPermutationTest.test_vectorized_nsamp_ptype_samplessq HHOOO + HHOOO + HHOOO + NN1 ( NN1 ( NN1 ( aDay H @ 44[A 11(;&(;(;T(BNN ""5"5t"<H4H$T)$.q:tMbffMMLtLtL '9 J t~~DIIF O))D DKKd;rTr)r"rg@rtg?g&_?g@gd~$I"@glu?)rPrY expected_lessexpected_2sidedexpected_Pr_gte_S_meanr expected_avg expected_std)okrGrrCjml) r1r4r2birDnr6rgd?gd?gYo5?g@@@g{C@gDmFA@casecUSnUSnUSnUSnUSnUSnUSnUSn S n [RR5n U R[S 5 [ X#4U [R S S 9n [ X#4U [R S S 9n SSS5 [W RX@RS9 [U RUSS9 [W RUSS9 [U RR5USS9 [U RR5U SS9 U RnU RR5n[U R5n[R"[R "U RU- 5[R "X- 5:5U- n[UU5 g!,(df  GN8=f)a. Results above from SAS PROC NPAR1WAY, e.g. DATA myData; INPUT X Y; CARDS; 1 1 1 2 1 3 1 4 2 1.5 2 2 2 2.5 ods graphics on; proc npar1way AB data=myData; class X; EXACT; run; ods graphics off; Note: SAS provides Pr >= |S-Mean|, which is different from our definition of a two-sided p-value. rPrYrr+r,r-r.r/cB[R"X5R$r]rrs r>r7TestPermutationTest.test_with_ties..statistic1d9rrTz$Ties preclude use of exact statisticr)r)rrNrg|=rgư>)r9testingrr) UserWarningrrrrrrrQr:rrhrr)rr8rPrYrr+r,r-r.r/rr-rbrRSr:ri Pr_gte_S_means r>test_with_ties"TestPermutationTest.test_with_tiess8 I I!"67_- 01!%&>!?N+ N+  0ZZ ) ) +s JJ{$J K"A6;BFF/57C#QFKRVV0;=D ,  '9 J M> _5A..335|$O..224lN MM$$))+ %% &rvvc&;&;D&@A"$&&. 12234 . >', +s A G Gzalternative, expected_pvalue))rg?)rgOO?)rg[R"UT5R$r])r spearmanrrrs r>rDTestPermutationTest.test_against_spearmanr_in_R..statistic1dks??1a(22 2rTr)rr)rrgvIh%<=rN)rr9rrrrr)rrrrPrrrbrYs @r>test_against_spearmanr_in_R/TestPermutationTest.test_against_spearmanr_in_RQs_( 1 0/ 3t[:+-66{L  '9 J O%@rTr-)rfrc[R"[SS9 [[R "/SQU55 SSS5 g!,(df  g=f)Nz`batch` must be positive.rr)r6r7r8rtr_batch_generator)rr-s r>test_batch_generator_iv+TestPermutationTest.test_batch_generator_ivts3 ]]:-H I --i? @J I Is #A ArrC)rrrzrr9ziterable, batch, expectedcP[[R"X55nXC:Xdegr])rtrrI)riterabler-rgots r>test_batch_generator(TestPermutationTest.test_batch_generator}s#;//@ArTc/SQn/SQnSn[R"X4USSS9nURURURpvnS[ R "XuS- :5-[U5- nXhs=:XaS :Xde eg) N)rrrr)rrr<rc4[R"X5S$rrrs r>rFTestPermutationTest.test_finite_precision_statistic..statistics>>!'* *rTFr)r rrr1gUUUUUU?)r rrrrQr9rrh) rrPrYrrbrrnull correct_ps r>rT3TestPermutationTest.test_finite_precision_statistics   +$$aVY56@B--S5J5J4t5y011CI= )c)))))rTr4N)-rVrWrXrYrrrr6rZr\r9rIrHrrr rrrrr]rrrrrrr[rrr rrr) tie_case_1 tie_case_2r@rFrJr@batch_generator_casesrOrarTrcr^rTr>rr s; D> 2Jl [[^bii.C.C.0ii.C.C.EF [[/CE"/EF"/J [[^bii.C.C.0ii.C.C.EF [[;)7+<+:??:?x [[ [[;>?A ? A< [[Wg.A/A$Ah2.#AhIy+AB#AhIy1a&+IJL [[8244 [[1**rTrc [R"/SQ[S9n[R"U5n[ 5nSn[ R "U5HUnUS- n[R"XASS5nUR[UVs/sHn[U5PM sn55 MW [R"[[U5S- 5Vs/sH3n[R"[!XS5[!XS-S55PM5 sn5n[#X75 [#[U5U5 gs snfs snf)N)rrrrrrrf)r9rintcumsumsetr_all_partitions_concatenatedsplitaddr; frozensetprodr@rhr rrr)rincall_partitionsrpartition_concatenated partitioningr_rs r> test_all_partitions_concatenatedrhs #&A 1BUNG"-"J"J1"M1 xx 63B@ 5 !E 1)A, !EFG#N ww!&s1vax2!0A c!B%j#a!g,?!023H#^$h/ "F2s D8:D=fun_namerrrcR^^[RRS5mTRSS9mU4SjnT4TSS.T4TSS.TUS .S .nS S 0n[[U5nX nUR U5 U4S jnU"SUSS.UD6 U"SUSS.UD6 SnU"SUSS.UD6 U"SUSS.UD6 g)Nlr:r0rc@>[RRUTS9$)Nr$)r rrrOs r>r&test_parameter_vectorized..rvsszz~~4c~::rTrA)rr5r/r)rr5r)rrrjr)rGc>URS:d[R"UT5(de[R"XS9$r)r r9rBr:)rPr$rs r>r,test_parameter_vectorized..statistics3vvzR^^Av6666wwq$$rT)rr TcRURS:Xde[R"U5$r)r r9r:rs r>rros vv{{wwqzrTFr^)r9rIrrJr r) rir fun_optionscommon_optionsr$optionsrr5rs @@r>test_parameter_vectorizedrts ))   /C ZZRZ F;+13+7!917 #"%8)888)888)889)99rTc\rSrSrSrSrg)TestMonteCarloMethodicSn[RRS5n[R"[ US9 [ R"URUS9 SSS5 g!,(df  g=f)Nz.Use of `rvs` and `rng` are mutually exclusive.ir)rr5)r9rIrr6r7r8r MonteCarloMethod)rr=r5s r>test_rvs_and_random_state.TestMonteCarloMethod.test_rvs_and_random_statesIBii##H- ]]:W 5  " "szzs ;6 5 5s A$$ A2r^N)rVrWrXrYryrcr^rTr>rvrvsrs JJ/);H LL--<8~#AB+;,C;"#AB>C>2#AB+D%=1(:2,C(:VAB#ABOCCO0#'  " +W]]_=C>C.(,,0%)+ 7 t12+399;="="J(+,/1+399;=+=+\G\#:;!Q(>)<>.EF#AB!Q(>)CG>,#AB (C (#AB6C6"2&/!T; < 2;<#AB;C=;.!&m]%CD-E-2+.+ ",/ "."+ , 4#AB>C>>45!>6!>HA.}8}8@ e3e3Ps *s *x0*NP:P:B<