e!hA ddlZddlZddlmZmZmZddlmZddl m Z ddl m Z m Z ddlmZmZddlmZmZddlmZdd lmZmZmZmZddlmcmZd Zej<j?d gd ej<j?d gddZ ej<j?d gd dZ!ej<j?d gd ej<j?d gdej<j?dddgdZ"ej<jFej<jIdej<j?d gd dZ%dddd Z&ej<j?de&jOdZ(dddd Z)dZ*dZ+ej<jFej<j?de)jOd Z,d!d"d#Z-ej<jFej<j?de-jOd$Z.ej<j?d d%d&gej<j?d dd'gd(Z/ej<jId)ej<j?d gd ej<j?d dd'gd*Z0ej<j?d gd d+Z1ej<j?d gd d,Z2d-Z3ej<j?d.dd/gd0Z4ej<jId1ej<j?d gd d2Z5d3Z6ej<j?d4d5d6gd7Z7ej<j?d8gd9ej<j?d gdd:Z8d;Z9ej<j?d gdd<Z:ej<jFej<j?d gd d=Z;ej<jFej<jyd>d?Z=d@Z>GdAdBZ?GdCdDZ@GdEdFZAdGZBej<j?dHgdIdJZCGdKdLZDy)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 d}tjt|5tdtj dddd}tjt|5tt tj dddd}tjt|5tgddgftj dddd}tjt|5tgdgdftj d dddd }tjt|5tdtj d dddd}tjt|5tgdftj ddddd}tjt|5tgdftj ddddd}tjt|5tgdftj ddddd}tjt|5tgdftj ddddd}tjt|5tgdftj ddddd}tjt|5tgdftj ddddd}tjt|5tgdftj d dddd}tjt|5tgdftj ddddd}tjt|5tgdftj d dddd!}tjt|5tgdftj d"#dddy#1swYxYw#1swYxYw#1swYaxYw#1swY&xYw#1swYxYw#1swYxYw#1swYxYw#1swYIxYw#1swYxYw#1swYxYw#1swYxYw#1swYixYw#1swY1xYw#1swYxYw#1swYyxYw)$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 i/var/www/html/diagnosisapp-backend/venv/lib/python3.12/site-packages/scipy/stats/tests/test_resampling.pytest_bootstrap_ivr?s5G z 1!RWW9G z 1$%'277#$OG z 1-9qc"BGG,-OG z 1C9l+RWWTBCAG z 11!RWW01+G z 139,c232G z 1@9,$?@>G z 1<9,U;<>G z 1=9,V<=9,"=>FG z 189,Q78=G y 089,Y788q$$--CC1133@@<<==667788>>8888sP #P& P%+#P23P?5 Q : Q? Q& Q3 R R  R R' R4" S PP"%P/2P<?Q  QQ#&Q03Q=R  RR$'R14R=S r/)basic percentileBCar$rrrc Ttjjdtjjddd}t |ftj d|d|d}tjjdt |ftj d||d}t |jj|jjt |jj|jjt |j|jy)Nrr0 d)r-r/ random_stater$r))r-r/r$r)) r9randomseedrandrr:rconfidence_intervallowhighstandard_error)r/r$xres1res2s r>test_bootstrap_batchrSPsIINN1 r2r"A aT277$v"#$C ADIINN1 aT277"VC 1D))--t/G/G/K/KL))..0H0H0M0MN$$d&9&9:ctjjdd}tjj|tjj|dddfd }tjt }t |f|d}t fdd}t|j|jt|j|jy) NrrGc2||z dzj|S)Nrr#r:rPyr$s r> my_statisticz+test_bootstrap_paired..my_statisticks1qD))rTc.|}|}||}|SN)ir$abresrZrPrYs r>my_paired_statisticz2test_bootstrap_paired..my_paired_statisticns& aD aD1a  rTr4T)rr5) r9rIrJrKarangelenrrrLrO) r/nrbr^rQrRrZrPrYs @@@r>test_bootstrap_pairedrhcsIINN1 A qA qA* #a&A aT.A 6D aV\$A >DD,,d.F.FGD'')<)<=rTrTFc tjjddd}d}||}tjj|}tjj|}tjj|}t |||f|||ddd} | j j | jj dzk(sJgd} || |<tj|j| |}tj|j| |}tj|j| |}t |||f|||d|d} t| jj| jjt| jj| jjt| j| jt|} | j|t!| jjj | t!| jjj | t!| jj | y) Nrcr|j||j|z|j|zSNr#rW)rPrYzr$s r>rZz/test_bootstrap_vectorized..my_statistics3vv4v 166t6#44qvv4v7HHHrT)r0rErFrG)rr/r5r$r))rG)rrrrc)r9rIrJrKrbootstrap_distributionshaperO broadcast_toreshaperrLrMrNlistpopr) r/r$rrZrn n_samplesrPrYrlrQrprR result_shapes r>test_bootstrap_vectorizedru}sIINN1I Ed I y!A y!A y!A aAY VF 5D  ' ' - -""((61 23 2GGDM  '*E2A  '*E2A  '*E2A aAY VF3 8DD,,00,,002D,,11,,113D'')<)<=;LT))--33\B))..44lC$$**L9rTz#MemoryError with BCa observed in CIctjjd}tjj ddd|}d}tj t|dz tj|tj|}|j| }|j}t|ftjd|| }tdi|d |i}t|j|d t|j |d|j#td|tdi||dd}t|jj$|j'|d |j#td|tdi||dd}t|jj(|j'd|z d y)Nl)ICbX r)locscalesizerHffffff?r)dfryrz) confidence)data statisticr)r/r5r&gMb@?rtolga2U0*3?atolrr)r1lessr& alternativegreaterr])r9rI default_rngr normrvstrfr:semintervalstddictrrrLrOupdaterNppfrM) r/r5ralphadistexpected_interval expected_seconfigras r>test_bootstrap_against_theoryrsc ))   3 4C ::>>aqt#> FD E 77c$ikrwwt}EIIdO LD  7((*K w"''tS*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, expectedctjgd}t|ftjd|d}t |j |dy)N)r0rF)@rg+@"2QYyi@Br)r)r/r5g{Gzt?r)r9arrayrr:rrL)r/expectedrPras r>test_bootstrap_against_RrsD" 5 6A QD"''wv CC++XEBrTiic`gd}gd}d}tjjd}tj||f|dd|}tj||f|dd|}tj||f|d d|}tj |j }tj |j }tj |j } d } || z t| z } || z t| z } | | z t| z } | d ksJ| d kDsJt| d ksJy)N)gbF?gp?g|Hܿgxy[e}u׿g*g_?ggi6g2?)gF?gK?g&F}/?g?g(iC?g}$ץ?ghu?gfj}!?chtj||}tj||}||z Srk)r skew)rPrYr$s1s2s r>rz1test_multisample_BCa_against_R..statistics+ ZZ % ZZ %BwrTltg 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_Rrs3  .A 9A  ))   2 3CA '&)s4I//1a&)L(+6Kooq!fi$'S2G 556I''+99:Kggg112GIi'Y7J)+S^;L)#S^3H    $   x=4  rTctjgd}tjgd}dd}||g}tj||}d}d}|j}d}t j ||||||\} } } t| dy) N) r0(.r4h)r&^cc`tj||tj||z Srkr9r:)rlrYr$s r>rz:test_BCa_acceleration_against_reference..statistic,s#wwqt$rwwqt'<<'test_BCa_acceleration_against_referencer"s 78A /0A= q6D //$ *C D E,,K E++D)T5,7@KAq%E/0rTc tjjdd}d}d}tjdd}|j }d}|j ||f }t|ftj ||d |d } | j} tj| d|k|| dkz} | |k(sJtj| ||j} | d kDsJy)NrrG?rwrryrzir{rrdrr&r)r-r/r$皙? r9rIrJr rr:rrrLsum binomtestpvalue) r/rrgr)r&r stat_truen_replicationsrracici_contains_truers r>#test_bootstrap_against_itself_1sampr;sIINN1 AK ::!1 %D IN 88.!,8 -D TG gg%5 +!  C Bvvr!uy0YA5FGH x '' '__-~-//5v C<my_statz4test_bootstrap_against_itself_2samp..my_statws,D)D)u}rTrrrrrrrcr)r/rn1n2r)r&rdist1dist2rrrrrarrrs r>#test_bootstrap_against_itself_2samprgs(IINN1 B BK JJ1A &E JJ3a (E uzz|+IN IINB/I 0E IINB/I 0E UEN%%5 +!  C Bvvr!uy0YA5FGH x '' '__-~-//5v C< zEtest_bootstrap_vectorized_3samp..statistic..86;;t$8rr$rs` r>rz2test_bootstrap_vectorized_3samp..statistic84888rTcD|D]}|jdk(rJ|ddiSNrr$rndimrrrs r> statistic_1dz5test_bootstrap_vectorized_3samp..statistic_1d4 $F;;!# ## $$'Q''rTrrwTrG)r r$r)r/r5Fr9rIrJrKrrrLrO) r/r$rrPrYrlrQrRrs @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 fd dfd}tjjdtjjdd}t |fd|dd|d}t |f|d |dd |d}t |j |j t |j|jy) Nrc&|j|SrkrWrPr$s r>rz2test_bootstrap_vectorized_1samp..statisticsvv4v  rTc:|jdk(sJ|dSNrrr#r)rPrs r>rz5test_bootstrap_vectorized_1samp..statistic_1ds vv{{##rTrrwTrG)r r$r)r-r/r5Fr0rr)r/r$rrPrQrRrs @r>test_bootstrap_vectorized_1sampr s!$ IINN1 q!A aT9D!$D D aTtest_bootstrap_degeneraters =D  [[ * HDCe99E H"''&AS44rvvrvv6FG H H H "''&9S,,.>?##Q'  H H H Hs$(D AC7D7D <DD c Vtjjd}tjdd}|j d|}|f}t |tj|dtjjd}t |tj|dtjjdd }t|j|jt|j|jt|jtjsJy) NlE<1 rrrrGr{rHi[%r/r)r5F)r/r)r5r ) r9rIrr rrrrrrLrO isinstancefloat64)r/r5rrrarefs r>test_bootstrap_gh15678rs ))   0C ::!1 %D 8838 /D 7D D%**V --d3 5C D%**V --d3 GCC++S-D-DEC&&(:(:; c(("** 55 5rTc tjjd}tjdd}|j d|}tj |}|f}t|tj ddtjjd }||jjk(sJttj| tjddtjjd }t|jj |jjt|jj |jjy) Nln!|'x`5rrrrGrrBifr)r9rIrr rrminrrLrMrmaxrrN)r5rrtrue_minrarRs r>test_bootstrap_minr#s ))   0 1C ::!1 %D 8838 /Dvvd|H 7D D"&&C --d3 5C s..22 22 2 bhhtn_bffU..t4 6DS,,000,,113S,,111,,002rTadditional_resamplesrc ztjjd}|jd}d}|}||z}tjjd}tj|ftj ||dd}tj|ftj ||d d | }tjjd}tj|ftj ||d d }t |j|jd t |j|jd y)Nl ]4MrGrrlmdpi4 r|rA)r)r5r&r/rrB)r)r5r&r/r1+=r) r9rIrr rr:rrOrL)r$r5rPrrn3rars r>test_re_bootstrapr(s ))   0 1C  A B B " "B ))   2 3C //1$RS+/  FC //1$RS+/+. 0C ))   2 3C //1$RS+/ ?CC&&(:(:GC++S-D-D rTzSensible to machine precisionctjjd}tjdd}|j d|f}t |tj|d}tjdi|d d i}|jt d | tjdi|d dd}tjdi|d dd}t|jj|jjdt|jj|jjdtj|jjsJtj|jjsJt!j"t$d5tjdi|ddidddy#1swYyxYw)Nl@QixtrrrrGrrd)rrr5r$r&rrrr|rrrr&rz`alternative` must be one ofrrz ekki-ekkir])r9rIrr rrrrrrrrLrNrMisneginfisposinfr6r7r8)r/r5rrrrlgs r>test_bootstrap_alternativer.sm ))   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 11 1 ;;q,,11 22 2 z)G H; :&:k:;;;s .GGcd}tjjdtjj|}t t j |}t|dD]<}|d|ddf}tj||d}tj||rtest_jackknife_resampler;1s EIINN1 A [ , ,Q /0A 59 -Q l99Q+~~c8,,, - [<>"a  rTrng_name RandomStatercttj|d}|tj|d|d}|d}d}d}tjj dtjj |}tj|||}t|D]?}|d|ddf} t|d|d|d} |d| f} tj| | r?Jy)Nz not available.rr0r0r4.rd) getattrr9rIr6skiprJrKr_bootstrap_resampler5rr7) r<r5rng1rng2r)rnrPrYr^r9jsrs r>test_bootstrap_resamplerEDs "))Xt ,C { xj01 q6D q6DK EIINN1 A'';DAA ; -Q l $59eBi 8S"W:~~c8,,,-rTscore)r?rcd}tjjdtjj|}t j ||d}d}|||ddz }t ||dy)N)r0rrdr#cNtjtj|||Sr\)r9apply_along_axisr percentileofscore)r_rFr$s r>vectorized_posz0test_percentile_of_score..vectorized_posfs""5#:#:D!UKKrTrGV瞯<)r9rIrJrKr_percentile_of_scorer)rFr$rnrPprNp2s r>test_percentile_of_scorerS^sf EIINN1 A((E;AL 5r *3 .BAr5!rTctd}tjjdtjj|}tjj|dddz}t j ||}t |dD]2}||}tj||||d}t||d4y)N)r0rIrrdrGr#rO) r9rIrJrKr_percentile_along_axisr5rAr)rnrPqrYr^rars r>test_percentile_along_axisrWns EIINN1 A cr #c)A**1a0A 58_.d==1qt"5Xu-.rTcdfd}tj|}tjj dtjj ddd}tjj ddd}tjj ddd}||||}|||||}t ||y) Nc,tfd|DS)Nc3@K|]}|jywr\rWrs r>rz>test_vectorize_statistic..statistic..rrrrs` r>rz+test_vectorize_statistic..statisticrrTcD|D]}|jdk(rJ|ddiSrrrs r>rz.test_vectorize_statistic..statistic_1drrTrrrwr1rr#)r_vectorize_statisticr9rIrJrKr) r$r statistic2rPrYrlrQrRrs @r>test_vectorize_statisticr^s9(11,?JIINN1 q!QA q!QA q!QA Q14 (D aAD )DD$rTctjjd}d}tjj |d|d}d}t |f||ddd }tj|jjj|k|jjj|kDzd }tj|d k\sJtj|d ksJ|jjjdk(sJ|jjjdk(sJ|jjdk(sJ|jjdk(sJy)Nls:)rrG)rGrGrctjtj||tj||dgS)Nr)ddof)r9asarrayr:r)rr$s r>rz/test_vector_valued_statistic..statistics5zz2774.66$1578 8rTrd')r/r$r)r-rr#ZrGrrG)rrGrc)r9rIrr rrrrrLrMTrNallrnrOrm)r/r5paramsrrracountss r>test_vector_valued_statisticrksA ))   +C F ZZ^^V*3 GF8 VIyb $C 1CVVS,,0022V;..3355>@F 66&B,   66&C-   " " & & , , 88 8  " " ' ' - - 99 9    # #x // /  % % + +~ == =rTzignore::RuntimeWarningctjjd}dfdfd}gdgdgdgdgd g}tj|j}t |||d }t ||d }t |jjd |jjd t |jjd |jjd y)Nl\Qc|j|}|j|}||dz ||dz zj|}d|z|j||j|z||z dzzz S)N).Nr)r:var)rPrYr$xmymcovs r> concordancez9test_vector_valued_statistic_gh17715..concordancest VVD\ VVD\ByM!a"Y-&78>>tDCAEE$K!%%+5bQFGGrTcR||z}||z}tj|||Sr\)r9 nan_to_num)tptnfpfnr$actualrrrs r>rz7test_vector_valued_statistic_gh17715..statistics-b7}}[4@AArTc>g||tjdfS)N.)r9newaxis)r$argsrs r>statistic_extradimz@test_vector_valued_statistic_gh17715..statistic_extradims#%$%%bjj#o66rT)rrrr)rrrr)rr1rr)rr1rr)rrrT)r5rrrOr) r9rIrrrgrrrLrMrN)r5r}rrarrrrs @@r>$test_vector_valued_statistic_gh17715rs ))   2 3CH B 7       D 88D>  D D,#d CC D)T :CC++//2++//e=C++003++00u>rTc*tjjd}|jd}|jd}d}tj||f|tj|j |j f|dd}t jt|5tj||dd d ff|dddt jt|5tj||dd d dff|dddt jt|5tj|j |j d ddd ff|ddddy#1swYxYw#1swYqxYw#1swYyxYw) NiF|)r0r)rErcFtj|||jSrkr ttest_indrrXs r>rz test_gh_20850..statisticsq!$/999rTrr#z-Ignoring the dimension specified by `axis`...rr0r) r9rIrr rrgr6r FutureWarning)r5rPrYrr=s r> test_gh_20850rsO ))   +C 7A 7A: OOQFI& OOQSS!##J 2>G m7 33 Acrc1fI 23 m7 35 Acrc1Q3hK()45 m7 3A acc!A#ss(m,ia@AA 3355AAs$( E1+"E=08F 1E:=F Fc feZdZdZdefdZdZedZedZ e jjedZ ee jjded d d Zee jjd d dZe jj$e jjd de jjdej&ddddZe jjdej,ej.fe jjd de jjdej&ddddZe jjdej2dd dZe jje jjdej&ddddZe jj$e jjdde jjdeddZd Zd!Zd"Ze jjAde jjCd#d$Z"y)%TestMonteCarloHypothesisTest皙?NcfdS)Nc:j|di|S)NrHdtype)rb)r|kwdsrrsrvs_inxps r>z6TestMonteCarloHypothesisTest.get_rvs..s*RZZ0V20VQU0V6;&0&=rTr])selfrrrrs ````r>get_rvsz$TestMonteCarloHypothesisTest.get_rvss = =rTcfd}|S)Ncj||}j||d}|j|}|||z dzz S)Nr#r)r$ correctionrG)r:rnrn)rPr$mvrgrs r>rz=TestMonteCarloHypothesisTest.get_statistic..statisticsH%Aqt2A A!cz> !rTr])rrrs ` r> get_statisticz*TestMonteCarloHypothesisTest.get_statistics " rTc jgd}d#fd }d}jdjdf}tjjtjjf}t j t|5t||dd dddd }t j t|5t|tjj|d dddd }t j t|5t|tjj|d dddd}t j t|5t|d|dddt j t|5jddgddgg}t|ddg|dddd}t j t|5jgdg}t|ddg|dddd}t j t|5t|tjjddddd}t j t|5t|tjj|ddddd}t j t|5t|tjj|ddddd}t j t|5t|tjj|ddddd}t j t|5t|tjj|ddddd}t j t|5t|tjj|d dddd!}fd"} t||jy#1swYYxYw#1swYxYw#1swYxYw#1swYxYw#1swYuxYw#1swY7xYw#1swYxYw#1swYxYw#1swYyxYw#1swY9xYw#1swYxYw#1swYxYw#t$r&}t|j|sJYd}~yd}~wwxYw)$N?@@c*j||SrkrW)rPr$rs r>statz@TestMonteCarloHypothesisTest.test_input_validation..stat s77147( (rTz/Array shapes are incompatible for broadcasting.rrw)rrwrcyNrr]rXs r>rzDTestMonteCarloHypothesisTest.test_input_validation..srTrdr#r!r"rr0`rvs` must be callable or sequence of callables.rrrg@c|Sr\r]rPs r>rzDTestMonteCarloHypothesisTest.test_input_validation..#arTIf `rvs` is a sequence...c|Sr\r]rs r>rzDTestMonteCarloHypothesisTest.test_input_validation..(rrTc|Sr\r]rs r>rzDTestMonteCarloHypothesisTest.test_input_validation..(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 statisticcbjtjj|S)Nr)rbr rr)r{rs r>rz?TestMonteCarloHypothesisTest.test_input_validation..rvsDs!::ejjnn$n78 8rTr\) rbzerosr rrr6r7r8rr<r:str startswith)rrrrr=tempres ` r>test_input_validationz2TestMonteCarloHypothesisTest.test_input_validationszz,' )D "((6"23zz~~uzz~~. ]]:W 5 G T3(<2 F G/ ]]:W 5 C T5::>>4c B CE ]]:W 5 I T5::>>4C H IE ]]9G 4 / T4 . / ]]9G 4 >::Bx"b23D TK#6 = >. ]]:W 5 E::|n-D TK#=t D E2 ]]9G 4 9 T5::>>4 8 9> ]]:W 5 L T5::>>4U K L> ]]:W 5 M T5::>>4V L M@ ]]:W 5 F T5::>>4u E F@ ]]:W 5 G T5::>>4v F G0 ]]:W 5 M T5::>>4V L M6 9 . T3 0i G G C C I I / / > >  E E  9 9 L L M M F F G G M M .q6$$W- -- .sO6(P(PP(P*&P7'&Q2(Q?(Q (Q+(Q8&(RR6PP PP'*P47QQQQ(+Q58RR  R?R::R?cfd}d}jgd}tr&t|tjj |yt jt|5t|tjj |dddt jt|5t|tjj jddddy#1swY_xYw#1swYyxYw)Nc&j|Sr\rW)rPrs r>non_vectorized_statisticzWTestMonteCarloHypothesisTest.test_input_validation_xp..non_vectorized_statisticMs771: rTz!`statistic` must be vectorized...rrFr) rbr rr rrr6r7r8r:)rrrr=rs ` r>test_input_validation_xpz5TestMonteCarloHypothesisTest.test_input_validation_xpKs 6L) B< VUZZ^^5M N  ]]:W 5 O VUZZ^^5M N O ]]:W 5 P VUZZ^^RWW O P P O O P Ps(&C,12C8,C58Dc tjjd}|j|j d}t |  fdd_d_|ddd}jtjjtjjd | |d <tdd d i|}tj dtjd jtjjtjjd | |d <d_tdd di|}tj dtjdjtjjtjjd | |d <d_tdi|}tj dtjdt|j|jt|j|jy)N)pr0rc|jdk(rdn|jd}t|j_xjdz c_j ||Sr )rrnr! batch_sizecounterr)rPr$rrrxp_tests r>rz:TestMonteCarloHypothesisTest.test_batch..statisticds\ffkqwwqzJ#&z93G3G#HI    " .4%%g.qt< test_batchz'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$rctjjd}gd}d||<|jdj}dt |vrdnd}t t|}t ||}|j||} tj| d | } |j| | } t| } |j| } |jtjj||| } t| | | d d|}|j| j | }|j| j"| }t%|j |t%|j"||j&y)Nit$rrrrGr64float32r)r{r)popmeanr$r)rrTi N)r r)r$r)r9rIrrbrrr?rr ttest_1samprrrrrrrrr r)rr$rr5r{ dtype_default dtype_strdtype_nprrPrrrrra ref_statistic ref_pvalues r>test_axis_dtypez,TestMonteCarloHypothesisTest.test_axis_dtypesD ii##G,T  2,, !%]);!;I) 2y)I&   T  :$$Q> JJqJ &!!$&&w/ ll5::>>3elCq#yT+0t= 8#5#5U C ZZuZ=  }5 JTYY?rTr) two-sidedrrc:tjjd}|jd}t j |d|}|j |}t|}|j|}|jtjj||}t||||} t| j|j |jt| j|j |j|j y)N)rJrrrrr)r9rIrrr rrbrrrrrrr rrr) rrrr5rPrrrrras r>test_alternativez-TestMonteCarloHypothesisTest.test_alternativesii##H-   R  (2;? JJqM!!$&&w/ ll5::>>32l6q#ykJ rzz#--'@A BJJszz$:KrT)rrr_grGrwcvtjjd}tjj |d|}tj |tjj}fd}|jtjj |}t|||dd}t|j|jd k(r-t|j|j|j yd k(r0td |jz |j|j yy) NrrJr_r{rHrcztj|tjjdjSNasymp)moder)r ks_1samprcdfr)rPrs r> statistic1dzGTestMonteCarloHypothesisTest.test_against_ks_1samp..statistic1ds+>>!UZZ^^'.9;;D9 ErTrFr)r rrrrr)r9rIrr skewnormrrrrrrrrrr) rrr_r5rPrrnorm_rvsras ` r>test_against_ks_1sampz2TestMonteCarloHypothesisTest.test_against_ks_1samps ii##H- NN  #  >>>!UZZ^^M E<< 4q(K+/E+68  x'9'9: ) # CJJdii H F " AcjjL(// J#rThypotestrrrr2rctjjd}tjj |d|}||}fd}|j tjj |}t|||d|} t| j|jt| j|j|jy) Ni rrc,||jSrk)r)rPr$rs r>rzLTestMonteCarloHypothesisTest.test_against_normality_tests..statisticsAD)33 3rTTr rr) r9rIrr rrrrrrrrr) rrrr_r5rPrrrras ` r>test_against_normality_testsz9TestMonteCarloHypothesisTest.test_against_normality_testss ii##H- NN  3  ?A;7 4<< 4q(I$+68  x'9'9: HOO$))DrTctjjd}tjj |d|}tj |}d}|jtjj |}t|||dd}t|j|jt|j|j|jy) Ni!MrrcDtj||jSrk)r normaltestrr s r>rzGTestMonteCarloHypothesisTest.test_against_normaltest..statistics##AD1;; ;rTTrrr)r9rIrr rrrrrrrrrr)rr_r5rPrrrras r>test_against_normaltestz4TestMonteCarloHypothesisTest.test_against_normaltestsii##H- NN  3  ?##A& <<< 4q(I$+46  x'9'9: HOO$))DrTctjjd}tjj |d|}tj |tjj}d}|jtjj |}t|||ddd}t|j|jt|j|j|j y) Ni rJrcrtj|tjjjSr\)r cramervonmisesrrrrs r>rzMTestMonteCarloHypothesisTest.test_against_cramervonmises..statistic1ds#''5::>>:DD DrTrFrrr)r9rIrr rrrrrrrrrrr)rr_r5rPrrrras r>test_against_cramervonmisesz8TestMonteCarloHypothesisTest.test_against_cramervonmisessii##I. NN  #  >''5::>>: E<< 4q(K+/E+46  x'9'9: HOO$))DrT dist_name)rlogisticr^c fd}t5}|jtt|d}dddjsJ|j d}t jjd}tjj|d|}tj|} | j} | jdz } fd} |jt!tj|} t5}|jtt#|| | dd d }dddt%j| t%|j&| d |j(z y#1swY@xYw#1swYRxYw)Nctjjd}tjj |d|}tj |}|j|jz S)NxrGr) r9rIrr tukeylambdarandersonrcritical_values)r_r5rPrrr^s r>funz?TestMonteCarloHypothesisTest.test_against_anderson..funs`))'' 2C!!%%ac%DA~~a3H%%(@(@(CC CrTr)x0rrGrcDtj|jSr\)r rr)rPrs r>rzGTestMonteCarloHypothesisTest.test_against_anderson..statistic1d/s>>!Y/99 9rTrFrrrr)rfilterRuntimeWarningrsuccessrPr9rIrr rrrrsignificance_levelrr?rrrr)rrr^rsupsolr_r5rPr expected_stat expected_prdist_rvsras `` r>test_against_andersonz2TestMonteCarloHypothesisTest.test_against_andersonsO D   "C JJ~ &sq/C "{{{ EE!Hii##I.    ! !!#C ! @>>!Y/ ** 003C7  :<<y 9 = =sC   LC JJ~ &"1h#.D.3LC L  }5 JQtyy[A5 " "& L Ls#E4 'F4E>F ctjjd}tjd}t ||jtj dd}|j dk(sJy)NlIL~QH(srGTrrg-C6?)r9rIrrrr:r)rr5rPras r>test_p_never_zeroz.TestMonteCarloHypothesisTest.test_p_never_zero<sVii##$78 HHSMq#**bgg*.FDzzV###rTctjjd}|jd|jdf}|j|jf}d}t j |||d}t j |d|d gd}t|j|jt|j|jd y) NF:fmPrrcFtj|||jSrkrrXs r>rzFTestMonteCarloHypothesisTest.test_against_ttest_ind..statisticJs??1ad3== =rTrdr#rr{Gz?r) r9rIrnormalr rrrrrrr5rrrrars r>test_against_ttest_indz3TestMonteCarloHypothesisTest.test_against_ttest_indEsii##$67zzvz&  (::jj#**$ >$$T3 Cood1gQyr: s}}5 CJJT:rTctjjd}|jd|jd|jd|jdf}|j|j|j|jf}d}t j |||dd }t j |d di}t|j|jt|j|jd y)Nr rfr)re)rf)rgc@tj|d|ijS)Nr$)r f_onewayr)r$r|s r>rzETestMonteCarloHypothesisTest.test_against_f_oneway..statisticYs>>43d3== =rTrdrr$rr${Gz?r) r9rIrrr rrrrrrs r>test_against_f_onewayz2TestMonteCarloHypothesisTest.test_against_f_onewayRsii##$67  )3::8:+D  )3::8:+DFjj#**cjj#**< >$$T3 1:<nnd,, s}}5 CJJT:rTz2Statistic may not depend on sample order on 32-bitctjjdd}fd}|d}tj||tj d|}tj |j|jk}tj |j|jdzk}||k7sJ|j|dz|dzz k(sJy) NlU5rzITestMonteCarloHypothesisTest.test_finite_precision_statistic..rvsks'%044$S4QQ QrTrGr)rr)g?r) r9rIrr rrnrnull_distributionrr)rr)rrPrac0c1r5s @r>test_finite_precision_statisticzVVC))S]]: ; VVC))S]]G-DD ERxxzzb1f{Q7777rT)#__name__ __module__ __qualname__rr9rrrrrr6markxslowr parametrizer5rrslowlinspacerr skewtest kurtosistestrrerrrrrr fail_slowxfail_on_32bitr!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 [[ TU8V8rTrceZdZdZej j dZej j dZdZ dZ y) TestPowerc tjjd}tj}|j |j f}d}d}t jt|5t|||ddddd}t jt|5t|d|dddt jt|5t||j df|dddd }t jt|5t||j f|dddt jt|5t||d dddd }t jt|5t|||d dddt jt|5t|||tjdd dddd}t jt|5t|||ddddd}t jt|5t||ddggdfdddt jt|5t||ddgddgfdgdidddd}t jt|5td||dddd}t jt|5t|||ddddt jt|5t|||ddddd}t jt|5t|||d dddt jt|5t|||d dddy#1swYxYw#1swYxYw#1swYxYw#1swYxYw#1swYixYw#1swYAxYw#1swYxYw#1swYxYw#1swYxYw#1swYxYw#1swY]xYw#1swY5xYw#1swYxYw#1swYxYw#1swYyxYw)!NlyB01)r0rFrrr"rrrr)r0z3`significance` must contain floats between 0 and 1.r significancerdrz`kwargs` must be a dictionaryr)kwargsz+shape mismatch: objects cannot be broadcastr0rE)rF rFr4rPz`test` must be callablez(`n_resamples` must be a positive integerir(g%@z"`batch` must be a positive integerr,) r9rIrr rrr6r7r8rr<r))rr5testrn_observationsr=s r>rzTestPower.test_input_validationysii##$78zz3::&!D ]]:W 5 = $^ < =E ]]9G 4 . $n - . ]]9G 4 > $V,n = >. ]]:W 5 7 $ ~ 6 7 ]]:W 5 $ $U # $H ]]:W 5 = $^! < = ]]:W 5 N $^"++b!:L M N2 ]]9G 4 ? $^I > ?@ ]]:W 5 7 $r2h 5 6 7 ]]:W 5 L $r2hR13 :J K L, ]]9G 4 - $^ , -= ]]:W 5 > $^ = > ]]:W 5 ? $^ > ?7 ]]:W 5 8 $^3 7 8 ]]:W 5 9 $^4 8 9 9Y = = . . > > 7 7 $ $ = = N N ? ? 7 7 L L - - > > ? ? 8 8 9 9s'N N N# N0N=;O .$O7O$,O1#O>"P PP%=P20P> NN #N-0N:=O OO!$O.1O;>P PP"%P/2P;>Qctjjd}fdd_d_t dd}tjjd}t d i||jdd}tjdtjdtjjd}d_t d i||jd d}tjd tjd tjjd}d_t d i||jdd}tjdtjdt|j |j t|j |j y) Nrc|jdk(rdn t|}t|j_xjdz c_t j |d|jSr )rrfr!rrr rr)rPr$rr6s r>r6z"TestPower.test_batch..testsSffks1vJ!*doo>DO LLA L$$Q5<< rzTestPower.test_batchsFii##K0 =  bdCii##K05t515T\\4(T__a(ii##K0 6t626T\\2&T__b)ii##K0 8t848T\\1%T__d+TZZ,TZZ,rTc  tjjd}tjddgdfd }tjgdddtj f}dgdi}tjgd}t j||j||| }g}|D]_}|D]X} |dD]N D]G d fd } t j| |j| | } |j| jIPZatj||jj}t|j|d d y)Nlx7'\r皙?c tjttd|jdz}t j ||||S)Nr)rr$)r9 expand_dimsr;r5rr r)rPrr$popmeans_expandedpopmeanss r>r6z*TestPower.test_vectorization..testsF "xuQQR ?S9T U $$Q(9{*.0 0rT)r0rrIrrGrr)rr皙?rr2r3c6tj||S)Nr)r r)rPr$ alternative_i popmean_is r>test2z+TestPower.test_vectorization..test2s!#(#4#4Q AN$PPrTr1r r)rrrc) r9rIrrrbr{r rrappendrprnr)rr5r6nxr3r2rarsignificance_inx_irEtmprCrDr?s @@@r>test_vectorizationzTestPower.test_vectorizations3ii##N388QH% 0ZZ- .q"**} =!ABzz":; kk$ B\!')* .N .%+M%:.M%-. P$kk%T7EG 399-.. . .jjciioo.  3T=rTctjjd}tj}|j ddd}|j |j f}tjgd}tj||||}tj|ddtjf|jj}t|j|d y) N&(r0rG)rr0r)rr@rr1rr) r9rIrr rintegersrrbrror{rnr)rr5r6r7rr2ras r>test_ttest_ind_nullzTestPower.test_ttest_ind_nullsii##O4b#G<jj#**$zz"34 kk$^,O|ArzzM'BCIIOOT  test_ttest_1samp_powerz TestPower.test_ttest_1samp_power sii##O45464467 _-$d| kk%++SZZ'3FD 3T2rTN) r"r#r$rr6r%r(rrKrOrQr]rTr>r/r/xsP69p [[!-!-F [[%>%>N <3rTr/c eZdZdZdZdZejjde jje jjgejjdgddZ ejjde jje jjgejjdd ejd d zfd d ej dd fgdZdZejj'dZdZejjdddZejjdddZejjdddZdZejj5ejjdddZejj&ejjdddZejjdddZdZejjdddZejj5ejjdd d!Z gd"gd#d$d%d&d'd(d)d*Z!gd+gd,d-d.d/d0d1d2d*Z"ejj5ejjd3e!e"fd4Z#ejj&ejjd5d6d7Z$ejjd8d9d:Z%e&d;d gfe&dd gd<gd=gfe&d>d gd<gd=dd?ggfgZ'ejjd@e'dAZ(ejjSd dBZ*yC)DTestPermutationTestr&cLtjjd|_y)Nl F)r9rIrr5)rs r> setup_methodz TestPermutationTest.setup_method$s99(()<=rTcd}d}tjt|5tgddgf|dddd}tjt|5td|dddtjt|5td|dddd}tjt|5tgdgdf|d dddd }tjt|5tgdgdf|d dddd}tjt|5tgdgdf|d dddd}tjt|5tgdgdf|ddddd}tjt|5tgdgdf|ddddd}tjt|5tgdgdf|ddddd}tjt|5tgdgdf|ddddd}tjt|5tgdgdf|d dddd}tjt|5tgdgdf|ddddy#1swYwxYw#1swYRxYw#1swY/xYw#1swYxYw#1swYxYw#1swYxYw#1swY{xYw#1swYNxYw#1swY!xYw#1swYxYw#1swYxYw#1swYyxYw)NcFtj||f|jSr\rrXs r>rz:TestPermutationTest.test_permutation_test_iv..stat+s??Aq640:: :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_ivz,TestPermutationTest.test_permutation_test_iv)s ;G ]]:W 5 5 i!-t 4 5K ]]:W 5 ) T4 ( ) ]]9G 4 & Q % &/ ]]:W 5 E i3T D E5 ]]:W 5 6 i3T.4 6 6E ]]:W 5 K i3Tc J K> ]]:W 5 N i3Tu M N> ]]:W 5 O i3Tv N O@ ]]:W 5 H i3T G H@ ]]:W 5 I i3T H I0 ]]:W 5 O i3Tv N OA ]]9G 4 J i3Ty I J JW 5 5 ) ) & & E E 6 6  K K N N O O H H I I O O J JsJ( J5 K:K4K.K)(K6"LLLL) L5(J25J?K KK&)K36LL LL&)L25L>rHrX)pairingssamples independentc |jjd}|jjd} fd d _d _d|dd}t j ||f fd|dd|}t jd t jdd _t j ||f fd |dd|}t jd t jd d _t j ||f fd|dd|}t jd t jdt |j|jt |j|jy) Nr0c|jdk(rdn t|}t|j_xjdz c_t j ||t j ||z SNrr#)rrfr!rrr9r:)rPrYr$rrs r>rz1TestPermutationTest.test_batch..statisticis_ffks1vJ#&z93G3G#HI    " 7714(27714+@@ @rTrrT)r)rXr r)r-rHrrrr)r5rIrrr rrr) rrXrHrPrYrrQrRrrs @r>rzTestPermutationTest.test_batch_sx HHOOB  HHOOB  A   #9I"$%%q!fiLq3??LFJLY&&-Y))1- %%q!fiLr3??LFJLY&&+Y))2. %%q!fiLt3??LFJLY&&*Y))40T[[$++.T[[$++.rTzpermutation_type, exact_sizerZrr)r[r~r\r1cv|jjd}|jjd}d}|dd}tj||f|fd|dd|}t |j j dtj||f|fi|}t |j j |y)Nrc`tj||tj||z SrkrrXs r>rz8TestPermutationTest.test_permutations..statistic#7714(27714+@@ @rTT)rXr r)r)rH)r5rIr rrrr{) rrX exact_sizerHrPrYrrras r>test_permutationsz%TestPermutationTest.test_permutationss HHOOA  HHOOA  A%5"$$$aVYKA2>q/KEIKS**//3$$aVY?$?S**//rzNTestPermutationTest.test_randomized_test_against_exact_both..statisticrbrTTr\rGr rXr-rr5r)rr) r binomr rrrr9infrrr rrr5rGny permutationsrPrYrrrrarRs r>'test_randomized_test_against_exact_bothz;TestPermutationTest.test_randomized_test_against_exact_boths% S*B }}R"Wb)L888 JJNNN # JJNNN #!t A# [FtY'L'!%'iLRVVLtL}}... DKKd;rTcd\}}d\}}}d|z|kDsJtjj|}tjj|}||f}d} ddd||d } t|| fd |i| } t|| fd tj i| } | j | j k(sJt| j| jd y) N)rN)rri}rrc6tj||z |SrkrrXs r>rzQTestPermutationTest.test_randomized_test_against_exact_samples..statistics771q5t, ,rTTr[rGrhr)rr) r rrrr9rjrrrrks r>*test_randomized_test_against_exact_samplesz>TestPermutationTest.test_randomized_test_against_exact_sampless + S,B "u|### JJNNN # JJNNN #!t -# [FtY'L'!%'iLRVVLtL}}... DKKd;rTc d|j}}d\}}}tj||kDsJtjj |}tjj | |g} fd}t j|} ddd||d} t|| fd |i| } t|| fd tji| } | j| jk(sJt| j| jd y) Nr)r~r~i@rc6tj|dSNrr pearsonrrPrYs r>rzTTestPermutationTest.test_randomized_test_against_exact_pairings..statistic1ds>>!Q'* *rTTr[rGrhr)rr)r5r factorialr rrrr\rr9rjrrr)rrr5rGrlrmrPrrrrrarRrYs @r>+test_randomized_test_against_exact_pairingsz?TestPermutationTest.test_randomized_test_against_exact_pairingss'S *B   $|333 JJNNN # JJNNN #s + 44[A " [FtY'L'!%'iLRVVLtL}}... DKKd;rTrrc|jjdd}|jjddd}tj||d}fd }t ||f|t j d |j }t|j|j|j t|j|j|j y) Nrrr{rzrwr)r{ryrzexact)rrcJtj||djSr)r ks_2samprrPrYrs r>rz>TestPermutationTest.test_against_ks_2samp..statistic1ds#>>!QW.9;;D9 ErTrr)rr5r) r5rr r~rr9rjrrrr)rrrPrYrrras ` r>test_against_ks_2sampz)TestPermutationTest.test_against_ks_2samps HHOO!O , HHOOO 3>>!QKgN E 1v{+4$((D  x'9'9 J HOO$))DrTc|jjdd}|jjdd}dddd }||}tj||| }d }t ||f|t j ||j }t|j|j|j t|j|j|j y)Nrrr{rwrrrrrrcBtj||jSr\r ansarirrws r>rz.statistic1d<<1%// /rTrr) r5rr rrr9rjrrrr) rrrPrYalternative_correspondencealternative_scipyrrras r>test_against_ansariz'TestPermutationTest.test_against_ansari s HHOO!O , HHOO!O ,/8173>&@"7{C<<12CD 01v{+6DHHF  x'9'9 J HOO$))DrTc tjjdd|j}tjjdd|j}tj||d|}d}t ||f|dt j|d|j }t|j|j|j t|j|j|j y) N)rrwrrr{ryrHr@rrcFtj|||jSrk)r mannwhitneyurrXs r>rz@TestPermutationTest.test_against_mannwhitneyu..statistic*s%%a6@@ @rTT)r r)rr$r5r) r uniformrr5rrr9rjrrrr)rrrPrYrrras r>test_against_mannwhitneyuz-TestPermutationTest.test_against_mannwhitneyu"s MM  9!$((  K MM  9$TXX  N%%a L A1vyT+-66{$%4885  x'9'9 J HOO$))DrTctjjdd|j}tjjddd|j}tj||d}d }t ||f|t jd |j }t|j|j|j t|j|j|j y) Nrr)r{rzrHrwr)r{ryrzrHr|r.cFtj||djSN asymptoticr.)r cramervonmises_2samprrws r>rz9TestPermutationTest.test_against_cvm..statistic1d;s#--a5ACCL9 MrTrrr) r rrr5rrr9rjrrrr)rrPrYrrras r>test_against_cvmz$TestPermutationTest.test_against_cvm4s JJNNN B JJNNqN I--a7C M 1v{+4$((D  x'9'9 J HOO$))DrTr$)rdrcXtjjd}|jd}|jd}|jd}|||f}d}d}tj|} tj|} tj |d} tj |d } tj |d } | | | | | }| | | | | }d |d d|j d}t||fdtji|}t||fddi|}t|j||jt|j|j|jt|j|dt|j|jdy)NlY}y!dhrr)rrr)rrrc:tj|jSr\)r kruskalrrs r>rzITestPermutationTest.test_vectorized_nsamp_ptype_both..statistic1dXs==$'11 1rTc:tj|jSr\)r rrrs r>pvalue1dzFTestPermutationTest.test_vectorized_nsamp_ptype_both..pvalue1d[s==$'.. .rT)rrr)rrrrr#Frr\r r$rrXr5r)rrgQ?rr) r9rIrrr\ror5rrjrrrr)rr$r5rPrYrlrrrrrx2r:z2expected_statisticexpected_pvaluerrarRs r> test_vectorized_nsamp_ptype_bothz4TestPermutationTest.test_vectorized_nsamp_ptype_bothFspii##$78 JJQJ  JJIJ & JJIJ &1ay 2 / 44[A 11(;__Q * __Q * __Q *&r2r= R$7$T)$1$((Dt[MbffMMkLtLtL  '9 J t~~DIIF O$? DKKd;rTctjjdd|j}tjjdd|j}d}d}fd}t j |}|||d }|d} |d} d dd |jt jd } t||z f|fi| } t||f|fi| } t| j| j|j dk7r"t| j| |j t| j| |j t| j| j|j y)N)rr1rrrr@cDtj|djSNrrr wilcoxonr)rls r>statistic_1samp_1dzETestPermutationTest.test_against_wilcoxon..statistic_1samp_1ds>>!8BB BrTcFtj||djSrrrws r>statistic_2samp_1dzETestPermutationTest.test_against_wilcoxon..statistic_2samp_1ds>>!QF;EE ErTc4tj||S)Nr)r rrs r>test_1dz:TestPermutationTest.test_against_wilcoxon..test_1ds>>!QK@ @rTrr#Fr[)r r$rrXr5r)rr) r rrr5rr\r9rjrrrrr)rrrPrYrrrr6rrrrrQrRs ` r>test_against_wilcoxonz)TestPermutationTest.test_against_wilcoxonvs7 MM  9!$((  K MM  9$TXX  N C F A//811%  a[ #Q{$-dhh!vv' 1(:CdCA(:CdC TYYG + % DNNM J Zdii@ T[[tyyArTc d|jjddd}d||dk(<d d}||dd}}}tj||||}tj|f|d d t j |j| }t|j|j|j y)Nrrr0rrdc6tj|dkD|S)Nrr#r9rr s r>rz=TestPermutationTest.test_against_binomtest..statistics66!a%d+ +rTrGrTr[)r rXr)r5rrr) r5rNr rrr9rjrrr) rrrPrkrgrQrras r>test_against_binomtestz*TestPermutationTest.test_against_binomtests HH  a  ,!q&  ,A,Ca1??1a D$$aT96?13TXX1<>  HOO$))DrTc|jjd}||jjdztj|d}fd}t |f|dt j |j}t|j|j|jt|j|j|jy) Nr1rr|r.cHtj|djSr)r kendalltaurrws r>rz@TestPermutationTest.test_against_kendalltau..statistic1ds##Aq>HH HrTrZ)rXr)r5r) r5rr rrr9rjrrrr)rrPrrrarYs @r>test_against_kendalltauz+TestPermutationTest.test_against_kendalltaus HHOOO # Q' '##Aq9 It[:+-66txxA  x'9'9 J HOO$))DrTcfd}tjjd}|jddkDjt}|jdd|zzdkDjtt j j|d}t|f|dtj||}t j|| }t|j|dy) Nc@tj|dk(dk(zSrrrws r>rz@TestPermutationTest.test_against_fisher_exact..statistics6616a1f-. .rTl (yM4r g333333?g?rrZ)rXr)rr5r) r9rIrastypefloatr contingencycrosstabrrj fisher_exactrr) rrrr5rPtabrarRrYs @r>test_against_fisher_exactz-TestPermutationTest.test_against_fisher_exacts /ii##$78 ZZ]S ( ( / ZZ]T!V #c ) 1 1% 8((A.q1tY+-66{#&(!!#;? DG,rT)r2rc|jjd}|jjd}|jjd}tj||}tj||}tj||}|d}|||f}d}d}t j |}t j |} |t j|d |i} | t j|d |i} d |d d dd } t||fdt ji| } t||fddi| }t| j| |jt| j|j|jt| j| |jt| j|jdy)N)rrrr)rrr)rrrr#rcFtj|ddjS)NTrrankedr/)r page_trend_testrrs r>rzLTestPermutationTest.test_vectorized_nsamp_ptype_samples..statistic1ds"((d0<>>Gi HrTcFtj|ddjS)NTr|r)r rrrs r>rzITestPermutationTest.test_vectorized_nsamp_ptype_samples..pvalue1ds"((d0799? @rTr$FrrZrr)rxrrr)r5rIr rankdatarr\r9broadcast_arraysrrjrrrr)rr$rPrYrlrrrrrrrrrarRs r>#test_vectorized_nsamp_ptype_samplesz7TestPermutationTest.test_vectorized_nsamp_ptype_sampless HHOOO + HHOOO + HHOOO + NN14 ( NN14 ( NN14 ( aD1ay H @ 44[A 11(;&(;(;T(BNN ""5"5t"<H4H$T)$.q:t[MbffMMkLtLtL '9 J t~~DIIF O$))D DKKd;rTr)r"rg@r;g?g&_?g@gd~$I"@glu?)rPrY expected_lessexpected_2sidedexpected_Pr_gte_S_meanr expected_avg expected_std)okrGrrjml) rrrbirnrrgd?gd?gYo5?g@@@g{C@gDmFA@casec|d}|d}|d}|d}|d}|d}|d}|d} d } tjj5} | jtd t ||f| tj d } t ||f| tj d } dddt j||jt| j|dt j|dt| jj|dt| jj| d| j}| jj}t| j}tjtj | j|z tj ||z k\|z }t||y#1swYAxYw)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. rPrYrrrrrrcBtj||jSr\rrws r>rz7TestPermutationTest.test_with_ties..statistic1d9rrTz$Ties preclude use of exact statisticr)r)rrNrg|=rgư>)r9testingrr UserWarningrrjrrrrrr:rrfrr)rrrPrYrrrrrrrrrarRSr:rg Pr_gte_S_means r>test_with_tiesz"TestPermutationTest.test_with_tiess8 I I!"67_- 01!%&>!?N+ N+  0ZZ ) ) + =s JJ{$J K"Aq6;BFF/57C#QFKRVV0;=D  =  '9 J M> _5A..335|$O..224lN MM$$))+ %% &rvvc&;&;D&@A"$&&4. 12234 . >' = =s AGG(zalternative, expected_pvalue))rg?)rgOO?)rgrzDTestPermutationTest.test_against_spearmanr_in_R..statistic1dks??1a(22 2rTrZ)rXr)rrgvIh%<=rN)rr9rjrrrr)rrrrPrrrarYs @r>test_against_spearmanr_in_Rz/TestPermutationTest.test_against_spearmanr_in_RQsa( 1 0/ 3t[:+-66{L  '9 J O%@rTr-)rdrctjtd5tt j gd|dddy#1swYyxYw)Nz`batch` must be positive.rr)r6r7r8rqr_batch_generator)rr-s r>test_batch_generator_ivz+TestPermutationTest.test_batch_generator_ivts? ]]:-H I A --i? @ A A As "AArrC)rrrwr~r ziterable, batch, expectedcPttj||}||k(sJyr\)rqrr)riterabler-rgots r>test_batch_generatorz(TestPermutationTest.test_batch_generator}s';//%@AhrTc gd}gd}d}tj||f|dd}|j|j|j}}}dt j ||dz k\zt|z }||cxk(rd k(sJJy) N)rrrr)rrr1r~c4tj||dSrtrurws r>rzFTestPermutationTest.test_finite_precision_statistic..statistics>>!Q'* *rTFrZ)r rXrr&gUUUUUU?)r rrrrr9rrf) rrPrYrrarrnull correct_ps r>r!z3TestPermutationTest.test_finite_precision_statistics   +$$aVY56@B--S5J5J46tq5y011CI= )c)))))rTN)+r"r#r$rrUrYr6r%r'r9rIr=rrr rxrirdrnr(rqryrrrrr&rrrrrr tie_case_1 tie_case_2rrrr5batch_generator_casesrr,r!r]rTr>rSrS 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**rTrSc vtjgdt}tj|}t }d}t j |D]Q}|dz }tj||dd}|jt|Dcgc] }t|c}Stjtt|dz Dcgc]3}tjt!||dt!||dzd5c}}t#||t#t||ycc}wcc}w)N)rrrrrrrd)r9rintcumsumsetr_all_partitions_concatenatedsplitaddr; frozensetprodr5rfr rirr)rgncall_partitionsrpartition_concatenated partitioningr^rs r> test_all_partitions_concatenatedrs #&A 1BUNG"-"J"J1"MH1 xx 63B@ 5 !E1)A,!EFGH ww!&s1vax2 c!AB%j#a!g,?23H(#^$h/ "F2s D18D6fun_namerrrcVtjjdjdfd}fddfdd|d d }d d i}tt|}||}|j |fd }|d|dd||d|dd|d}|d|dd||d|dd|y)Nlr:r0rcFtjj|S)Nr)r rrrs r>rz&test_parameter_vectorized..rvsszz~~4c~::rTrA)rr5r/r[)rr5rX)rrrr)rGc|jdkDstj|sJtj||Sr_)rr9r7r:)rPr$rs r>rz,test_parameter_vectorized..statistics2vvzR^^Av666wwqt$$rT)rr TcN|jdk(sJtj|Sr)rr9r:rs r>rz,test_parameter_vectorized..statistics vv{{wwqzrTFr])r9rIrr?r r) rr fun_optionscommon_optionsroptionsrr5rs @@r>test_parameter_vectorizedrs ))   /C ZZRZ F;+13+7!917 #"%8)888)888)889)99rTceZdZdZy)TestMonteCarloMethodcd}tjjd}tjt |5t j|j|dddy#1swYyxYw)Nz.Use of `rvs` and `rng` are mutually exclusive.ir)rr5)r9rIrr6r7r8r MonteCarloMethod)rr=r5s r>test_rvs_and_random_statez.TestMonteCarloMethod.test_rvs_and_random_statesRBii##H- ]]:W 5 <  " "szzs ; < < r r srs 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<