o ֭jj8@sZddlmZddlZddlZddlZddlZddlmZm Z m Z m Z ddlm Z ddlmZddlmZddlmZmZmZmZmZmZmZddlmZmZmZdd lmZdd lm Z dd l!m"Z"m#Z#Gd d d Z$GdddZ%GdddZ&Gddde Z'GdddZ(GdddZ)GdddZ*GdddZ+GdddZ,GdddZ-dS) )productN)assert_ assert_equalassert_allcloseassert_almost_equal)raises) distributions)epps_singleton_2sampcramervonmises_cdf_cvmcramervonmises_2samp_pval_cvm_2samp_exact barnard_exactboschloo_exact) mannwhitneyu _mwu_state_MWU)check_named_results)_TestPythranFunc)SmallSampleWarningtoo_small_1d_not_omitc@s<eZdZddZddZddZddZd d Zd d Zd S)TestEppsSingletoncCsJtgd}tgd}t||\}}t|dddt|ddddS)N) gffffffֿgffffff@gGz?\(\?ffffff?gQ@gq= ףp?gGzgGz׿gp= #(@) gffffffg333333ÿgףp= @g @gGz @g)\(@g@g(\@g(\ @333333!@gHzG.@rdecimalgQ,r?)nparrayr rselfxywpr']/var/www/html/CropPilot/venv/lib/python3.10/site-packages/scipy/stats/tests/test_hypotests.pytest_statistic_1s z"TestEppsSingleton.test_statistic_1cCsBtd}td}t||\}}t|dddt|ddddS) N)rrr*r*r*rrrrr,r,r, r.r.r.)r.r+rr,r.r.rr,r-r.rrr/rrr,r0r.g!@MbP?atolg&J?rr)rr r rrr!r'r'r(test_statistic_2%s  z"TestEppsSingleton.test_statistic_2c Cstjdtdtd}}tt|t|\}}tt|t|\}}t||\}}t||ko8|knt||koI|kdSdS)N)rrandomseedaranger listtupler) r"r#r$w1p1w2p2w3p3r'r'r(test_epps_singleton_array_like/s &z0TestEppsSingleton.test_epps_singleton_array_likecCsjdtd}}tjttdt||}t|jtj t|j tj WddS1s.wYdS)Nrr*rr+r.match) rr:pytestwarnsrrr r statisticnanpvaluer"r#r$resr'r'r(test_epps_singleton_size:s  "z*TestEppsSingleton.test_epps_singleton_sizecCs0dddddtjftd}}ttt||dS)Nrr*rr+r,r.)rinfr: assert_raises ValueErrorr r"r#r$r'r'r(test_epps_singleton_nonfiniteBsz/TestEppsSingleton.test_epps_singleton_nonfinitecCs2tdtd}}t||}d}t||dS)Nr6)rIrK)rr:r r)r"r#r$rM attributesr'r'r( test_namesGs zTestEppsSingleton.test_namesN) __name__ __module__ __qualname__r)r4rCrNrSrVr'r'r'r(rs    rc@sxeZdZddZddZddZddZd d Zd d Zd dZ ddZ e j dddggddZddZddZdS)TestCvmcC ttgddgddddS)N)gy;i?g#^ ?gE>?g D )?r+{Gz?皙??g+?-C6?r2rr r"r'r'r( test_cdf_4R   zTestCvm.test_cdf_4cCr[)N)g8*5?g@߾?gHm?g%1 ?r.)r]r^r_g333333?r`r2rarbr'r'r( test_cdf_10XrdzTestCvm.test_cdf_10cCr[)N)g}tg?g`?gI5o?gׁsF?r\r`r2rarbr'r'r( test_cdf_1000^rdzTestCvm.test_cdf_1000cCsttgdgddddS)N)a+e?+?&pn?+MJA?r\r`r2rarbr'r'r( test_cdf_infds   zTestCvm.test_cdf_infcCs4ttddgdddgttddgdddgdS) NgX(~$?gUUUUU5f@irrgaah?g"@)rr rbr'r'r(test_cdf_supportjszTestCvm.test_cdf_supportcCs$ttgddtgddddS)N)rhrirjrkd'r`r2rarbr'r'r(test_cdf_large_nos    zTestCvm.test_cdf_large_ncCsLtdtddko dkntdtdko dkdSdS)NgwJ?gt@rf?)rr rbr'r'r( test_large_xvs"*zTestCvm.test_large_xcCs<d}tt|dd}tt|j|dkt|jddS)N 皙?normrrr)r ronesrr rIrrK)r"nrMr'r'r( test_low_pszTestCvm.test_low_pr#r'?cCsZtjttdt|d}t|jtjt|j tjWddS1s&wYdS)NrErv) rGrHrrr rrIrrJrKr"r#rMr'r'r(test_invalid_inputs  "zTestCvm.test_invalid_inputcCstgdd}t|jdddt|jdddtgddd}t|jdddt|jd ddtgd d }t|jd ddt|jd dddS)N)g333333r*rg?r+皙?333333?rvgZ %q?ư>r2gE ж?)rrzg!O!W*?gz"W`?) rr*r,ffffff?gQ? ?@exponge.?gnz\(r?)r rrIrK)r"rMr'r'r( test_values_RszTestCvm.test_values_RcCs|tdd}}t|tjj}t|d}t|j|jf|j|jft|tj j|}t|d|}t|j|jf|j|jfdS)Nr,)rffffff?rbeta) rr:r rrcdfrrIrKr)r"r#argsr1r2r'r'r(test_callable_cdfs  zTestCvm.test_callable_cdfN)rWrXrYrcrergrlrnrqrsryrGmark parametrizer|rrr'r'r'r(rZNs   rZc@speZdZejddgidgiggdgddZddZd d Zgd Z gd Z d dddgddddgddddgd dddgddddgddddggZ ejde ddZ d dddgddddgdddd gd d!dd"gdd!dd#gdd!dd$ggZ ejde d%d&Zd'd(Zgd)gd*gd+d,Zgd-gd.gd/gd0d1Zgd2gd3gd4gd5gd6d7Zgd8gd9gd:gd;gd<gd=d>Zd?d@ZdAdBZdCdDZd dddEgddddFgddddGgd dddEgddddFgddddEggZejdHedIdJZdKdLZejdMddgdNdOZdPdQZgd1dRdSdTdUejdRdVdWdXdXdYg dZd[fgd1dRdSdTdUejejdVdWdXdXdYg d\d]fdWdVejdXgdRdSdTdUejdRdVdWdXdXdYg d^d_fdWdVejdXgdRdSdTdUejejdVdWdXdXdYg d`dafdWejejdXgdRdSdTdUejejdVdWdXdXdYg dbdcfgZejddededfZgdggdhgdigdjgdkgdlgdmgdngdog Z ejdpe dqdrZ!dsdtZ"gd,dudvgddwggd,dudvgddwggd,dudvgd dxggd,dVgddyggd,dVgddyggd,dVgd dzgdWdVgdWdVgdd{gdWdVgdWdVgdd{gdWdVgdWdVgd d|gg Z#ejgd}e#d~dZ$ddZ%ejdgdddZ&dS)TestMannWhitneyU kwargs_updater#r$r#r$cCstddg}tddg}t||d}||tjttdtdi|}t |j tj t |j tj WddS1sAwYdS)Nrr*rr+rrEr') rr dictupdaterGrHrrrrrIrJrK)r"rr#r$kwargsrMr'r'r( test_emptys  "zTestMannWhitneyU.test_emptycCs tddg}tddg}ttddt||ddWdn1s&wYttd dt||dd Wdn1sCwYttd dt||d d Wdn1s`wYttddt||ddWddS1s~wYdS)Nrr*rr+z`use_continuity` must be onerEekki)use_continuityz`alternative` must be one of alternativez`axis` must be an integerrzaxisz`method` must be one ofmethod)rr rPrQrrRr'r'r(test_input_validations"z&TestMannWhitneyU.test_input_validationcCstjdd}tj|d}tj|d}t||}t||dd}t||dd}|j|jks3J|j|jks;Jtj|d}tj|d}t||}t||dd}t||dd}|j|jksfJ|j|jksnJt||}t||dd}t||dd}|j|jksJ|j|jksJtj|d}tj|d}t||}t||dd}t||dd}|j|jksJ|j|jksJtj|d}tj|d}|d|d<t||}t||dd}t||dd}|j|jksJ|j|jksJdS)Nrr0 asymptoticrexactr)rr8r9randrrK)r"rxr#r$autorrr'r'r( test_autosH       zTestMannWhitneyU.test_auto)gm9Aj@g+H3[@gi>s@) g#hA{@glz@gcDf@gdz*h@gZA@gI9^YQa@g`@g՞p@g:q@g&@gZ |@g`r@gMc3g@ two-sidedrrr) +?less)r +?greater)r缌%c?r)rg9:?)rg9:?)rg*::?)kwdsexpectedcCs$t|j|jfi|}t||dSN)rr#r$rr"rrrMr'r'r( test_basic szTestMannWhitneyU.test_basicT)rr)r)rr)rrF)rgl,KNh?)rgiژ?)rgl,KNh?cCs(t|j|jfddi|}t||dS)Nrr)rr$r#rrr'r'r(test_continuity2s z TestMannWhitneyU.test_continuityc Csgd}tgd}tgdd}tgdd}|d||||||||||dg}t||ddd}gd }gd }t|j|t|j|dS) NrDrr*rr+r,)rrrrrr])rrrrrr)rr)r. !@r0rr/r-)rg]U?g[?gi\ ?gZX<_?gx.?g ?)rr rrrIrrK) r"r#y0dydy2r$rM U_expected p_expectedr'r'r(test_tie_correct@s* z!TestMannWhitneyU.test_tie_correct)g?r_g?)r}皙?皙?r~)r^r}rrr_g?rr*r)rrr~)gx&?g/$?gJ +?rr~)y&1?v/?gv/?rgjt?~jt?ʡE?) gy&1?gV-?rr}gS?gv?gʡE?g'1Z?gm?rD)gK7A`?gZd;O?r_gMbX?)Mb?RQ?RQ?M?rr) g;On?;On?V-?g?gJ +?rgx&?r_gCl?) Mb?Mb?Mb?gy&1?rM?g|?5^?gn?g\(\?!rh?K7?)Mbp?rrr~jt?g333333?g"~j?ףp= ?gzG?K7?gGz?gl?r_gI +?r)rrg1Zd?r)rrr1Zd?g%C?rr) g~jt?g~jt?rgsh|??gS㥛?rrg +?rr) g{Gzt?r]g~jt?gL7A`?rgjt?gPn?gI +?gX9v?gQ?gMb?gsh|??gK7A`?)Mb`?rg;On?gQ?g9v?gˡE?gT㥛 ?gbX9ȶ?grh|?gQ?rgx&?gv/?gMbX?g(\?gQ?)r1rrrg9v?g/$?rrgL7A`?g ףp= ?gQ?rgK7?g`"?g7A`?rgV-?gjt?gˡE?)rr*rr+r,r-c Csttdtdd|j|j|j|jd}|D]n\}}|D]e\}}t dt |}tj ||t tj j|d|ddt d||d}t tj j|dtj j|dtj j|ddtj j|d}t ||dddtj ||tj j|d} t || q qdS) Nsr)rr+r,r-kr1r2rr)setattrrrpn3pn4pm5pm6itemsrr:lenr set_shapesrrsfpmf) r"p_tablesrxtablemr&uu2rpmf2r'r'r(test_exact_distributionps,   z(TestMannWhitneyU.test_exact_distributioncCstjdtjd}tjd}t||dd}t||dd}|j|jks(Jt|j|jdks5Jtjd}tjd}t||dd}t||dd}|j|jksWJt|j|jdksdJdS) Nrr,rrrr](r1)rr8r9rrrIabsrK)r"r#r$res1res2r'r'r(test_asymptotic_behaviors     z)TestMannWhitneyU.test_asymptotic_behaviorcCsrtgdddgddd}tgdddgddd}t|j|j|jdks&Jtgdddgd dd}t|d dS) Nrrz@rrrrr_r)rr)rrrK)r"res_lres_grMr'r'r(test_exact_U_equals_meansz)TestMannWhitneyU.test_exact_U_equals_meanrr)rr_)rg郡E?)rresultcCsttdi||dS)Nrr*rr*)rr)r"rrr'r'r(test_scalar_datasz!TestMannWhitneyU.test_scalar_datacCsHttdddddttdddddttddddddtjfdS) Nrrr)r_rrF)rrr_)rrrrJrbr'r'r(test_equal_scalar_datas z'TestMannWhitneyU.test_equal_scalar_datarcCsztjdd}d\}}tj|dd}tjd|ddd}t||||d }d }|jj|ks1J|jj|ks9Jt||d t||d }}|d }|j |j ksTJt |||f}t |||f}|jdd |ksqJ|jdd |ks|Jt |} t |} t d d|DD]} || } || } t| | |d}|j| | <|j| | <qtj |j| tj |j| dS)Nr)r/r.rr0r-rr})rr)r-rr0r)N.cSsg|]}t|qSr')range).0ir'r'r( sz8TestMannWhitneyU.test_gh_12837_11113..r)rr8r9rrrKshaperImoveaxisndim broadcast_tozerosrtestingr)r"rrrrxr#r$rMr statisticspvaluesindicesxiyitempr'r'r(test_gh_12837_11113s4     z$TestMannWhitneyU.test_gh_12837_11113cCs~gd}gd}t||}tj|d<t||}t|j|jt|j|jtj|d<t||}t|jtjt|jtjdS)NrD) rr-r/r0rrr*rr+r+r,r+)rrrOrrIrKrJ)r"r#r$rrres3r'r'r( test_gh_11355s     zTestMannWhitneyU.test_gh_11355rr-r/r0r*rr+r,r.g+zQ?rg}$k\?g1@g!˛G*?rg,s?8@gFHQ?)r#r$rIrKcCs2t||dd}t|j|ddt|j|dddS)Nrr-q=r2)rrrIrK)r"r#r$rIrKrMr'r'r(test_gh_11355b szTestMannWhitneyU.test_gh_11355b)Trrg&?)TrrgO?)TrrgO?)Frrg9@VN!x?)Frrg9M>?)Frrg9M>?)Trrg?UV?)TrrgߺVJH?)TrrgVJH?)rrr pvalue_expc Cs:d}d}d}t|||||d}t|j|t|j|dS)N#) rug(\?g= ףp=?gp= ף?g333333?gGz?g(\?g= ףp=?rg\(\?)gffffff?g)\(?rgGz?g\(\?rrr)rrrIrrK) r"rrrr statistic_expr#r$rMr'r'r( test_gh_91841s zTestMannWhitneyU.test_gh_9184cCsfttjtjtjtjtjg}ttjtjtjtjtjg}t||}t|jtjt|jtjdSr)rr rJrrrIrK)r"abrMr'r'r( test_gh_4067Ps  zTestMannWhitneyU.test_gh_4067rzr)rga׀}?)rrr)rzg?h?)rzr)r*g5&#\?)r*r)r#r$rrcCs$t||d|dd}t||dddS)NTrrrrtol)rr)r"r#r$rrrMr'r'r( test_gh_2118hs zTestMannWhitneyU.test_gh_2118cCstjd}d\}}|j|d}|j|d}ttdtddtjtj ||dd}tjj j }|dt |j |||j d ksDJtj ||dd|tjj j ksUJtjtj |d|dd d tjj j }|dd ksrJtj d||dd d |tjj j ksJdS) Ng>mjK)r,rsizerrrrrrr)rr)rr8 default_rngrrrrresetstatsrconfigurationsrminrI)r"rngrrxr#r$rMrr'r'r(test_gh19692_smaller_tableps"     &  z+TestMannWhitneyU.test_gh19692_smaller_tabler)rrrcCsxtjd}|jdd}|jdd}tj||t|dd}tj||d|dd}t|j|jdd t|j|jdd dS) Nr)r*r,r)r*r-r)rrrrV瞯.rr) ALL_INTEGER ALL_FLOATdtypesrr:r/ functoolspartialr somersd partialfuncr)r" input_arrayr'rbr( setup_methods  zTestSomersD.setup_methodcGs6|j|}t|j|jjddt|j|jjdddS)Nr%r2)r7rrIrrKr"rrMr'r'r( pythranfuncs zTestSomersD.pythranfunccCsfgdgdgdg}t|}|tj}tj|fi|}t|j|jddt|j|jdddS)N)rmr/r)r/r=rrt)rrr*r/r%r2)r r6get_optional_argsrrIrK)r"rr optional_argsrr'r'r(test_pythranfunc_keywordss   z%TestSomersD.test_pythranfunc_keywordscCsngd}gd}d}t||}t|j|dddt|j|dddgd}gd }d}t||}t|j|dddt|j|dddgd }gd }d }t||}t|j|dddt|j|dddtd }td }d}t||}t|j|dddt|j|dddtd }tgd}d}t||}t|j|dddt|j|dddtd }td ddd}d}t||}t|j|dddt|j|dddtd }tgd}d}t||}t|j|dddt|j|dddgd}gd}d}t||}t|j|dddt|j|dddtgdgd}t|jtjt|jtjtgdgd}t|jtjt|jtjtgdgd}t|jtjt|jtjtdgdg}t|jtjt|jtjtgg}t|jtjt|jtjtd}td}t t tj||dS)N)r,r*rrr-r+r/r0)r,r*r-rrr0r/r+rrrr%r2r) rr,r*rrr-r+r/r0) r,r*rr-rrr0r/r+)r,r*rrr-r+r/)r,r*r-rrr/r+)g+$I$I¿g=/3n+?r.rrr) rr*rrr+r-r,r/r0r)gs'}'?rDr)gr) rr/r0r-r,rr+r*rr)g}'}'rD)rtr*rrtr*)rr+r/rr)g.ʂ?)r*r*r*)r*rr*g$@g4@) r r6rrIrKrr:r rJrPrQ)r"r#r$rrMx1x2r'r'r(test_like_kendalltaus                 z TestSomersD.test_like_kendalltaucCsgd}gd}d}d}d}t||}t|j|ddt|j|ddt|jjd t||}t|j|ddt|j|ddt|jjd dS) N)rrrr*r*r*r*r*rrrr*r*r*r*r*r*r*rrrrrr)rrrrrrrrrrr*r*r*r*r*r*r*r*r*r*r*r*r*r*gCE]t?g^_?gO((Ƿ?r%r2r`)rr*r*r)r r6rrIrKrrr)r"r#r$d_crd_rcr&rMr'r'r(test_asymmetry)s  zTestSomersD.test_asymmetrycCstddgddgddgddgddgg}|j}d}tt|j|td d gd d gd d gg}d \}}tt|j|tt|jj|td d gd d gd d gg}d}tt|jj|dS)Nr0r*r-r,rr+rgHHHHHH?r<rUr6)gM&w?rrgtE]t)rr Trr r6rI)r"rdyxdxyr'r'r(test_somers_originalBs(z TestSomersD.test_somers_originalc Cs6d}d}t|}tjdtjj|t||d|}t |}tj |dt |ddd}t |}tj |dt |ddd}t |} tj |dt |dddd} t | } t |j dd d t |j |j t |j | j t |j | j t |jd d d t |j|jt |j| jt |j| jdS) Nror+r-rr&r*rrgayr%r2gPj$?)rprodr8r9r  multinomialrvsrwreshaper6insertrrrIrK) r"NrrrrMs2rs3r s4res4r'r'r(*test_contingency_table_with_zero_rows_colsXs(        z6TestSomersD.test_contingency_table_with_zero_rows_colsc Csd}d}t|}tjdtjj|t||d|}|d}d}t t |d t |Wdn1s;wY|d}d }t t |d t |Wdn1s\wYd }t t |dt ggWdn1szwYt t |dt d ggWdn1swYt d }t t |d t |Wdn1swYd |d <t t |dt |WddS1swYdS)NrorSrrTr*z:All elements of the contingency table must be non-negativerEr]z5All elements of the contingency table must be integerz?At least two elements of the contingency table must be nonzero.r)rrr) rrUr8r9r rVrWrwrXrPrQr6r) r"rZrrrs5messages6s7r'r'r(test_invalid_contingency_tablesws<        "z+TestSomersD.test_invalid_contingency_tablescCsbgd}ddtjg}gd}ddtj g}t||}t||}t|j|jt|j|jdS)Nrrg@)rr*rrrF)rrOr r6rrIrK)r"r#rHr$y2rMrr'r'r(test_only_ranks_matters   z"TestSomersD.test_only_ranks_mattercCs6td}td}t||}t|jtddS)Nr.)rr:r r6rreyerLr'r'r(test_contingency_table_returns   z)TestSomersD.test_contingency_table_returncCs`gd}gd}tj||dd}|jdksJtj||dd}t|j|jt|jd|jdtj||d d}t|j|jt|j|jd|tj||dd}|jdks\Jtj||d d}t|j|jt|jd|jdtj||dd}t|j|jt|j|jdtjt d d tj||d dWddS1swYdS) Nr)r,r-r/r0r/rrrrrr*rz`alternative` must be...rE ekki-ekki) r r6rIrrrKreverserGrrQ)r"rGrHrrMr'r'r(test_somersd_alternatives,"z$TestSomersD.test_somersd_alternativepositive_correlation)FTcCstd}|r |nt|}|rdnd}tj||dd}|j|ks#J|jdks*Jtj||dd}|j|ks9J|j|r?dndksDJtj||dd}|j|ksSJ|j|rYdndks^JdS) Nr.rrrrrrr)rr:flipr r6rIrK)r"rlrGrHexpected_statisticrMr'r'r( test_somersd_perfect_correlations  z,TestSomersD.test_somersd_perfect_correlationcCsVddg}d}tdtj||d}tj||d}d}t||j}t||dddS) Nrr*@Bi_rgHz Yr%r2)r8r9choicesr r6rIr)r"classes n_samplesr#r$ val_sklearn val_scipyr'r'r(!test_somersd_large_inputs_gh18132s z-TestSomersD.test_somersd_large_inputs_gh18132N)rWrXrYr9r;rBrIrMrRr_rdrfrhrkrGrrrorvr'r'r'r(r.s o#  )  r.c@seZdZdZejdddgddggdfdd gd d ggd fd d gdd ggdfd dgddggdfd dgddggdfd dgddggdfdd gddggdfd dgddggdfddgdd ggdfdd gddggdfd d gdd ggdfg d d!Zejdddgddggd"fdd gd d ggd#fd d gdd ggd$fd dgddggd%fd dgddggd&fd dgddggd'fdd gddggd(fd dgddggd)fddgdd ggd*fdd gddggd+fd d gdd ggd$fg d,d-Zd.d/Z ejdddgddggd0fgd1d2Z ejddd gddggd3e j ffd dgddggd3e j ffgd4d5Z ejdd d gdd ggd6fd d7gd8dggd9fd:d;gdd?d@gdAdBZdCS)DTestBarnardExactz8Some tests to show that barnard_exact() works correctly.input_sample,expected+rr.')gXyq @g{2s&Q7?ror*rfr,)gllgEA]0K?r/r0)*)1%g_?r)g_c1?g=?rT)g5PyQgQ@2?rr<)ggJ"?)g_c1gwݝل?rr+)g7@g?r)g~t,?3O?r-)gr?~CY7?cCs(t|}|j|j}}t||g|dS)zThe expected values have been generated by R, using a resolution for the nuisance parameter of 1e-6 : ```R library(Barnard) options(digits=10) barnard.test(43, 40, 10, 39, dp=1e-6, pooled=TRUE) ``` NrrIrKrr" input_samplerrMrIrKr'r'r( test_precise szTestBarnardExact.test_precise)g7\@gA2?)gXS;gh?)g>!Ɏg6?)gSy@?g^F?)g-gXI#?)gaЍgo?)gb]?gFugH ?)g 6ҭ@g?)gi( r})gNXzr~cCs,t|dd}|j|j}}t||g|dS)zThe expected values have been generated by R, using a resolution for the nuisance parameter of 1e-6 : ```R library(Barnard) options(digits=10) barnard.test(43, 40, 10, 39, dp=1e-6, pooled=FALSE) ``` F)pooledNrrr'r'r(test_pooled_param's z"TestBarnardExact.test_pooled_paramcCd}tt|dtddgddggddWdn1swYd }tt|dttd ddWdn1sBwYd }tt|dtd dgddggWdn1sdwYd }tt|dtddgddggdWddS1swYdS)N7Number of points `n` must be strictly positive, found 0rErr*rr+rrx,The input `table` must be of shape \(2, 2\).r-*All values in `table` must be nonnegative.rzI`alternative` should be one of {'two-sided', 'less', 'greater'}, found .* not-correct)rPrQrrr:rXr" error_msgr'r'r( test_raisesD$"zTestBarnardExact.test_raisesrEcC6t|}|j|j}}t||dt||ddSNrrrrIrKrrr'r'r(test_edge_cases^sz TestBarnardExact.test_edge_casesrrcCrrrrr'r'r(test_row_or_col_zerojz%TestBarnardExact.test_row_or_col_zero)r{gE\/??,)ggQ5rDr7i)g&X}>rDrrrc Csf|\}}|dkrt|dddddf}| }t||d}|j|j}}t||g||gdddS)a "The expected values have been generated by R, using a resolution for the nuisance parameter of 1e-6 : ```R library(Barnard) options(digits=10) a = barnard.test(2, 7, 8, 2, dp=1e-6, pooled=TRUE) a$p.value[1] ``` In this test, we are using the "one-sided" return value `a$p.value[1]` to test our pvalue. rNrrHz>r2)rr rrIrKr) r"rrr expected_statless_pvalue_expectrMrIrKr'r'r(test_less_greaterws  z"TestBarnardExact.test_less_greaterN)rWrXrY__doc__rGrrrrrrrrJrrr'r'r'r(rwsr      rwc@seZdZdZdZejdddgddggdfdd gd d ggd fdd gd dggdfd dgd d ggdfddgd dggdfdd gddggdfddgddggdfddgddggdfd dgddggdfg ddZejdddgd dggd fddgddggd!fdd gd d ggd"fdd#gd d ggd$fdd gd dggd%fddgd dggd&fdd gddggdfddgd'dggdfddgddggd!fddgddggd(fd dgddggd)fg d*d+Z ejdddgd dggd,fddgddggd-fdd gd d ggd.fdd gd dggd/fddgd dggd0fdd gddggd1fddgddggd-fddgddggd2fgd3d4Z d5d6Z ejdddgdd gge j e j ffddgd dgge j e j ffgd7d8Zd9d:Zejd;d<d=d>Zd?S)@TestBoschlooExactz9Some tests to show that boschloo_exact() works correctly.rrxr*r/r0)<vB\?g/??r,rr.)gM?gA>?rrTr<)_VѶ?g֭?)u%?gc'?rr+rrr)r_g?rt)+f?gXc}v?%)gZыD?ggi]?cC2t|dd}|j|j}}t||g||jddS)aThe expected values have been generated by R, using a resolution for the nuisance parameter of 1e-8 : ```R library(Exact) options(digits=10) data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE) a = exact.test(data, method="Boschloo", alternative="less", tsmethod="central", np.interval=TRUE, beta=1e-8) ``` rrr2NrrIrKrATOLrr'r'r( test_lesss zTestBoschlooExact.test_lessryrrz) k\2?g0,%?)gKv?gN3?)rg'&5?r|)gw@_?g7?)gi{?gɑ)z?)օa?g1|?r-)gY<;?gND?)ge?gG`?cCr)aThe expected values have been generated by R, using a resolution for the nuisance parameter of 1e-8 : ```R library(Exact) options(digits=10) data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE) a = exact.test(data, method="Boschloo", alternative="greater", tsmethod="central", np.interval=TRUE, beta=1e-8) ``` rrr2Nrrr'r'r( test_greaters zTestBoschlooExact.test_greater)rgqQS,5?)rgG?/??)rgKE`?)rghr1ֽ?)rgrfb?)r_g?)rgP:pRv?cCs4t|ddd}|j|j}}t||g||jddS)aThe expected values have been generated by R, using a resolution for the nuisance parameter of 1e-8 : ```R library(Exact) options(digits=10) data <- matrix(c(43, 10, 40, 39), 2, 2, byrow=TRUE) a = exact.test(data, method="Boschloo", alternative="two.sided", tsmethod="central", np.interval=TRUE, beta=1e-8) ``` r@)rrxr2Nrrr'r'r(test_two_sidedsz TestBoschlooExact.test_two_sidedcCr)NrrErr*rr+rrrr-rrzK`alternative` should be one of \('two-sided', 'less', 'greater'\), found .*r)rPrQrrr:rXrr'r'r(rrzTestBoschlooExact.test_raisescCrr)rrIrKrrr'r'r(rrz&TestBoschlooExact.test_row_or_col_zerocCs`ddgddgg}t|ddj}t|ddj}dt||dks!Jt|ddj}|d ks.JdS) Nrrrtrrrr*rrr)rrKr")r"tblplpgptr'r'r(test_two_sided_gt_1s z%TestBoschlooExact.test_two_sided_gt_1r)rrcCs>ddgddgg}t||dj}tj||dd}t||dS)Nr*r/r0rr)rrIr  fisher_exactr)r"rr boschloo_statfisher_pr'r'r(test_against_fisher_exact)sz+TestBoschlooExact.test_against_fisher_exactN)rWrXrYrrrGrrrrrrrrJrrrr'r'r'r(rsr     rc@seZdZejdgedfeddgfgddZddZ dd Z d d Z ejd gd ddZ ejj ddZddZddZddZdS) TestCvm_2samprr,rcCsXtjttdt|}t|jtjt|j tjWddS1s%wYdS)NrE) rGrHrrr rrIrrJrKr:r'r'r(test_too_small_input4s "z"TestCvm_2samp.test_too_small_inputcCsNtd}d}tjt|dt||dWddS1s wYdS)Nr,z/method must be either auto, exact or asymptoticrExyz)rr:rGrrQr )r"r$msgr'r'r(r|<s "z TestCvm_2samp.test_invalid_inputcCsNgd}gd}t||}tt|t|}t|j|jf|j|jfdS)N)r*rr+r/r-)rrrtr>)r rr rrIrKr"r#r$rrr'r'r(test_list_inputBs  zTestCvm_2samp.test_list_inputcCs>gd}gd}t||}t|jdddt|jddddS)N) gffffff@g @rgffffff!@皙"@g#@g333333$@g333333%@gffffff&@)g@g@g@@333333@gffffff @g333333"@g#@g%@g&@g'@g(@g)@g*@g333333-@gS㥛?r1r2g ףp= ?r])r rrIrKr"r#r$rr'r'r(test_example_conoverIs  z"TestCvm_2samp.test_example_conoverzstatistic, m, n, pval))ir,r-gcj`?)iir/r/gtE]t?)i@r+r-g88?)ir-r/gXwS?cCstt||||dSr)rr )r"rIrrxpvalr'r'r(test_exact_pvalueSs zTestCvm_2samp.test_exact_pvaluecCstjdtjjdd}tjjdd}t||}td|jko$dknt||d}td|jko=dkdSdS)Nirpri rrr}) rr8r9rrvrWr rrKrr'r'r(test_large_sample^s  (zTestCvm_2samp.test_large_samplecCsdtjdtjd}tjd}t||dd}t||dd}t|j|jt|j|jdddS) Nrr/r0rrrr]r2) rr8r9rr rrIrrKrr'r'r(test_exact_vs_asymptoticjs   z&TestCvm_2samp.test_exact_vs_asymptoticcCsttd}gd}t||dd}t||dd}t|j|jtd}t||dd}t||dd}t|j|jdS)NrT)r_g@g333333*@rrrrr)rr:r rrKrr'r'r(test_method_autoss  zTestCvm_2samp.test_method_autocCsVtd}t||}t|j|jfdt|dd|dd}t|j|jfddS)Nr|rCr+)rr:r rrIrKr{r'r'r(test_same_inputs  zTestCvm_2samp.test_same_inputN)rWrXrYrGrrrr:rr|rrrxslowrrrrr'r'r'r(r3s"      rc@s.eZdZgdgdgdfZgdgdgdfZgdgdgdfZdZdZd Ze j j d eed feed feed ffgd dddZ dZ dZe j j d ee dfeedffddgdddZddZddZddZdd Zd!d"Zd#d$Ze j d%d&d'd(Ze j d)gd*d+d,Zd-d.Zd/S)0 TestTukeyHSD)r 7@ffffff:@皙;@fffff=@)ffffff<@皙A@=@皙@@皙>@)g:@gL<@gL8@g333333:@g;@)r rgHzG:@rrrrr)r rr) rrrrrrrrrraK Comparison LowerCL Difference UpperCL Significance 2 - 3 0.6908830568 4.34 7.989116943 1 2 - 1 0.9508830568 4.6 8.249116943 1 3 - 2 -7.989116943 -4.34 -0.6908830568 1 3 - 1 -3.389116943 0.26 3.909116943 0 1 - 2 -8.249116943 -4.6 -0.9508830568 1 1 - 3 -3.909116943 -0.26 3.389116943 0 aS Comparison LowerCL Difference UpperCL Significance 2 - 1 0.2679292645 3.645 7.022070736 1 2 - 3 0.5934764007 4.34 8.086523599 1 1 - 2 -7.022070736 -3.645 -0.2679292645 1 1 - 3 -2.682070736 0.695 4.072070736 0 3 - 2 -8.086523599 -4.34 -0.5934764007 1 3 - 1 -4.072070736 -0.695 2.682070736 0 aS Comparison LowerCL Difference UpperCL Significance 2 - 3 1.561605075 4.34 7.118394925 1 2 - 1 2.740784879 6.08 9.419215121 1 3 - 2 -7.118394925 -4.34 -1.561605075 1 3 - 1 -1.964526566 1.74 5.444526566 0 1 - 2 -9.419215121 -6.08 -2.740784879 1 1 - 3 -5.444526566 -1.74 1.964526566 0 zdata,res_expect_str,atolr`g|=)equal size samplezunequal sample sizezextreme sample size differences)idsc Cstj|ddddtdd}tj|}|}|D]G\}}} } } } t |dt |d}}t |j ||f| |dt |j ||f| |dt |j ||f| |dt |j||fd k| dkqdS) a SAS code used to generate results for each sample: DATA ACHE; INPUT BRAND RELIEF; CARDS; 1 24.5 ... 3 27.8 ; ods graphics on; ODS RTF;ODS LISTING CLOSE; PROC ANOVA DATA=ACHE; CLASS BRAND; MODEL RELIEF=BRAND; MEANS BRAND/TUKEY CLDIFF; TITLE 'COMPARE RELIEF ACROSS MEDICINES - ANOVA EXAMPLE'; ods output CLDiffs =tc; proc print data=tc; format LowerCL 17.16 UpperCL 17.16 Difference 17.16; title "Output with many digits"; RUN; QUIT; ODS RTF close; ODS LISTING;  -  r,Ndtype)r-r-rr2r^rasarrayreplacesplitfloatrXr  tukey_hsdconfidence_intervalintrlowrIhighrK) r"datares_expect_strr3 res_expect res_tukeyconfrjlrhsigr'r'r(test_compare_sass! zTestTukeyHSD.test_compare_sasz 1 2 -8.2491590248597 -4.6 -0.9508409751403 0.0144483269098 1 3 -3.9091590248597 -0.26 3.3891590248597 0.9803107240900 2 3 0.6908409751403 4.34 7.9891590248597 0.0203311368795 z 1 2 -7.02207069748501 -3.645 -0.26792930251500 0.03371498443080 1 3 -2.68207069748500 0.695 4.07207069748500 0.85572267328807 2 3 0.59347644287720 4.34 8.08652355712281 0.02259047020620 rrrzunequal size samplec Cstj|tdd}tj|}|}|D]E\}}} } } } t|dt|d}}t |j ||f| |dt |j ||f| |dt |j ||f| |dt |j ||f| |dqdS)an vals = [24.5, 23.5, 26.4, 27.1, 29.9, 28.4, 34.2, 29.5, 32.2, 30.1, 26.1, 28.3, 24.3, 26.2, 27.8] names = {'zero', 'zero', 'zero', 'zero', 'zero', 'one', 'one', 'one', 'one', 'one', 'two', 'two', 'two', 'two', 'two'} [p,t,stats] = anova1(vals,names,"off"); [c,m,h,nms] = multcompare(stats, "CType","hsd"); rrr-rr2N)rrrrrXr rrrrrrIrrK) r"rrr3rrrrrrrrr&r'r'r(test_compare_matlabs  z TestTukeyHSD.test_compare_matlabc Csd}tj|ddddtdd}gdgd gd f}tj|}|}|D]E\}}}} } } t |d t |d }}t |j ||f| d d t |j ||f|dd t |j ||f| dd t |j||f| d d q,dS)a+ Testing against results and p-values from R: from: https://www.rdocumentation.org/packages/stats/versions/3.6.2/ topics/TukeyHSD > require(graphics) > summary(fm1 <- aov(breaks ~ tension, data = warpbreaks)) > TukeyHSD(fm1, "tension", ordered = TRUE) > plot(TukeyHSD(fm1, "tension")) Tukey multiple comparisons of means 95% family-wise confidence level factor levels have been ordered Fit: aov(formula = breaks ~ tension, data = warpbreaks) $tension z diff lwr upr p adj 2 - 3 4.722222 -4.8376022 14.28205 0.4630831 1 - 3 14.722222 5.1623978 24.28205 0.0014315 1 - 2 10.000000 0.4401756 19.55982 0.0384598 rrr,Nrr)r66r<F43rCrmr=r)rT,)r>rrr?rtr>rr6$*rrrrzr7rrzr)rrrr>r.ryr7r|rrTrrr?rr|r|rr7rrr2rgh㈵>r) r"str_resrrrrrrrrrr&r'r'r(test_compare_r s& zTestTukeyHSD.test_compare_rc Csgd}gd}gd}gd}t||||}|}tgdgdgdgdg}tgd gd gd gdg}d D]$\} } t|j| | f|| | fd dt|j| | f|| | fd dq@dS)zp Example sourced from: https://www.itl.nist.gov/div898/handbook/prc/section4/prc471.htm )皙@g@333333@gffffff@g@)g @rg333333@gffffff"@r)g @g%@g333333 @rr)rgffffff@gffffff@gffffff@g@)rrrg)g(\?rgq= ףpgp= ף?)gGz?rrg ףp= ?)rrrr)rrrgzG?)gzG@rg?g= ףp=@)g= ףp=@rrg@)rr)r*r)rrrrJr]r2N)r rrrrrrr) r"group1group2group3group4rMrlowerupperrrr'r'r(test_engineering_stat_handbook2s,  "z+TestTukeyHSD.test_engineering_stat_handbookcCstjdtjdd}tj|}|}t|j|j j tt |j |j dtt |j|jdt|j |j j tt |j dt|j |j j tt |j ddS)Nr5rrorrrr)rr8r9rr rrrrrrOdiagonalrIrK)r"rrMrr'r'r(test_rand_symmLs  zTestTukeyHSD.test_rand_symmcCsNttddtgddtjggdWddS1s wYdS)Nz...must be finite.rErr*)r-r/r)rPrQr rrrOrbr'r'r( test_no_inf_s"zTestTukeyHSD.test_no_infcCsTttddtddgddggddggdWddS1s#wYdS)Nz...must be one-dimensionalrErr*rr,)r,rr-rPrQr rrbr'r'r( test_is_1dcs$"zTestTukeyHSD.test_is_1dcCsHttddtgddggdWddS1swYdS)Nz...must be greater than onerEr*r,)r+r,r-r rbr'r'r( test_no_emptygs"zTestTukeyHSD.test_no_emptynargsrcCsFttddtjgdg|WddS1swYdS)Nz...more than 1 treatment.rErr/rr )r"rr'r'r(test_not_enough_treatmentsks"z'TestTukeyHSD.test_not_enough_treatmentscl)rFrrr*cCsVttddtgdddgddg}||WddS1s$wYdS)Nzmust be between 0 and 1rErrr+r)rPrQr rr)r"rrr'r'r(test_conf_level_invalidps "z$TestTukeyHSD.test_conf_level_invalidcCsPtj|jdd}tj|jdd}t|j|jdt|j|jddS)Nr*rr)r rdata_diff_size ttest_indrrK)r"r res_ttestr'r'r(test_2_args_ttestvszTestTukeyHSD.test_2_args_ttestN)rWrXrYdata_same_sizer extreme_size sas_same_size sas_diff_size sas_extremerGrrr matlab_sm_sizmatlab_diff_szrrrr r rrrrrr'r'r'r(rs\    % )    rc @seZdZejdgdgdfddZejdgdgdgd gd gd gd gd gdfddZddZddZ ddZ dS)TestPoissonMeansTestzc1, n1, c2, n2, p_expect)rrorrogea?)r*ror-rog c?cCs$t||||}t|j|dddS)Nr`r2r poisson_means_testrrK)r"c1n1c2n2p_expectrMr'r'r(test_paper_examplessz(TestPoissonMeansTest.test_paper_examplesz c1, n1, c2, n2, p_expect, alt, d)rTr.rTr.g{}?rr)r.r.r.r.goPF?rr)2r|rrgae?rr^)rrorTrg/V-=?rr)rrtr+rTg")?rr)r+rTrrog_'Qm~?rr)r+rTrr.g|?rr)rrr)r|g 0ݷ?rrc Cs,tj||||||d}t|j|ddddS)N)rdiffg>gؗҜ<r3rr!) r"r#r$r%r&r'altdrMr'r'r(test_fortran_authorssz)TestPoissonMeansTest.test_fortran_authorscCs0d\}}d\}}t||||}t|jddS)N)rprprr!r"count1count2nobs1nobs2rMr'r'r(test_different_resultssz+TestPoissonMeansTest.test_different_resultscCs0d\}}d\}}t||||}t|jddS)Nr rrr!r/r'r'r(test_less_than_zero_lambda_hat2sz4TestPoissonMeansTest.test_less_than_zero_lambda_hat2cCsd\}}d\}}d}tt|dtd|||Wdn1s#wYtt|dt||d|Wdn1sAwYd}tt|dtd|||Wdn1sawYtt|dt||d|Wdn1swYd}tt|dt|d||Wdn1swYtt|dt|||dWdn1swYd }tt|dtj||||dd Wdn1swYd }tt|dtjd d d d ddWddS1swYdS)Nr rz`k1` and `k2` must be integers.rErz1`k1` and `k2` must be greater than or equal to 0.rz%`n1` and `n2` must be greater than 0.z(diff must be greater than or equal to 0.)r*zAlternative must be one of ...rr*errorr)rP TypeErrorr r"rQ)r"r0r1r2r3rar'r'r(rs>$z*TestPoissonMeansTest.test_input_validationN) rWrXrYrGrrr(r.r4r5rr'r'r'r(r ~s(   r c@s`eZdZddZddZejdgdddZejdgd d d Z d d Z ddZ dS) TestBWSTestcCsltjd}|jdd\}}d}tjt|dt||g||gWdn1s,wYd}tjt|dttjg|Wdn1sMwYd}tjt|dt|gWdn1slwYd}tjt|dtj||d d Wdn1swYd }tjt|dtj||d d WddS1swYdS)Ng(?cCsRtjd}|jdd\}}tj|||d}t|j|ddt|j|ddd dS) Nr9r:rrvIh%<=rr]r}r+rr8rr r;rrIrKr"rrIrKr#r#r$rMr'r'r(test_against_Rs zTestBWSTest.test_against_R))rgD5H?gdԕ?)r`ч?gȲךX?)rrCg4)?cCsZtjd}|jdd}|jdd}tj|||d}t|j|ddt|j|dd d dS) Nl.sZrrr0rr?rr]r}r+r@rAr'r'r(test_against_R_imbalanced s  z%TestBWSTest.test_against_R_imbalancedcCstjd}|jdd\}}tjd}tjd|d}tj|||d}t|jdks,Jtjd}tjd|d}tj|||d}t|j|jtjd}tjd|d}tj|||d}t |j|jrgJdS)N/HN()r*r.rr.) n_resamplesr#rlVC A) rr8rr r&r;rnull_distributionrallclose)r"r#r#r$rrrr r'r'r( test_methods    zTestBWSTest.test_methodcCstjd}|jdd}|d}tj||dd}|jdksJt|jdt|j tj||dd}|jdks9Jt|jdtj||dd}|jdksNJt|jdt|j tj||dd}|jdkshJt|jddS) NrEr,rrrrrr) rr8rr r;rIrrKrrG)r"r#r#r$rMr'r'r(test_directions2s   zTestBWSTest.test_directionsN) rWrXrYr<r=rGrrrBrDrIrJr'r'r'r(r8s     r8). itertoolsrnumpyrr8r4rG numpy.testingrrrrrrP scipy.statsr rscipy.stats._hypotestsr r r r r rrscipy.stats._mannwhitneyurrr common_testsrscipy._lib._testutilsrscipy.stats._axis_nan_policyrrrrZrr.rwrrrr r8r'r'r'r(s@    $  7]roXtT