5[gZddlZddlZddlmZmZmZmZmZmZm Z m Z m Z m Z m Z mZmZmZddlZddlmZmZmZmZddlm ZddlmZddlmZej<j?dgdd Z ej<j?dd d gd Z!d Z"dZ#dZ$dZ%dZ&dZ'dZ(dZ)ej<j?dejTe+e,fgdzdZ-dZ.dZ/dZ0dZ1dZ2dZ3dZ4ej<j?dgd d!Z5d"Z6d#Z7d$Z8d%Z9d&Z:Gd'd(Z;Gd)d*ZGd0d1Z?Gd2d3Z@y)4N) betabinom betanbinom hypergeom nhypergeom bernoulli boltzmannskellamzipfzipfianbinomnbinomnchypergeom_fishernchypergeom_walleniusrandint)assert_almost_equal assert_equalassert_allclosesuppress_warnings)r ) root_scalar)quadzk, M, n, N, expected, rtol)) g<C0u=r/ctd}d}d}tjgdgdg|||}t|gdgdgd y) Nrrr)rr9r)r9rZrr)r9rrr9)rrr9rrRr#)rrJr)r(r)rKrJs r-test_nhypergeom_r0r[[s; A A A .., 5q!Q ?CC, 5EBr/cdtjdgddgdggd}|jdk(sJy)N)r:rG r7rFrr9rZrsize)rrvsshape)xs r-test_nhypergeom_rvs_shapereds2 r9tbTlFA 77l "" "r/c:tjjdtjdddd}tjjdtjj d}tj |ddd}t||y)Nrr]r: dr`)rUrandomseedrrbuniformppfr)rdr,ys r-test_nhypergeom_accuracyrnlsmIINN1r1bs+AIINN1 s#Aq"a$AAr/cVtjdd}d}tj|d|}|dk(}t ||tj d}d}tj|||}gd}t ||d tj|||}gd }t ||d y) Nrr9grh|?rrZr)rrr$I$I?rSg$I$I?rrrRr#)rrrrqg۶m۶m?r9r9r9)rUrVrrJrlogrr%)r'r*r,r+lamcs r-test_boltzmann_upper_boundruws "aA A a"AAvHH &&)C A aa A-HAxe, aa A+HAxe,r/cd}tj|dz}t|ddj|}tjd|dzz |dz}t ||y)Nr9)rUrVrrJrepeatr)r)r'r,r+s r-test_betabinom_a_and_b_unityrysZ A !a%A!Qq!Ayya!ea!e,H8$r/dtypesrc|\}}}|d|d|d}}}ttj|||ddy)NrrZrr'momentsgaa)rrstats)rzn_typea_typeb_typer)abs r--test_betabinom_stats_a_and_b_integers_gh18026rsC$FFFRj&)VAY!qAIOOAq!S9;NOr/cd}d}tjd}td||j|}t |||zz j|}t ||y)Ngffffff@g)\(?rZr9)rUrVrrJrr)rrr'r,r+s r-test_betabinom_bernoullirsY A A ! A!Qq!Aa!e%))!,H8$r/cTd\}}}ttj|||dy)N)g?rr9 confidencer)r,rr)rr intervalalphar)r,s r-test_issue_10317rs$KE1aEQ!>"-q#6:1a%r*1a%q)r/cRd}tjjdd}d}tj|||}t |dtj ddd}tj|||}t |dd}tj|||}t ||y) Nrrhrr`r{Gz?gGz?r9)rUrirr rlrlinspace)r,r)rdrls r-test_issue_5122rs A #B'A A ))Aq! Cb D$#A ))Aq! Ca A ))Aq! Car/c vttdtjddj ddy)Nrpirr)rr rUlogspacerlrr/r-test_issue_1603rs)tR[[T2377=qAr/cd}tjddd}ttj|d|z|ddy)Nrrr7rZr)atol)rUrrr r%)r,rds r-test_issue_5503rs6 A Ar2AEIIa1a(#D9r/zx, n, p, cdf_desired) )i,r333333?g24?)rrrgU7i(?)i0uirgέ̑ ?)ii@Brg8@ ?)i-逖rgf(G?)irgwtg?)irrgV}ֈ?)iprrgz̍?)i0rrg nC:cFttj||||yNrr r%)rdr)r, cdf_desireds r-test_issue_5503pt2rsEIIaA& 4r/cFttjddddy)NrZlJ)g-q=gB.+n?rrr/r-test_issue_5503pt3rsEIIa13IJr/cFttjddddy)Nr6gAA?g_[Cˮi8)rr r3rr/r-test_issue_6682rsFIIc2w/1FGr/cvtjgddd}dtj dg}t ||y)N)rrr9rrgG޿g1D&+)r logcdfrUinfr)result references r-test_issue_19747rs0 ]]:q# .Frvvg{3IFI&r/cRd}d}d}ttj|||dy)Nrri)rr rJ)r)r,r's r-%test_boost_divide_by_zero_issue_15101rs( A A AEIIaA&,r/cZgd}tjd||}gd}t||y)N) r9rrhriiiipr) g>l?gzo?gVrRs?g$?gVWCW?g4yB?g4 .?g.Lr?gD?)r r%r)mur% cdf_expecteds r-test_skellam_gh11474rs- 9B ++aR C PLC&r/c4eZdZdZdZdZejjdejejdddejjdd dfjZ ejj!d e d Zd Zy ) TestZipfiancd}d}tjdd}ttj|||t j||ttj |||t j ||ttj|||t j||ttj||dt j|dy)Ng@rr9msvkr|) rUrVrr rJr r%r3r~)selfrr*r's r-test_zipfian_asymptoticz#TestZipfian.test_zipfian_asymptotics   IIa  Aq!,dhhq!n= Aq!,dhhq!n= 1a+TWWQ]; aF; 1f5 7r/cd\}}d}tjd|dz}ttj|||tj|||dttj |||tj |||dttj |||tj |||dttj||dtj||ddy)N)gG?g1?r9gƠ>r#rr|)rUrVrr rJr%r3r~)ralt1agt1r*r's r-test_zipfian_continuityz#TestZipfian.test_zipfian_continuitys, d  IIaQ  AtQ/Qa1H! # AtQ/Qa1H! # 1dA. 1dA0F! # dAv> dAv>T Kr/ctjjdtjjddd}tjj ddzdz}tjjddd}gd}gd}t t j|||dd|ddd t t j|||dd|ddd y) Nrr9rwrr`rh) gĹԨ?gءk> ]x[?g>g W4?g -?g$fhn=g9[.>g`֔>gnm?) gaj?gZk=?rg?g90 ?gɸU?rTgJ_9?g2h2-?gB0?ư>r#g-C6 ?) rUrirjrrrr rJr%)rr'rr)rJr%s r-test_zipfian_RzTestZipfian.test_zipfian_R(s q II  a"  - IINN2 r !A % II  a2  .+F  Aq!,QR0#ab'E Aq!,QR0#ab'Er/rr9rrZ(za, nc tjd tj fd}tj|dz}||||}tj|}tj||}tj||z dz|}|dz} tj||z | z dz|} tj||z | z dz|dz } t t j||||t t j||||t t j||d ||| | gy) Nc\dtjd|dz|zz jS)z$Naive implementation of harmonic sumr9)rUrVsum)r)ss r-Hnsz+TestZipfian.test_zipfian_naive..HnsHs,bii1Q3'**//1 1r/c>|dks||kDryd||zz ||z S)z#Naive implementation of zipfian pmfr9rr)r'rr)rs r-pzipz,TestZipfian.test_zipfian_naive..pzipMs.1uA1a4x#a)++r/r9)weightsrZrrrmvskr|) rU vectorizerVcumsumaveragerr rJr%r~) rrr)rr'rJr%meanvarstdskewkurtosisrs @r-test_zipfian_naivezTestZipfian.test_zipfian_naiveDs!  2  2  ,  , IIacN1amiinzz!S)jj!d(Q43hzzAdFC.pmf_mean_var..fs1"2q)"2q1u-BwA~%r/cJtfdtdzDS)Nc3:K|]}||zzywrr).0rmrr's r- zYTestNCH.test_nchypergeom_fisher_naive..pmf_mean_var..P..s @.?1Q419.?sr9)rrange)r'rxlxus`r-PzFTestNCH.test_nchypergeom_fisher_naive..pmf_mean_var..Ps@eBQ.?@@@r/r9rZrUmaximumminimum)rdr*rr)rrP0P1P2rJrrrrr r s ``` @@@@r- pmf_mean_varz;TestNCH.test_nchypergeom_fisher_naive..pmf_mean_varsRBAqt$BAr"B &  A1B1B1BA$)C7Dr'R"WqL(Cc> !r/mr|v) rdr*rr)rrUrrrrJr~) rrdr*rr)rrrJrrs r-test_nchypergeom_fisher_naivez%TestNCH.test_nchypergeom_fisher_naives FFDFFDGGTVVTYYF1b!T  "  ",&aB48T3*..q!RDA3G*00B4M *00B4M r/c tjjdd}d}tjjd||}tjjd||}||z}tjd||j }tj d||z }tj||}tj|||j } tjj| j dz} dtjfdt5} | jtd ttj|||| |||| d dddtjfd } t5} | jtd ttj |||| d| |||| d dddtjfd} | | |||| }tj"| |||| }d\}}tj$||z ||tj$|zzk}|j'tj(|dz kDsJt+||||||| |D]\}}}} ||z }|||| \}}tj,||dz} | | |||| j'dksJttj"| |||| j'dy#1swYxYw#1swY]xYw)NrZrrhr9r`rcr||z }tjd||z }tj||}||fS)Nrr)r*rr)rrr r s r-rWz9TestNCH.test_nchypergeom_wallenius_naive..supports8RBAqt$BAr"Br6Mr/ct|z |\}}fd}t|||fjS)Nc2|z d|z z z zzdz SNr9r)urrr)rs r-funzCTestNCH.test_nchypergeom_wallenius_naive..mean..funs(tqAaC8|a//!33r/)bracket)rroot) r*rr)rr r rrrWs ``` @r-rz6TestNCH.test_nchypergeom_wallenius_naive..means>RBQAq)FB 4sRH5:: :r/z!invalid value encountered in mean)messageg{Gz?r#c||z }||||}|||z z}||z ||z|z z}||z|z|dz ||z||zzzz Srr) r*rr)rrrrrrs r-variancez:TestNCH.test_nchypergeom_wallenius_naive..variancesgRBQAq!AR!V A1q2vz"AQ3q5QqSRTBqD[12 2r/rr|g?ch |z  |\}} fd fd}|S)Ncjz zz z z}d||z zz zd|d|z zz z zz}|Srr)tDresrrr)rrds r- integrandzHTestNCH.test_nchypergeom_wallenius_naive..pmf..integrandsQrAvJ"!*-QqSzAo1qs8qs(;; r/cpt|}t|z }tdddd}||z|dzS)Nrr9gؗҜ<)epsrelepsabs)rr)rdrr the_integralr*rrr)s r-rz@TestNCH.test_nchypergeom_wallenius_naive..pmf..fsH"2q)"2q1u-#Iq!+0 @ Bwa00r/r) rdr*rr)rr r rr*rrWs ` ``` @@r-rJz5TestNCH.test_nchypergeom_wallenius_naive..pmfs:RBQAq)FB    1Q4Kr/)rrr)rUrirjrrbrcrrrrrfilterRuntimeWarningrrrr~rJabsrprodziprV)rrcmax_mrrr*r)r r rdrsupr$rJpmf0pmf1rr$irrWs @@r- test_nchypergeom_wallenius_naivez(TestNCH.test_nchypergeom_wallenius_naives q YY  q%e  4 YY  q%e  4 G KK1177 + ZZ1R4  ZZ2  KKBRXX . IINNAGG $Q &   ;  ; C JJ~B  D 166q"aC B1-D :!  3  3 C JJ~B  D %++Ar1aEB1% !    $1aQ"$((Ar1a8 d FF4$; $bffTl):": :uuw!++++q!ub!faeQrU;KAr1aRBQAq)FB "b1f%Aq!RA&**,r1 11 155aB1EIIKQ OM c Pd}d}d}d}tjd}tjgd}d}d}tt j ||||||d d tt j |||||d tt j|||||d y) Nr6rrwg@r)g9T;grgyztj>gt>g+k?gMeGD?gqNZVɒo?gc2n?g~ZN?gE}?Y?g?g·aQ?g]??gr)?gp{?g39 y?g m>?gJ򕴝-@gih@rR)r$rr#gdy=)rUrVarrayrrrJrr) rr(r)r*rr5rJrrs r-test_wallenius_against_mpmathz%TestNCH.test_wallenius_against_mpmaths   iimhh/0*"  -11#q!QEs" 0-221aDA5 *-11!Q4@% )r/cttd}||}|jdddgdgggdd}|jdk(sJy) Nrr6rrrw)rrTg@r_r`)rrrbrc)rrrrrds r-test_rvs_shapezTestNCH.test_rvs_shape-sN (:*?AY HHRrdRD\?H Nww,&&&r/N)rrrrUrirjrcr4rrrr*rbr)rr rr rdrrrrrrrr9r<r>rr/r-rrksDIINN1 E E   1e%  0B   1e%  0B RA Aqqww'A Aqt B Ar B B*A 299>>177 #A %D [[[24KLN4N4 DZPx()T [[[24KLN'N'r/rzmu, q, expected)rxg@ #)irg"qVUc^d}||||zz }}ttj||||y)Nrw)rr r)rqr+rar)r,s r-test_nbinom_11465rB:s3 D tBwqAFMM!Q*H5r/ctjddd}d}t|j|}|dzdk(}t |d|t |dd|z t ||dy)Nrr9rgg?r)rUrrrJrr)rdr,rJr8s r- test_gh_17146rDGsl Aq"A A A,  1 C Q!ACGQCFAaC aR!r/ceZdZejj dgdgdgdgdZejj dgdgdgd gd Zy ) TestBetaNBinomzx, n, a, b, ref)rgSArrwg:[Y))rhr6rrwgdg?)rrrrwgL&[>cLttj|||||dy)NrIr#)rrrJ)rrdr)rrrefs r-test_betanbinom_pmfz"TestBetaNBinom.test_betanbinom_pmfTs  q!Q2CeDr/z n, a, b, ref)rrr6g:o?)rr^r^g>["@)rhrrg >[?cNttj|||d|dy)Nr'r|g [n% )r/N)rrrrrrrIrKrr/r-rFrFSso [[.EBHJK E K E [[^C:?AB) B)r/rFceZdZdZy)TestZipfctjdd}|jtj}t d}|j |}|j |}t ||y)Nrrr^)rUrVrrr rJrrs r- test_gh20692zTestZipf.test_gh20692sR IIa ((288$Awhhqkhhw' S+&r/N)rrrrOrr/r-rMrMs'r/rM)Ar itertools scipy.statsrrrrrrr r r r r rrrnumpyrU numpy.testingrrrr scipy.specialrscipy.optimizerscipy.integraterrrr.r4rCrNrXr[rernruryproductintfloatrrrrrrrrrrrrrrrrrBrDrFrMrr/r-rZs """" 1& 5:;,;,56:;< , < , H4) >C#-&%#49#4#4U|nQ6F#GHPIP%F E *"B: / 2 5 5K H'- '['['|L'L'^*;79:6:6 *)*)Z''r/