+ LV-jeRt^RIt^RIt^RIt^RIt^RIt^RIHt^RIt^RI t ^RI t ^RI H t H t HtHtHtHtHtHt^RIt^RIHt^RIt^RIHtHt^RIHt^RIHt^RIH t ^R I!H"t"H#t#H$t$H%t%^RI&H't'^R I(H)t)^R I*H+t+^RI,t^R I-H.t.H/t/H0t0^R I1H2t2H3t3^RI4H5t5H6t6^RI7H8t8H9t9^RI:H;t;HH?t?]PP^8tB]PR8H;'d] P!4R8HtDRR.tEERbRltFRtGRtH!RR4tIERcRltJRtK]PPR.ERdO4R4tN!RR 4tO!R!R"4tP!R#R$4tQ!R%R&4tR!R'R(4tS!R)R*4tT!R+R,4tU!R-R.4tV!R/R04tW!R1R24tX!R3R44tY!R5R64tZ!R7R84t[!R9R:4t\!R;R<4t]!R=R>4t^!R?R@4t_!RARB4t`!RCRD4ta!RERF4tb!RGRH4tc!RIRJ4td!RKRL4te!RMRN4tf!RORP4tg!RQRR4th!RSRT4ti!RURV4tj!RWRX4tk!RYRZ4tl!R[R\4tm!R]R^4tn!R_R`4to!RaRb4tp!RcRd4tq!ReRf4tr!RgRh4ts!RiRj4tt!RkRl4tu!RmRn4tv!RoRp4tw!RqRr4tx!RsRt4ty!RuRv4tz!RwRx4t{!RyRz4t|!R{R|4t}!R}R~4t~!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t^^]EP !^^]EP ,, 4,^ERe]EP !^^]EP ,, 4,^ERf]EP !^^]EP ,, 4,^ERg]EP !^^]EP ,, 4,^iERh]EP !^^]EP ,, 4,. 3R]EP !^^e]EP ,, 4^ERi]EP !^^e]EP ,, 4,^R]EP !^^e]EP ,, 4,R, ^R]EP !^^e]EP ,, 4,R, ^iR]EP !^^e]EP ,, 4,R, . 3ERjERj]EP !^]EP ,4, ^ERk^ ]EP !^]EP ,4,, ^ERl^d]EP !^]EP ,4,, ^ERm^]EP !^]EP ,4,, ^iERnR]EP !^]EP ,4,, RERoR]EP !^]EP ,4,, RERpR]EP !^]EP ,4,, R.3.t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4tRt!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4tRt!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4t!RR4tRt]PEPm]BRR7R4tRtRtRtRtRtRtRtRt!RR4tRtRtERtERtERtERtERtERt]PPERER]P4!ERER .43ER ]P4!ER ER .43ER ]P4!ERqERr.43ER]P4!ERER.43.4ER4tERtERtERtERt]PPER.ERsO4ER4t]PPER.ERtO4ER4t]PPER.ERuO4ER4tERt]PPER.ERvO4ER4tERt!ERER ]'EP4t!ER!ER"]'EP4t!ER#ER$]'EP4t!ER%ER&]'EP4t!ER'ER(4t!ER)ER*4t]PEPm]BER+R7ER,4tER-tER.tER/tER0tER1tER2tER3tER4tER5tER6t]PPER.ERwO4ER74t]PPER8ERxERy.4ER94t]PPER:^^.ERzO3^^ .ERzO3^^.ER{O3.4ER;4tER<tER=tER>t!ER?ER@4tERAt!ERBERC4t!ERDERE4t!ERFERG4t!ERHERI4tERJt]PPERK]6]3,4ERL4tERMtERNt]PPERER|ER}.4ERO4t]PPER.ER~O4ERP4tERQtERRtERSt!ERTERU4t!ERVERW4t!ERXERY4t!ERZER[4t!ER\ER]4t!ER^ER_4Et]PPER`.ERO4ERa4EtR#(z! Test functions for stats module N)Path) assert_equalassert_array_equalassert_almost_equalassert_array_almost_equalassert_allcloseassert_assert_array_lessassert_array_max_ulp)raises) typecodesarray)rec_append_fields)special)check_random_state)IntegrationWarningquad trapezoidcumulative_trapezoid) argsreduce)_XMAX)xlogy polygammaentr)distcont invdistcont) distdiscreteinvdistdiscrete) FitDataError _argus_phi)rootfmindifferential_evolution)productdarwinx86_64 tukeylambdapearson3cHVfV RV 2p\\W4VR7R#)Nz does not have attribute msg)rhasattr)abr*s&&&u/Volumes/fast/ai/experiments/ui-tars-smoke/.venv/lib/python3.14/site-packages/scipy/stats/tests/test_distributions.py_assert_hasattrr/5s& {,QC0 GAMs#cN\\PPR4R#)f_genN)r/scipystats distributionsr0r.test_api_regressionr7;sEKK--w7r0cB\\PPP4p\ Uu.uF q^,NK pp\ Uu.uF q^,NK pp.ROpW#,V,p\\RV44pW8XgQhR#uupiuupi)c@\V4PR4'*#)<)str startswithss&r..test_distributions_submodule..JsA(9(9#(>$>r0N) rv_discrete rv_continuous rv_histogramentropy)setr3r4r5__all__rrfilter)actualdist continuousdiscreteotherexpecteds r.test_distributions_submodulerO@s **22 3F&./hdq''hJ/$01LDQLH1 E$u,H6>IJH   01s BBc a]tRt^OtoRt]P PR.RO4]P PR^^]P^ ^d.4R44t Rt Rt ]P PR.RO4R4t ]P PR .R O4R 4tR t]P PR .R!O4R 4tRt]P P$]P PR^^.4]P PR.R"O4]P PRRR.4]P PRRR.4R44444t]P P(R4t]P PRR#^.4R4tRtRtRtRtVtR#)$ TestVonMisescN\PPR4VnR#)lExqNnprandom default_rngrngselfs&r. setup_methodTestVonMises.setup_methodP99((4r0kxcRpRpV!V^^V4V!V^^V4V!V^^ V4V!V^^V4V!V^^V4V!V^^ V4R#)c\P!WVR7p\VPV4VPV^\P ,V,,44R#locscaleN)r4vonmisesrpdfrTpir]Lr?r^vms&&&& r.check_vonmises_pdf_periodicHTestVonMises.test_vonmises_periodic..check_vonmises_pdf_periodicVs>2B q 266!q255y1}2E+F Gr0c\P!WVR7p\VPV4^,VPV^\P ,V,,4^,4R#ra)r4rercdfrTrgrhs&&&& r.check_vonmises_cdf_periodicHTestVonMises.test_vonmises_periodic..check_vonmises_cdf_periodicZsK2B q A "qAIM': ;a ? Ar0Nr6)rYr]r^rkros&&& r.test_vonmises_periodic#TestVonMises.test_vonmises_periodicSsd H A $Aq!Q/#Aq!Q/#Aq"a0#Aq!Q/#Aq!Q/#Aq"a0r0c\\PP\P )4\\PP \P 4R#N)rr4 vonmises_liner,rTrgr-rXs&r.test_vonmises_line_support'TestVonMises.test_vonmises_line_supportgs8U((**RUUF3U((**BEE2r0ch\P!R4p\VP^4R4R#) ?N)r4rerrn)rYrjs& r.test_vonmises_numerical$TestVonMises.test_vonmises_numericalks" ^^C BFF1Is+r0zx, kappa, expected_pdfc^\PPW4p\WCRR7R#V瞯<rtolN)r4rerfr)rYr^kappa expected_pdfrfs&&&& r.test_vonmises_pdfTestVonMises.test_vonmises_pdfys"nn  *6r0zkappa, expected_entropyc^\PPV4p\W2RR7R#vIh%<=rN)r4rerEr)rYrexpected_entropyrEs&&& r.test_vonmises_entropy"TestVonMises.test_vonmises_entropys"..((/>r0cpRp\PPV4p\PPV4p\PPV4p\P!^^^R7P VR7p\P!^^\P ,^R7P VR7p\P!^^^\P ,\V4, ^,R7P VR7p\WVRR7\WWRR7R#)i@}rb random_stateratolN) rTrUrVr4rervsrgabsr)rYseedrng1rng2rng3rvs1rvs2rvs3s& r.test_vonmises_rvs_gh4598%TestVonMises.test_vonmises_rvs_gh4598syy$$T*yy$$T*yy$$T*~~aQa044$4G~~aQruuWA6:::M~~aQ%&ruuWSY%6q%8;;>3D3;Q //r0zx, kappa, expected_logpdfc^\PPW4p\WCRR7R#r~)r4relogpdfr)rYr^rexpected_logpdfrs&&&& r.test_vonmises_logpdf!TestVonMises.test_vonmises_logpdfs"&&q0e.Ar0c^#rr6rs&r.r@rrr0c<\P!RV,4#)y?)rTexprs&r.r@rsBFF2a4Lr0 complex_funcN) rTrUrVr4reexpectr issubdtypedtypefloatingrganglecomplexfloating)rYrWrcrlbresboundss& r.test_vonmises_expect!TestVonMises.test_vonmises_expectsNii##$78A+Bnn299+FQ}}SYY 4444!bee)^#nn299+OOQ}}SYY 4444!bee)^#nn299:PS;ASPQS sag7}}SYY(:(:;;;;r0rvs_loc rvs_shapefix_locTF fix_shapecvV'dV'd\P!R4\PP R4p\ P PVRVVR7pR^/pV'dWR&V'dW'R&\\ P V\ P P3/VBR#)zNothing to fit.rsizercrfscaleflocf0N) pytestskiprTrUrVr4rer_assert_less_or_close_loglikennlf)rYrrrrrWdatakwdss&&&&& r.test_fit_MLE_comp_optimizer(TestVonMises.test_fit_MLE_comp_optimizers  KK) *ii##$78~~!!)$G/2"4!} "L "J%ennd&+nn&9&9 C=A Cr0cp.ROp\P!V4pR\P,p\PP WR7wr4pV\P !\4P8XgQh\\PV\PP^VR7R#)g!*K?r)rrN)g!*Kg4,%OԿgZ۵?g ge_@g-8gPghn?grNg d?)r3r4rerrTrgrXs&r.test_vonmises_fit_bounds%TestVonMises.test_vonmises_fit_boundssT    !WA 6    "%%'9"5!6Q ?r0rWN)皙?re))r{Gz?g|65?)r9@g?U?)rryg"?)@rgDfI?)rrg1<)rry))rgA ?)g,eޞ?)dgBb9d22)rg G)g ))rrg EJ?)rrgoι?)rryg0,)rrgdty)rrgʹw[A)rryg)rrחA)__name__ __module__ __qualname____firstlineno__rZrmark parametrizerTrgrqrvr{rrrrrxslowrslowrrrrr__static_attributes____classdictcell__ __classdict__s@r.rQrQOs5 [[S-0 [[S1aC"891:11$3, [[5/0707 [[6;< ? < ? 0  [[8BC=C=<4 [[ [[YA/ [[[-8 [[Yu 6 [[[4-8C9790 C$ [[ O O [[Vb!W- (. ( Q9 @ @r0rQc Vf VPpVP!V3/VBp\\V4V`!V3/VBpV'g!\P !WV8g4'gQhV!WQ4pV!Wa4pWx8:g!\P !WxRR7'gQhRV9dVR,VR,8XgQhRV9dVR,VR,8XgQhR#R#)aF This utility function checks that the negative log-likelihood function (or `func`) of the result computed using dist.fit() is less than or equal to the result computed using the generic fit method. Because of normal numerical imprecision, the "equality" check is made using `np.allclose` with a relative tolerance of 1e-15. Nrrrrr)rrsupertyperTanyallclose) rJrfuncmaybe_identicalrmle_analytical numerical_optll_mle_analyticalll_numerical_opts &&&&, r.rrs |yyXXd+d+N$t*d/==M vvn56666^2M0  1 KK)% H HJ I ~b!T&\1114b!T(^333r0c RR.pVP'd3\VPPR44pV.R ORV, p\\ V\ P !\V4444p.R Op\P!\RR7;_uu_4VP!V3/VBRRR4\P!\RR7;_uu_4VP\ P.4RRR4\P!\RR7;_uu_4VP\ P.4RRR4\P!\RR7;_uu_4VPV^R7RRR4\P!\R R7;_uu_4VP!V.^.\V4^, ,O5!RRR4R# +'giELF;i +'giEL ;i +'giL;i +'giL;i +'giR#;i) rr,N3All parameters fixed. There is nothing to optimize.matchz#The data contains non-finite valueszUnknown keyword arguments:) extra_keywordzToo many positional arguments.)rf1f2r)shapeslensplitdictziprTarangerr RuntimeErrorr ValueErrornaninf TypeError)rJparamnshapes all_fixedrs& r.assert_fit_warningsr0>s X E {{{dkk'',- #HW--S #e* 567I D |& ' ' ## ' zB D D "&& D zB D D "&& D y(D E E Q' F y(H I I -SZ!^,- J I ' ' ' D D D D D F E I I Is<#G4#!H0!H=H/=-I4 H  H  H, / H?  I rJc\\4pW,p\\V4p\ VP !VP .VO5!^4\VP!VP .VO5!\P)4\ VP !VP.VO5!^4\VP!VP.VO5!\P)4R#)zgh-6235N) r%rgetattrr4rrfr,rrrTr+r-)rJdctargss& r. test_supportr5Us x.C 9D 5$ D/$/3TVV+d+bffW5/$/3TVV+d+bffW5r0c8a]tRtRtoRtRtRtRtRtVt R#) TestRandIntigcN\PPR4VnR#)l{2NrSrXs&r.rZTestRandInt.setup_methodhr\r0c\PP^^^dVPR7p\ \ P !V^84\ P !V^84,4\ \V4^d8H4\PP^^RVPR7p\ \ P!V4R8H4\ VPP\R,94\PP^^.VPR7p\ V^8V^.8,4\ \V\ P4\\V44R7\P!^^.4P^VPR7p\ VPP\R,94R#)rr AllIntegerrr)Nr 2)r4randintrrWrrTrr#shapercharr isinstance ScalarTypereprrrYvalsvals& r.test_rvsTestRandInt.test_rvsks:}}  BStxx Htby!BFF419$556D S !}}  BW488 L')* 9\#::;mmBTXX>sRx() 3 .DcODmmB#'''A )L"99:r0c\PR,p\P!V^8V^8,R^4p\PP V^^4p\ W24R#):r9$NN{Gz?)rTr_wherer4r>pmfr)rYr]outrEs& r.test_pdfTestRandInt.test_pdfxsN EE$KhhQ1r6*J:}}  Ar*!$,r0c,\P!^^$^d4p\P!V4p\P!V^8V^8.RVR, ^,R, .^4p\P P V^^4p\WC^ R7R#)r9?@decimalNr)rTlinspacefloorselectr4r>rnr)rYr^r]rOrEs& r.test_cdfTestRandInt.test_cdf~st KK2s # HHQKiib!q&)C!C%'F1C+DaH}}  Ar*!$R8r0rN) rrrrrZrGrPrZrrr s@r.r7r7gs5 ;- 99r0r7cVa]tRtRtoRtRtRtRtRtRt Rt R t R t R t VtR #) TestBinomicN\PPR4VnR#)i&8jNrSrXs&r.rZTestBinom.setup_methodr\r0c \PP^ RRVPR7p\ \ P !V^84\ P !V^ 8*4,4\ \ P!V4R8H4\ VPP\R,94\PP^ RVPR7p\ \V\44\P!^ R4P^VPR7p\ \V\ P44\ VPP\R,94R# rrr;rNr<)r4binomrrWrrTrr?rr@r rAintndarrayrDs& r.rGTestBinom.test_rvss{{r4gDHHMtqy!BFF42:$667')* 9\#::;kkoob$TXXo> 3$%kk"d#'''A 3 +, )L"99:r0c\PP^d^d^4p\PP^^d^4p\VRR^R7\VRR^R7R#)rrSrrrNr4rcrNr)rYvals1vals2s& r.test_pmfTestBinom.test_pmfsJ S!, 3*sQ7sQ7r0c\P!^R4p\P!.RO4p\ \ W"44)pVP 4p\WC4\P!^R4pVP 4p\VR4\P!^R4pVP 4p\VR4R#)r rzrrSN)?rzro) r4rcrTr sumrrErr)rYr- expected_p expected_hhs& r. test_entropyTestBinom.test_entropys KK3 XX/0 % 788 IIK& KK3  IIKQ KK3  IIKQr0c p\P!4;_uu_4\P!R\4\ \ P !^^R7P4^4\ \ P !^^R7P4^4RRR4R# +'giR#;i)errornpN) warningscatch_warnings simplefilterRuntimeWarningrr4rcrstdrXs&r. test_warns_p0TestBinom.test_warns_p0sl  $ $ & &  ! !'> : qA.335q 9 qA.224a 8' & & &s A 3$%ood#'''A 3 +, )L"99:r0c\P!R4pR\P!R4,R\P!R4,, pVP 4p\ W24\P!R4pVP 4p\ VR4\P!R4pVP 4p\ VR4R#)rorrrSNп)r4rrTlogrErr)rYr-rrrss& r.rtTestBernoulli.test_entropys OOD !266$<'$rvvd|*;; IIK& OOC  IIKQ OOC  IIKQr0rN) rrrrrZrGrtrrr s@r.rrs5 ;  r0rc&a]tRtRtoRtRtVtR#) TestBradfordicRp\P!RR4p\PP W!4p\PP W14p\ W$4R#)rN)rTlogspacer4bradfordrnrr)rYcr^rxxs& r. test_cdf_ppfTestBradford.test_cdf_ppfsG  KKR  NN  q $ ^^   %r0r6N)rrrrrrrr s@r.rrsr0rca]tRtRtoRt]P PRR R RRRRRRRR]R3. 4R4t ]P PR.RO4]P PR] PP^3] PPR3.4R44t]P PR.RO4R4t]P PR.RO4]P PR] PP ^3] PP"R3.4R 44tR tVtR #) TestCauchyicV\PPR4pVR8XgQhR#)ZbtirN)r4cauchyrfrs& r.test_pdf_no_overflow_warning'TestCauchy.test_pdf_no_overflow_warnings# LL  U #Cxxr0x, refcf\PPW).4p\W2V.RR7R#r~)r4rrr)rYr^rlogps&&& r. test_logpdfTestCauchy.test_logpdf s*||""Ar7+Cju5r0z method, sgnc<V!WA,4p\WRRR7R#r~r)rYr^rmethodsgnrzs&&&&& r. test_cdf_sfTestCauchy.test_cdf_sfs 35MU+r0c\PPV4p\W2RR7\PP V)4p\WBRR7R#+gHP)Mb@?gP)rg*6y)?g8@(J)gNgmCg'K@$[U)rgonьgs 32r))g7y1gu lY<)g^?)rro)r9rz)rr)rg)nE?)g7y1C?))gRyXsg!;v+͋)(,* EgjLҳ)$@g )r9B.)g$gCį )g :@sgg>Z)g 5g*\B))#B ;g8)& .>g /)ro)rzr)rrS) !?gqmA)g?gqJRB)rrrrrrrrrrr4rrnrrrrisfrrrr s@r.rrsd [[ # & % ' $ # # # # $ $ % ' (6 (6  [[ %  [[ ,,  A  ,,//2   ,  ,  [[X;<0<0 [[ /  [[ ,,  A  ,,  B  !", " ,r0rca]tRtRtoRtRtRtRtRtRt ] PPRR ]3R.4R 4t ] PPR.RO4R 4tR tVtR #)TestChiiTgnT9;g2 >h?@ct\PP^ ^4p\WPRR7R#)rbrrN)r4chirr CHI_SF_10_4rYr?s& r.test_sfTestChi.test_sf]s& IILLQ ++%8r0cv\PPVP^4p\ V^ RR7R#)rrrN)r4rrrrrYr^s& r.test_isfTestChi.test_isfas( IIMM$**A .2E*r0chRp^p\PPW4p\VRRR7R#)rrrNga )r4rrr)rYr^dfrs& r. test_logcdfTestChi.test_logcdfes-  !!!( 6UCr0chRp^p\PPW4p\VRRR7R#)rrrN@y`)r4rrr)rYr^rrs& r. test_logsfTestChi.test_logsfls+   &5EBr0df, ref@@c\\\PPV4VRR7R#-q=rN)rr4rrrYrrs&&&r. test_meanTestChi.test_meanys  r*Ce?gG6ƴ?g9b?|^?߈܁r?gT C?rrgm?gFZۤ?g HB?gTXek?g?bȓE?gLZ;x?gJ>)g2c}`?g>t\f?g˂ q?gg'T){?g mw?g>t\?h'T)?7-h J?u8n?A#?g:Q?rrrrgZrԢ?g߄q?g)aL?gr'O)g3c}?g>t\?g˂ ?f'T)?#n?5-h J??u8n?vL1??#?||P?g:Q?rrrrrrrr)rTrWr4 crystalballrfr rrYX calculatedrNs& r.rPTestCrystalBall.test_pdfs KKc2 &s +&&**13#*> 88 K L 59&&**13#*> 88 J K 59&&**13#3c*R 88 G H 59r0c \P!R R^4RR p\PP VRRR7p\P !.R O4p\ W2RR7\PP VRRR7p\P !.RO4p\ W2RR7\PP VRRR RR 7p\P !.RO4p\ W2RR7R#)rrTNrSrrrrr rzr rr)ghQ| )?gj,,?g d}9?gH6?g|^?gEv?glQ| )?g d}9?g|^?go[?gsX?g& Ռ?g~M+?gZR c?gtE_s?gi?gB?gEg?gvP ?gx훟?)gpr?gt\v?ge}?gݠn^?gCȘU?at\?ݠn^?~[n|D?gS4s/?2s?`qfhG?IgL?荩s?rA?IH?guk?ghi"?gP ?g.q?gs?)go?r#ge?r$gPU,?r%gjݗ ?r&g?r'g_!Y?r(ga?r)gR?r*g( ?r+g܉x?r,)rTrWr4rrnr rrs& r.rZTestCrystalBall.test_cdfs KKc2 &s +&&**13#*> 88 E F 59&&**13#*> 88 E F 59&&**13#3c*R 88 E F 59r0zx, beta, m, rootrefcb\PPWVR7p\WTRR7R#)rrrN)r4rrr)rYr^r r rootrefrs&&&&& r.rTestCrystalBall.test_sfs)   ! !!! ! 4%0r0c \P!.RO4p\P!.RO4p\P!.R O4p\PP ^W4p\ W4RR7\P!.R O4p\P!R R \P R R.4pWe, p\PP ^W4p\ WxRR7\P!\P \P \P RR.4pWe, p \PP ^W4p \ WRR7\P!\P \P \P \P R.4pWe, p \PP ^W4p \ WRR7\P!\P \P \P \P R.4pWe, p \PP ^W4p\ WRR7\P!\P \P \P \P R.4pWe, p\PP ^W4p\ VVRR7R#)z[ All values are calculated using the pdf formula and the integrate function of Mathematica MbP?rg @gn(@gMg'#@N)rrSr rr )r r r@"@)rSrSrSrSrS)gԚh@g [&@gͪ@g2dQ@gn6-@g?ޫV&̿gBgQCRgQF\jg?Җ)g('UH)rTr r4r_munprr+)rYr r expected_0th_momentcalculated_0th_momentnormr,expected_1th_momentcalculated_1th_momentexpected_2th_momentcalculated_2th_momentexpected_3th_momentcalculated_3th_momentexpected_4th_momentcalculated_4th_momentexpected_5th_momentcalculated_5th_moments& r. test_momentsTestCrystalBall.test_moments"s xx12 HH. /!hh'@A % 1 1 7 74 C+O xxHI HHh"&&)YG Hh % 1 1 7 74 C+O HHbffbffbfffh? @h % 1 1 7 74 C+O HHbffbffbffbffjA Bh % 1 1 7 74 C+O HHbffbffbffbffg> ?h % 1 1 7 74 C+O HHbffbffbffbffh? @h % 1 1 7 74 C+-BOr0c\P!^^4pVP4pR^RrTp\P!W4V4p\ \ VPV44V4p\W'RR7R#)r @ Hz>rNi) r4rrErTrWrrrfr)rYcbres1lohiNr^res2s& r.rtTestCrystalBall.test_entropyOs\   q! $zz|B KK "bffQi!,.r0r6N))(@rSrgI9)r4r3?gy <)rgjt?g +,)grr gGA*?)g>rzrTg^v'?)gd~Qr皙?g眭?)rrrrrPrZrrrrrCrtrrr s@r.rrsM<:|<:D [[ 011+PZ//r0rc8a]tRtRtoRtRtRtRtRtVt R#) TestNBinomiZcN\PPR4VnR#)l\:NrSrXs&r.rZTestNBinom.setup_method[r\r0c\PP^ RRVPR7p\ \ P !V^844\ \ P!V4R8H4\ VPP\R,94\PP^ RVPR7p\ \V\44\P!^ R4P^VPR7p\ \V\ P44\ VPP\R,94R#ra)r4nbinomrrWrrTrr?rr@r rArdrerDs& r.rGTestNBinom.test_rvs^s||DwTXXNtqy!"')* 9\#::;llr4dhh? 3$%ll2t$(((B 3 +, )L"99:r0c 0\\P!\PP RRR44\PP RRR44\PPP ^^^4p\V^4R#)iigp= ף?N) rrTrr4rXlogpmfrNr3r)rYrFs& r.rlTestNBinom.test_pmfishu||223TBC ((c48 :kk  ''1a0S!r0c\PP.RORRR7p\P!\PP .RORRR74p\ W4R#)r9333333@?rxN)r9rr9r)r4rXrrTrrnrrYrErs& r.test_logcdf_gh16159TestNBinom.test_logcdf_gh16159qsI||""<3$"?ffU\\%%lcT%BC"r0rN) rrrrrZrGrlrarrr s@r.rTrTZs5 ;##r0rTca]tRtRtoRt]P PR4t]P PR4t ]P PR4t ]P PR4t Rt Rt R tR tR tVtR #) TestGenInvGaussixcN\PPR4VnR#)lnNrSrXs&r.rZTestGenInvGauss.setup_methodyr\r0c\P!RR4p\P!VPRVPR7VP 4wr#\ VR8R4R#)ffffff@rr皙?TNr4 geninvgausskstestrrWrnrrYgigrrzs& r.test_rvs_with_mode_shift(TestGenInvGauss.test_rvs_with_mode_shift|KS)||CGGDHHGEswwOQXt$r0c\P!RR4p\P!VPRVPR7VP 4wr#\ VR8R4R#)rrrirrjTNrkrns& r.test_rvs_without_mode_shift+TestGenInvGauss.test_rvs_without_mode_shiftsKT*||CGGDHHGEswwOQXt$r0c\P!RR4p\P!VPRVPR7VP 4wr#\ VR8R4R#)r皙?rirrjTNrkrns& r.test_rvs_new_method#TestGenInvGauss.test_rvs_new_methodrrr0caV3Rlp\V!^R4R4\V!^R4R4\V!^R4R4R#)c<\P!W4pVPRSPR7p\P!W2P 4^,R8#)rirrj)r4rlrrWrmrn)rzr-rorrYs&& r. my_ks_check4TestGenInvGauss.test_rvs_p_zero..my_ks_checksF##A)C''t$((';C<<WW-a047 7r0rwTrrN)r)rYr|sf r.test_rvs_p_zeroTestGenInvGauss.test_rvs_p_zeros= 8 [C($/[C($/[C($/r0c \\P!R^4P^ RR7^\P!R^4P^ RR7, 4R#)rrN)rr4rlrrXs&r.test_rvs_negative_p#TestGenInvGauss.test_rvs_negative_psP!!$*..BT.JE%%c1-11r1MM Or0c>\PPRR ^RR7p\\P!VR^.R7^,R8R4^d\ P !R^^ 4r2\PPVR ^V, VR7p\V\P!V4PV44\PPVR ^V, VR7p\V\P!V4PV44R #) riiFPW)rrzr-rinvgaussr4333333?Tr)rzr-rdNr) r4rlrrrmrTrWrfrrrn)rYigmur^pdf_igcdf_igs& r. test_invgaussTestGenInvGauss.test_invgausss    " " " SU\\"js;A>EtLR[[q"-A""&&qDAF"&Er 2 6 6q 9:""&&qDAF"&Er 2 6 6q 9:r0c\P!.RO4p\P!R^^ 4p\V\P P VR^44R#);rrzN) rg8y8?g\}Z?ghܲ0|?gRq?)rTr rWrr4rlrf)rYvals_Rr^s& r. test_pdf_RTestGenInvGauss.test_pdf_RsJ,- KKa $ 1 1 5 5aa @Ar0c\\PP^RR4^4\\PPR^2^4^4R#)r9rzg>AN)rr4rlrfrXs&r. test_pdf_zeroTestGenInvGauss.test_pdf_zerosBU&&**1c37; U&&**3A6:r0rN)rrrrrZrrrrprtrxr~rrrrrrr s@r.rdrdxs5 [[%%  [[%%  [[%%  [[00O ; B;;r0rdca]tRtRtoRtRt]PPR.RO4R4t ]PPR.RO4R4t Rt Rt R t R tR tR tR tVtR#)TestGenHyperbolicic\P!.RO4p^^^rCpRRreW#V,WF,3p\P!VRVRV/p\P!R ^ ^ 4p \ VP V 4V^RR7R#) F^T=rzrrcrdrrrN) rg{X'=gNunƢsy>g7!t?gɫzI?g̸?g?gq 6.?gmPZz?g:PgN?)rTr r4 genhyperbolicrWrrf rYrlmbdaalphar rdeltar4ghr^s & r. test_pdf_rTestGenHyperbolic.test_pdf_rs} 1dEU{DJ/  $ r0c\P!.RO4p^^^rCpRRreW#V,WF,3p\P!VRVRV/p\P!R ^ ^ 4p \ VP V 4V^RR7R#) =R0W<=rzrrcrdư>rN) rgEB=g#Nb>g_PG>g! 1Or?gUM?g u./?gfy?gl5\?gkV?r)rTr r4rrWrrnrs & r. test_cdf_rTestGenHyperbolic.test_cdf_rs} 1dEU{DJ/  $ ),rbrrorrg:;)rrrrr9rgw=)Krrrr9rg3Z0&K<)rrrrrrrrrrrrrGrrrrrrr s@r.rrs?,>4 [[% < . . [[% @ - -=. %$"*r0rca]tRtRto]P P RR R .4R4t]P P RR R .4R4tRt Rt Vt R#) TestHypSecantiz x, referencec^\PPV4p\W2RR7R#r)r4 hypsecantrr)rYr^ referencers&&& r.rTestHypSecant.test_sfs"__   "E2r0z p, referencec^\PPV4p\W2RR7R#r)r4rrr)rYrzrr^s&&& r.rTestHypSecant.test_isfs" OO   "51r0cRpRp\PPV4p\W2RR7\PP V)4p\WBRR7R#)I@rrNg!b)r4rrrrrs& r.rTestHypSecant.test_logcdf_logsfsJ %''*%0%%qb)/r0r6N)gA0=)r=g!b;)rg9bT*@)rgTW-;@) rrrrrrrrrrrrr s@r.rrso  [[^9:<=3=3 [[^89;<2<200r0rca]tRtRtoRtRt]PPR.R O4R4t Rt Rt Rt R t R tVtR #) TestNormInvGaussic\P!.RO4p\P!.RO4p\PP V^RR7p\ W1RR7R#)p~٧>rzr,r-rrN)rg @X>g}e?g{)߳?gwH?irr9)rTr r4rrnr)rYr_cdfx_testvals_cdfs& r. test_cdf_RTestNormInvGauss.test_cdf_RsM <=,-%%))&A)=d3r0c\P!.RO4p\P!.RO4p\PP V^RR7p\ W1RR7R#)p϶>rzrrrN)rg\H55$?geСt&?go NiH?g]u>r)rTr r4rrfr)rYr_pdfrvals_pdfs& r.rTestNormInvGauss.test_pdf_RsK<=,-%%))&A)=d3r0zx, a, b, sf, rtolc\PPWV4p\WdVR7\PP WBV4p\WqVR7R#rNr4rrrr)rYr^r,r-rrr?rs&&&&&& r.test_sf_isf_mpmath#TestNormInvGauss.test_sf_isf_mpmathsH    ! !! *D)    " "2! ,4(r0cR ^.p^^.p^pRR.p\PPWV4p\WTRRR7\PP WBV4p\WaRR7R#) r!7??=rrrhrrNrr)rYr^r,r-rr?rs& r.test_sf_isf_mpmath_vectorized.TestNormInvGauss.test_sf_isf_mpmath_vectorizedsj H F  "8 9    ! !! *E6    " "2! ,4(r0c\P!^^4p\P!^^^4pVP V4pVP V4p\ WB4R#rN)r4rrTr'rrr)rYdstr^rrs& r. test_gh8718TestNormInvGauss.test_gh8718sH  A& IIaQ  VVAYggbkr0c ^Rr!\P!V^,V^,, 4pW#, V^,V^,, RV,V\P!V4,, R^^V^,,V^,, ,,V, 3p\V\PPWRR74R#)rrzr mvskmomentsN)rTrrr4r)rYr,r-gammav_statss& r. test_statsTestNormInvGauss.test_statss#11q!t $9adUAXosQw!bggen:L/M!a!Q$hAo-.68We0066qV6LMr0c^Rr!\P!.RO4p\PP W1V4p\ V\PP WAV44R#)rrzNr2rzr)rTr r4rrrrn)rYr,r-rrEs& r.test_ppfTestNormInvGauss.test_ppfsM#1-.!!%%f3 2 2 6 6t BCr0r6N))rrr9rr)rrr9rr)rrrgy5`?r)rbrrgMLC9r)rrrrrrrrrrrr rrrrr s@r.rrsY 44 [[0JK ) K ) ) NDDr0rcba]tRtRtoRtRtRtRtRtRt Rt R t R t R t R tR tVtR#)TestGeomicN\PPR4VnR#)lyNrSrXs&r.rZTestGeom.setup_methodr\r0c\PPRRVPR7p\ \ P !V^844\ \ P!V4R8H4\ VPP\R,94\PPRVPR7p\ \V\44\P!R4P^VPR7p\ \V\ P44\ VPP\R,94R#r)r4geomrrWrrTrr?rr@r rArdrerDs& r.rGTestGeom.test_rvsszz~~dtxx~Htqy!"')* 9\#::;jjnnTn9 3$%jj""1488"< 3 +, )L"99:r0cb\PP\P!R4^VP R7pVP \P8XgQh\P!V\P!\P4P84'gQhR#)#rNi) r4rrrTrrWrint64riinfoint32maxrYrs& r. test_rvs_9313TestGeom.test_rvs_9313skjjnnRVVC[qtxxnHyyBHH$$$vvcBHHRXX.2223333r0cf\PP.ROR4p\V.RO4R#)rrzNr)rzror)r4rrNrrYrEs& r.rlTestGeom.test_pmfs!zz~~i-!$(:;r0c*\P!\PP .ROR44p\PP .ROR4p\ WR^R7\PP ^^4p\VR4R#)rrzrrhrNr)rTrr4rrNr[rr)rYrjrkrFs& r. test_logpmfTestGeom.test_logpmfscuzz~~i56 !!)S15q9jj1%S#r0c\PP.ROR4p\PP.ROR4p\ .RO4p\ W4\ V^V, 4R#rrzNrrzr?)r4rrnrr rrYrEvals_sfrNs& r.rTestGeom.test_cdf_sfsKzz~~i-**-- 3/+,!$1!'1X:6r0c*\PP.ROR4p\PP.ROR4p\ .RO4p\ V\ P!V44\ V\ P!V)44R#r1) r4rrrr rrTrlog1pr4s& r.rTestGeom.test_logcdf_logsfsazz  C0**""9c2+,!$x(89!'288XI+>?r0cz\PP.ROR4p\.RO4p\ W4R#)rzNr2)rSrr )r4rrr rrYrErNs& r.rTestGeom.test_ppfs)zz~~0#6)!$1r0c^\\PPRR4RRR7R#)rrSr rN)rr4rrrXs&r.test_ppf_underflowTestGeom.test_ppf_underflow#s ue4cFr0cj\P!R4P4p\VRRR7R#)g<,Ԛ?g۩@rrN)r4rrEr)rYrss& r.test_entropy_gh18226TestGeom.test_entropy_gh18226's( JJv  & & (,59r0c\PPR4p\PP R^ VR7^8P 4'gQhR#)iKH9rN)rTrU RandomStater4rrr)rYrs& r.test_rvs_gh18372TestGeom.test_rvs_gh18372-sEyy,,Y7  u2LIAMRRTTTTr0rN)rrrrrZrGr(rlr.rrrr>rArFrrr s@r.rrsE5 ;4<7@2 G: UUr0rc,a]tRtRtoRtRtRtVtR#) TestPlancki5cz\PP.ROR4p\.RO4p\ W4R#)rrTNr)g|dy?g8'\>ggrU!>)r4planckrr rr;s& r.rTestPlanck.test_sf7s/||y"-23 "$1r0cz\PP.ROR4p\.RO4p\ W4R#)rN)r@@p@)gP.gh>gTF)r4rKrr rr;s& r.rTestPlanck.test_logsf>s,||!!"7?:;!$1r0r6N)rrrrrrrrr s@r.rIrI5s222r0rIc8a]tRtRtoRtRtRtRtRtVt R#) TestGennormiDc.ROp\PPV^4p\PPV4p\ W#4R#rNr)r4gennormrfrrrYpointspdf1pdf2s& r. test_laplaceTestGennorm.test_laplaceEs9}}  +}}  (D'r0c.ROp\PPV^4p\PPVRR7p\ W#4R#rrdNrg;f?)r4rUrfr8rrVs& r. test_normTestGennorm.test_normLs;}}  +zz~~fF~3D'r0c\P!R4p\PP R4pVP RVR7p\P !W1P4PR8gQh\P!^4pVP RVR7p\PP RVR7p\P!W44PR8gQh\P!^4pW!n VP RVR7p\PP RRVR7p\P!W54PR8gQhR#)rzrrrrdrrNg;f?) r4rUrTrUrVrrmrnpvaluerks_2samprr8)rYrJrWr rvs_laplacervs_norms& r.rGTestGennorm.test_rvsSs}}S!ii##J/hhDsh3||C*11C777}}QhhDsh3mm''T'D ~~c/66<<<}}QhhDsh3::>>t#>N~~c,33c999r0c\P!RR.RR..4p\PP R4pVP .ROVR7p\P !VR ,\P!R4P4^,R8gQh\P !VR ,\P!R4P4^,R8gQh\P !VR ,\P!R4P4^,R8gQh\P !VR ,\P!R4P4^,R8gQhR#) rzrSrrTrbrrN)rr r )rr9r9)rr9r)rrr9)rrr)r4rUrTrUrVrrmrn)rYrJrWrs& r.test_rvs_broadcasting!TestGennorm.test_rvs_broadcastinges}}sBi"b23ii##J/hhLsh;||CL%--*<*@*@A!DsJJJ||CL%--*<*@*@A!DsJJJ||CL%--*<*@*@A!DsJJJ||CL%--*<*@*@A!DsJJJr0r6N) rrrrrZr_rGrjrrr s@r.rRrRDs ((:$KKr0rRcha]tRtRto]P P RRR.4R4tRtVt R#) TestGibratiox, sfxc\\PPV4VRR7\\PP V4VRR7R#+=rN)rr4gibratrrrYr^sfxs&&&r. test_sf_isfTestGibrat.test_sf_isf}s8  *Ce< ((-qu=r0r6N)rgD.?)gJ]b<) rrrrrrrrurrr s@r.rmrmos8 [[X(A(E(GH>H>r0rmca]tRtRtoRt]P PR.R O4R4tRt Rt ]P PR.R O4R4t R t Vt R #) TestGompertzic\PP\PPR^4^4p\ VR4R#0.++N)r4gompertzrrnrrs& r.test_gompertz_accuracy#TestGompertz.test_gompertz_accuracys3 NN  u~~11&!g:0N<?)rrg:?) rrrrrrrrZrrr s@r.rrs: [[[+N+M+OPCPCr0rca]tRtRto]P P RRR.4R4t]P P R.RO4R4t]P P RRR.4]P P R .RO4]P P R R R .4]P P R R R .4R4444t Rt Rt Vt R#) TestHalfNormirnc\\PPV4VRR7\\PP V4VRR7R#r)rr4halfnormrrrss&&&r.ruTestHalfNorm.test_sf_isfs: ))!,c>**3/?r0rc\\\PPV4VRR7R#r~)rr4rrnrYr^rs&&&r.rZTestHalfNorm.test_cdfs **1-s?r0rh㈵>r rvs_scalerTF fix_scalec\PPR4p\PP WRVR7pV'd[V'dSRp\ P!\VR7;_uu_4\PPWaVR7RRR4R#/pV'dWR&V'dW(R&\\PV3/VBR R /BR# +'giR#;i rrrcrdrrrrrrNrrrT) rTrUrVr4rrrr r( halflogisticrr rYrrrrrWr error_msgrs &&&&& r.r(TestHalfNorm.test_fit_MLE_comp_optimizersii##$78~~!!gT/2"4 y%I|9==""&&t)&L>  "L &N &ennd = s 8"C C0 c\P!\4;_uu_4\PP .RO^R7RRR4R# +'giR#;irrNr)rr rr4rrrXs&r.test_fit_errorTestHalfNorm.test_fit_error s7 ]]< ( ( NN  yq  1) ( ( ( $A A" r6N)r |N?)rb Oul2;))Ww'&l7g;G<^͛7)C]r2  "L &N%e&6&6EE>= s 8"C C, c\P!\4;_uu_4\PP .RO^R7RRR4R# +'giR#;ir)rr rr4rrrXs&r.rTestHalfCauchy.test_fit_error+s9 ]]< ( (     3) ( ( (rr6N) rrrrrrrrrrrr s@r.rrs [[Yt 5 [[[4+6 [[Yu 6 [[[4-8F9776F,44r0rca]tRtRto]P P RRR.4R4t]P P R.RO4R4tRt Rt ]P P RR R .4]P P R .RO4]P P R R R.4]P P RR R.4R4444t Rt Rt VtR#)TestHalfLogistici1rc\\\PPV4VRR7R#r~)rr4rrrs&&&r.rTestHalfLogistic.test_sf9s! **--a0#EBr0q, refc\\\PPV4VRR7R#r~)rr4rrrYrrs&&&r.rTestHalfLogistic.test_isfEs! **..q13UCr0cfRpRp\PPV4p\W2RR7R#)g>@rrNguVJ)r4rrrrYr^rrs& r.rTestHalfLogistic.test_logcdfLs, $##**1-%0r0cfRpRp\PPV4p\W2RR7R#)rqrrNg+)r4rrrrYr^rrs& r.rTestHalfLogistic.test_logsfSs, $""((+/r0rrrrrTFrc\PPR4p\PP WRVR7p/pV'd[V'dSRp\ P!\VR7;_uu_4\PPWaVR7RRR4R#V'dWR&V'dW'R&\\PV3/VBR R /BR# +'giR#;ir) rTrUrVr4rrrr r(rr) rYrrrrrWrrrs &&&&& r.r,TestHalfLogistic.test_fit_MLE_comp_optimizerZsii##$78!!%%'36&8 y%I|9==""&&t)&L>  "L &N &e&8&8$ <$ <6: <>= s :"C C0 cRp\\VR7;_uu_4\PP .RO^R7RRR4R# +'giR#;i)z; Maximum likelihood estimation with 'halflogistic' requiresrrN)r9r r) assert_raisesrr4rrrYr*s& r.test_fit_bad_floc"TestHalfLogistic.test_fit_bad_flocws<L )g?gD??&&u48D//$$U6C??&&u48DHH\2M88$89LHH%ABMQw ")E F  4  2  4 }Q'7 a8H'H I #c  +HH]+E??&&u48DHHc2773< +A4%HIM  4Qw#EOO$6$6t$$G$68#EOO$7$74$H$68$EOO$7$74$H$68#EOO$6$6t$$G$68//%%d6C s+;'; roNir!r")r9pGC@FfD??&&u48D//$$U6C??&&u48DHH\2M88$89LHH%KLMQw "*=!>  4  2  4 }Q'7 a8H'H I #c  +HH]+E??&&u48DHHc2773< +CT%JKM  4??&&t$7D  ,Qw#EOO$6$6t$$G$68#EOO$7$74$H$9;$EOO$7$74$H$68#EOO$6$6t$$G$9;//%%d6C s+;'; EG fg^K[X)rTr r+rr4rrr)rYr,r-r^rNs& r.test_logcdf_gh17064!TestTruncnorm.test_logcdf_gh17064s HHrvvgwRVVGR8 9 HHbffbffaRVV4 5 HH* ++ -44Q7BQB/55qb98Dr0cl\P!R^4P^4pRp\W4R#)r gl).0S?Nr )r4rrr)rYrrs& r.test_moments_gh18634"TestTruncnorm.test_moments_gh18634s.oob!$++A. !r0rN))r9rr)333333?ffffff?gܲ l)rg>g+))rrghiɮ9?)r|rr)g0.+rr)}Ô%IrZ_2?)rrrSrr)rrbrrSrgLFu)r r!rgMFmz%?rgſ)r r rghI}?rgHNMgQ63Egr1"D)rr!gU*?bn?gnA??c3TXʿ)r rgU*ҿregnArf)rg"$7"gOxr?gp7`g'>K@)rrrrrZrrrrtrr rrrr1r9r>rTr+_truncnorm_stats_datar rCrKrNrRrUr[r^rrr s@r.rrsP5 [[[BCH CH  [[[DE H E H2 2 00&!MF-P^,  1 = < BFF       &&!            ;"FHH%:; [[V%:;D<D + 0 /H E""r0rca]tRtRto]P P R.RO4R4t]P P R.RO4R4t]P P RRR.4R4t ]P P RRR.4R 4t ]P P RRR.4R 4t ]P P RRR.4R 4t R t VtR #)TestGenLogistici x, expectedcbRp\PPW4p\WBRR7R#)rrrN)r4 genlogisticrr)rYr^rNrrs&&& r.rTestGenLogistic.test_logpdfs)   ''-U3r0rc\\\PPV4VRR7R#r)rr4rmrErs&&&r.rtTestGenLogistic.test_entropys! ))11!4cFr0rc\\\PPW4VRR7R#r)rr4rmrrs&&&&r.rTestGenLogistic.test_sfs! )),,Q2CeDr0 q, c, refc\\\PPW4VRR7R#V瞯rg@N)r!ffffff@333333@rgffffff@@r)r9gpR(g㘍R2gJޱp9g?Mf= G?r)g@rg@g@)rrg) 9g"1)rTr r4rrr) rYorangespears fruits_eatenquantileeatenrrN quantilesrM expected2s & r.test_precision2TestHypergeom.test_precision2*sxx =>D (*(5!!(eOWL( *889:AD90 !!)u_guM> ad;*s0CcD\P!^^^4pVP4p\P!RR.4p\P !\ W344)p\W$4\P!^^^4pVP4p\VR4R#)rrrorN) r4rrErTr rprrr)rYhgrsrqrrs& r.rtTestHypergeom.test_entropy=sv __Q1 % JJLXXtTl+ ffU::;; & __Q1 % JJLQr0cRpRpRpRp\PPWW44pR p\WV^R7^pRpRpRp\PPWW44pR p\WV^R7R#) rcA.Aj@rU@XrNgog_7 j$)r4rrrrYr]rryrLresultrNs& r.rTestHypergeom.test_logsfIsv    &&qQ2Fa8    &&qQ2 Fb9r0c^pRpRpRp\PPWW44pR p\WV^R7^(pRp^2pRp\PPWW44pR p\WV^R7^}pRp^pRp\PPWW44pR p\WV^R7\P !.R O4pRp^2pRp\PPWW44p\P !^R 4p\WV^R7R#) rrrrrUrrrNg)\Ug*@Vgר)rrr)r4rrrrTr fullrs& r.rTestHypergeom.test_logcdfas    ''a3Fa8    ''a3&Fb9    ''a3 Fb9 HH\ "   ''a377112Fb9r0cRpRpRp\PPWV4pW!, V,p\WE4R#)ipiPi.N)r4rrr)rYrryrLhmrms& r.test_mean_gh18511TestHypergeom.test_mean_gh18511s<    __ ! !! * UQYr0cJ^ pRp\P!^^4pRV,p\PP V^, WBV4p\P !V^84'gQh\P !\P !V4^84'gQhR#)rbrNr)rTr'r4rrrdiff)rYryrLrpopulation_sizerzs& r.test_sf_gh18506TestHypergeom.test_sf_gh18506sz   IIa q& OO  q1uo! <vva!e}}}vvbggaj1n%%%%r0rN)rrrrrZrGrrrrrtrrrrrrrrrr s@r.rrsX5 9>FL<& :0,:\   [[ & &r0rc^a]tRtRtoRt]P PR.RO4R4t]P PR.RO4R4t R t R t R t R t ]P PR RR.4R4t]P PR.RO4R4tRtVtR#) TestLoggammaicN\PPR4VnR#)lLNrSrXs&r.rZTestLoggamma.setup_methodr\r0z x, c, cdfr2rc\PPW4p\WCRR7\PP W24p\WQRR7R#r)r4loggammarnrr)rYr^rrnrzys&&&& r.rTestLoggamma.test_cdf_ppfs@ NN  q $U+ NN  s &5)r0zx, c, sfc\PPW4p\WCRR7\PP W24p\WQRR7R#r)r4rrrr)rYr^rrr?rs&&&& r.ruTestLoggamma.test_sf_isfs@ NN  a #E* NN  r %5)r0cb\PPR^4p\VRRR7R#)rr rN g@)r4rrr)rYlps& r.rTestLoggamma.test_logpdfs& ^^ " "4 +G%0r0cjRpRp\PPW4pRp\W4RR7R#)r3@rrNgI.)r4rrr)rYr^rrrs& r.rTestLoggamma.test_logcdfs/  &&q,%%0r0cjRpRp\PPW4pRp\W4RR7R#)r @rrNg9g!{1+ҷ)r4rrr)rYr^rrrs& r.rTestLoggamma.test_logsfs/  $$Q*#/r0c\P!.RO4PR^4pVF8wr#rEp\PPVRR7p\ WsWEV.^R7K: R#)rzmsvkrrUN)rzg"~jg46<@goŏr3rSg h"lxg Q?gZd;333333@rOg{/L@g??gгYҿgh|?5?r)rTr reshaper4rr)rYtablerrvarskewkurtcomputeds& r.rTestLoggamma.test_statssg wr1~  ). $AS~~++Av+>H %hs$0G./ 1).r0rc\PPVRVPR7p\P !V4P 4'gQh\PPV4p\P!\P!W#84\V44pVPRR7pVPRu;8dVP8gQhQhR#)rrr)confidence_levelrzN)r4rrrWrTisfinitermedian binomtest count_nonzeror# proportion_cirr)rYrr^medbtestcis&& r.rGTestLoggamma.test_rvss NN  qvDHH  E{{1~!!####nn##A& 0 0 93q6B  % 8vv%bgg%%%%%r0rc\\\PPV4VRR7R#r)rr4rrErs&&&r.rtTestLoggamma.test_entropy s ..q13UCr0rN))rr g;Y%?)rgRٲa6_<)gHr2gi|e?)r2g6l ?)i+rg֬CCf9)irgֺl ))rrgQOu3;)rg\Z{0)rr2geɐ?)ig{Gzd?gr0T?))rg3H໱k3@)rr~)rg8,[H^} )rg90$)rgtum\)rrrrrZrrrrrurrrrrGrtrrr s@r.rrs5 [[[EF*F* [[ZBC * C * 110 1 [[S3,/ &0 & [[X<= D = Dr0rcla]tRtRto.ROt]P PR]4R4tRt Vt R#) TestJohnsonsui r@cVwr#rEp\\PPW#V4VRR7\\PP WSV4W&R7R#vIh%,=rN)rr4 johnsonsurrrYr@r^r,r-rtols&& r.ruTestJohnsonsu.test_sf_isf# FaS**13ReD++B15qCr0r6N))rrrg?r)rrrgH';r)rrrg?x7r rrrrcasesrrrrurrr s@r.rr s3 ;E [[VU+D,Dr0rcla]tRtRto.ROt]P PR]4R4tRt Vt R#) TestJohnsonbi* r@cVwr#rEp\\PPW#V4VRR7\\PP WSV4W&R7R#)rrrN)rr4 johnsonsbrrrs&& r.ruTestJohnsonb.test_sf_isf9 rr0r6N))rrrrrG)r rrg0jA:r)rrrgQbF5rrr s@r.rr* s3 =E [[VU+D,Dr0rca]tRtRtoRtRtRtRtRtRt ] PPRRR.4R 4t R t] PPR R R .4R4tRtRtVtR#) TestLogistici@ cN\PPR4VnR#)l}$NNrSrXs&r.rZTestLogistic.setup_methodA r\r0c\P!R^4p\PP V4p\PP V4p\ W4R#rQNr)rTrWr4logisticrnrrrYr^rrs& r.rTestLogistic.test_cdf_ppfE sB KKR  NN  q ! ^^   "r0c\P!R^4p\PP V4p\PP V4p\ W4R#r)rTrWr4rrrrrs& r.ruTestLogistic.test_sf_isfK sB KKR  NN  a  ^^   "r0cRpRp\\PP^V, 4V4\\PP V4V4R#)g.funcm sDAD AA~.RVVTXN335689A>B$(a&&$(a014&&$(a014678:;$(#&(TAX./ &ennd;r0rN)g9x&?gi '?)gB"?gs^&?)rrrrrZrrur rrrrrr r%r,r6rrr s@r.rr@ s5  84)  [[03J3K3MNCNC.F  [[\D%=9 1: 1<a]tRtRtoRtRtRtRtRtRt Vt R#) TestLogseri cN\PPR4VnR#)i5- NrSrXs&r.rZTestLogser.setup_method 99((3r0c\PPRRVPR7p\P !V^84'gQh\P !V4R8XgQhVPP\R,9gQh\PPRVPR7p\V\4'gQh\P!R4P^VPR7p\V\P4'gQhVPP\R,9gQhR#r) r4logserrrWrTrr?rr@r rArdrerDs& r.rGTestLogser.test_rvs s||7Jvvdai    xx~(((zz)L"9999llt$((;#s####ll4 $$QTXX$>#rzz****yy~~rNrrYr s& r.test_pmf_small_pTestLogser.test_pmf_small_p s$ LL  Q & 7#r0c\\PPR4p\VR4R#)rgW?N)r4r>rrrAs& r.test_mean_small_pTestLogser.test_mean_small_p s" LL  d # ;'r0cR.R.R..p.ROp.RO.RO.R O.p\PPW!4p\PP WCRR7R#) rzrgYnrNwJ?)rrbrr)g5q?gD,t"?g02L9g@Y )g%Lӷ>gU5r9r9)g0gu8?g+?g㳒?g,nq?)r4r>rrTtestingr)rYrzr]rrs& r.rTestLogser.test_sf sYUTFXJ ' A@89 llooa# ""3&"9r0rN) rrrrrZrGrBrErrrr s@r.r9r9 s#4 9 $( : :r0r9c a]tRtRto]P !RR7R4t]PPR] P] P.4]PPR.RO4]PPR.RO4]PPRR R .R R .34R 4444t ]PPR ] P^3] PR3.4R 4t RtVtR#)TestGumbel_r_li functionscopec@\PPR4#rrTrUrVrXs&r.rWTestGumbel_r_l.rng yy$$T**r0rJrrzfix_loc, fix_scaleTFcVP^dW#VR7p\4pV'd V^,VR&V'd V^,VR&\W3/VBR#)rrrrN)rr%r) rYrJrrrrrWrrs &&&&&&& r.r%&TestGumbel_r_l.test_fit_comp_optimizer sXxxSg%(*v "Q;DL &]DN &d9D9r0z dist, sgncV\P!.RO4,pVPV4wrE\WBR,4\VRRR7R#)r!g@g3qtw>rrN)r!r!r!r!r!r!r!gW@)rTr rr)rYrJrzrcrds&&& r.r TestGumbel_r_l.test_fit sD :; ;XXa[  !3345DAr0r6Nrrr9rrrr)rrrrrfixturerWrrr4gumbel_rgumbel_lr%r rrr s@r.rLrL s ^^*%+&+ [[Vennenn%EF [[Y 3 [[[*5 [[1#UmeT];=:=64G :  [[[ENNA+>,1NNB+?+ABBBBr0rLc a]tRtRtoRtRt]PPR4]PPR^^.4]PPR^^.4]PPR^^.4R4444t ]PPRR ^.4]PPR^^.4]PPR^^.4]PPR ] !R R .^R 7Uu.uFpR V9gK VNK up4] P!RR7R44444t] P!RR7R4tRtRtRtVtR#uupi) TestParetoi c j \P!4;_uu_4\P!R\4\P P RRR7wrr4\ V\P4\ V\P4\ V\P4\ V\P4\P P RRR7wrr4\ V\P4\ V\P4\ V\P4\ V\P4\P P RRR7wrr4\ VR4\ V\P4\ V\P4\ V\P4\P P RRR7wrr4\ VR4\ V\P4\ V\P4\ V\P4\P P RRR7wrr4\VR4\VR4\ V\P4\ V\P4\P P RRR7wrr4\VR4\VR 4\ V\P4\ V\P4\P P R RR7wrr4\VR4\VR4\VR \P!R4,4\ V\P4\P P R RR7wrr4\VR4\VR4\VR\P!R4,4\ V\P4\P P R RR7wrr4\VR4\VR4\VR\P!R4,4\VR4RRR4R# +'giR#;i)rwrzrrrSrr rrrrrr3g2@Ng?grq@gffffff?gK~?g۶m۶m?UUUUUU?gqq?rg%I$I?gQ?gUUUUUU@grq?g88Nb@) r{r|r}r~r4paretorrTr+r*rrrYr rGr?r]s& r.rTestPareto.test_stats s $ $ & &  ! !'> :++C+@JA! BFF # BFF # BFF # BFF #++C+@JA! BFF # BFF # BFF # BFF #++C+@JA! C BFF # BFF # BFF #++C+@JA! C BFF # BFF # BFF #++C+@JA! Ay ) A2 3 BFF # BFF #++C+@JA! As # At $ BFF # BFF #++C+@JA! Ay ) A2 3 A 2777+;; < BFF #++C+@JA! Ay ) Az * A0277;3GG H BFF #++C+@JA! Ay ) A2 3 A RWWW-== > AL Mo' & & &s Q9R!! R2 cRp^pRp\PPW^VR7pW1, V,p\WE4R#)eArrbN)r4rcrr)rYr^r-rdrzrNs& r.rTestPareto.test_sf4 s<   LLOOAauO 5Ga<$r0z2ignore:invalid value encountered in double_scalarsrrrc \PPR4p\PP ^dWW$R7p\PP V^RR7^,p\PP V^RR7^,p\PP V^RR7^,pYgu;8XdTu;8Xd R8XgQhQh\PP ^dWV^,VR7p\PP V^R7wrp \V ^,VP44V^, p V P^,p \V V \P!\P!WP4, 44, 4\V ^4R#)rrr-rdrcrp= ף?)rr)rfix_b)rfbrN) rTrUrVr4rcrrrrr?rpr)rYrrrrWrshape_mle_analytical1shape_mle_analytical2shape_mle_analytical3 shape_mle_a loc_mle_a scale_mle_a data_shiftndatas&&&& r.r TestPareto.test_fit< sv ii##D)||SI$+ ?!& 0 0A$ 0 G J % 0 0AT 0 J1 M % 0 0A$ 0 G J%.%.)-. /. /.||SI%,q[ E.3ll.>.>t!.>.L+ [1_dhhj1AX   #[RVVBFF:nn6F+F$GHH JY"r0rfix_shape, fix_loc, fix_scaleTFrepeatignoreinvalidc \PPR4p\PP ^dWW'R7p/p V'dWR&V'dW)R&V'dW9R&\ \PV3/V BR#)rrjrrrN)rTrUrVr4rcrr rYrrrrrrrWrrs &&&&&&& r.r&TestPareto.test_fit_MLE_comp_optimizerY spii##D)||SI$+ ? "J "L &N%ellDADAr0c ^^^r2p\PPWV^d\PP R4R7p\ \PV4R#)rið&rN)r4rcrrTrUrVr)rYr?locationrdrs& r.test_fit_known_bad_seed"TestPareto.test_fit_known_bad_seedp sO "#Aq||S-/YY-B-B7-K M%ellD9r0c\\P4\\\PP .RO^R7\\\PP .RO^^R7R#)rrrNr)rr r!)r0r4rcrrrrXs&r.test_fit_warningsTestPareto.test_fit_warningsz sCELL)lELL$4$4iaHlELL$4$4ia r0c\PPR4p\PP R^^dVR7p\ V^4\PPV4pR#)r)rcr-rrNi~)rTrUrVr4rcrr r)rYrWrrs& r.test_negative_dataTestPareto.test_negative_data sSii##D)||DACcJ$" LL  T "r0r6N)rrrrrrrrfilterwarningsrr r#rTerrstaterrrrrrrzr s0@r.r`r` s:Nx% [[!12 [[[1a&1 [[YA/ [[[1a&1#2022 #0 [[[2q'2 [[YA/ [[[1a&1 [[<)0$q)I-)IA!&! a)I-.[["B#.203B [[":#: ##K-s  E7E7r`c @a]tRtRtoRtRtRtRtRtRt Rt R t R t ] PPR ^.RO3RRR^ ]P"!^4, ]P$.3RRRR]P"!^4,R.3R.RO3.4R 4tR tRtVtR#) TestGenparetoi chRFep\P!V4p\PP V4wr#\ VR4\ \P!V44Kg \P!R4p\PP V4wr#\W#.RR.4R#)rSrrzNrSr) rTrr4 genpareto _get_supportrrisposinfr)rYrr,r-s& r.test_abTestGenpareto.test_ab s~A 1 A??//2DA B  BKKN #  JJsO++A.S *r0cT\P!RR7p\P!^R^4p\ VP V4\P P V44\ VPV4\P PV44\ VPV4\P PV44\P!RR^ 4p\ VPV4\P PV44R#)rrrrSN) r4rrTrWrrfrrnrrrYrvr^rs& r.test_c0TestGenpareto.test_c0 s__r " KK3 #q 5;;??1#56q 5;;??1#56a%++.."34 KKB #q 5;;??1#56r0c\P!RR7p\P!^R^4p\ VP V4\P P V44\ VPV4\P PV44\ VPV4\P PV44\P!RR^ 4p\ VPV4\P PV44\ VP^4^4R#)rSrrrNr) r4rrTrWrrfuniformrnrrrrs& r.test_cm1TestGenpareto.test_cm1 s __s # KK3 #q 5==#4#4Q#78q 5==#4#4Q#78a%--"2"21"56 KKB #q 5==#4#4Q#78  ! a(r0c\P!RR7p\VP\P 4VP \P 4.RR.4\\P!VP\P 444\P!RR7p\VP\P 4VP \P 4.RR.4\\P!VP\P 444\P!RR7p\VP\P 4VP \P 4.RR.4\\P!VP\P 444R#)rrrrSNr) r4rrrfrTr+rnrisneginfrrYrs& r. test_x_infTestGenpareto.test_x_inf s __s #82r(C BIIbff-./ __r "82r(C BIIbff-./ __s #82r(C BIIbff-./r0c\P!^^ ^4pRFp\PP W4pRF6p\PP WV,4p\ W5RR7K8 \PP W4pRF6p\PP WV,4p\ WgRR7K8 K R#)r9rrN)r9rr g+)r r )rTrWr4rrfrrn)rYr^rpdf0dcpdfccdf0cdfcs& r.test_c_continuityTestGenpareto.test_c_continuity s KK2r "A??&&q,D%**1"f57&??&&q,D$**1"f57%r0c \P\P!RRRR7\P!R^^RR7R\P!RRRR7, 3,pR F^p\P P W4pR F6p\P P WV,4p\W5RR7K8 K` R# rrrbaseFendpointrSrN)rrr)rTrLrrWr4rrr)rYrrppf0rppfcs& r.test_c_continuity_ppf#TestGenpareto.test_c_continuity_ppf EE"++eT4++dArE:r{{5$S99: ;A??&&q,D%**1"f57&r0c \P\P!RRRR7\P!R^^RR7R\P!RRRR7, 3,pR F^p\P P W4pR F6p\P P WV,4p\W5RR7K8 K` R#r)rTrLrrWr4rrr)rYrrisf0risfcs& r.test_c_continuity_isf#TestGenpareto.test_c_continuity_isf rr0c x\P\P!RRRR7\P!R^^RR7R\P!RRRR7, 3,pR FLp\ \ P P\ P PW4V4VRR7KN R #) rrrrFrrSrrN)rgC]r2rgV瞯Ҽ) rTrLrrWrr4rrnr)rYrrs& r.test_cdf_ppf_roundtrip$TestGenpareto.test_cdf_ppf_roundtrip s EE"++eT4++dArE:r{{5$S99: ;/A EOO//0C0CA0I1ME +/r0cb\PPRR^^4p\VR4R#)rrNgpEȜ)r4rrrrs& r.rTestGenpareto.test_logsf s'$$T31512r0zc, expected_statscd\PPVRR7p\W2RRR7R#)rrrrrhN)r4rr)rYrexpected_statsrs&&& r.rTestGenpareto.test_stats s)&&q&&9UGr0c`\PPR4p\VRRR7R#)rgvǼ ?rrN)r4rrr)rYrGs& r.test_varTestGenpareto.test_var s# OO   %,59r0r6N)rrr rrorbgqq @gqq?g?g$I$ ?rq?gll0@r)rzUUUUUU?r9333333)rrrrrrrrrrrrrrrrrTrr*rrrrr s@r.rr s + 7 ) 0 888+3  [[ \  T2bggaj="&&1 2 ]T2771:$5v> ? " # %& H & H::r0rcPa]tRtRtoRtRtRtRtRtRt Rt R t R t Vt R #) TestPearson3i cN\PPR4VnR#)lTJNrSrXs&r.rZTestPearson3.setup_method r\r0c\PPRRVPR7p\P !V4R8XgQhVP P\R,9gQh\PPRVPR7p\V\4'gQh\P!R4P^VPR7p\V\P4'gQhVP P\R,9gQh\V4^8XgQhR#)rrAllFloatrzrNr<) r4r'rrWrTr?rr@r rArrer#rDs& r.rGTestPearson3.test_rvs s~~!!#G$((!Kxx~(((zz)J"7777nn  488 <#u%%%%nnS!%%adhh%?#rzz****yy~~:!66663x1}}r0c\PP^.RO4p\V\P !.RO4RR7\PPRR4p\V\P !R.4RR7\PP.R OR4p\V\P !.R O4RR7R#) r rrrND}i?Nrrrw)gtT?g+q?g?r r r rr9r)rgj?gdC?gs?g'M?)r4r'rfrrTr r+s& r.rPTestPearson3.test_pdf s~~!!!_5bhh'KL! #~~!!"c*bhh |44@~~!!"4c:bhh(@AGK Mr0cT\PP^.RO4p\V\P !.RO4RR7\PPRR4p\VR.RR7\PP.R OR4p\V.R ORR7R#) r rrrgE&"J?Nr)grE?gTp]0?g~?r r)giZ*J?gh+8w?gcK?g|\PP^^4p\V^4\ V\ P 4'gQh\PP^R4p\V^4\ V\ P 4'gQhR#)rrN)r4r'rrrArTnumber)rYrs& r.test_return_array_bug_11746(TestPearson3.test_return_array_bug_117462 sr&&q!,VQ&")),,,,&&q(3VQ&")),,,,r0c.ROpRp\PP\PPW!4V4p\ W24\ P !R.R..4p\ P!R^4p\ \PPWT4\PPV)V)44\ \PPWT4\PPV)V)44\ \PPWT4\PPWT)4)4R#)r!rzrNrrr ) r4r'rrnrrTr rWrfrr)rYrrrrr^s& r.test_ppf_bug_17050TestPearson3.test_ppf_bug_17050= s!nn  !3!3F!BEJ$xx$#( KKA **13**A2u5 7**13))1"te4 6**13++Au55 7r0c.ROp.ROp.ROp\\PPW!4VRR7\\PPV^4\PPV4RR7R#)rrqrN)rrzrS皙)rTrrr)g7>g~׼=gp'+7g;)rr4r'rr8)rYrr^rs& r.rTestPearson3.test_sfP sZ% ">))!2CeD))!Q/q1ANr0rN)rrrrrZrGrPrZrrrrrrr s@r.rr s65 MK ( -7& O Or0rcJa]tRtRtoRtRtRtRtRtRt Rt R t Vt R #) TestKappa4i\ c.ROpRpRFNp\PPWV4p\PPW)4p\ WE4KP R#)rrSN)rrrwrz) gffffffrr皙ɿrrrwrzrSffffff?)r4kappa4rnrrrYr^rsr]rE vals_comps& r.test_cdf_genparetoTestKappa4.test_cdf_genpareto] sL  A<<##A!,D++Ar2I D , r0c\P!R^^ 4pRp\P!R^^ 4p\PP WV4p\P P W4p\ WE4R#)rrNrr )rTrWr4rrn genextremerrs& r.test_cdf_genextremeTestKappa4.test_cdf_genextremeh s_ KKAr "  KKAr "||a($$((. (r0c\P!^^ ^ 4pRpRp\PP WV4p\P P V4p\ WE4R#r)rTrWr4rrnrrrs& r.test_cdf_exponTestKappa4.test_cdf_exponq sN KK2r "  ||a(KKOOA& (r0c\P!R^^ 4pRpRp\PP WV4p\P P V4p\ WE4R#)rrNr)rTrWr4rrnr]rrs& r.test_cdf_gumbel_rTestKappa4.test_cdf_gumbel_rz sP KKAr "  ||a(NN&&q) (r0c\P!R^^ 4pRpRp\PP WV4p\P P V4p\ WE4R#)rrNrr)rTrWr4rrnrrrs& r.test_cdf_logisticTestKappa4.test_cdf_logistic sP KKAr "  ||a(NN&&q) (r0c\P!R^^ 4pRpRp\PP WV4p\P P V4p\ WE4R#)rrSNr)rTrWr4rrnrrrs& r.test_cdf_uniformTestKappa4.test_cdf_uniform sP KKAr "  ||a(MM%%a( (r0c4\P!^^4R#r )r4rrXs&r.test_integers_ctorTestKappa4.test_integers_ctor s  Qr0r6N) rrrrrrrrrrrrrr s@r.rr\ s- -)))))r0rc>a]tRtRtoRtRtRtRtRtRt Vt R#) TestPoissoni cN\PPR4VnR#)l, NrSrXs&r.rZTestPoisson.setup_method r\r0c\P!^4p\PP .ROV4pRV^, V^,^, .p\ W#4R#)r rzNr9rr )rTrr4rrNr)rYln2rErNs& r.test_pmf_basicTestPoisson.test_pmf_basic sEffQi}}  C0QQq)'r0c\PP.RO^4p.ROp\W4\PP R^4p\ VR4R#)r9ffffff?Nr)rr9r9r9r9)r4rrNrintervalr)rYrErNrs& r.test_mu0TestPoisson.test_mu0 sE}}  A.4*==))$2Xv&r0c\PPRRVPR7p\P !V^84'gQh\P !V4R8XgQhVPP\R,9gQh\PPRVPR7p\V\4'gQh\P!R4P^VPR7p\V\P4'gQhVPP\R,9gQhR#)rzrr;rNr<) r4rrrWrTrr?rr@r rArdrerDs& r.rGTestPoisson.test_rvs s}}  7 Jvvdai    xx~(((zz)L"9999mm$((;#s####mmC $$QTXX$>#rzz****yy~~>#q#9!9!$QrT3!45HkknnQ*G G .*r0c\P!^^^ 4R,p.ROpVFapV\P!V4, p\PP WB4p\P !V4p\VR4Kc R#)r9rNNrN)r=rrirr)rTrWrr4r#rnrr )rYr^ns_xxnprobsdiffss& r.test_cdf_sqrtnTestKSTwo.test_cdf_sqrtn se KK1b !" % ,Bbggbk!BKKOOB+EGGENE eT * r0c\P!^^^ 4pRFTp\PP W4p\PP W4p\ V^V, 4KV R#r9Nr )rTrWr4r#rnrr)rYr^ryrr5s& r.rTestKSTwo.test_cdf_sf sP KK1b !)A{{q,HkknnQ*G %hG <*r0c$\P!^^^ 4pRFqpV\P!V4, p\PP W24p\PP W24p\V^V, 4Ks R#r5)rTrWrr4r#rnrr)rYr^ryr/rr5s& r.test_cdf_sf_sqrtnTestKSTwo.test_cdf_sf_sqrtn sa KK1b !)ARWWQZB{{r-HkknnR+G %hG < *r0c:\P!^^^ 4pRF|pWRV, 8,p\PP W24p^V8VR8,p\PP WB4p\ We,W5,RR7K~ R#)r9rzrrrNr )rTrWr4r#rnrrrYr^ryr/rcondrEs& r.test_ppf_of_cdfTestKSTwo.test_ppf_of_cdf ss KK1b !)As1u9B{{r-HLX_5D;;??8/D DJt < *r0c:\P!^^^ 4pRF|pWRV, 8,p\PP W24p^V8VR8,p\PP WB4p\ We,W5,RR7K~ R#)r9rzrSrrNr )rTrWr4r#rrr)rYr^ryr/vals_isfr<rEs& r.test_isf_of_sfTestKSTwo.test_isf_of_sf ss KK1b !)As1u9B{{r-HLX^4D;;>>(.D DJt < *r0cn\P!^^^ 4pRFpV\P!V4, VRV, 8,p\PP W24p^V8VR8,p\PP WB4p\We,W5,4K R#r))rTrWrr4r#rnrrr;s& r.test_ppf_of_cdf_sqrtnTestKSTwo.test_ppf_of_cdf_sqrtn s KK1b !)Abggaj.!c!e),B{{r-HLX^4D;;??8/D DJ 1 *r0cn\P!^^^ 4pRFpV\P!V4, VRV, 8,p\PP W24p^V8VR8,p\PP WB4p\We,W5,4K R#)r9rzrNr )rTrWrr4r#rrr)rYr^ryr/r5r<rEs& r.test_isf_of_sf_sqrtnTestKSTwo.test_isf_of_sf_sqrtn" s KK1b !)Abggaj.!c!e),BkknnR+GKGdN3D;;??7.D DJ 1 *r0c\P!^^^ 4R,pRFLp\PP W4p\PP W24p\ WA4KN R#)r9r,Nr )rTrWr4r#rrnr)rYr0ryr/rs& r.rTestKSTwo.test_ppf, sO Aq"%b))A*B{{r-H H ,*r0c .ROp\P!.RO4p\P!.RO.RO.RO.R O.R O.R O.4p\V4FwrEV\P!^4,\P!\P ^, V, 4,p\ PPWe4p\WsV,RR7K R#) rbrrN)rbr=rrrr)rogUUUUUU?rzrr r!)gRT>gr ?gK5 ?gn "5?gU ?gk?)g&6#>gb E>gI5?g\ ?g;?g%1?)g6>gL>g<?g䓀^?gO|6$?g qO?)g}+:>gwQ9>gH{?g]qv (?g)~/v?gmJ?)gV27J=gGk(>gWX?g Y?g{0`?gZ H?)g_R=g4>g̪?gT_@?g<e?g2?) rTr enumeraterrrgr4r#rnr)rYr-ratiosrNidxryr^rs& r.test_simard_lecuyer_table1$TestKSTwo.test_simard_lecuyer_table13 s+8988 - - - - - -    mFC"RWWRUU1WQY%77A{{q,H Hsm$ ?$r0r6N)rrrrrZrr2rr8r=rArDrGrrOrrr s@r.rr sC0( / +====22-@@r0rc2a]tRtRtoRtRtRtRtVtR#)TestZipfiU cN\PPR4VnR#)lp\#NrSrXs&r.rZTestZipf.setup_methodV r\r0c\PPRRVPR7p\P !V^84'gQh\P !V4R8XgQhVPP\R,9gQh\PPRVPR7p\V\4'gQh\P!R4P^VPR7p\V\P4'gQhVPP\R,9gQhR#)rrr;rNr<) r4zipfrrWrTrr?rr@r rArdrerDs& r.rGTestZipf.test_rvsY szz~~cdhh~Gvvdai    xx~(((zz)L"9999jjnnStxxn8#s####jjo!!!$((!;#rzz****yy~~ \0 &enndIvvt{+q01111%ennd1M &enndC &ennd-1XX__Q-?-B Dr0c\\P4\P!\ 4;_uu_4\PP .RO^R7RRR4R# +'giR#;ir)r0r4rrr rrrXs&r.test_fit_raise_errors"TestInvgauss.test_fit_raise_errors sCENN+ ]]< ( ( NN  yq  1) ( ( (s $A** A; c.R Op.ROp\PPRVR7p\W#4\PPRRR7p\ VR4\PP ^nRR7p\ VR4\PPRR4p\ VR 4\PPR R4p\ VR 4R #)4wT;?皙?rr2?gfFgr^:g9̗?rgx_;:g!?gٷ ??N)rgjdV&}?g{i\0>gDV_h?gqac X?)rrrrr)r4rrnrrr)rYrrNrI cdf_actual sf_actuals& r.rTestInvgauss.test_cdf_sf s."##CB#/X&^^''$'7  $9:NN%%cd%3  #78 ##GV4 34##Hf5 12r0cj\PPRRR7p\VR4\PP^nR4p\VR4\PP RRR7p\VR4\PP ^nR4p\VR4R#) rrrr2NgiJ.gr^ߺgpgFgpQ^L)r4rrrr)rYrrs& r.rTestInvgauss.test_logcdf_logsf s&&v$&7 12&&sD1 56$$Ut$456$$S$/01r0zmu, refc\\\PPV4VRR7R#r)rr4rrE)rYrrs&&&r.rtTestInvgauss.test_entropy s ..r2CeDr0c&\P!\PR7p\P!RRR7pRRrC\ VP V4VP V44\ VPV4VPV44\ VPV4VPV44\ VPV4VPV44\ VPV4VPV44\ VPV4VPV44\ VPV4VPV44\ VPV4VPV44R#)rrzr^rSN)r4rrTr+invgammarrrfrrnrrrr)rYrJdist0r^rzs& r.test_mu_inf_gh13666 TestInvgauss.test_mu_inf_gh13666! s~~(s#.31 A Q8 UYYq\2 A Q8 UYYq\2 1 u{{1~6 EHHQK0 UYYq\2 UYYq\2r0r)r r9r)g ףp=@g r@gʡE6@)gX9v>@g @gGz@))rg,9)r2g% !)rg40)rgK' @)rg`p @)rrrrrZrrrr rrrrrtrrr __classcell__rr s@@r.rprp s5 [[7')>?A>A>: [[7')>?A DA DD2 !3F2. [[Y)DE E E E33r0rpcna]tRtRto]P P R.RO4R4tRtRt Vt R#) TestLandaui2 namec"VR9d\P!RR^43pM\P!R^^ 43p\\PV4p\\P V4pV!V!pV!.VO^N^N5!p\ WVRR7R#)rrrr rN>rrr )rTrWr2r4landau levy_stabler)rYrr^ landau_method levy_methodrrs&& r.test_landau_levy_agreement%TestLandau.test_landau_levy_agreement3 s > ! Ca()A B2&'A d3 e//6 Q#1#a##u-r0c\\PPRR7\P3^,4\\PP ^4\P4R#)rrN)rr4rrTr*rrXs&r.rCTestLandau.test_momentsC sEU\\'''7"&&1EU\\((+RVV4r0r6N)rfrnrrr) rrrrrrrrrCrrr s@r.rr2 s4 [[V%GH .I .55r0rcaa]tRtRto]P P R.R O4]P P R.RO4R44t]P P R.RO4V3Rl4tRt Rt R t R t R t R tVtV;t#) TestLaplaceiI rrc\PPR4p\PP ^dWVR7p\P !V4p\P!\P!WE, 44\V4, p\PPV4wrx\WuRRR7\WRRR7\PPWER7wrx\WRRR7\PPWFR7wrx\Wu4V^,p\P!\P!WG, 44\V4, p\PPWGR7wrx\Wh4\PPWFR7wrx\WW4\\\PPWEVR7\\\PP\P .4\\\PP\P".4R#)rrrrrr#rN)rTrUrEr4rrrrprr#rrrrr(r)r*r+) rYrrrWrloc_mlerzrcrds &&& r.r TestLaplace.test_fitJ s ii##D)}}  cw.1!3))D/FF266$.12SY> ]]&&t, 5u=u5A]]&&t&: u5A]]&&t&> %kFF266$*-.T: ]]&&t&6 Y&]]&&t&> W" lEMM$5$5t& ( j%--"3"3bffX>j%--"3"3bffX>r0rvs_loc,rvs_scalec<\PPR4p\PP RWVR7pRp\PP V4wrg\\\P4\P`V4wrV!WgV4p V!WV4p W8g$\P!WRRR7'gQhR#R#)rrrc R\V4)\P!^V,4,^V, \P!\P!W , 44,, ,#)rr)r#rTrrpr)rcrdrs&&&r.ll3TestLaplace.test_fit_MLE_comp_optimizer..ll sMCI%8E'266"&&*<#==>? ?r0rrN) rTrUrEr4rrrr rr) rYrrrWrrrcrdloc_opt scale_optll_mlell_optrs &&& r.r'TestLaplace.test_fit_MLE_comp_optimizery sii##D)}}  d.1!3 ? ]]&&t, "4 #6#(==66:) C%G-"++f38u#F#F FF#Fr0c\P!.RO4p\PP V^R7wr#\ V^RRR7\PP V^R7wr#\ V^RRR7R#)rSrrrr#N)rSrSr rTrr=)rTr r4rrr)rYrrcrds& r.test_fit_simple_non_random_data+TestLaplace.test_fit_simple_non_random_data scxx78]]&&t!&4 qu59]]&&tA&6 QU7r0cNRp\PPV)4pVR8XgQh\PPV4pVR8XgQh\PPV4pVR8XgQh\PPV)4pVR8XgQhR#)rrrSN)r4rrnr)rYr^p0p1s& r.test_sf_cdf_extremes TestLaplace.test_sf_cdf_extremes s  ]]  r "Syy]]  q !Syy ]]  a Syy]]  qb !Syyr0c^p\PPV4p\V\P !V)4^, RR7R#)rrrN)r4rrrrTr)rYr^rzs& r.rTestLaplace.test_sf s5  MM  Q 2661":a66"++e&6&6&<& r6Nrr r!r)rzrSr)rg @) rrrrrrrrrrrrrr s@r.rr sDE [[S,/O0O [[S/2 [[\L>: [[WtUm4 [[[4-8F95;3Fr0rc a]tRtRto]P P RRR.4R4t]P!RR7R4t ]P P R.RO4]P P R.RO4]P P R .RO4]P P R ] !R R .^R 7Uu.uFpR V9gK VNK up4R4444t Rt Rt RtRtVtR#uupi) TestPowerlawi zx, a, sfc\\\PPW4VRR7R#r~)rr4powerlawr)rYr^r,rs&&&&r.rTestPowerlaw.test_sf s ))!/%@r0rMrNc@\PPR4#rQrRrXs&r.rWTestPowerlaw.rngrTr0rrrrwTFrxcJ\PP^WW7R7p\4p V'dWR&V'd8\P !VP 4\P)4V R&V'dW9R&\\PV3/V BRR/BR#))rr,rcrdrrrrrTN) r4rrr%rT nextafterrr+rr~s &&&&&&&& r.r(TestPowerlaw.test_fit_MLE_comp_optimizer s~~!!si(1"Ev "J << RVVG)rrzrrr rZr[)rrrrrrrrr\rWr#rrrrrrrs0@r.rr s [[Z1HJKAKA ^^*%+&+ [[[*=> [[Y 3 [[[*5 [[<)0$q)I-)IA!&! a)I-.<.64? <$ D4666{-s  C>C>rc@a]tRtRto]P P R.R O4R4t]P P RR R .4R4t]P P R.R O4R4t ]P P R.RO4R4t Rt Vt R #)TestPowerLogNormiSz x, c, s, refc^\\PPWV4VRR7R#r)rr4 powerlognormrrYr^rr?rs&&&&&r.rTestPowerLogNorm.test_sfds# **--aA6%Hr0z q, c, s, refc^\\PPWV4VRR7R#)g|=rN)rr4rr)rYrrr?rs&&&&&r.rTestPowerLogNorm.test_isfns# **..qQ75Ir0c^\\PPWV4VRR7R#)gt =rN)rr4rrnrs&&&&&r.rZTestPowerLogNorm.test_cdfts# **..qQ75Ir0c^\\PPWV4VRR7R#)gAfc=rN)rr4rrfrs&&&&&r.rPTestPowerLogNorm.test_pdfs# **..qQ75Ir0r6N))rrQrga!R()r2rQrg;?)r2rrgE?)MDrrg=)g?rrr)gR4 -rQrr))rrrg?)rrrgv6)rrrg+ 94|7)rrrg?))rrrgS^Ռe 9)@xDr2rg 0g;)g\)c=Hr2rg7H7) rrrrrrrrrrZrPrrr s@r.rrSs" [[^EF I F I  [[^ACEFJFJ [[^AB J B J [[^GHJ HJr0rca]tRtRto]P P R.R O4R4t]P P R.R O4R4t]P P R.R O4R4t Rt Vt R#) TestPowerNormirc\\\PPW4VRR7R#r)rr4 powernormrrs&&&&r.rTestPowerNorm.test_sfs **10#EBr0rsc\\\PPW4VRR7R#rrN)rr4rrrws&&&&r.rTestPowerNorm.test_isfs ++A13UCr0c\\\PPW4VRR7R#r)rr4rrnrs&&&&r.rZTestPowerNorm.test_cdfs ++A13UCr0r6N))rZrg|<)rQr ge()rrg%!@6)rrgT;-))rrQgĿ)rHrg!X)r rg缂U)rjrgͫ1@)rr gD@)Jz5rgU!zYr@))irZg.u T9)r rZ?)rrZgsK-)rrgk8>)rrg^u1-) rrrrrrrrrrZrrr s@r.rrs [[[CD C D C [[[=>D>D [[[BC D C Dr0rca]tRtRtoRtRtRtRtRt] PPRR R .4R4t R t VtR #) TestInvGammaic*\P!4;_uu_4\P!R\4\P P RRR7p.R Op\ W4.R Op\P P VRR7p.R O\PRR.\PRR.\P\PR.3p\W4FwrE\WE4K R R R 4R # +'giR #;i) rwg(\O3@rr\g^Yb?g,r?gOD@g2*Z@gFN8@N)g~?gCX&?g*7gS?gBqFq@)rR@gffffff@)rg Ey?gў3oz?) r{r|r}r~r4rrrTr+r*r&r)rYrrNr,r^rs& r.test_invgamma_inf_gh_1866&TestInvGamma.test_invgamma_inf_gh_1866s $ $ & &  ! !'> :>>''%'@D%H D +A>>''!V'!5I9>!#/>>s )B8 BcR\PP^ RRR7wrr4\\P!V44\\P!V44\\P!V44\\P!V4'*4R#)r@rrN)r4r$rrTrrds& r.test_f_momentsTestF.test_f_momentssaWW]]2sF]; a A A ABKKN"#r0c\P!4;_uu_4\P!R\4\P P ^ .^,.RORR7RRR4R# +'giR#;i)rwrdfndfdrN)r rrr)r{r|r}r~r4r$rXs&r.test_moments_warningsTestF.test_moments_warningssL  $ $ & &  ! !'> : GGMMrd1f,M G' & & &s AA00 B c\P!^.^ ..4p\P!^ ^ .4p\PPWRR7wr4rVW"^, , .^,p\ W74^V^,,W,^, ,V, V^, ^,, V^, , p\ WH4^V,V,^, \P !^V^, ,4,V^, \P !WV,^, ,4,, p \ WY4^ V^V,^, ,W,^, ,V^, V^, ^,,,,p W^, ,V^, ,W,^, ,p W, p \ Wl4R#)r!rr1N)rTr r4r$rr) rYr2r3r rGr?r]rhv2s2k2numk2denk2s & r.test_stats_broadcastTestF.test_stats_broadcast"sYhhbT{#hhBx WW]]sV]D a1Wo q  aZ39q= )C /37Q, >#' Jus{Q"''!S1W+"66a2773c A #6779cQsURZ(CIM:Qw37Q,./0Qw37+sy1}= ]r0r6N) rrrrrr.r4r<rrr s@r.r"r" s0$H r0r"c a]tRtRtoRtRtRt]PPRRR.4R4t ]PPR.RO4]PPR^^.^^.^^..RO.RO.^^.^^..^.^...4R 44t ]PPR^^.^^.^^..RO.RO.^^.^^..^.^...4R 4t R t ]PPR .RO4R 4tRtVtR#) TestStudentTi4c`\\PP^^.4RR.4R#)rg?g?N)rr4rrrXs&r.test_rvgeneric_stdTestStudentT.test_rvgeneric_std5s"!%''++q!f"5 J7OPr0c\\PP^RR7\P\P \P \P 34\\PPRRR7R\P\P \P 34\\PP^RR7R\P\P \P 34\\PPRRR7RR \P \P34\\PP^RR7\P \P34\\PPRRR7R\P34\\PP^RR7R\P34\ \PPRRR7R R R 7R #)rr)rr)\(?rgGz@r[gGz@g ףp= @r rNg i@)rgq@)rr4rrTr+r*rrXs&r.test_moments_tTestStudentT.test_moments_t9s]UWW]]a]8VVRVVRVVRVV4 6UWW]]dF];"&&"&&"&&1 3UWW]]a]8"&&"&&"&&1 3UWW]]dF];/266266: <UWW]]a]68HIUWW]]dD]9C=IUWW]]a]6bff F t <>U" $r0cl.ROp.ROp\\PPV4VRR7R#)rrrN)rr rr)g?@gTM]?ggY?U?rr4rrE)rYrrNs& r.test_t_entropyTestStudentT.test_t_entropyIs' ;+XEBr0zv, refc\\\PPV4VRR7R#rrI)rYrGrs&&&r.test_t_extreme_entropy#TestStudentT.test_t_extreme_entropyPs *Ceggoo!o4ZZ''AQ'7 ww"[/qJ S_j1S%|4r0c.ROp.ROp.ROp.ROp\\PPW4VRR7\\PP W4VRR7R#)rr rN)rrrbr)rd~QJrr)gZ_2g 9g!ǟuIgdg)g73?r9gմFhAW;gm0_?)rr4rrrf)rYr^r logpdf_refpdf_refs& r.test_logpdf_pdfTestStudentT.test_logpdf_pdfsO  #@ >q-zF A*G%@r0 x, df, refc\PPW4p\WCRR7\PP V)V4p\WSRR7R#r)r4rrrr)rYr^rrrrs&&&& r.rTestStudentT.test_logcdf_logsfs@ &%0 qb"%/r0r6N)rrH)rrc)rfrrnrrr)r9rr9)rrr))gRrguSaG)r9rr)R@rgGn`?)rrrrrArErJrrrrMr^rdrjrrrr s@r.r?r?4sWQ$ C [[X68:; =; = [[Z*<= [[\QFQFQF-6 ,B./VaV,<./S1#J,89C 9= C  [[\QFQFQF-6 ,B./VaV,<./S1#J,895 95& A [[\CD0 D0r0r?cha]tRtRtoRtRtRtRtRtRt Rt R t R t R t R tR tRtVtR#)TestRvDiscreteicN\PPR4VnR#)i}NrSrXs&r.rZTestRvDiscrete.setup_methodr<r0c.ROp.R OpRp\P!RW3R7pVPW0PR7p\ V\ P 4'gQh\W4F9wrg\\WV8H4\V4, V, 4R8dK9Qh VPVPR7p\ P!\V4\ P4'gQhR#) rrsample)rvaluesrrjrN)rr9rr r!r)rrrrrr)r4rBrrWrArTrer&rrprrrinteger)rYstates probabilitysamplesrr^r?rzs& r.rGTestRvDiscrete.test_rvss$4    8V4I J EEwXXE 6!RZZ((((,DAs16{5>1A56= ==- EEtxxE (}}T!Wbjj1111r0c>\P!.RO4p\P!.ROV3R7p\ \ W44)pVP 4p\WC4\P!.RO.RO3R7pVP 4p\VR4R#)rorvrN)ror_rr)rSr9r9) rTr r4rBrprrErr)rYr.rzrrrss& r.rtTestRvDiscrete.test_entropysu*+   i%7 8%-.. IIK&   i%= > IIKQr0c.R Op.R Op\P!W3R7pRR.R^..p\VPV4RR.RR..R R 7R #)rrzrrwr}rSr3r rr rNrrzrrw)r4rBrrN)rYxkpkrr^s& r.rlTestRvDiscrete.test_pmfs[     rh /"X!W q sc$*/ 1r0c.ROp.ROp\P!W3R7p.ROp.ROp\VPV4VRR7\VUu.uFqcPV4NK upVRR7R#uupi) rr}r rNrr)r rSrRrrr rr)r9rzrzrzrrrr)r4rBrrn)rYrrrx_valuesrNrs& r.rZTestRvDiscrete.test_cdfsf     rh /55x((? h7hh7 u .7A4c.ROp.ROp\P!W3R7p.ROp.ROp\VPV4VRR7\VUu.uFqcPV4NK upVRR7R#uupi) rr}r rNrr)rrzr`rrrS)rrr r rr)r4rBrr)rYrrrq_valuesrNrs& r.rTestRvDiscrete.test_ppfsf     rh /0%x((? H5HqH5 u .5rc .RO.RO3p\P!VR7p\VPVP VP RR4R,4VP R,4R#)rr}Nrr,)rr rr)rrwrrrr)r4rBrrrnr)rYrErs& r.test_cdf_ppf_next TestRvDiscrete.test_cdf_ppf_nextsS!:;   d +266"&&s"4t";<559 &r0cT\P!^ 4PR4p\P!.RO.RO.RO.4p\P !W3R7p\ VP4\P!VPVP,4RR7R#)rr}r rN)r!r)rrrrj)rrrjrj) rTr'rr r4rBrrrprrrYrrrs& r.test_multidimension"TestRvDiscrete.test_multidimensionsq YYr] " "6 * XX---/0  rh / RVVBEEBEEM%:Gr0c.ROpRR.p\\\P3/\ W3R7B.ROp\\\P3/\ W3R7B.ROp.ROp\\\P3/\ W3R7B.ROp.ROp\\\P3/\ W3R7B^^.pRR.p\\\P3/\ W3R7BR#)rrzr}Nr)rz333333?gffffff濩rr r!rr)rrrrr)rr)r4rBr%rYrrs& r.test_bad_inputTestRvDiscrete.test_bad_inputs 3Zj%"3"3MtB87LM j%"3"3MtB87LM  j%"3"3MtB87LM  'j%"3"3MtB87LMV3Zj%"3"3MtB87LMr0cr\P!^4PR4\P!RR4r!\ \ \ P3/\W3R7B\P!^4PR4\P!RR4r!\ \ \ P3/\W3R7B\P!^4PR4\P!RR4r!\\ P!W3R7P^4R4R#)rr}Nr r r r!gUUUUUU?r!r ) rTr'rrrr)r4rBr%rrNrs& r.test_shape_rv_sample#TestRvDiscrete.test_shape_rv_samples1%%f-rwwvs/CBj%"3"3MtB87LM1%%f-rwwvs/CBj%"3"3MtB87LM1%%f-rwwvs/CBU&&rh7;;A>Dr0c.ROp.ROp\P!W3R7p\VP4\P !VP VP,4RR7R#)rr}r rN)rr rrrr)rrwrwrwrwr)r4rBrrrTrprrrs& r. test_expect1TestRvDiscrete.test_expect1sE  +   rh / RVVBEEBEEM%:Gr0c .RNRNRNRNRNRNRNRNRNR NR NR NR NR NRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNRNR NR!NR"NR#NR$NR%NR&NR'NR(NR)NR*NR+NR,NR-NR.Np.R/NR0NR1NR2NR0NR0NR3NR4NR0NR0NR0NR5NR6NR7NR8NR9NR:NR;NRNR?NR@NR0NRANRBNRCNRDNRENRFNRGNRHNR0NR0NRINRJNR6NRKNR0NRLNR0NRMNR0NR0NRNNRONR0Np\P!W3RP7p\VP4VP 4RQRR7\VP4\ RS\ W444RQRR7\VPRT4\ RU\ W444RQRR7RV#)Wgi@gr@gy@g@@g@g@@g @rg0@g@gP@g@gp@g@g@g @g@rNgh@g0@g@g@g@gP@g@g@g@rOg8@g@gȩ@g@gX@g @g@g@gx@g@@g@gh@g̰@g0@g@g@g\@g@g-C6:?rgF%uk?g:vz?g.nr?g|гY?g@6*?g4U0*C?g@ +?gC6z?gnr?g ?gOjM?g n?g1w-!?g:pΈ?gףp= ?g_Le?g^)p?g~jtx?gJ4q?gGzt?gJ4q?gMb@?g@Psׂ?g@^)ˀ?gN@s?g1w-!_?gൄ|г?gH}}?g F%u?g^)ˀ?gv?r}r rc36"TFwrW,xK R#5irtr6.0rGws& r. .TestRvDiscrete.test_expect2..Es9jdaAEEjscV^,#r r6rs&r.r@-TestRvDiscrete.test_expect2..HAqDr0c3F"TFwrV^,V,xK R#5i)r Nr6rs& r.rrIs<AqD1HHs!N)r4rBrrrrpr&)rYrpyrs& r. test_expect2TestRvDiscrete.test_expect2&s =U =E =5 =% = = =u =e =V =  = =# =%+ =-3 =5; ==C =EK =  = =# =%+ =-3 =5; ==C =EK = = =$ =&, =.4 =6< =>D =FL =  =  =$  =&,  =.4  =6<  =>D  =FL  =   =   = $  = &,  = .4  = 6<  =?f?c?6?#7??c?"?$*?,/?14?69?$?&;?#?%9?$ ?&: ? " ? $7 ? 9M ? # ? %8 ? :O ??)?+??$?&;?$?&;?#?%9???.?0E?$?&9?;>?#?%(?*>?@C?EH?#?%9?;>?   qg .  RWWYU; 9c!j99 G  .1<Q<<5 Jr0rN)rrrrrZrGrtrlrZrrrrrrrrrr s@r.rqrqsL4 2  1 . .&HN( EH#J#Jr0rqc,a]tRtRtoRtRtRtVtR#)TestSkewCauchyiLc\P!R^^d4p\\PP V^R7\P P V44\\PPV^R7\P PV44\\PPV^R7\P PV44R#rrZNr) rTrWrr4 skewcauchyrfrrnrrs& r. test_cauchyTestSkewCauchy.test_cauchyMs KKAs #!%"2"2"6"6qA"6">"',,"2"21"5 7!%"2"2"6"6qA"6">"',,"2"21"5 7!%"2"2"6"6qA"6">"',,"2"21"5 7r0c\PP^4pVP^ 4^,^, pVP^ 4^ ,^, p.ROp.ROp\ \ P PW24V4\ \ P PW24V4\ \ P PWR4V4R#)r9N) g&45?gQʶ?g-4[u?g a2?g Rn?gѦ}cي?) gN@Q?gPPz@ ?g?g; 7?glg?g!?%?gQҙ?gC"?go?g@Sg?) rTrUrErandrr4rrfrnr)rYrWr,r^rfrns& r.test_skewcauchy_R TestSkewCauchy.test_skewcauchy_RVs&ii##A& HHRL1 q  HHRL2  !%$ ((,,Q2C8((,,Q2C8((,,S4a8r0r6N)rrrrrrrrr s@r.rrLs7 9 9r0rca]tRtRtoRt]P R4t]PPR.RO4R4t Rt Vt R#) TestJFSkewTiyc^;rV^,p.ROp.ROp\P!W4p\P!V4p\VP V4VP V44\VP V4VP V44\VP V4VP V44\VPR4VPR44R#)rrN)rrrSr)rrrorrrS)r4 jf_skew_trrrfrnr)rYr,r-rr^rjfrs& r.test_compare_tTestJFSkewT.test_compare_tzs  U ! - __Q " GGBKq 1558,q 1558,q 1558,(!''&/:r0c \P!\\4PR, 4p\P P VRR7#)atSample data points computed using the `ST5` distribution from the GAMLSS package in R. The pdf has been calculated for (a,b)=(2,3), (a,b)=(8,4), and (a,b)=(12,13) for x in `np.linspace(-10, 10, 41)`. N.B. the `ST5` distribution in R uses an alternative parameterization in terms of nu and tau, where: - nu = (a - b) / (a * b * (a + b)) ** 0.5 - tau = 2 / (a + b) z"data/jf_skew_t_gamlss_pdf_data.npyz x,pdf,a,bnames)rTloadr__file__parentrec fromarraysr$s& r.gamlss_pdf_dataTestJFSkewT.gamlss_pdf_datasAww N ! !$H H vv  [ 99r0za,bc VVR,V8HVR,V8H,,pVR,VR,re\V\P!W#4PV4RR7R#)zCompare the pdf with a table of reference values. The table of reference values was produced using R, where the Jones and Faddy skew t distribution is available in the GAMLSS package as `ST5`. r,r-r^rfrrN)rr4rrf)rYrr,r-rr^rfs&&&& r.test_compare_with_gamlss_r&TestJFSkewT.test_compare_with_gamlss_rs[  S !Q &?3+?1+D E cDK3U__Q266q9Fr0r6N)r)rr)r ) rrrrrrr\rrrrrrr s@r.rrysK;  ^^ : : [[U$>? G@ Gr0rriQi'i0ybilHO1Zi3rii Ni(rFica]tRtRtoRtRtRtRtRtRt Rt ] PPR ]4R 4tR tR tR tRtVtR#) TestSkewNormic&\R4VnR#)rN)rrWrXs&r.rZTestSkewNorm.setup_methods%d+r0c\P!R^^d4p\\PP V^R7\P P V44R#r)rTrWrr4skewnormrfr8rs& r. test_normalTestSkewNorm.test_normalsB KKAs #!%.."4"4Q!"4"<"'**.."3 5r0c\R4pRp\PPRW!R7p\ W#P 4\PPRW!R7p\ W#P 4R#)rr)r,rrN)r!rrr )rr4rrrr?)rYrWr?r^s& r.rGTestSkewNorm.test_rvss[ & NN  E  DUGG$ NN  %  BUGG$r0c\R4p\PP^\ R4^^VR7p\ P !V4\ P!V4\P!V4\P!V4.p\PP^^^RR7p\WC^R7\PPR\ R4^^VR7p\ P !V4\ P!V4\P!V4\P!V4.p\PPR^^RR7p\WC^R7R#)rr)r,rrcrdrr)r,rcrdrrUNr) rr4rrrdrTrrrkurtosisr)rYrWrrNrs& r.rCTestSkewNorm.test_momentss & NN  Sq,/  1GGAJq 5::a=%..:KL>>''!!V'L!(a@ NN  #c(,/  1GGAJq 5::a=%..:KL>>''"!1f'M!(a@r0c.RO.RO.RO.RO.RO.R O.R O.pVF2wr#p\PPW#4p\WTRR7K4 R#) rrrN)rrgEAV)rrgV|@$a)rr9g0rY)rrzZk)r5rr)r rļ)r gcr)r4rrr)rY logpdfvalsr^r, logpdfvalrs& r.test_pdf_large_xTestSkewNorm.test_pdf_large_xsN 2 3 0 2 4 5 5  *OA)>>((.D D$ 7 *r0c\PP.ROR4p\V\P !^4RR7\PP^R4p\VRRR7R#)rbr rrrSN)rbrQrr)r4rrnrrTonesrs& r.test_cdf_large_xTestSkewNorm.test_cdf_large_xsP NN  |R 02771:E2 NN  r3 '3U+r0c.RO.RO.RO.RO.RO.pVFawr#p\PPW#4p\WTRR7\PP V)V)4p\WTRR7Kc R#) rrrN)rXrgzbLe9)rr gn'/2;)r rg::)rgrg|<)rrr)r4rrnrr)rYcdfvalsr^r,cdfvalrzs& r.test_cdf_sf_small_values%TestSkewNorm.test_cdf_sf_small_values sl . + , . /  $LA&""1(A AD 1!!1"qb)A AD 1 $r0z a, momentsc\V^R7F1wr4\PPW14p\ WTRR7K3 R#)r)startr rN)rLr4rrr)rYr,rorderrNrs&&& r.test_noncentral_moments$TestSkewNorm.test_noncentral_momentss3(:OE..''1C C 6 ;r0c \PPR4pR RRrCp\P!W#V4pVP ^dVR7p\PP VR ^R7wrxp \PP VR^R7wrp Yu;8Xd ^8XgQhQhWz8wgQh\PP V^^RR7wrpV^8XgQh\P!WV4pVPRR7p\P!V4\P!V43p\VV4\PP ^^dVR7p\PP V4p\P!\P!V44'gQh\PP VRR 7wppp\P!V4'gQh\P!V4\P!V4pp\VVV\P!^\P , 4,,4\VV^,^^\P , , ,4\PP VR R 7wppp\PP V)R R 7wppp\VVV.V)V)V.4\PP VRR 7wpp p!\PP V)RR 7wp"p#p$\V"V#V$.V)V )V!.4R #)l# kerrzrrmmrrmsrrmleNr rg)rTrUrVr4rrrrrrrcrrisinfrrrg)%rYrWr,rcrdrJra2loc2scale2a3loc3scale3a4loc4scale4dist4rrra5loc5scale5r rGa6ploc6pscale6pa6mloc6mscale6ma7ploc7pscale7pa7mloc7mscale7ms%& r.r TestSkewNorm.test_fit#sii##$78C~~ae,hhCch2!>>--c4a-@& >>--c4a-@& q     xx!>>--c1Qt-L&{{r0kk$k'ggclEJJsO+S!||#>nn  &vvbkk#&''''!>>--d4-@D&xx|||wwt}bffTl14&2771RUU7+;";;<619AI 67$nn00U0CUG#nn00#e0DUGeW-ufg/FG#nn00T0BUG#nn00#d0CUGeW-ufg/FGr0c~a \P!.R O^ ^.,,^.,4p\PP V4p\PP W!4p\PP VRR7p\PP WA4pW5R, 8gQh\P PR4o RRRRRR/pV 3Rlp\P!\PWVR 7p\PPW(PR R 7R #)rT)superfitrzly NDr,rcrdc<\WSR7#)r)r")funrrWs&&r. optimizer0TestSkewNorm.test_fit_gh19332..optimizerhs)#3? ?r0)rrrN)rrr)rr)rrb)rrb) rTr r4rrrrUrVrIrparams) rYr^rr params_superrrr fit_resultrWs & @r.test_fit_gh19332TestSkewNorm.test_fit_gh19332Ss HH*R1#X5; <##A&nn!!&, ~~))!d); nn!!,23Yii##$78wy';G @YYu~~qIN  ""6+<+<4"Hr0c^\\PPRR4RRR7R#)rrgi}#%?rrN)rr4rrrXs&r.rTestSkewNorm.test_ppfns" **457HuUr0rN)rrrrrZrrGrCrrrrrr_skewnorm_noncentral_momentsrr rrrrr s@r.rrsi,5 % A8 ,2" [[\+GH7I7 .H`I6VVr0rc8a]tRtRtoRtRtRtRtRtVt R#) TestExponiucX\\PP^4^4R#r)rr4rrfrXs&r. test_zeroTestExpon.test_zerovsU[[__Q'+r0c\\PPR4R4\\PP \PP ^(44^(4R#)rN)rr4rrnrrrXs&r. test_tailTestExpon.test_tailys=U[[__U+U3U[[__U[[^^B%78"=r0c\P!RRRRR\P.4p\\\ P PV4R#o_?46@yX5ͻ@j+?UX @N)rTr r*rr)r4rrrs& r.test_nan_raises_errorTestExpon.test_nan_raises_error}5 HHfffffbffE Fj%++//15r0c\P!RRRRR\P.4p\\\ P PV4R#r,)rTr r+rr)r4rrrs& r.test_inf_raises_errorTestExpon.test_inf_raises_errorr4r0r6N) rrrrr&r)r2r6rrr s@r.r$r$us,>6 66r0r$cza]tRtRtoRtRtRt]PPR^^.4R4t Rt Vt R#) TestNormic\P!RRRRR\P.4p\\\ P PV4R#r,)rTr r*rr)r4r8rrs& r.r2TestNorm.test_nan_raises_error5 HHfffffbffE Fj%**..!4r0c\P!RRRRR\P.4p\\\ P PV4R#r,)rTr r+rr)r4r8rrs& r.r6TestNorm.test_inf_raises_errorr<r0cd.ROp\\\PPVRR7R#)rshrimp)plateNr)rr,r4r8rrs& r.test_bad_keyword_argTestNorm.test_bad_keyword_args i(Cr0rccRp\PP^ V,^ V,VR7p\W2RR7\PP^ V,)^ V,)V)R7p\W2RR7R#)gg {F.:rcrrN)r4r8 _delta_cdfr)rYrcrNrs&& r.test_delta_cdfTestNorm.test_delta_cdfsi) %%bfbf#%>e4 %%3i2c6%Ee4r0r6N) rrrrr2r6rBrrrrGrrr s@r.r9r9s>5 5 D [[UQF+ 5, 5r0r9c0a]tRtRtoRtRtRtRtVtR#) TestUniformigh-10300c\P!RRRRR\P.4p\\\ P PV4R#r,)rTr r*rr)r4rrrs& r.r2!TestUniform.test_nan_raises_error7 HHfffffbffE Fj%--"3"3Q7r0c\P!RRRRR\P.4p\\\ P PV4R#r,)rTr r+rr)r4rrrs& r.r6!TestUniform.test_inf_raises_errorrNr0r6N) rrrr__doc__r2r6rrr s@r.rJrJs8 88r0rJc Ja]tRtRtoRtRtRtRt]PPR.R O4R4t ]PPR.RO.RO.RO.RO.RO.4R 4t ]PPR.RO.RO.RO.RO.RO.RO.4R 4t R tVtR #) TestExponNormicLRp^^^rCpRWC,, p\PPWRVRR7p\Wa!WCV44R^RrCpRWC,, p\PPWRVRR7p\Wa!WCV44^^^rCpRWC,, p\PPWRVRR7p\Wa!WCV44R^ RrCpRWC,, p\PPWRVRR7p\Wa!WCV44R#) c&R^W,^,, ,p^W,^,, VR,,pR^W,^,,R,,pV^V, ,W,RW,, ,WE.#)rS@rr r6)lamsigropK2exp_skewexp_kurts&&& r.get_moms,TestExponNorm.test_moments..get_momssncg\))DCI>)D4L8Ha39q.0B77H3J#sw- 7L Lr0rSrrarrNr r)r4 exponnormr)rYr\rrXrWKstss& r.rCTestExponNorm.test_momentss M!Q 39 oo##AS&#IC#B!781c 39 oo##AS&#IC#B!78!Q 39 oo##AS&#IC#B!782s 39 oo##AS&#IC#B!78r0c\P!RRRRR\P.4p\\\ P PV^^R7R#r-r.r/r0r1rN)rTr r*rr)r4r^rrs& r.r2#TestExponNorm.test_nan_raises_error; HHfffffbffE Fj%//"5"5qqKr0c\P!RRRRR\P.4p\\\ P PV^^R7R#rc)rTr r+rr)r4r^rrs& r.r6#TestExponNorm.test_inf_raises_errorrer0cV\\PPR^4R4\\PPR^4R4\\PPRR4R4\\PPRR4R4R#)rrNi|)rr4r^rfrXs&r.test_extremes_xTestExponNorm.test_extremes_xsnEOO//a8#>EOO//a8#>EOO//d;SAEOO//d;SAr0zx, K, expectedc\\\PPW4VRR7R#r)rr4r^rf)rYr^r_rNs&&&&r. test_std_pdfTestExponNorm.test_std_pdfs ++A18%Hr0zx, K, scale, expectedc\PPWVR7pVR8Xd VR8XgQhR#\WTRR7R#)r^rrrN)r4r^rnrrYr^r_rdrNrzs&&&&& r.test_cdf_small_KTestExponNorm.test_cdf_small_Ks9 OO  E  2 s?8O8 Ae 4r0c\PPWVR7pVR8Xd VR8XgQhR#\WTRR7R#)r^r-a=rN)r4r^rrrps&&&&& r.test_sf_small_KTestExponNorm.test_sf_small_Ks9 OO  q5  1 s?8O8 Ae 4r0r6N))rQrg6 N-)rrg՜n+H?)rrgi?)rrgbJI-)rbrgI8?)rbrgQ3|-0?)r9rrgVAҤ?)r{Gzt?rg'^ >)r+rrr)r+rrgbr ;)rr2rgVMe?)rbrrgh G}$;)r rwrg ]fJ?)rrwrzg31"g#;)rbrwrzgf +-)rQrwrzr)r r2rga9S?)rrrrrCr2r6rjrrrrmrqrurrr s@r.rSrSs96L L B$ [[-BCICI [[4>D4G? AB 5 B 5& [[4B@DF3@ BC5C5r0rSca]tRtRtoRtRt]PPR.R O4R4t ]PPR.R O4R4t Rt Vt R#) TestGenExponi(c^RIHp\PP \ P !^^ R4RRR4p\V!VRR7^^4R#)r9)simpsonrrzr)dxN)scipy.integrater{r4genexponrfrTr'r)rYr{rzs& r.test_pdf_unity_area TestGenExpon.test_pdf_unity_area)s?+ NN  ryyB5sC EGA$/A6r0c\PP\P!^^ R4RRR4p\P !^V8*V^8*,4'gQhR#)r9rrzrN)r4r~rnrTr'r)rYrns& r.test_cdf_boundsTestGenExpon.test_cdf_bounds/sKnn  1b$!7c3GvvqCxC1H-....r0z x, p, a, b, cc\PPWWE4p\WbRR7\PP W#WE4p\WqRR7R#rp)r4r~rrr)rYr^rzr,r-rrrs&&&&&& r.ruTestGenExpon.test_sf_isf?sD^^  qQ *E*nn  q,U+r0c\PPWWE4p\WbRR7\PP W#WE4p\WqRR7R#rp)r4r~rnrr)rYr^rzr,r-rrnrs&&&&&& r.rTestGenExpon.test_cdf_ppfMsDnn  q,U+nn  q,U+r0r6N))rgM((++A13WE((,,S4aeDr0r6N)rrrrrurrr s@r.rrZsEEr0rc&a]tRtRtoRtRtVtR#) TestExponpowifc\\PPRR4R4\\PP \PP ^R4R4^4R#)rrrrN)rr4exponpowrnrrrXs&r.r)TestExponpow.test_tailgsLENN..ub95AENN..u~~/@/@B/GL r0r6N)rrrrr)rrr s@r.rrfsr0rc2a]tRtRtoRtRtRtRtVtR#) TestSkellamimc\P!R^4p^ ^r2\P!.RO4p\\P P WV4V^R7R#)rbrUNr)gYjP'?g$S?gQ`s2?gC/qF?gD<]Y?g/Xj?g86y?g]a¨?g:?g?gQ>?g?g:?g]a¨?g96?g0X?gA<]?gB/q?gQ`s?g$S?gYjP'?gX_?g+rx?gfSr?gJXx~?)rTr'r rr4skellamrN)rYr]mu1mu2skpmfRs& r.rlTestSkellam.test_pmfnsK IIc2 qS - . EMM--ac:FBOr0c\P!R^4p^ ^r2\P!.RO4p\\P P WV4V^R7R#)rbrUNr)gƃQ?gS `X'?gVy)>?g: S?gei.f?gȼJ]x?g)?g{^'+i?gLIg>?gnZ ?g?rvͪ?gA?g0?g#J?gqSm . ?g~B?gd?gFBD?gQ9?gmF^?gcMy?gGJq?gC^c?g.){?gñݸ?)rTr'r rr4rrn)rYr]rrskcdfRs& r.rZTestSkellam.test_cdfsK IIc2 qS - . EMM--ac:FANr0c^^Rr2p\\PPWV4^RR7\\PP WV4^RR7R#)r9gs1CrrN)rr4rrNrn)rYr^rrs& r.test_extreme_mu2TestSkellam.test_extreme_mu2sH. ))!#6F ))!#6Fr0r6N) rrrrrlrZrrrr s@r.rrmsP*O*GGr0rc a]tRtRtoRtRt]P!RR7R4t]PPRR^.4]PPR .RO4]PPR .RO4]PPR ] !R^R 7Uu.uFpR V9gK VNK up4] P!RR7R44444tRtRtVtR#uupi) TestLognormic\P!4;_uu_4\P!R\4\P P .RO^4p\V.RO4RRR4R# +'giR#;i)rwNr9rzr)rg~r?ge3E?)r{r|r}r~r4lognormrfrrYrfs& r.rPTestLognorm.test_pdfsS $ $ & &  ! !'> :--##K3C %c+H I' & & &s A A44 B c R^Rr2p\\PPW, VR7\PP\ P !W, 4V, 44\\PPW, VR7\PP\ P !W, 4V, 44R#)g(\5i@gOn?r>N)rr4rrr8rTrr)rYrrsigmas& r.rTestLognorm.test_logcdfsU ((%(8 bffRUmE&9: < ++BEU+; ((u)<= ?r0rMrNc@\PPR4#rQrRrXs&r.rWTestLognorm.rngrTr0rrrrrwFrxrzr{c\PP^dWW'R7p/p V'dWR&V'dW)R&V'dW9R&\\PV3/V BRR/BR#)r)rr?rdrcrrrrrTN)r4rrrr~s &&&&&&&& r.r'TestLognorm.test_fit_MLE_comp_optimizersk}}  cY%,!@ "J "L &N &emmT a"&&!"56r0c h.ROp\\\PPV^^RR7R#)rr@)rrrANrrzr`)rr,r4r rrs& r.test_fit_bad_keyword_args"TestBeta.test_fit_bad_keyword_argss$ i1$ &r0cf.ROp\\\PPVRRR7R#)rrz)fafix_aNr)rr)r4r rrs& r.#test_fit_duplicated_fixed_parameter,TestBeta.test_fit_duplicated_fixed_parameters! j%**..!3Gr0zOverflow, see gh-14901reasoncdRRRr2p\\PPWV4R4R#)g?rog0xAgx)>N)rr4r r)rYrzr,r-s& r.test_issue_12635TestBeta.test_issue_12635s)%dJa qQ/1FGr0cN\P!.RO4p\P!.RO4pRp\PP W2^,RV, 4p\ WA4\PP WB^,RV, 4p\ WS4R#)@3@?rrN)rg^?gNC?)rbrr)rTr r4r rrr)rYinv_R count_listrzinvrs& r.test_issue_12794TestBeta.test_issue_12794sv12XXo. jjnnQQ0CD#jjmmCa*1DEr0cRp\P!^^4pRp^V, V^,W2, rep\PP WEV4p\PP WuV4p\ V^V, 4R#)h㈵>rN)rTr'r4r rrnr) rYalpha_2count_nobsrr,r-rrs & r.test_issue_12796TestBeta.test_issue_12796si1b!g+vz4=ajjnnQ1%jjnnSQ'Q[)r0cN^Rr!\\PP^W4\P 4R^r!\\PP^W4\P 4^^r!\\PP^W4^4\\PPRW4^4^^r!\\PP^W4^4\\PPRW4^4R#)rrzrw+phNrS)rr4r rfrTr+)rYr,r-s& r.rTestBeta.test_endpointss#1UZZ^^Aq,bff5A1UZZ^^Aq,bff5!1UZZ^^Aq,a0UZZ^^FA115!1UZZ^^Aq,a0UZZ^^Ha3Q7r0cPRRRr2p\PPWV4R#)gףp= ?gvH7B@0BN)r4r r)rYrr,r-s& r.test_boost_eval_issue_14606$TestBeta.test_boost_eval_issue_146064sa qQr0ra, bcBRpV!WBV4R# \dR#i;i)rN) OverflowError)rYrr,r-rzs&&&& r. test_beta_ppf_with_subnormal_a_b)TestBeta.test_beta_ppf_with_subnormal_a_b8s+   1O   s  z x, a, b, refc`\PPWV4p\WTRR7R#r)r4r rrrYr^r,r-rrs&&&&& r.rTestBeta.test_logcdfQs$""1+%0r0c^\PPWV4p\WTR4R#)rN)r4r rrrYr^r,r-rrs&&&&& r.rTestBeta.test_logsfYs$   q)E*r0rcb\\P!W4P4V4R#rtrr4r rErs&&&&r.rtTestBeta.test_entropyjs  1(002C8r0za, b, ref, tolcd\\P!W4P4W4R7R#rr)rYr,r-rrs&&&&&r.rTestBeta.test_extreme_entropyrs!H  1(002CBr0c F\P!\P4p\P!.RO4p\P!.RO4pV!W#R7P 4p\P!V!V^,V^,R7P 4V!V^,V^,R7P 4V!V^,V^,R7P 4V!V^,V^,R7P 4.4p\ WE4R#)SArN)rrrgrb)rrgrrQ)r4make_distributionr rTrrEr)rYBetar,r-rrs& r.test_entropy_broadcasting"TestBeta.test_entropy_broadcastings&&uzz2 JJ* + JJ* +Qn$$&jj$11.66811.66811.66811.668:; !r0r6N)r)@)rr)rrrg>׵NlIq)& .!>rlrg'rO)rrrgp}^U)rrlrg 1*))rzrzgο)r2rg)rrg&׾k )rrg*)M))rrbgCpr )rbrQg8EOr)gNAgNAg(4Fr)rgSAg8TAUr)rg _Bg[VVuD&r)rgrggR#Lr)r rg?r5r)r r\r )rbrgnQc}4r)rbrg{CEr )rbrg%\r )rr rr )rrrr )rgrbrr )rrrrrrrrrrskipif MACOS_INTELrrrrrrr4r rrrrrrtrrrrr s@r.rrs,7& H  [[ ,DE HF H [[ ,DE F " [[ ,DE *F *8.  [[X  'GH [[Vnn%EFGI. [[^GCEF1F1  [[^GEGH+H+ [[[DE 9 E 9 [[ . C/. C " "r0rca]tRtRto.ROtRtRt]PPR.RO4R4t ]PPR]4R4t Rt ]PPR ].RO,4R 4t ]PPR RR.4R 4tR tR tRtRtRtRtVtR#) TestBetaPrimeic\RRr!\P!.RO4p\P!W4p\ \P !VP V44P44\VPV4\P!VP V444R#)rrNr) rTr r4 betaprimerrrrrrfrrs& r.rTestBetaPrime.test_logpdfsj4t HH_ % OOE ( AHHQK(,,./a"&&!"56r0c \PP^RR4p\VR4RRr2\P !.R O4p\PPWV4p\ \P!V4P44\PPpVUu.uFqe!\PWbV4NK pp\WG^RR7R#uupi) r9rwrrrr-=rNr) r4r rnrrTr rrrrC _cdf_singler)rYr^rr r(gen_cdfrFcdfs_gs& r.rZTestBetaPrime.test_cdfs OO  3 ,Q4t HH_ %""1T2 D!%%'(%%11HIJ'%//3t<J159Ks4 C&zp, a, b, expectedc`\PPWV4p\WTRR7R#r)r4r rr)rYrzr,r-rNr^s&&&&& r.rTestBetaPrime.test_ppfs$ OO  a (%0r0z x, a, b, pc^\\PPWBV4VRR7R#rp)rr4r r)rYr^r,r-rzs&&&&&r.test_ppf_gh_17631TestBetaPrime.test_ppf_gh_17631s++A!4aeDr0c\P!R4p\P!R4p\P!R4p\\PP W1V4RRR7R#)rSrzrvrN)rTr rr4r _ppf)rYr,r-rzs& r. test__ppfTestBetaPrime.test__ppfsF HHSM HHSM HHSM,,Q15sGr0zx, a, b, expectedc^\\PPWV4VRR7R#r)rr4r rn)rYr^r,r-rNs&&&&&r.test_cdf_gh_17631TestBetaPrime.test_cdf_gh_17631s! ++A!4hUKr0cr\PPWV4pVR8gQh\WTRR7R#)rSgh㈵>rN)r4r rnr)rYr^r,r-rNrs&&&&& r.test_cdf_extreme_tails$TestBetaPrime.test_cdf_extreme_tails s/ OO  a (3ww$/r0c.ROp.ROp.ROp.ROp\PPW1V4p\WTRR7R#)rrrN)rrr rjrjrjrjY@r rjrjrjrr)r!r rrrrrrjrjr r r rr)rrꌠ9Y>)Frr^ 9^;r|rrrr r|rr)g/:g U/7gKH9g<![a?gw 4?g)J?g5?z ^J?$5Є*?m嗽?D>?g%?gsM )r|rrgrNs: )r|rjr gK>))rrPrgV)?)rrPrgWJB9*=)rrPrgOGJ<)rBro@gΣ01)?ror. g8ǘ`?)>q?rjrgMD))rrrr )rrjrgV[?)rrjrr0 )rgrjrg%?)rgr rjgoXr?)rrrrcdf_valsrrZrrrrr r r r rrrr) r, rrr s@r.r r sH 7:> [[ B11 [[\84E5EH [[   L L [[ . / 120 204*100 " "r0r c fa]tRtRtoRtRtRtRtRtRt ] PPRR R .4R 4t ] PPR .RO4R 4t] PPR.RO4] PPR.RO4] PPR.RO4] PPRRR.4] PPRRR.4] PPRRR.4R444444tRtVtR#) TestGammaiOc\PP^ZRRR7p\VR4\PP^^ RR7p\VR4R#)Zir^g CTb?gBΖ?Nrw)r4rrfrrs& r.rPTestGamma.test_pdfPsFkkoob#To2C-kkooa4o0C+r0c^\PP^^4p\V^4R#r)r4rrrrs& r.rTestGamma.test_logpdfXs$##Aq)FA&r0cf.ROp\\\PPV^RR7R#)rr@)rrANr)rr,r4rrrs& r.r#TestGamma.test_fit_bad_keyword_args^s i!!8Lr0c\P!\PP R^4RRR7'gQh\P!\PP R^d4RRR7'gQhR#)rBgplC@r rrg1}t@rN)rTiscloser4rrrXs&r.rTestGamma.test_isfbsf"zz%++//%3+%99 99zz%++//%5+%99 99r0cj^Pp^pRp\PPW4p\WCRR7R#)PrrNgޟ򑁺)r4rrrrs& r.rTestGamma.test_logcdfxs/  $##A)%0r0cjRpRpRp\PPW4p\WCRR7R#)r2r rrNgDm)r4rrrrs& r.rTestGamma.test_logsfs/  % !!!'/r0rdrSrTc\PPV^,V^,^VR7p\VRRR7R#)r^g>S*rrN)r4rrFrrYrdrs&& r.rGTestGamma.test_delta_cdfs6 &&uSy%)Qe&L6UCr0z a, ref, rtolcZ\\PPV4W#R7R#r)rr4rrE)rYr,rrs&&&&r.rtTestGamma.test_entropys  ++A.?r0r,rcrTFrrc\PPR4p\PP WV^dVR7p/p V'dWR&V'dW)R&V'dW9R&^\ V 4, p V ^8XdURp \P!\V R7;_uu_4\PP!V3RR/V BR R R 4R #\PP!V3RR/V Bp \P!V !p V ^8d/\V P4\P!V44V ^8d7\V P^4\P!V^,44V ^8d9\V P^4\P!V^,44R #R # +'giR #;i) rrrrrrrrrN)rTrUrVr4rrr#rr r)rrrr)rYr,rcrdrrrrWrrrrthetarJs&&&&&&& r. test_fit_mmTestGamma.test_fit_mmsPii##$78{{qS,/1 J L "NCI  A:MIz;; :T:T:<  :T:T:{{E" A: DIIK 7 A: DKKNBGGD!G,< = A: DKKNBGGD!G,< = <; s ($G G r6N))rgt.r)r r~r)rgX>k-־ @r)rgwM7@r)gNgmCg%og$6@r)roܘ#]@r)rrr )rr9r )rrrrrPrrrrrrrrrGrtrL rrr s@r.r3 r3 Os,,' M9,10 [[WsCj1 D2 D [[^BC@C@ [[S.1 [[UN3 [[Wn5 [[WtUm4 [[Yu 6 [[[4-8>975642 >r0r3 cRp\P!RV^,, ^V^,, ^^R7p\VPR4R4R#)r2rbrSgV{q@Ng;f?)r4rrrf)jitterZs r.test_pdf_overflow_gh19616rR sCF 6FAI-q619}!1MAAEE#J 12r0ca]tRtRtoRt]P PR.R O4R4t]P PR.R O4R4t ]P PR.R O4R4t Rt R t Vt R #) TestDgammaic\PPR4p^ pVP^ VR7pVP ^ VR7p\ P PW44p\ PP\P!V4V4^, p\WV4\ P !V4p\VPV4VRR7R#)lOP($rdr)rrrvrN) rTrUrVnormalrr4dgammarfrrr)rYrWrr^r,rrrJs& r.rPTestDgamma.test_pdfsii##$78 JJRdJ + KKRdK +llq$kkoobffQi+a/!||A Su5r0zx, a, expectedcl\PPW4p\WCRR7\PP W24p\WQRR7\PP V)V4p\WcRR7\PP W24p\Wq)RR7R#r)r4rX rnrrrr)rYr^r,rNrnrrrs&&&& r.test_cdf_ppf_sf_isf_tail#TestDgamma.test_cdf_ppf_sf_isf_tailsllq$E2llx+U+ \\__aR #51llx+Re,r0rc\\\PPV4VRR7R#rrr4rX rErs&&&r.rtTestDgamma.test_entropys  ,,Q/5Ar0c\\\PPV4VRR7R#rr^ rs&&&r.test_entropy_entreme_values&TestDgamma.test_entropy_entreme_valuess"  ,,Q/5Ar0c\P!.RO4p\PP V4p\ \ V44F5pW#,\PP W,48XdK5Qh R#)rN)rrrr)rTr r4rX rErr#)rYr^rrs& r.test_entropy_array_input#TestDgamma.test_entropy_array_input!sX HH' ( LL  #s1vA45<<//55 55r0r6N))rrggrU>)r5rg$IC<)irg2Wj<)rg~)A.=)rrgf?))rgmqn@)?g?)rRgx?))r|rb)rg>_)rg5i)rgo@)rg4Fa3@)rgrO]@)rrrrrPrrrr[ rtra rd rrr s@r.rT rT s 6, [[-:; - ; - [[X89 B 9 B [[X:; B; B66r0rT c a]tRtRtoRt]P PR.R O4R4t]P PR.R O4R4t Rt ]P PR.RO4R4t R t R t VtR #)TestChi2i(c\\PPRR4R^R7\\PP^d^d4R^R7R#)rgޞw1D?rUg6:֜?Nrr4chi2rfrXs&r.rTestChi2.test_precision+s@EJJNN468N$& (EJJNN346J$& (r0rlc^\PPW4p\WCRR7R#r)r4rl rr)rYr^rrrs&&&& r.rTestChi2.test_logcdf2s"""1)%0r0c^\PPW4p\WCRR7R#r)r4rl rr)rYr^rrrs&&&& r.rTestChi2.test_logsf=s"   '/r0c~Rp\PPRV4p\VRRR7\PPRV4p\VRRR7^ p\PPRV4p\VRRR7\PPR V4p\VR RR7R #) r^gλm:=6g'bd5)rgƯʉ?)rgn>@)g͔v0@)rg9:3@)rrrrrrrrrrrrtrw rrr s@r.ri ri (s( [[ + 1  1  [[ + 0  0 F0 [[Y:; C ; C==r0ri ca]tRtRtoRtRtRtRt]PPRR ^.4R4t Rt Vt R #) TestGumbelLimcN\PPR4VnR#)l>SNrSrXs&r.rZTestGumbelL.setup_methodnr\r0c\P!RR4p\PP V4p\PP V4p\ W4R#rNrr)rTrWr4r^rnrrrs& r.rTestGumbelL.test_cdf_ppfrsB KKb ! NN  q ! ^^   "r0c \P!RR4p\PP V4p\PP V4p\P !V4p\P!V4)p\WE4R#r ) rTrWr4r^rrrrexpm1r)rYr^rrXurGs& r.rTestGumbelL.test_logcdf_logsfxs` KKb ! NN ! !! $ NN  # FF1I ]]1  r0c\P!R^4p\PP V4p\PP V4p\ W4R#r)rTrWr4r^rrrrs& r.ruTestGumbelL.test_sf_isfsB KKQ  NN  a  ^^   "r0rcc\PP^dWPR7p\PP W!R7wr4\ W14R#)rrrN)r4r^rrWrr)rYrcr fitted_locrs&& r.test_fit_fixed_param TestGumbelL.test_fit_fixed_paramsB~~!!s((!K**4*: Z%r0rNr)rrrrrZrrrurrrr rrr s@r.r{ r{ msB5   [[URG,&-&r0r{ c,a]tRtRtoRtRtRtVtR#) TestGumbelRic\\\PP^24RRR7R#)r=?~T}%m;r rN)rr4r]rrXs&r.rTestGumbelR.test_sfs# ))"-/E" $r0c\\\PPR4RRR7R#)rBr< r rN)rr4r]rrXs&r.rTestGumbelR.test_isfs# **513D" $r0r6N)rrrrrrrrr s@r.r r s$$$r0r c a]tRtRtoRt]P !RR7R4t]P R4t]P R4t ]P R4t ]PP]PPR ]P!^24]P!R ]PPR 7.4]PPR R R.4]PPR.R:O4]PPRR;R<.4R44444t]PP ]PPRR^.4R44tRt]PP'RR7]PPRR=R>.4]PPRR R.4R444tRtRtRt]PP]PPR]P!.R?O.R@O.RBO4]P!.RCO.RDO.REO]PPR 7]P!.RFO]P2!R ^^4]P2!RA^^4]PP R 7.4R!44t]PPR]P!.R?O.R@O.RBO4]P!.RCO.RDO.REO]PPR 7]P!.RFO]P2!R ^^4]P2!RA^^4]PP R 7.4R"4t]PPR#^^.4]PPR$R%R&.4R'44t]PPR(R)^.R*^..4R+4t]PPR,RG^]P<]P>]P>3.RHRI..4R-4t ]PPR..RJO4]PPR/]!PDPFR]P2!RK^^ 4R03]!PDPHR]P2!RK^^ 4R03]!PDPFRL]P2!^^^ 4R3]!PDPHRL]P2!^^^ 4R03.4R144t%]PPR2.RMO4R34t&]PP&]PPR2.RNO4R444t']P R54t(R6t)R7t*R8t+Vt,R9#)OTestLevyStableicN\PPR4VnR#)lYNrSrXs&r.rZTestLevyStable.setup_methodr\r0T)autousec R\PnR\PnR\Pn\P P \PnR#)z2Setup default parameters for levy_stable generatorS1 piecewiseN)r4rparameterizationcdf_default_methodpdf_default_method _levy_stable _QUAD_EPSquad_epsrXs&r.reset_levy_stable_params'TestLevyStable.reset_levy_stable_paramssH.2*/:,/:,%*%7%7%A%A"r0c \P!\\4PR, 4p\P P VPRR7pV#)a{Sample data points for pdf computed with Nolan's stablec See - http://fs2.american.edu/jpnolan/www/stable/stable.html There's a known limitation of Nolan's executable for alpha < 0.2. The data table loaded below is generated from Nolan's stablec with the following parameter space: alpha = 0.1, 0.2, ..., 2.0 beta = -1.0, -0.9, ..., 1.0 p = 0.01, 0.05, 0.1, 0.25, 0.35, 0.5, and the equivalent for the right tail Typically inputs for stablec: stablec.exe << 1 # pdf 1 # Nolan S equivalent to S0 in scipy .25,2,.25 # alpha -1,-1,0 # beta -10,10,1 # x 1,0 # gamma, delta 2 # output file z.data/levy_stable/stable-Z1-pdf-sample-data.npyx,p,alpha,beta,pctrrTrrrrrrTr$s& r.nolan_pdf_sample_data$TestLevyStable.nolan_pdf_sample_datasL6ww N ! ! < = vv  /C D r0c \P!\\4PR, 4p\P P VPRR7pV#)aSSample data points for cdf computed with Nolan's stablec See - http://fs2.american.edu/jpnolan/www/stable/stable.html There's a known limitation of Nolan's executable for alpha < 0.2. The data table loaded below is generated from Nolan's stablec with the following parameter space: alpha = 0.1, 0.2, ..., 2.0 beta = -1.0, -0.9, ..., 1.0 p = 0.01, 0.05, 0.1, 0.25, 0.35, 0.5, and the equivalent for the right tail Ideally, Nolan's output for CDF values should match the percentile from where they have been sampled from. Even more so as we extract percentile x positions from stablec too. However, we note at places Nolan's stablec will produce absolute errors in order of 1e-5. We compare against his calculations here. In future, once we less reliant on Nolan's paper we might switch to comparing directly at percentiles (those x values being produced from some alternative means). Typically inputs for stablec: stablec.exe << 2 # cdf 1 # Nolan S equivalent to S0 in scipy .25,2,.25 # alpha -1,-1,0 # beta -10,10,1 # x 1,0 # gamma, delta 2 # output file z.data/levy_stable/stable-Z1-cdf-sample-data.npyr rr r$s& r.nolan_cdf_sample_data$TestLevyStable.nolan_cdf_sample_datasMJww N ! ! < = vv  /C D r0c n\P!\\4PR, 4pV#)zSample data where loc, scale are different from 0, 1 Data extracted in similar way to pdf/cdf above using Nolan's stablec but set to an arbitrary location scale of (2, 3) for various important parameters alpha, beta and for parameterisations S0 and S1. z1data/levy_stable/stable-loc-scale-sample-data.npy)rTrrrrr$s& r.nolan_loc_scale_sample_data*TestLevyStable.nolan_loc_scale_sample_datas/ww N ! ! ? @  r0 sample_sizerimarksr S0r z alpha,betarSrzz gamma,deltacV\Pn\P!W#WER7p\P!VP W`P R7VP 4wrV R8gQhR#))rr rdrcrrjN)r4rr rmrrWrn) rYr rr rrr lsrrzs &&&&&&& r.rGTestLevyStable.test_rvs s^(.>*   % || FF ((F ;RVV 4xxr0r c  RpRpRp\PPW!W4RVPR7p\P!VRW!W43R7wrgVR8gQhR #) z3Additional test cases for rvs for alpha equal to 1.rSrzrrwrrrrN)r4rrrWrm)rYr rrcrdr^statrzs&& r.test_rvs_alpha1TestLevyStable.test_rvs_alpha1*sc    ! !%3'+$(( " D,,q-%*#$=?4xxr0c|.R Op\PPV4wr#rE\VR^RR7\ VR ^4\ VR^4\ VR^4V.R O,p\PPV4wrxr\ V^4\ VR 4\ V R^4\ V R^4R#) gwg?67',?Gz?rrhgۧ1?g9]c?gv?N)wgr rrrr羼u?r r r r r r r r r tA}˜.?r r r r )\(̿)r r r r r r)r4rrrrr rYr^alpha1beta1loc1scale1ralpha2beta2rrs & r.r TestLevyStable.test_fit7s  ',&7&7&A&A!&D#t148E4+FGQ/ '1 9 9&+&7&7&A&A"&E#tVQUBFFA.D&!,r0zUnknown problem with fitstart.rzalpha,beta,delta,gammaparametrizationc  V\Pn\PPWW4RVPR7p\PP V4pVwrr\ WW4.WW.RR7R#)z7Test that fit agrees with rvs for each parametrization.rrrrN)r4rr rrWrr) rYrr rrr rr alpha_obsbeta_obs delta_obs gamma_obss &&&&&& r. test_fit_rvsTestLevyStable.test_fit_rvsOsw-<)  $$ Ue$((% ))$/471 Y % ' ) 7 r0cF\P!.RO4p\PP V4wr#rE\PP V)4wrgr\ V^4V^8wgQh\ Wb4\ Ws)4\ W)4\ W4R#rN) rrr!r!rbrbrbrrrr)rTr r4rrrr) rYr^r r r r r r rrs & r.test_fit_beta_flip!TestLevyStable.test_fit_beta_fliphs HH? @&+&7&7&A&A!&D#t&+&7&7&A&A1"&E#tUAqyyF+E6*D%(F+r0c:^p\P!.RO4p\PP V)4wr4rV\PP V)V,4wrxr\ Ws4\ W4\ WV,4\ W4R#r )rTr r4rrr) rYSHIFTr^r r r r r r rrs & r.test_fit_delta_shift#TestLevyStable.test_fit_delta_shifttsx HH? @&+&7&7&A&A1"&E#t&+&7&7&A&A1"u*&M#tF+E)D,/F+r0cz.ROp\PPV4wr#rEV^8g QRV 24hV\V48gQR\V4 RV 24h.ROp\PPV4wrxrV^8g QRV 24hV \ V48gQR\ V4 RV 24hR#) rzExpected alpha < 1, got zExpected loc < z, got zExpected alpha > 1, got zExpected loc > N) rrr!r!rbrbrbrrr ) rrr!r!rbrbrbrrr )r4rrrr&r s & r.test_fit_loc_extrap"TestLevyStable.test_fit_loc_extraps 6&+&7&7&A&A!&D#tz>5fX>>zc!f}DAxvdVDD} 7&+&7&7&A&A"&E#tz>5fX>>zc"g~FR vFF~r0z pct_range,alpha_range,beta_rangerc JaaaaTp\P!4pVPR8H;'dVPR8HoRP VPVPVP .4pRRVVV3Rl.RRVVV3R l.RRVVVV3R l.RR VVV3R l.R RVVV3Rl.R RVVV3Rl..p\ V4EFwp wrp V \Pn V eW\!V4,MTp Rp\P!4;_uu_4\P!RV\4\PPV R,V R,V R,^^R7p\P !RR7;_uu_4\#V .R#OV\P$!WR,, 4\P$!WR,, 4\P$!V R,4, .4pRRR4XVR,V 8\P&!V4,,pRV RV RV RVP(P* R V 2 p\-VV R,V VR!R"7RRR4EK R# +'giL;i +'giEK;i)$z2Test pdf values against Nolan's stablec.exe outputLinuxi686/dnirGc > <\P!VR,S4\P!VR,S4,\P!VR,S4,VR,^8HVR,R8H,VR,R8VR,R8,VR,R8H,,VR,R8*\P!VR,RR.4,,VR,R 8*\P!VR,R R .4,,VR,R 8*\P!VR,R R.4,,VR,R 8H\P!VR,RR.4,\P!\P!VR,4.RO4,,VR,R 8H\P!VR,R.4,\P!\P!VR,4R .4,,VR,R 8H\P!VR,RR.4,\P!\P!VR,4.RO4,,VR,R 8HVR,R8H,VR,R8H,,VR,R 8HVR,R8H,VR,R8H,,VR,R 8HVR,R 8H,VR,R8H,,VR,R8H\P!VR,R.4,\P!\P!VR,4.RO4,,VR,R8H\P!VR,RR.4,\P!\P!VR,4.RO4,,VR,R8H\P!VR,RR.4,\P!\P!VR,4RR .4,,VR,R8,(,#)pctrr rzrr( rrrrrjrrwrrorffffff??rSrR)rzr`ra333333ӿ)gٿr rrrzrrwrr)rrrSr)rTisinrr alpha_range beta_range pct_ranges&r.r@7TestLevyStable.test_pdf_nolan_samples..sGGAeHi0GGAgJ 45GGAfIz23 vY!^uX_.vY#-wZ3.0uX_. wZ3.GGAeHsCj9: wZ3.GGAeHsCj9:$wZ3.GGAeHr2h78%,wZ3.GGAeHsCj9:GGBFF1V9$5|DE-6wZ3.GGAeHrd34GGBFF1V9$5t<=7@wZ3.GGAeHsCj9:GGBFF1V9$57MNOA#JwZ3.vY#-/uX-/K(TwZ3.vY$.0uX-/U-^wZ3.vY#-/uX-/_2hwZ2-GGAeHrd34GGBFF1V9$57GHIi7rwZ2-GGAeHsCj9:GGBFF1V9$5|DEs<|wZ2-GGAeHsCj9:GGBFF1V9$5RyAB}AF7s*GC,EHr0r rc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8,VR,R8g,#)r rr rwrSrTr r s&r.r@r siGGAeHi0GGAgJ 45GGAfIz23wZ#%'wZ2% 'r0c<VR,R8HS'*,\P!VR,S4,RS9,\P!VR,S4,#)rrSr r r )rr r is_linux_32r s&r.r@r sXwZ2%$_&GGAeHi01;&(GGAfIz2 3r0g& .=c<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8*,#)r rr rwr r s&r.r@r sVGGAeHi0GGAgJ 45GGAfIz23wZ3&(r0 fft-simpsonrc<VR,R8\P!VR,S4,\P!VR,S4,\P!VR,S4,#)rrr r r r s&r.r@r sVwZ3&GGAeHi01GGAgJ 45GGAfIz23r0rc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,^8,VR,R8,#)r rr rr r s&r.r@r (shGGAeHi0GGAgJ 45GGAfIz23wZ!^%wZ#% 'r0Nz0Density calculations experimental for FFT methodrzr^rr rdrcrrelerrrzz pdf test  failed with method 'z ' [platform: z]  Ferr_msgverbosecalcabserrr )platformunamesystemmachinejoin processorrLr4rr r{r|rr~rfrTrrrisnanrrr)rYr r r r rr platform_desctestsixdefault_methodr filter_funcsubdatar*rzsubdata2failuresmessager s&&fff @r.test_pdf_nolan_samples%TestLevyStable.test_pdf_nolan_sampless4% llg-II%--62I  \\5==%// :< tJL ^U% U% W' t& t& kC L"+5!1 B3AE   0 + 7;t,=A DC((**''#~F%%))CLG$FO *[[h//04FF1s|#34FF1s|#34rvvgcl7KK H0$H%-((1+  t#88HI""/HNN4H4H3IH:W CL#! 9+* "20/+**s'A5JA7I= .~sGGAeHi0GGAgJ 45GGAfIz23 wZ2-GGAfI/ABCuX-/ wZ2-GGAfI?@uX-/  r0rjc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8H\P!VR,.RO4,VR,R8H,,VR,R8H\P!VR,.RO4,VR,R8H,,#r r r s&r.r@r sGGAeHi0GGAgJ 45GGAfIz237r)& +=>?5T)+ 7r)& ?;<5T)+ r0r rc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8,#)r rr 333333?r r s&r.r@r sVGGAeHi0GGAgJ 45GGAfIz23wZ#%'r0rc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8,VR,R8*,#)r rr rr r r s&r.r@r iGGAeHi0GGAgJ 45GGAfIz23wZ#%'wZ3& (r0r2c<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8,VR,R8*,#)r rr rg rr r s&r.r@r r r0rc<\P!VR,S4\P!VR,S4,\P!VR,S4,VR,R8,VR,R8*,#)r rr rSrg r r s&r.r@r r r0NrzzZCumulative density calculations experimental for FFT method. Use piecewise method instead.r^rr r rr rzz cdf test r z' r Fr r )rLr4rr r{r|rr~rnrTrrrr rrr)rYr r r r rr r r rr r rzr r r s&&fff r.test_cdf_nolan_samples%TestLevyStable.test_cdf_nolan_samplesas 2%U% ,T$ &t& t& t& t& GL \"+5!1 B3AE   0 + 7;t,=A ((**''>" $ %%))CLG$FO *[[h//04FF1s|#34FF1s|#34rvvgcl7KK H0$H%-((1+'rd*??Os&nn2232hZACL#! ;+* "2"0/+**s'A5HA7H A/H HH H. r-r@rfrnc \P!4pVPR8H;'d R\P!4^,9pV'dVR8Xd\P !R4TpR\ PnR\ Pn WfR,V8H,pRV 2\ Pn VR9gQhVR8Xd\ PPM\ PPpV!VR,VR ,VR ,^^R 7p \WV,R 4R #)zGTests for pdf and cdf where loc, scale are different from 0, 1 r 32bitrfz4Test unstable on some platforms; see gh-17839, 17859r r-Sr^rr r rN)rfrn)r r r architecturerrr4rr r r rfrnr) rYr r-r@r r rr rMr$s &&&& r.test_location_scale"TestLevyStable.test_location_scales llg-WW'X=R=R=TUV=W2W  45= KKN O*0;,/:,G}-./0[*~%%%%)U]E   ! !8I8I8M8M  CL''*GFO1!  DM40r0zmethod,decimal_placesr r c  D\P!.RO4p\P!.RO4p\P!.R O4p\P!RR7;_uu_4\P!4;_uu_4\P !R\ RR7V\Pn \PPV^V^^R7p\WdW!4RRR4RRR4R# +'giL;i +'giR#;i) zwsample points extracted from Tables and Graphs of Stable Probability Density Functions - Donald R Holt - 1973 - p 187. rzrzDensity calculation unstable.*)categoryr r N)r9r9r9r9rrrrr r r r r!r!r!r!rrrr)gtF_?gܵ|?gQ?g@?gec]?gHP?gPkw?g!rh?gTN?g ?gDioɴ?g r?gH?g+?gHPs?gX?g&?g 0*?gv/?ga2U0*?)r9rorzrr9rorzrr9rorzrr9rorzrr9rorzr) rTr rr{r|rr~r4rr rfr)rYrdecimal_placesxsdensitybetasrfs&&& r.'test_pdf_alpha_equals_one_beta_non_zero6TestLevyStable.test_pdf_alpha_equals_one_beta_non_zerosXX H ((     [[X & &(?(?(A(A  # #H'8 4:E   0##''AuA1'EC n )B & &(A(A & & &s%(DA#C; )D; D D D zparams,expectedc\PPV^,V^,V^,V^,RR7p\W24R#)r9rraN)r4rr)rYrrNobserveds&&& r.rTestLevyStable.test_stats;sE$$** 1IvayfQivay+  H/r0rzfunction,beta,points,expectedrc ^Tu;8d ^8gQhQh\V!WAVR7\P!\V4V44R#)a4Ensure the pdf/cdf routines do not return nan outside support. This distribution's support becomes truncated in a few special cases: support is [mu, infty) if alpha < 1 and beta = 1 support is (-infty, mu] if alpha < 1 and beta = -1 Otherwise, the support is all reals. Here, mu is zero by default. rr N)rrTrr#)rYrrMr rWrNs&&&&&&r.!test_distribution_outside_support0TestLevyStable.test_distribution_outside_supportIs>P5}1}}} Vt 4 GGCK * r0zx,alpha,beta,expectedc R\Pn\\PP WVR7V4R#)aTest pdf for x equal to zeta. With S1 parametrization: x0 = x + zeta if alpha != 1 So, for x = 0, x0 will be close to zeta. When case "x equal zeta" is not handled properly and quad_eps is not low enough: - pdf may be less than 0 - logpdf is nan The points from the parametrize block are found randomly so that PDF is less than 0. Reference values taken from MATLAB https://www.mathworks.com/help/stats/stable-distribution.html Afc=r3 Nr4rr rrfrYr^rr rNs&&&&&r.test_x_equal_zeta TestLevyStable.test_x_equal_zetaws7Z&-"    ! !!t ! <  r0c R\Pn\\PP WVR7V4R#)aTest pdf for x near zeta. With S1 parametrization: x0 = x + zeta if alpha != 1 So, for x = 0, x0 will be close to zeta. When case "x near zeta" is not handled properly and quad_eps is not low enough: - pdf may be less than 0 - logpdf is nan The points from the parametrize block are found randomly so that PDF is less than 0. Reference values taken from MATLAB https://www.mathworks.com/help/stats/stable-distribution.html r7 r3 Nr8 r9 s&&&&&r.test_x_near_zetaTestLevyStable.test_x_near_zetas7D&-"    ! !!t ! <  r0c,\P!4#rt) threadingLockrXs&r.levy_stable_lockTestLevyStable.levy_stable_locks~~r0c r\PPp\RRRRR7p\P p\P !R /VBpT;_uu_4RVnRVnVP!R /VBR^ RV!R 4/BpVP^ V!R 4R 7p\P!Wv8H4'dQhRVnRVnVP!R /VBR^ RV!R 4/BpVP^ V!R 4R 7p \W4\W4R R R 4R # +'giR #;i) rrrrSrr rcrdr r rr rNr6) rTrUrVr%r4rr rrr rYrB rWr"unfrozenfrozen unfrozen_afrozen_a unfrozen_bfrozen_bs && r.$test_frozen_parameterization_gh208213TestLevyStable.test_frozen_parameterization_gh20821sii##Ccs#>$$"",V,  (,H %&*F #!UURUc)nUJzzrIzGHvvh455 55(,H %&*F #!UURUc)nUJzzrIzGH  .  .   s CD%% D6 c B\PPp\RRRRR7p\P pV;_uu_4RVn\P !R /VBpVP!R /VBR^ RV!R4/BpVP^ V!R4R 7p\Wv4R Vn\P !R /VBpVP!R /VBR^ RV!R4/BpVP^ V!R4R 7p \W4R R R 4R # +'giR #;i) rrrrSrE r rrrF rr Nr6) rTrUrVr%r4rr rrrG s && r.%test_frozen_parameterization_gh20821b4TestLevyStable.test_frozen_parameterization_gh20821bsii##Ccs#>$$  (,H %&&00F!UURUc)nUJzzrIzGH  .(,H %&&00F!UURUc)nUJzzrIzGH  .   s B?D  D rN))rSr9)rSr)rr9)rrzrr9r)rrr r!)rSrr r!)rrzr)rrr r)rr9r)rrjrzrr)rrzrrr )grr9rr`r) rrjrror rzr rrrr)r r r9r)r rrbr)rbrr9r9)rorzrrf r))r9/iZ?s0?gz5F?)r9~H?|/mg>m?)r9.d5 ?L ,ٿg깃\?)r9>?ЀjoؿgA]4?)r9pa2?Յ пgqQ?) )rrT rU gz¶?)rrV rW gV?)rrX rY gRE)3\?)rrZ r[ gB?)rr\ r] g/^?)-C6rT rU gH?)r^ rV rW gn?)r^ rX rY gx\?)r^ rZ r[ gVM%?)r^ r\ r] gm?)-rrrrrZrr\r r r r rrrr-rGrr r xfailr r r r rTrWr r r$ r- r+r*rr4rrnrfr4 r: r= rB rN rQ rrr s@r.r r s5 ^^D!B"B ^^B ^^))V ^^   [[ [[ LL fll4v{{7G7GH   [[/$> [[C [[]VV,<=>? $ [[ [[Vc1X. / -0 [[>? [[    [[D$<  @   , , G [[ [[* LL  LL(#(kk&&   LLD CB' B2&kk''  - *~+,~@ [[* LL  LL(#(kk&&   LLD CB' B2&kk''  - *z+*zx [[Wq!f- [[VeU^415.1@ [[ AJ !  @ [[ !RVVRVVRVV!< =  / 00 [[W&78 [['!!%% CB'  !!%% CB'  !!%% Ar2&  !!%% Ar2&  ' : ;9<  [[& )8 98 0 [[ [[   !" 0 ^^  /.//r0r c,a]tRtRtoRtRtRtVtR#)TestArrayArgumenticN\PPR4VnR#)l4cNrSrXs&r.rZTestArrayArgument.setup_methodr\r0c\PP\P!^4\P !^4RVP R7p\VPR4R#)rrNr) r4r8rrTr'rrWrr?r's& r.test_noexception"TestArrayArgument.test_noexceptionsEjjnn"))A,rwwqz")BSYY(r0rN)rrrrrZre rrr s@r.ra ra s5))r0ra c,a]tRtRtoRtRtRtVtR#) TestDocstringicJ\PPe5\R\PPP 494\P Pe7\R\P PP 494R#R#)Nrayleighr)r4rj rQrlowerrrXs&r.test_docstringsTestDocstring.test_docstringssg >> ! ! - J%.."8"8">">"@@ A ?? " " . K5??#:#:#@#@#BB C /r0cZ\P!4\P!4R#rt)r4rCrBrXs&r.test_no_name_argTestDocstring.test_no_name_args  r0r6N)rrrrrl ro rrr s@r.rh rh sDr0rh c\.RO4p\V^8V^4wr\V.RO4\V^.4\^^8V^4wr\W4\V^.\P!V4,4\V^8V^4wr\W4\V^.\P!V4,4R#)rN)rr!r rr r!r!)r!r r r!r!)r rrrTrr,r-rs r.test_args_reducers s #$A a!eQ "DAq/*q1# a!eQ "DAqq1# *+ a!eQ "DAqq1# *+r0cfa]tRtRto.ROtRt.ROt]PPR] 4R4t Rt Rt RtRtR tR tR t]PP&]PPR R R.4R44t]PPR R R.4R4tRtVtR#) TestFitMethodi,cN\PPR4VnR#)l=NrSrXs&r.rZTestFitMethod.setup_method0r\r0z dist,argsc WP9d\P!V R24\P!RRRRR\P .4p\P!RRRRR\P .4p\\V4p\\VPV^R7\\VPV^R7R#) rKz fit known to fail or deprecatedr-r.r/r0r1r#N) fitSkipNonFiniterrrTr r*r+r2r4rr)r)rYrJr4r^rdistfuncs&&& r.!test_fit_w_non_finite_data_values/TestFitMethod.test_fit_w_non_finite_data_values6s (( ( KK4& @A B HHfffffbffE F HHfffffbffE F5$'j(,,!<j(,,!.mlefuncsAAB!"BKK&Eufw{{1~566ufw{{1~5668DKr0rzrrr)rrrr)rrr)rrr)rrr)rrrrN)rrorzr) rTr r4r rrrrrr))rYr r^r,r-rcrddadbrrb2rrrs& r.test_fix_fit_betaTestFitMethod.test_fix_fit_beta~s  HH' ( ::>>!!A>>cS!UAa(1a&t< HH' ( ::>>!!>DcQS!UAq!AD)U$zz~~bQQq~IRT1VQ$AD)B" j%**..!#aH HHZ j%**..!!AFj%**..!!A!Lj%**..!!A!L j%**..!aQ (r0c\P!.RO4p\PP V4wr#\ V^4\ V^4\PP V^R7wr#\ V^4\ V^4\PP V^R7wr#\ V^4\ V^4R#)r r#rN)r r rrrrrr)rTr r4rrrr s& r.test_expon_fitTestFitMethod.test_expon_fits HH- .[[__Q' S!UA[[__Qq_1 S!UA[[__QQ_/ S!UAr0c ^\P!.RO4p\P!V^, 4p\PP V^R7wr4p\ W2P4RR7\V^4\ V\P!VP44RR7\PP V^^R7wr4p\ V\P!V\P!^4, ^,P44RR7\V^4\V^4\PP V^RR7wr4p\VR4\V^4\ V\P!VP44RR7R#) rrrrrr)rfix_sN)rr!rbr;) rTr rr4rrrrrrrr)rYr^lnxm1r?rcrds& r.test_lognorm_fitTestFitMethod.test_lognorm_fits+ HH- .q1u !MM--aa-8Eyy{7S!rvvejjl3%@!MM--aa-BErww):Q(>'D'D'FG" $S!UA!MM--aat-DEUD!S!rvvejjl3%@r0c\P!.RO4p\PP V4wr#\ W!P 44\ V\P!V44\PP V^R7wr#\ V^4\ W1P44\PP V^ R7wr#\ V^4\ V^ 4\\\PPVRR7\\\PPVRR7R#)rSrr#rrTN)rSrRrr4) rTr r4rrrrrr&rr)r s& r.test_uniform_fitTestFitMethod.test_uniform_fits HH) *]]&&q) S%%'"UBFF1I&]]&&qq&1 S!UEEG$]]&&q&4 S!UBj%--"3"3QSAj%--"3"3QsCr0rMLEMMc PRRr2\PPW#^dVPR7p\PP VRVR7p\PP VRVR7p\ WVRRR7\PP VRVR7p\ WVRRR7\PP VRVR7p\PP VRVR 7p\ WxRRR7\PP VRVR 7p\ WxRRR7\ \\PPV^^VR 7\ \\PPV^^^^VR 7\PP VR^^VR 7p V wrr\WV ..RO4Rp\PPV^dVPR7p\PP WVR7wrp \W4R#)r r3r)rr)rrrr)rr)rr)rmr)rlr)rrr)rrrrr)rrrrN)r r9r) r4r rrWrrrr)rr)rYrr,r-r^res_1res_2res_3res_4res_5aabbrssrs&& r. test_fshapesTestFitMethod.test_fshapess 21 JJNN1cN A qR7 qR75u= q6:5u= qR7 qR75u= q6:5u= j%**..!aO j%**..!aQv 7  qRa&Ibb\:. {{qsB[[__T_? Rr0c\PpVP^^dVPR7p\ RR7p\ \ VPV3/VBRV/BR#)r )r_rr) enikibenikirNi)r4r^rrWr%rr,r)rYrrJrr3s&& r.test_extra_paramsTestFitMethod.test_extra_paramssIxx!#DHHx=t$i4F3FvFr0rN)ncfr"r# irwinhall)rr8rr )rrrrrrZry rrrrr{ r r r r r r r rr r rrr s@r.ru ru ,s 1D5A [[[(3=4=@ $3L5(n A(D$ [[ [[Xt}5$6$L [[Xt}5G6Gr0ru cPa]tRtRtoRtRtRtRtRtRt Rt R t R t Vt R #) TestFrozenic\Pp\P!RRR7pVPR4pVPRRRR7p\W44VP R4pVP RRRR7p\W44VP R4pVP RRRR7p\W44VP R4pVP RRRR7p\W44VPR4pVPRRRR7p\W44VP4pVPRRR7p\W44VP4pVPRRR7p\W44VP4pVPRRR7p\W44VP4pVPRRR7p\W44VP4pVPRRR7p\W44VP^4pVP^RRR7p\W44\VPVP4\VPVP4R#)rr rbr~ roN)r4r8rfrrnrrrrrrrrErr,r-)rYrJrI result_frs& r.r_TestFrozen.test_norm!szzC0::d#$D4X&::d#$D4X&::d#$D4X&::d#$D4X&99T?4s3X&==?S1X&;;=t3/X&::<d#.X&::<d#.X&>>#$c2X&==#QD4X&VXXtvv&VXXtvv&r0cRp\Pp\P!V4pVPR4pVPRV4p\WE4VP R4pVP RV4p\WE4VP R4pVP RV4p\WE4VP R4pVP RV4p\WE4VPR4pVPRV4p\WE4VP4pVPV4p\WE4VP4pVPV4p\WE4VP4pVPV4p\WE4VP4pVPV4p\WE4VP4pVPV4p\WE4VP^4pVP^V4p\WE4\VPVPP4\VP VPP 4R#)rr~ rorN)r4rrfrrnrrrrrrrrErr,rJr-)rYr,rJrI r rs& r. test_gammaTestFrozen.test_gammaTs {{Q::d#$"X&::d#$"X&::d#$"X&::d#$"X&99T?q!X&==?QX&;;=1X&::<!X&::<!X&>>#aX&==#Q"X&VXXv{{}}-VXXv{{}}-r0c\P!^4pVP^4pVPRR7VP^4p\W#4R#)rrrN)r4rrr)rYrI m1rhs& r.test_regression_ticket_1293&TestFrozen.test_regression_ticket_1293sFq! ]]1  V $]]1  Rr0cRp\P!VR7pVPPV4wr4\ W4.RR.4Rp\PP ^VR7\ VPPV4^\ P.4Rp\P!VR7pVPPV4wr4\ W4.RR.4Rp\PP ^V4\ VPPVPP3\PPV44\P!RR7p\VPVPJ4RFp\ P!V4p\P!VR7pVPVPrC\ VR4\\ P!V44\ P!R4p\PPV4wr4\W4.RR.4K R#) rrrrrzNrrr)r4rrJrrrfrTr+r,r-rrrr)rYrrr,r-rv1s& r.rTestFrozen.test_abs  __q !ww##A&aVb#Y'  A#RWW))!,q"&&k:  __q !ww##A&aVb#Y'  Aq!bggii+U__-I-I!-LMoo$'(A 1 A1%B44q B  BKKN # 3A??//2DA QFRI .r0cL\\\PR44R#) rv_frozenN)rr+r4r5rXs&r.test_rv_frozen_in_namespace&TestFrozen.test_rv_frozen_in_namespaces++[9:r0cj\P!4p\\VR44^*Vn\ VPP 4\PP^*4P 44\PPR4pVP^VR7R#)rrrN) r4r8rr+rr get_staterTrUrEr)rYrI rndms& r.test_random_stateTestFrozen.test_random_states/0!V((224YY**2.88: <yy$$T*  -r0c\P!RR4p\P!R4p\P!.RO.R O3R7pWV3FpRVnVP ^R7\ P!V4pVP ^R7p\ P!V4pVP ^R7p\Wh4VPR4VPR4.p \V ^,V ^,4\VPV ^,4VPV ^,44K R#) gdsz@g@e?r r}rrTrzN)r9rr r!r ) r4r rrBrrpickledumpsloadsrrrn) rYr poissrudistfnr?r0 unpickledr1medianss & r. test_picklingTestFrozen.test_picklingszz,.AB b!""<+?+ABF+F"&F  JJAJ  V$A#B QIA&B  zz# c(:;G WQZ 0 GAJ/"wqz2 4,r0c Rp\P!^^^R7p\P!RRR7;_uu_4VP V^^RR7p\PP VR ^^^^RR7pRRR4\ XX4\P !^^R7pVP V4p\P P VR ^R 7p\ Wg4R# +'giLt;i) cV#rtr6rs&r.r$TestFrozen.test_expect..funcsHr0)r,rcrdrz)r|divideT)rub conditionalr4rcrdrr r NrEr4rcrr!)r4rrTrrrr)rYrgmgm_val gamma_valrzp_val poisson_vals& r. test_expectTestFrozen.test_expects [[1!1 - [[( ; ;YYtaTYBF **4d./A4+II<  * MM! #mm**4d*B +< ;s .s1qs r0rbrUcV#rtr6rs&r.r@r sr0rjrc^#rr6rs&r.r@r Qr0rcrdrr rc^#rr6rs&r.r@r ar0Trcrdrr r rSN)r4r8rrr)rYrGr rr prob90prob90cs& r.r_TestExpect.test_norms JJ  3!  DAq"- JJ  m!  <Aq"- ZZ^^Daq^ 1 ZZ^^Daq^ 1"";AQ2"MFC4**##KQaB04$6GR4r0c .\PPRR^^R7p\VR^ R7\PPRRRRR7p\VR^ R7\PP R^ ^ ^^R7p\PP R ^ ^ ^^R7p\PPR RRRWCR R 7p\VR ^ R7\PPRR^^WCRR 7p\VR^ R7R#)c0VR, VR, ,#)UUUUUU@r6rs&r.r@&TestExpect.test_beta..s51U7(;r0)r4rcrdrSrUcV#rtr6rs&r.r@r sr0rTrrrbrjcR#)rSr6rs&r.r@r sRr0Fr rc^#rr6rs&r.r@r $r r0TNrqq?r )rbrb)r4r rrr)rYrGr r rr r s& r. test_betaTestExpect.test_betas JJ  ;'"#1  .Avr2 JJ  kRr  JAub1 ZZ^^D"baq^ 9 ZZ^^D"baq^ 9"".-r r0r rUc"VR, ^,#r4r6rs&r.r@r 0s adQYr0c"VR, ^,#r r6rs&r.r@r 5s QrTAIr0)r4rcrr c^#rr6rs&r.r@r <sqr0c^#rr6rs&r.r@r As1r0T)r4rcrr r c^#rr6rs&r.r@r Fs!r0)r4rr N)rQrbr)r4rrrrNrp) rYm_truev_truer rGv_bounds prob_true prob_boundsprob_bcprob_bs & r.test_hypergeomTestExpect.test_hypergeom(sn..r2qb.A OO " ";[b " IAr2 OO " "#6[') # +Ar2??))*=/:.0Q2*? Hb9eoo))1b'2r1!)DHHJJ oo,,[{13b-B KB?//((;B,-"$)HGQ3'' ++,(4FAr2r0c\PPRR^RR7p^\PP^^4, p\ W^R7\PPRR^RR7p\ V^^R7R#)c^#rr6rs&r.r@)TestExpect.test_poisson..LsQr0F)r4rr rUc^#rr6rs&r.r@r Qsr0TNr)r4rrrnr)rYr prob_b_trueprob_lbs& r. test_poissonTestExpect.test_poissonJswmm**;Ta7<+>  ))!Q// KbA--&&{!37'9GQ3r0c\PpVPRR7pVPRR7VPRR7p\W#^R7R#)rrrUN)r)rz)r4genhalflogisticrr)rYhalflogrIrMs& r.test_genhalflogisticTestExpect.test_genhalflogisticUsG''~~6~*F#~~6~*D3r0c \\P!\PP RR444\\P!\PP RRR744\\P!\PP RRR744\\P!\PP RRR744R#)rGz?c^#rr6rs&r.r@/TestExpect.test_rice_overflow..dr0rc^#rr6rs&r.r@r) er* r0c^#r r6rs&r.r@r) fr* r0N)r' )rrTrr4ricerfrrXs&r.test_rice_overflowTestExpect.test_rice_overflow_s  EJJNN3567 EJJ--k-HIJ EJJ--k-HIJ EJJ--k-HIJr0c<R^r!\PPRV3R7p\VWR, , \P !RV, 4, RR7\PPRV3VR7p\WCV,RR7R#) rcV#rtr6r]s&r.r@(TestExpect.test_logser..kr r0rrSrrcV#rtr6r2 s&r.r@r3 qr r0r N)r4r>rrrTr)rYrzrcres_0res_ls& r. test_logserTestExpect.test_logserhsxa3 ##Kqd#;V rvvb1f~5E C ##Kqd#Ds{7r0c^^r!\PPRW3R7p\PPRW3R7p\W1V, RR7\WC^,, W,RR7R#)cV#rtr6rs&r.r@)TestExpect.test_skellam..xAr0rcV^,#rr6rs&r.r@r< yrr0rrN)r4rrr)rYrp2r rhs& r. test_skellamTestExpect.test_skellamtscRB ]] ! !+RH ! = ]] ! !.x ! @G%0U BG%8r0c ^^qr!\PPRW34p\V\ R\ W444W!, , RR7R#)r9cV#rtr6rs&r.r@)TestExpect.test_randint..r r0c3$"TFqxK R#5irtr6rrs& r.r*TestExpect.test_randint..s5}!A}srrN)r4r>rrrpr)rYrJrKrs& r. test_randintTestExpect.test_randint}sICBmm"";95uR}55A Or0c\P!\4;_uu_4\PP RR4RRR4R# +'giR#;i)cV^,#rr6rs&r.r@&TestExpect.test_zipf..s1r0Nr)rwarnsr~r4rVrrXs&r. test_zipfTestExpect.test_zipfs4 \\. ) ) JJ  nd 3* ) ) )s "A A c\PPRRR7p\PPRRR^ RR7p\W^R7R#) c^#rr6rs&r.r@/TestExpect.test_discrete_kwds..r= r0rc^#rr6rs&r.r@rR r= r0rr)r4maxcount chunksize tolerancerUNr)r4rrr)rYn0n1s& r.test_discrete_kwdsTestExpect.test_discrete_kwdssJ ]] ! !+D ! 9 ]] ! !+D+/2"OBB/r0cRpRF6p\PP^V4p\W1!V4RR7K8 R#)cV^,^ V^,,,^V^,,,^V^,,,V,#r r6rs&r. poiss_moment5-TestExpect.test_moment..poiss_moment5s8q52b!e8#bQh.BE9B> >r0rrN)rr)r4rrr)rYr] rm5s& r. test_momentTestExpect.test_moments: ?B%%a,B B b 1 >r0c\\PP^$RR7^$4\\PP^(RR7^(4\\PP^ RR7^ 4\\PPRR7^4\\P P^UR7^U4R#)rJrSrbrrrEN))rr4r8rrrrXs&r.test_challenging_cases_gh8928(TestExpect.test_challenging_cases_gh8928s  ))b)>$2!IIOOD!4 #56H!I(*C" $)" $99>>&!,IIOOFA. ,.@A"3%@ =r0c n\\PP^^\ R44^4R#)r Ni)rr4rl rnrrXs&r.test_cdf_large_ncTestNct.test_cdf_large_ncs! aE%L91=r0zx, df, nc, expectedc^\\PPWV4VRR7R#r)rr4rl rf)rYr^rrw rNs&&&&&r.test_pdf_large_ncTestNct.test_pdf_large_ncs  aR0(Gr0r6N)rrbrgAs69)rrrgV1)rrrrrm rq rt rx r r rrrr rrr s@r.rj rj sWC2$  >,>6 [[2.+5  H  Hr0rj c8a]tRtRtoRtRtRtRtRtVt R#)TestRecipInvGaussicX\PP^R4pVR8XgQhR#)r9r`rN)r4 recipinvgaussrfrs& r.test_pdf_endpoint#TestRecipInvGauss.test_pdf_endpoints%    # #As +Cxxr0cv\PP^R4pV\P)8XgQhR#)r9r`N)r4r rrTr+rs& r.test_logpdf_endpoint&TestRecipInvGauss.test_logpdf_endpoint s,""))!S1wr0cd\PPRR4pRp\WRR7R#)rjrzg9V9s;r rN)r4r rnrrYrzrNs& r.test_cdf_small_x"TestRecipInvGauss.test_cdf_small_x$s+    # #D# .(%0r0cb\PP^PR4pRp\WR4R#)r@ rzglH EJJNN1N3488:; EJJ--a2-67;;=>  EJJNN1N3488:;zz62 D!%%'(  q!,ejjnnQ.BQ (r0c\PPp\V!RVPR7P ^4\V!RRVPR7P R4R#)r )r-r)r-rrN)r!r)r4r- rrrWrr?r's& r. test_rice_rvsTestRice.test_rice_rvsXsGjjnnS2DHH5::A>S2FBHH&Qr0c\PP\P!^ ^^ 4\P!^ ^^ 44p.ROp\ W4\P!R^R4p\PP VR^R7p.ROp\ WE4\PPP R\P!^ ^^ 44p.ROp\ WE4R#) rbrr^rzN)g/ Ƹ?g[jNJ\?g ?g9LJ?g#O?gXx?glO?g٥a%?gVzӯ?gTQ+xQ?g?`J?go?gbu=?gIo|?g"Z6?g@_@) g<~@gwd @gC?@gb 0@g` A@@g'NxP@g33-a@g4pv@g$F@)gJ-$@g+f4@gce*D>@g$D@gnHGI@grN@g< uQ@giffT@g9r[V@g3bXQY@gq xJ[@g`*DD^@g@`@gdAa@)r4r- rnrTr'rrr3)rYrncdf_exp probabilitiesrppf_exps& r.test_rice_gh9836TestRice.test_rice_gh9836]sjjnnRYYr33RYYr35KL O % #q#. jjnn]En; L %kk""3 "c2(>? ; %r0rN) rrrrrZr r r rrr s@r.r r ?s4(*R '&'&r0r c2a]tRtRtoRtRtRtRtVtR#) TestErlangicN\PPR4VnR#)lILNrSrXs&r.rZTestErlang.setup_methodr\r0c \P!4;_uu_4\P!R\4\ \\ P PR^^^VPR7.ROp\ P PV^R7p\ PPV^R7p\W#RR7RRR4R# +'giR#;i)rwrg rrwr2rN)rzrSrr3) r{r|r}r~rr4erlangrrWrrr)rYr result_erlang result_gammas& r.test_erlang_runtimewarning%TestErlang.test_erlang_runtimewarnings $ $ & &  ! !'> : .%,,*:*:CQ ! B (D!LL,,Ta,8M ;;??4A?6L Md C' & & &s B'C C c\\PPR^R.R7\PPR^R.R74R#)rzrZNr)rr4r rfrrXs&r.test_gh_pr_10949_argcheck$TestErlang.test_gh_pr_10949_argchecks<U\\%%caW%5[[__SQG_4 6r0rN) rrrrrZr r rrr s@r.r r s5D&66r0r ca]tRtRtoRtRtRt]PPRRR.4R4t ]PPRRR.R R ..4R 4t R t R t RtVtR#) TestRayleighicN\PPR4VnR#)l@XNrSrXs&r.rZTestRayleigh.setup_methodr\r0c\\PP^24p\VR4R#)r=Ng)Zx)r4rj rrrYrs& r.rTestRayleigh.test_logpdfs! NN ! !" %./r0c\\PP^24p\VR4R#)r=Ni)r4rj rrr s& r.rTestRayleigh.test_logsfs NN  $5!r0rc\PP^VW PR7pRpV!W14p\PP W1R7wrg\ Wa4\ Wu4\PP VRR7wrg\ VR4\PP V4wrg\ Wt!W644R#)rrc\P!W, ^,4^\V4,, R,#r rz)rTrpr#)rrs&&r.rz(TestRayleigh.test_fit..scale_mles*FFDKA-.!c$i-@RG Gr0rr`r#N)r4rj rrWrr)rYrrrrz scale_expectrcrds&&& r.r TestRayleigh.test_fits~~!!s(1"J H!/ ^^'''; S"U)^^''R'8 UB ^^''- UId01r0r' rg̡=E?g@1d?c\PP^VW PR7p\ \PV4R#)rrN)r4rj rrWr)rYrrrs&&& r. test_fit_comparison_super_method-TestRayleigh.test_fit_comparison_super_methods7 ~~!!s(1"J%ennd;r0c8\\P4R#rt)r0r4rj rXs&r.rTestRayleigh.test_fit_warningssENN+r0c\PPR4p^2RRrCp\PP W#WAR7p\PP V4wrgV\P!V48gQh\PP WSR7wrhV\P!V48gQhW8XgQhR#)irrrr#N)rTrUrVr4rj rrr) rYrWrcrdrrrrrs & r.test_fit_gh17088TestRayleigh.test_fit_gh17088sii##C(sCDnn  $ I^^'', $$$"^^///B$$$!!!r0rN)gB)Q?got|?)gS2P?g3?)rrrrrZrrrrrr r rr rrr s@r.r r s50" [[03K3K3MN2N20 [[0D$<4> 3K3MN<N<, " "r0r cRa]tRtRtoRtRtRt]PPR.R O4R4t ]PPR.RO4R4t ]PPRRR.4R 4t ]PPRRR.4R 4t R tVtR #) TestExponWeibicRpRpRp\PPWV4p\PPWV4p\ WE.RR.4R#)rrSr g:23W+Ng~`9k)r4 exponweibrfrr)rYr^r,rrzrs& r.test_pdf_logpdfTestExponWeib.test_pdf_logpdfsS    OO  a (%%aA. /1DE Gr0cf\P!RR^4p^p^dp\PP WV4p\P P W4p\ WE4\PPWV4p\P PW4p\ We4R#)rNrr)rTrr4r rf weibull_minrrrYr^r,rrzrNrs& r. test_a_is_1TestExponWeib.test_a_is_1s KKB "   OO  a ($$((.$%%aA.$$++A1'r0cf\P!R^^ 4p^p^p\PP WV4p\P P V4p\ WE4\PPWV4p\P PV4p\ We4R#)rNrX)rTrr4r rfrrrr s& r.test_a_is_1_c_is_1 TestExponWeib.test_a_is_1_c_is_1s KKAr "   OO  a (;;??1%$%%aA.;;%%a('r0z x, a, c, refc`\PPWV4p\WTRR7R#r)r4r rr)rYr^r,rrrs&&&&& r.rTestExponWeib.test_sfs$__  a (e,r0z p, a, c, refc`\PPWV4p\WTRR7R#r)r4r rr)rYrzr,rrrs&&&&& r.rTestExponWeib.test_isf>s$oo!!!*u-r0c`\PPWV4p\WTRR7R#r)r4r rr)rYr^r,rrrs&&&&& r.rTestExponWeib.test_logcdfKs$''a0%0r0c`\PPWV4p\WTRR7R#r)r4r rr)rYr^r,rrrs&&&&& r.rTestExponWeib.test_logsfSs$%%aA./r0r6N))rrrg=J ?)r=rrg$+PR>)}rrg ><)rrrggX.X:)r?rgRy9)r!r r g'vR0=)g_QPrrgQ:<))rorrg;y$(@)gMg [[ 5 - -* [[ 4..  [[^AGIJ1J1  [[^IFHI0I0r0r c,a]tRtRtoRtRtRtVtR#)TestFatigueLifei[cb\PPRR4p\VRRR7R#)rrX9rrN)r4 fatigueliferrrs& r. test_sf_tailTestFatigueLife.test_sf_tail]s)     ,0u=r0cfRp\PPVR4p\VRRR7R#)r rrrrN)r4r rr)rYrzrs& r. test_isf_tailTestFatigueLife.test_isf_tailks+ !    ! !!S )5u-r0r6N)rrrrr r rrr s@r.r r [s >..r0r ca]tRtRtoRtRt]PPRRR.4R4t Rt ]PPR R R.4R 4t R t Vt R #) TestWeibullirc^\PP^^4p\V^4R#r)r4r rrr s& r.rTestWeibull.test_logpdfts$    $ $Q *Qr0cRpRpRp\PPWVR7p\V\P !R4^, 4\PP WVR7p\VR\P!^4, 4\PPWVR7p\V\P!R4)4\PPWVR7p\V\P!\P!R4)44\PPWVR7p\V\P !R44\PPWVR7p \V R4\PP^^^R7p\V\P !R44\PP^^^R7p \V R4Rp\PPWVR7p\V\P !R4^, 4\PP WVR7p\VR\P!^4, 4\PPWVR7p\V\P !R44\PPWVR7p\VR4\PPWVR7p\V\P!R4)4\PPWVR7p \V \P!\P!R4)44\PPR^^R7p\V\P!R 4)4\PPR^^R7p \V \P!\P!R 4)44R#) rrr r^Nrrrg& .gCe)r4r rfrrTrrrrnrr rrr weibull_max) rYr^r,r-rzrrlcr?r s & r.test_with_maxima_distrib$TestWeibull.test_with_maxima_distribys'.       ! !!a ! 0266%=?+    % %a! % 4EBFF1I-.    ! !!a ! 0GMM%001    % %a! % 4BFFGMM%$8#89:    Q /266%=)    $ $Q $ 3E"    Qa 0266$<(    $ $R! $ 4D!     ! !!a ! 0266%=?+    % %a! % 4EBFF1I-.    ! !!a ! 0266%=)    % %a! % 4E"    Q /GMM%001    $ $Q $ 3BFFGMM%$8#89:     3GMM*@AAB    $ $UAQ $ 7BFFGMM2H$I#IJKr0rdrSrc\PPVR,V^,^VR7p\VR4R#)rYr^g>N)r4r rFrrF s&& r.rGTestWeibull.test_delta_cdfs;"!!,,U3Ya38-:56r0c\PPR4p^RRrCp\P!W#V4pVP ^dVR7p\PP VR^R7wrxp \PP VR^R7wrp Yu;8Xd ^8XgQhQhWz8wgQh\PP V^^RR7wrpV^8XgQh\P!WV4pVPR R 7p\P!V4\P!V43p\VV4R #) l>rrzrrrr( rrrrN) rTrUrVr4r rrrrr)rYrWrrcrdrJrc2rrc3rrc4rrrrrs& r. test_fit_minTestWeibull.test_fit_mins$ii##$783  /hhCch2!,,00c0B& ,,00c0B& q     xx!,,00a$0O&{{!!"F3kk$k'ggclEJJsO+S!r0rc\\PPW4VRR7\\PP W24VRR7R#r)rr4r rrrs&&&&r.ruTestWeibull.test_sf_isfs> )),,Q2CeD))--c5quEr0r6N)r=rr )rrgq뮕qH))rrrrrr rrrrGr rurrr s@r.r r rsx RLh [[WsCj1727("< [[[+J,D+EFFFFr0r c,a]tRtRtoRtRtRtVtR#) TestDweibullicD\PPR4p^ VP^d^ R7,p\P P V4p\PP V4\P!R4, p\W4RR7R#lq>~).rV rzrrN) rTrUrVrW r4dweibullrEr rr)rYrWrrrs& r.rtTestDweibull.test_entropy sqii##$78  2 . .nn$$Q'''*RVVC[8u-r0cJ\PPR4p^ VP^^ R7,p^ VP 4,p\ P PW24pR\ PPW24,p\WERR7R#r ) rTrUrVrW rr4r rr r)rYrWrr^rrs& r.rTestDweibull.test_sfsyii##$78   , ,  nn%E%%((..u-r0r6N)rrrrrtrrrr s@r.r r s ...r0r cVa]tRtRtoRtRtRtRtRtRt Rt R t R t R t VtR #) TestTruncWeibulli!cj\PPRR.RRR4p\VRR.4R#)rrg)\(?gףp= ?rN)r4truncweibull_minrfrr s& r.test_pdf_bounds TestTruncWeibull.test_pdf_bounds#s0  " " & &Sz3d CQc #r0c\PPRRR\P4p\ VR4\PPRRRR4p\ VR4R#)rrSrr3g~?N)r4r rrTr+rrr s& r.rTestTruncWeibull.test_logpdf(sU  " " ) )#sC @Q  " " ) )#sC =./r0cj\PPRR.RRR4p\VRR.4R#)rrSrrN)r4r rrr s& r.test_ppf_bounds TestTruncWeibull.test_ppf_bounds0s0  " " & &Sz3S AQc #r0c.ROp\PPVRRR4p\PPVRRR4p\ W4R#rrr N)rrrorzrrrS)r4r rrnrrYrr^q_outs& r.test_cdf_to_ppf TestTruncWeibull.test_cdf_to_ppf5sG 0  " " & &q"b" 5&&**1b"b9!r0c.ROp\PPVRRR4p\PPVRRR4p\ W4R#r )r4r rrrr s& r.test_sf_to_isfTestTruncWeibull.test_sf_to_isf;sG 0  " " & &q"b" 5&&))!RR8!r0ca a a aRo Ro Ro V V V 3Rlo\PP^S S S 4p\VR4\PP^S S S 4p\ V3RlS S 4wr4\ W#4\PP^S S S 4p\ V3RlS S 4wrd\ WV4\PP^S S S 4p\ V3RlS S 4wr\ Wx4\PP^S S S 4p \ V3RlS S 4wr\ W4R#) rrSr cb<W,\PPVSSS4,#rt)r4r rf)r^ryr,r-rs&&r.xnpdf)TestTruncWeibull.test_munp..xnpdfFs'4..221aA>> >r0c<S!V^4#rr6r^r s&r.r@,TestTruncWeibull.test_munp..M a r0c<S!V^4#rr6r" s&r.r@r# Qr$ r0c<S!V^4#r r6r" s&r.r@r# Ur$ r0c<S!V^4#)rr6r" s&r.r@r# Yr$ r0N)r4r rrrr)rYrCr m1_expectedrrh m2_expectedm3 m3_expectedri m4_expectedr,r-rr s& @@@@r. test_munpTestTruncWeibull.test_munpAs    ? # # * *1aA 6R  # # * *1aA 63Q: (  # # * *1aA 63Q: (  # # * *1aA 63Q: (  # # * *1aA 63Q: (r0c FRpRpRp\P!^\P!R\P!R\P!R4,)4,4, 4p\P P WCW4p\VR4\P PWCW4p\V\P!R4)4\P PRW1V4p\Wt4\P PWCW4p\VR4\P PWCW4p \V \P!R4)4\P PRW1V4p \W4R#)rSr rrzrN) rTrrrr4r rnrrrrrr) rYr,r-rx_medrnr rrr rs & r.test_reference_values&TestTruncWeibull.test_reference_values\s*   BFF3"rvvbz/0B)C#CDDE$$((18S!  # # * *5Q :RVVBZK($$((aA6#  # # & &u 6C  # # ) )%A 9RVVBZK($$((aA6#r0cRpRpRp\PpRp\PP WVR7p\P P WW4VR7p\ Wg4\PPWVR7p\P PWW4VR7p \ W4\PPWVR7p \P PWW4VR7p \ W4\PPWVR7p \P PWW4VR7p \ W4\PPWVR7p\P PWW4VR7p\ W4\PPWVR7p\P PWW4VR7p\ VV4\P P^^W4^R7p\ V\P!R44\P P^^W4^R7p\ VR4R#)rrrr r^Nr) rTr+r4r rfr rrrnrrrr)rYr^rr,r-rdrzp_truncrlp_truncrn cdf_truncr lc_truncr?s_truncr ls_truncs& r.test_compare_weibull_min)TestTruncWeibull.test_compare_weibull_mints    FF    ! !!e ! 4((,,Q1u,E#    % %a% % 8))00q50I%##A#6**..qQ.G '    % %a% % 8))00q50I%    U 3((++A!e+D#    $ $Q $ 7))//a%/HH%  " " % %b!Q % ;266$<(  # # ) )"aQ ) ?D!r0cRRRr2p\P!W#^d4p\PP WAW#4p\PP WAW#4p\P P W14\P P W!4, p\P P WA4V, p\P P WA4\P P W!4, V, p \PPWX4\PPWi4R#)rrorPN) rTrWr4r rfrnr rIr) rYrr,r-r^rXrr8rYrs & r.test_compare_weibull_min2*TestTruncWeibull.test_compare_weibull_min2stTa KKc "%%))!5%%))!5  $$Q*U->->-B-B1-HH  $$Q*T1!!%%a+e.?.?.C.CA.II4O ""4. ""4.r0r6N)rrrrr rr r r r- r1 r: r= rrr s@r.r r !s8$ 0$ " " )6$0'"R//r0r cta]tRtRtoRtRt]PPR.RO4R4t Rt Vt R#) TestRdistic\Pp.ROp\VPVP VR4R4V^R7R#)r2g@rUNr)r4rdistrrnr)rYr rvs& r.test_rdist_cdf_gh1285TestRdist.test_rdist_cdf_gh1285s5$FJJvzz&%'@%H"A /r0c&\P!RR^ 4pRp\R\P!V^, V^, 4P V^,^, 4,\P !V4P V44R#)rg@rzNgGz)rTrWrr4r rfrB )rYr^rs& r.test_rdist_betaTestRdist.test_rdist_betasd KKtR ( C 1Q3! 4 8 8!a% CC!KKN..q1 3r0rc\\\PPW4VRR7R#r)rr4rB rrs&&&&r. test_rdist_sfTestRdist.test_rdist_sfs  q,c>r0r6N))rig?)rg;|bZ?)rzig79)riUg޻S/) rrrrrC rF rrrrI rrr s@r.r@ r@ s:/3 [[ ??r0r@ cDa]tRtRtoRtRtRtRtRtRt Rt Vt R #) TestTrapezoidicF.ROpVFp^V^.p\\PPW2V4\PPW244\\PP W2V4\PP W244K R#)r9N)r9rrzr)rr4rrftriangrn)rYmodesmoder^s& r.test_reduces_to_triang$TestTrapezoid.test_reduces_to_triangsv DD! A  3 3AT B % 0 0 9 ;  3 3AT B % 0 0 9 ; r0cV\P!^^^ 4p\\PP V^^4\P P V44\\PPV^^4\P PV44R#r)rTrWrr4rrfrrnrs& r.test_reduces_to_uniform%TestTrapezoid.test_reduces_to_uniformsl KK1b !EOO//1a8%--:K:KA:NOEOO//1a8%--:K:KA:NOr0c\\PP^^^4^4\\PP^^^4^4\\PPR^R4R4\\PPRRR4R4\\PPRRR4R4\\PPRRR4R4\\PPR RR4R4\\PP RRR4R 4\\PP RRR4R 4\\PP RRR4R4\\PP R RR4R 4\\PP RRR4R4R #)r9rzrrrwrSrg?rPrr rg?N)rr4rrfrnrXs&r. test_casesTestTrapezoid.test_casessfEOO//1a8!<EOO//1a8!<EOO//Q</ 1EOO//S#>/ 1 EOO//S#>FEOO//S#>EEOO//S#>FEOO//S#>HEOO//S#>FEOO//S#>DEOO//S#>HEOO//S#>Dr0ca a a a a Rwo o o o S S , S S , , S S , S S , , S S S , 3wrr4^S S ,S , S , , o V V V V V 3RlpV!^4pV!^4V^,, pRS S , S ,S , ,S S ,S , S , , \P!RS S ,S , S , ,4,p\\PP WW44V^ R7\\PP WW44V^ R7\\PPWW44V^ R7\\PP ^^R^4R^ R7\\PP ^^R^4^^ R7\\PP ^^R^4^^ R7R#)r!c <SSV^,,SV^,,, SS, , SV^,,SV^,,, SS, , , ,V^,, V^,, #rr6)ryr,r-rdrss&r.r6TestTrapezoid.test_moments_and_entropy..momentsn!ac(Q1X-!A#6!Hq1Q3x/AaC89:qSqS" #r0rzrUNr r)r rr r!)rTrrr4rrrrE)rYrr? rcrdrrrrEr,r-rr\ rss& @@@@@r.test_moments_and_entropy&TestTrapezoid.test_moments_and_entropys" 1acac]QqSQqSM1acA 1QqM # # ayQi$'!1Qq/QqSU1W-sac!eAg0GGEOO00D " .EOO//C -EOO33BCG#R 1 EOO00Ar1=r2NEOO00Ar1=q"MEOO//1b!r0c\PPR4pRp\P!V!pVP RVR7p\PP V4pVR,wrg\P!V^RR7'dQh\P!V^RR7'dQh\PPW$4p\PPWT4p W8gQhR#)l\"MWrrNr NrrN)rrarr ) rTrUrVr4rrrrr) rYrW true_args true_distr^ fitted_argsfitted_cfitted_dtrue_llf fitted_llfs & r.&test_trapezoid_fit_convergence_gh235034TestTrapezoid.test_trapezoid_fit_convergence_gh23503)s ii##$89# OOY/ MM$SM 1oo))!, )_;;x5555;;x5555??'' 5__))+9 $$$r0r6N) rrrrrR rU rX r^ ri rt rrr s@r.rM rM s+;P E(M6?4%%r0rM c&a]tRtRtoRtRtVtR#) TestTriangiBcn\P!RR7;_uu_4\\PP ^^4R4\\PP R^4R4\\PP ^^4R4\\PP ^^4^4\\PP R^4R4\\PP ^^4^4\\PP RR4R4\\PP RR4R4\\PP RR4R4\\PP RR4R4\\PP RR4R4\\PP RR4^4RRR4R# +'giR#;i) raiser' rrzrSrrroN)rTrrr4rO rfrnrXs&r.test_edge_casesTestTriang.test_edge_casesCsd [[W % % ))!Q/ 4 ))#q12 6 ))!Q/ 4 ))!Q/ 3 ))#q12 6 ))!Q/ 3 ))"b12 6 ))#r2D 9 ))#r2C 8 ))"b12 6 ))#r2D 9 ))"b11 5& % % %s G9H## H4 r6N)rrrrrz rrr s@r.rw rw Bs66r0rw ca]tRtRto]P P RR R .4R4t]P P R.R O4R4tRt Rt Rt Vt R #) TestMaxwelliVrc\\\PPV4VRR7R#r)rr4maxwellrrs&&&r.rTestMaxwell.test_sf[s  ((+Su=r0rc\\\PPV4VRR7R#r~)rr4r rrs&&&r.rTestMaxwell.test_isfds  ))!,c>r0cbRp\PP^ 4p\W!RR7R#)gi6u=>  [[X:;? ;?1 00r0r} c,a]tRtRtoRtRtRtVtR#) TestMielkeixcRRr!\\P!W4P^4\P 4\\P!VR4P^4\P 4\ \P!\P!VR4P^444R#)g rh@gv?rSrDN)rr4mielkerrTr+rr)rYr]r?s& r.rCTestMielke.test_momentsysve1U\\!'..q1266:U\\!S)003RVV< ELLD188;<=r0c\P!R^d^24pRRr2\\PP WW#, 4\P P WV44R#)rg@gHzG@N)rTrWrr4burrrfr )rYr^r]r?s& r.test_burr_equivalence TestMielke.test_burr_equivalencesG KKc2 &T1 qQS15<<3C3CA!3LMr0r6N)rrrrrCr rrr s@r.r r xs>NNr0r c8a]tRtRtoRtRtRtRtRtVt R#)TestBurric \PR^.\PR^.\PR^.\PR^.\PR R.\PRR.\PR R..pVUUUu.uF#wr#qBP!VP .VO5!NK% ppppVUUUu.uFwr#qfNK pppp\ WW4VUUUu.uF#wr#qBP!VP .VO5!NK% ppppVUUUu.uFwr#p\P!V4NK pppp\ WW4R#uupppiuupppiuupppiuupppi) rrzrSrNr)rzr )rrr )rrz)rr ) r4fiskr burr12rfr,rrrTr)rYrr%r&rr(r)r*s& r.test_endpoints_7491TestBurr.test_endpoints_7491sZZq ! ZZ1 % ZZ # ZZ1 % \\8S ) \\63 ' \\63 ')9== 1vvbdd#U#=9=>!5I9>!#/;?@4<2ayy&&4@AEF)=I266)$F!#/ >>AFs)E: E )E"EcR^r!\P!W4P4wr4RRre\W54\WF4R#)rTg4U?g&MI?N)r4r r)rYrr\ rvariancemean_hc variance_hcs& r.test_burr_stats_9544TestBurr.test_burr_stats_9544s<A1A)//1 23F&.r0c2R^r!\P!W4P4wr4\\P!V44\\P!V44R^r!\P!W4P4wr4\\P !V44\\P!V44R^r!\PP \P!.RO4W4wrVrx\\P!V44\\P!V44\\P!V44\\P!V44R^r!\PP .ROW4wrVrx\\P !V44\\P!V44\\P!V44\\P!V44R^r!\PP .ROW4wrVrx\\P !V44\\P !V44\\P!V44\\P!V44R^r!\PP .ROW4wrVrx\\P !V44\\P !V44\\P !V44\\P!V44R^r!\PP .ROW4wrVrx\\P !V44\\P !V44\\P !V44\\P !V44R#)rzrrrrNr)r4r rrTr rr5r ) rYrr\ rr e1e2e3e4s & r.test_burr_nan_mean_var_9544$TestBurr.test_burr_nan_mean_var_9544ssA1A)//1"#A1A)//1 D!""#A1))"((<*@!G    A1)),= B    A1)),= B  B   A1)),= B  B  B  A1)),= B  B  B  B r0ctRRr!.ROp.ROp\\PPW1V4VRR7R#)rTr r rN)rrrr)g+&?gI$_@gKT@gؗVA)rr4r r)rYrr\ rrs& r. test_burr_isfTestBurr.test_burr_isfs2C1 &# qQ/5Ar0r6N) rrrrr r r r rrr s@r.r r s!0( /(!TBBr0r ca]tRtRto]P P RR R .4R4tRt]P P R.R O4R4t Rt Vt R#) TestBurr12izscale, expectedch\PPRR^^VR7p\W2RR7R#)gjAgjAr^rrN)r4r rFr)rYrdrNrs&&& r.rGTestBurr12.test_delta_cdfs-( ''S!Qe'De4r0ch^^r!\P^, p^\P^,^, , p\P^,^ VR,,, p\PpW4WV.p\P!W4PR4p\ WRR7R#)r rrr rN)rTrgr*r4r r) rYrr\ rrrrrrs & r.test_moments_edgeTestBurr12.test_moments_edges} !1uuQw"%%(2+ouuaxCH%66$)ll1 &&v.u-r0z p, c, d, refc`\PPWV4p\WTRR7R#r)r4r rr)rYrzrr\ rr^s&&&&& r.test_isf_near_zeroTestBurr12.test_isf_near_zero s$ LL  Q1 %U+r0r6N)rSgd~)rgR톇T ))rrQrzggx4/@)gҶOɃ;rQrzgI8aS@)rror"gj\/@)g#aMR/ror"g!GA) rrrrrrrrGr r rrr s@r.r r sd [[.<;=>5>5( .  [[ 0,,r0r c Xa]tRtRto.ROt.ROt.ROt]P!]]]34t .ROt ^^^ ^^x]P.t .ROt ]!]!]!] ] ] 4] 44tR R!^^]PR3^^]PR3^^ ]PR3.t]P(P*R4t]P(P.R4t]P4P7]4tRt]!]P4PA]]44;_uu_4t!]"PF!]!4t$RRR4]P(PKR ]$R ,4R 4t&]P(PKR ]$R ,4R 4t']P(P.]P(PQR4]P(PKR ]$R,4R444t)]P(P*R4t*]P(P.R4t+]P(PKR]4R4t,]P(P.]P(PQR4R44t-]P(P.R4t.Rt/Rt0Rt1Rt2Vt3R# +'giEL;i)"TestStudentizedRangei gv"?gd暕?ga-A?cVPF7wrVwr4p\PPW%V4p\ WcRR7K9 R#)rrN)rr4studentized_rangernr)rYpvkr p_expectedrGr]res_ps& r.test_cdf_against_tables,TestStudentizedRange.test_cdf_against_tables; s>iiFC" J1++//a8E ED 9 r0cVPF7wrVwr4p\PPW5V4p\ WbRR7K9 R#)rrN)rr4r rr)rYr q_expectedrzrGr]res_qs& r.test_ppf_against_tables,TestStudentizedRange.test_ppf_against_tablesB s=#yyOCGA!++//a8E ED 9 )r0z&data/studentized_range_mpmath_ref.jsonN case_resultcdf_datacVR,pVR,pVR,VR,VR,3p\PP!V!p\WSVR,VR,R7R# src_case mp_resultrr]rG expected_atol expected_rtolrN)r4r rnrrYr r r qkvrs&& r.test_cdf_against_mp(TestStudentizedRange.test_cdf_against_mpN ^z* , smXc]HSM9%%))3/%o6%o6 8r0pdf_datacVR,pVR,pVR,VR,VR,3p\PP!V!p\WSVR,VR,R7R#r )r4r rfrr s&& r.test_pdf_against_mp(TestStudentizedRange.test_pdf_against_mpY r r0z+intermittent RuntimeWarning: invalid value. moment_datacJVR,pVR,pVR,VR,VR,3p\P!RR7;_uu_4\PP!V!pRRR4\ XVVR,VR ,R 7R# +'giL/;i) r r r r]rGrzr{Nr r r)rTrr4r rr)rYr r r mkvrs&& r.test_moment_against_mp+TestStudentizedRange.test_moment_against_mpd sz* , smXc]HSM9[[ * *))00#6C+ Y%o6%o6 8+ *s B B" c^^ r!\\PP^\P W3R7p\ V^,^4R#)r!rN)rr4r rfrTr+r)rYr]rGrs& r.test_pdf_integration)TestStudentizedRange.test_pdf_integrationu s9"15**..266GA"r0c^^ r!\P!^^ RR7p\PP W1V4R,p\PP W1V4p\ WS4p\WdRR7R#)r!r)stepr,rrN)rTr'r4r rnrfrr)rYr]rGr^y_cdf y_pdf_rawy_pdf_cumulatives& r.test_pdf_against_cdf)TestStudentizedRange.test_pdf_against_cdf| sj"1 IIa$ '''++A!4R8++//a8 / = (d;r0 r_case_resultcVwr#rE\P!RR7;_uu_4\PP W#V4pRRR4\ XV4R# +'giL;i)rzr{N)rTrr4r rnr)rYr rr]rGr_resrs&& r.test_cdf_against_r'TestStudentizedRange.test_cdf_against_r sL'a [[ * *))--aA6C+U#+ *s !A A, c\P!RR7;_uu_4\PP ^^.^^.^ ^ .4pRRR4\ XP R4\P!\RR7;_uu_4\PP ^^^..RO4RRR4R# +'giLz;i +'giR#;i)rzr{Nz...could not be broadcast...rr)rbrr) rTrr4r r5rr?rr r)rAs& r.test_moment_vectorization.TestStudentizedRange.test_moment_vectorization s[[ * *''--q!fq!fr2hGA+ & ]]:-K L L  # # ) )!aV\ BM L + * M L Ls(B: &C : C C c \P!4;_uu_4\P!RR7;_uu_4\P!R\ 4\ PP.RO4wrpRRR4RRR4\\ PPXX44R# +'giLC;i +'giLN;i)rzr{Nr) r{r|rTrr}rr4r rr _argcheck)rYr]rrs& r.test_fitstart_valid(TestStudentizedRange.test_fitstart_valid s  $ $ & & H(E(E  ! !(,> ?11;;IFKA1a)F & ''11!R89 )F(E & &s# C AB: ?C : C C  C c~\PP^^ \P4p\PP^^ R4p\ WRRR7\PP ^^ \P4p\PP ^^ R4p\ WRRR7R#)r!韆rrN)r4r rfrTr+rrn)rYr res_finites& r.test_infinite_df%TestStudentizedRange.test_infinite_df s%%))!R8,,00B> d>%%))!R8,,00B> d>r0c *\PP^^ R4p\PP^^ R4p\PP^^ R4p\\\ WRRR7\ W#RRR7\PP ^^ R4p\PP ^^ R4p\PP ^^ R4p\\\ WRRR7\ W#RRR7R#)r!rr irrN)r4r rfrAssertionErrorrrn)rYrr res_sanitys& r.test_df_cutoff#TestStudentizedRange.test_df_cutoff s%%))!R8,,00B> ,,00B> nosd , TE%%))!R8,,00B> ,,00B> nosd , TEr0c~R^Rr2p\PPWV4p\V^RR7V^8gQhR#)gaRA@iSrrN)r4r rr)rYrr]rGrzs& r. test_clipping"TestStudentizedRange.test_clipping s>&q#a  # # & &qQ /15)Av vr0r6)gQ1@g33333F@g ףp=*K@gHzM@g/$@g~j!@g33333$@g{Gz&@gh|?5 @gQ8@g&1@g|?5@g@gOn@g!rhm@gB`"@rYgʡEs@g|?5^@g/$@gV-@g/$@gF@g?5^I @)gRV@gffffffl@gp@gr@gʡ @gHzG/@gQ82@gQ3@g!rh@g@g +@gʡEs @gjt@gB`"[@g@gnJ@gV- @gFx@gsh|?@g+N@g/$ @g @g@gGz@)gfffff"@ii i gHzG2@g(\A@gQC@gfffffE@g @gNbX9"@gHz$@g(\&@g-@gn@@gw/@g= ףp @gJ +@g-'@gn@gHz@gV-@gnJ@g$@g/ݤ@)rrr)r rrrQ)rr!i)#gupm1f?)rrbrg?<-=A?)4rrrrq05q01q001rTrqspsr+vskslistr&r#rr_datarrrr rr ospathdirnamer path_prefix relative_pathopenr filejsonrpregenerated_datarr r xfail_on_32bitr r r r r r r r r rrr s@r.r r  s 'C  'C  (D c4( )B B QBRVV $B B GBB', -D *( Arvv() Arvv() B)* F [[::  [[:: ''//(+Kg}XEQ=gѭEG,,RVVS9:;..rvvs;<=++BFFC89:,,RVVS9:;//<=>,,RVVS9:;,,RVVS9:;# ' & & &s R1S S, c \PPR4p\P!.RO4p\P P V^,V^,V^,^dVR7p\P!RR7;_uu_4\P!\P PVRR74pRRR4\XV^R 7\P!\P PVR^R 74p\W1^R 7\P!\P PVR^R 74p\W1^R 7\P!\P PVR^R 74p\W1^R 7\PPR4p^pR p\PP VR^dVR7p\P!\PPW%R74p\P!V\P!W%, ^,P44.4p\W6^R 7R# +'giEL;i)i.rorrzrrz)r rNrU)rrc)rrcr rr)rorrz) rTrUrVr r4rrrrrr8rr)rWtruer^rrcrrNs r.test_frozen_fit_ticket_1536r8 !s ))   %C 88O $D $q'47DG#CPA H % %%--++AB+78 &a0 XXemm''#1'= >Fa0 XXemm''d': ;Fa0 XXemm''d'; H!4) & % %s 6I55 J c\PPR4p\PP ^dVR7p\PP V4pRp\W#^R7R#)i rrUN)g ףp= ?gʡE?)rTrUrVr4rrrr)rWrrrNs r.test_regression_ticket_1530r: !sR ))   'C ,,  #  6C \\  c "FH!4r0c\PPR4pVP^*4pRFEp\P P W,4wr4\W2RR7\VRRR7KG R#)rrSrr`N)rg6A)rTrUrErandnr4rrr)rWr^offsetrcrds r.test_gh_pr_4806r> !s[ ))   %C " A'\\%%aj1 #.s-(r0c\\P!\PP R^444R#)riN)rrTrr4rr[r6r0r.test_poisson_logpmf_ticket_1436r@ !s$ BKK ,,T37 89r0c|RR.pVF1wr\PPVRR7p\W24K3 R#)aTest the powerlaw stats function. This unit test is also a regression test for ticket 1548. The exact values are: mean: mu = a / (a + 1) variance: sigma**2 = a / ((a + 2) * (a + 1) ** 2) skewness: One formula (see https://en.wikipedia.org/wiki/Skewness) is gamma_1 = (E[X**3] - 3*mu*E[X**2] + 2*mu**3) / sigma**3 A short calculation shows that E[X**k] is a / (a + k), so gamma_1 can be implemented as n = a/(a+3) - 3*(a/(a+1))*a/(a+2) + 2*(a/(a+1))**3 d = sqrt(a/((a+2)*(a+1)**2)) ** 3 gamma_1 = n/d Either by simplifying, or by a direct calculation of mu_3 / sigma**3, one gets the more concise formula: gamma_1 = -2.0 * ((a - 1) / (a + 3)) * sqrt((a + 2) / a) kurtosis: (See https://en.wikipedia.org/wiki/Kurtosis) The excess kurtosis is gamma_2 = mu_4 / sigma**4 - 3 A bit of calculus and algebra (sympy helps) shows that mu_4 = 3*a*(3*a**2 - a + 2) / ((a+1)**4 * (a+2) * (a+3) * (a+4)) so gamma_2 = 3*(3*a**2 - a + 2) * (a+2) / (a*(a+3)*(a+4)) - 3 which can be rearranged to gamma_2 = 6 * (a**3 - a**2 - 6*a + 2) / (a*(a+3)*(a+4)) rrN)rS)rzrrr)r)gUUUUUU?r g^cQg333333)r4rr)rr, exact_mvskrs r.test_powerlaw_statsrC !s>>, = ?E ~~##Av#6!$3r0c^\PP^^4p\VR4R#)r9rN)r4rrrrzs r.test_powerlaw_edgerF !s" a#ACr0cx\PP^^4p\VR4\PP ^.RO4p\V\ P RR.4\PP^.RO4p\V\ P R\ P ).4R#)r9rrSN)rorSr)r4rrrrfrTr+rE s r.test_exponpow_edgerH !s a#AC 1./ARVVS#&' a!12ARVVS266'*+r0ca]tRtRtoRt]P PR.R O4R4tRt Rt Rt Rt Vt R #) TestGenGammai!c`\PP^^^4p\VR4R#)r9rSN)r4gengammarfrrs& r.test_gengamma_edgeTestGenGamma.test_gengamma_edge!s$ NN  q!Q 'Qr0za, c, ref, tolcZ\\PPW4W4R7R#r)rr4rL rE)rYr,rrrs&&&&&r.test_gengamma_extreme_entropy*TestGenGamma.test_gengamma_extreme_entropy!s" ..q4cDr0c\PP^^R4pVR8XgQh\PP^^R4pV\P )8XgQhR#)r9rNr)r4rL rfrrTr+)rYrzrs& r.!test_gengamma_endpoint_with_neg_c.TestGenGamma.test_gengamma_endpoint_with_neg_c "sO NN  q!R (Cxx~~$$Q2.wr0c\PPR^R4p\VR4\PPR^ R4p\VR4R#)r rSNr g>gqq?)r4rL r5rrs& r.test_gengamma_munpTestGenGamma.test_gengamma_munp"sF NN S" -Az* NN R ,Av&r0c\\PP.RO^R4\P )RR.4R#)r9N)r9rrr)rr4rL rrTr+rXs&r.)test_gengamma_logpdf_broadcasting_gh245746TestGenGamma.test_gengamma_logpdf_broadcasting_gh24574"s,--iB?266'2rARSr0r6N))g`6Arg2j!@r)r rggA@r)rrrN r)rOrg#@r)g`FArbgI~mr)gl$W}Lr|g.I#Ur)rrrrrM rrrrP rS rV rY rrr s@r.rJ rJ !sP  [[-NO EO E 'TTr0rJ c \P!.R!NRNRNRNR"NR#NRNRNRNRNRNR$NR NR%NR&NR NR NR NR NRNRNRNRNRNRNRNRNRNR'NR(NR)NR*NRNR+NRNR,NR-NR.NRNR/NR0NR1NRNR2NR3NR4NRNR5N4p\P!RR7;_uu_4\P!4;_uu_4\P !RR\ 4\P !RR\4\PPV4R R R 4R R R 4R # +'giL;i +'giR #;i)6g>X*?gPix#?g=?g&^?gO_ZQ?g\?g~?g)/?gRܶ?gH߃[g?g7k?g#BF?g KU= ?gT?gZg[Qs?g2?g6Hwƭ8?g*- ?ga]?g'?g޻?g?J_?g~sw?g&zM?g F?gTzy?g 't?g,?rzr{z:The maximum number of subdivisions .50. has been achieved.z-floating point number truncated to an integerNg>X*ȿgoϿg뵏˿g핚οg S"gm_ пg0FfͿgg"@ӿg?%:˿g?Gg~,g HgW2ҿgwm ǿg@gW=ֿg".v1glg.Hꮿg*˿) rTr rr{r|rrr~r4r"r)r\ s r.test_ksone_fit_freezer\ "sZ  / /j /* /j /) /  / /!+ /-7 /9B /  /  /", /.9 /;F /  / /!+ /-7 /9C /  / /!+ /-7 /9B /  /  / ", / .9 / ;F /  / " / $- / /: / B""1biiS&91= #  BKK  # # %&  " "g~6UZZ^^E1b115 ))%B79KL # " # " # " "sBE>A2F> F  F" zmethod, expectedrng u%>g ;Q=rfggRz\R>rrgfHu@ge:@cj\\PV4!R^^.^ R7p\W!RR7R#)r)rw rrrN)r2r4rr r)rrNrs&& r.test_ncx2_zero_ncry "s+ UZZ (!QB ?FF51r0c\PP^ ^^R7p\PP^ ^R7p\ WRR7R#)rb)rrw r)rrrrN)r4rr rrl r)rrNs r.test_ncx2_zero_nc_rvsr{ "s<ZZ^^raa^ 8Fzz~~!~4HF51r0c^ \P!^^ 4,p\\PP R^VR7^4R#)rbrrrw N)rTr'rr4rr rn)rw s r.test_ncx2_gh12731r~ "s4 RYYq"  Bb115r0c\P!.RO4p^Rr!\PP WVR7p.ROp\ W4RR7R#)ю] @gh?8@r} rrN) r gj%@gb y/7@gH@gIZ@g Xl@g'/2?8@gݰmА@gU @gy]`'@g@g*\@g48c@) rSrSrSrSrSg?gy[qD?rrrrrr)rTr r4rr rr)r^nurWr sf_expecteds r.test_ncx2_gh8665r "sH " #A $ qC (B'K B%0r0c RpRp\P!\PP RW4\PP RW4RR7p\PP W V4p\P P W V,\P!^V,^V,,44p\W4RR7R#) i,ir2rrrn rrN) rTrWr4rr rrfr8rr)rrw r^ncx2_pdf gauss_approxs r.test_ncx2_gh11777r "s B B EJJNN5"1JJNN5"1u >Azz~~aR(H::>>!"Wbgga"fq2vo.FGLH6r0zx, c, expectedc\\PPW4p\W2R4R#r)N)r4 foldcauchyrrr^rrNrs&&& r.test_foldcauchy_sfr #s$     Q "BB%(r0rkc\\PPV4p\W!R4R#r )r4rrr)r^rNrs&& r.test_halfcauchy_sfr #s$     Q BB%(r0z p, expectedcZ\PPV4p\W!4R#rt)r4rrr)rzrNr^s&& r.test_halfcauchy_isfr *#s" QAA r0cl\P!^^R7p\VP^4^4R#)r9r^N)r4rrrn)rs r.test_foldnorm_zeror 6#s$  #BAr0c\\PPW4p\W2R4R#)r N)r4rrrr s&&& r.test_foldnorm_sfr J#s"   1 BB%(r0c\PP.RO^R4p\PPRR.^R4p\;QJd.RV4FNK 5M !RV44p\W 4\PP.RO4p\PP^R.4p\;QJd.RV4FNK 5M !RV44p\W 4R#) rrzrSc3p"TF,p\P\PV3,xK. R#5irtrTrLr*rF s& r.r-test_stats_shapes_argcheck..Y#s$8Cq"%% **C46rc3p"TF,p\PV\P3,xK. R#5irtr rF s& r.rr _#s$8Cq"%%266 **Cr N)rrzrS)r rr)r4rtuplerr)mv3mv2 mv2_augmenteds r.test_stats_shapes_argcheckr T#s ..  3 7C ..  Sz1c 2CE8C8EE8C88M$ --  l +C --  q#h 'CE8C8EE8C88M$r0c&a]tRtRtoRtRtVtR#) _distr_genii#c^*#*r6rYr^r,s&&&r._pdf_distr_gen._pdfj#sr0r6Nrrrrr rrr s@r.r r i#sr0r c&a]tRtRtoRtRtVtR#) _distr2_genin#c"^*V,V,#r r6r s&&&r._cdf_distr2_gen._cdfo#Avzr0r6N)rrrrr rrr s@r.r r n#sr0r c,a]tRtRtoRtRtRtVtR#) _distr3_genis#cW#,#rtr6rYr^r,r-s&&&&r.r _distr3_gen._pdft#s u r0c"^*V,V,#r r6r s&&&r.r _distr3_gen._cdfw#sAvzr0r6Nrrrrr r rrr s@r.r r s#sr0r c,a]tRtRtoRtRtRtVtR#) _distr6_geni}#c W!,V,#rtr6r s&&&&r.r _distr6_gen._pdf#s sQwr0c"^*V,V,#r r6r s&&&&r.r _distr6_gen._cdf#r r0r6Nr r s@r.r r }#sr0r ca]tRtRtoRtRtRtRtRtRt Rt R t R t R t R tR tRtRtRtRtRtRtVtR#)TestSubclassingExplicitShapesi#cZ\RRR7p\VP^^R7^*4R#)dummyr,rr"rZNr rrfrY dummy_distrs& r.test_correct_shapes1TestSubclassingExplicitShapes.test_correct_shapes#s% gc: [__Q!_,b1r0cl\RRR7p\\VP^3/\ ^R7BR#)r Ar rZN)r rr,rfr%r s& r.test_wrong_shapes_11TestSubclassingExplicitShapes.test_wrong_shapes_1#s' gc: i!AtayAr0ct\RRR7p\^^^R7p\\VP^3/VBR#)r za, b, cr rr N)r r%rr,rf)rYr r3s& r.test_wrong_shapes_21TestSubclassingExplicitShapes.test_wrong_shapes_2#s1 gi@ Q!q!i!;s;r0cJ\R^*R7p\\\3/VBR#)r r N)r%rr,r rYr3s& r.test_shapes_string0TestSubclassingExplicitShapes.test_shapes_string#s+i3s3r0cJ\RRR7p\\\3/VBR#)r z(!)r Nr%r SyntaxErrorr r s& r.test_shapes_identifiers_17TestSubclassingExplicitShapes.test_shapes_identifiers_1#s.k:55r0cJ\RRR7p\\\3/VBR#)r 4chanr Nr r s& r.test_shapes_identifiers_27TestSubclassingExplicitShapes.test_shapes_identifiers_2#0k:55r0cJ\RRR7p\\\3/VBR#)r zm(fti)r Nr r s& r.test_shapes_identifiers_37TestSubclassingExplicitShapes.test_shapes_identifiers_3#s1k:55r0cJ\RRR7p\\\3/VBR#)r za=2r Nr r s& r."test_shapes_identifiers_nodefaults@TestSubclassingExplicitShapes.test_shapes_identifiers_nodefaults#s.k:55r0cJ\RRR7p\\\3/VBR#)r z*argsr Nr r s& r.test_shapes_args.TestSubclassingExplicitShapes.test_shapes_args#r r0cJ\RRR7p\\\3/VBR#)r z**kwargsr Nr r s& r.test_shapes_kwargs0TestSubclassingExplicitShapes.test_shapes_kwargs#s 3k:55r0cJ\RRR7p\\\3/VBR#)r za, b, c, lambdar Nr r s& r.test_shapes_keywords2TestSubclassingExplicitShapes.test_shapes_keywords#s(9:k:55r0c!RR\P4pV!RR7p\VPR^R7\PPR4^,4R#)c&a]tRtRtoRtRtVtR#)FTestSubclassingExplicitShapes.test_shapes_signature.._dist_geni#cN\PPV4V,#rtr4r8r r s&&&r.r KTestSubclassingExplicitShapes.test_shapes_signature.._dist_gen._pdf#zzq)A--r0r6Nr r s@r. _dist_genr # . .r0r r,r"rzrZNr4rCrrfr8rYr rJs& r.test_shapes_signature3TestSubclassingExplicitShapes.test_shapes_signature#sI .++ .$TXXcQX')?r0c !RR\P4pV!RR7p\\VPR3/\ ^^R7BR#)c&a]tRtRtoRtRtVtR#)STestSubclassingExplicitShapes.test_shapes_signature_inconsistent.._dist_geni#cN\PPV4V,#rtr r s&&&r.r XTestSubclassingExplicitShapes.test_shapes_signature_inconsistent.._dist_gen._pdf#r r0r6Nr r s@r.r r #r r0r rr rzrN)r4rCrr,rfr%r s& r."test_shapes_signature_inconsistent@TestSubclassingExplicitShapes.test_shapes_signature_inconsistent#s; .++ .'i3A$a.Ar0c!RR\P4pV!RR7p\VPR^!R7\PPR4^!,4\VPR^!4\PPR4^!,4\ \ VPR3/\^!R7BR#)c&a]tRtRtoRtRtVtR#)?TestSubclassingExplicitShapes.test_star_args.._dist_geni#c`V^,p\PPV4V,#r9r )rYr^r4 extra_kwargs&&* r.r DTestSubclassingExplicitShapes.test_star_args.._dist_gen._pdf#s#"1g zzq)K77r0r6Nr r s@r.r r#s 8 8r0r rr rzr)xxxN)r4rCrrfr8rr,r%r s& r.test_star_args,TestSubclassingExplicitShapes.test_star_args#s 8++ 8  .TXXcrX2EJJNN34G4JKTXXc2& s(;B(>?i3?$2,?r0cl!RR\P4pV!RR7p\VPR^o^!R7\PPR4^!,^o,4\VPR^o^!4\PPR4^!,^o,4R#)c&a]tRtRtoRtRtVtR#)ATestSubclassingExplicitShapes.test_star_args_2.._dist_geni#cnV^,p\PPV4V,V,#rr )rYr^r= r4rs&&&* r.r FTestSubclassingExplicitShapes.test_star_args_2.._dist_gen._pdf#s("1g zzq)K7&@@r0r6Nr r s@r.r r #s A Ar0r zoffset, extra_kwargr rz)r= rNr r s& r.test_star_args_2.TestSubclassingExplicitShapes.test_star_args_2#s A++ A  56TXXc#2X>ZZ^^C(+c1 3TXXc3+ZZ^^C(+c1 3r0c!RR\P4pV!RR7p\VP^^R7\PP^44R#)c&a]tRtRtoRtRtVtR#)BTestSubclassingExplicitShapes.test_extra_kwarg.._distr_geni#crVPR^4p\PPV4V,#r)popr4r8r )rYr^r4kwargsrs&&*, r.r GTestSubclassingExplicitShapes.test_extra_kwarg.._distr_gen._pdf#s,%jj: zzq)K77r0r6Nr r s@r.r r#s 8 8r0r rr rNr )rYr rJs& r.test_extra_kwarg.TestSubclassingExplicitShapes.test_extra_kwarg#sD 8,, 8/TXXaQX/1BCr0c!RR\P4pV!RR7p\VPR4\PPR44R#)c&a]tRtRtoRtRtVtR#)ITestSubclassingExplicitShapes.test_shapes_empty_string.._dist_geni#c@\PPV4#rt)r4r8rfrs&&r.r NTestSubclassingExplicitShapes.test_shapes_empty_string.._dist_gen._pdf#szz~~a((r0r6Nr r s@r.r r#s ) )r0r r rzNr r s& r.test_shapes_empty_string6TestSubclassingExplicitShapes.test_shapes_empty_string#s? )++ )#TXXc]EJJNN3$78r0r6N)rrrrr r r r r r r r r r r r r r rrr!rrr s@r.r r #sd2B< 4 6 666666 @B @ 3 D99r0r ca]tRtRtoRtRt]PP] RR7R4t ]PP] RR7R4t Rt R t R tR tR tVtR #)TestSubclassingNoShapesi#cX\RR7p\VP^^R7^*4R#r rrZNr r s& r.test_only__pdf&TestSubclassingNoShapes.test_only__pdf$s# g. [__Q!_,b1r0cX\RR7p\VP^^R7^4R#r&)r rrfr s& r.test_only__cdf&TestSubclassingNoShapes.test_only__cdf$s#!w/ KOOAO3Q7r0docstring strippedrc\RR7p\VP^4\VPR4\P !RVP 4p\\V4^8H4R#)r r'r,zlogpdf\(x, a, loc=0, scale=1\)N) r rnumargsr"refindallrQrr#rYr rs& r.test_signature_inspection1TestSubclassingNoShapes.test_signature_inspection $sX!g. [((!,[''-jj:$,,.CA r0c\RR7p\VP^4\VPR4\P !RVP 4p\\V4^8H4R#)r r'rz!logpdf\(x, a, b, loc=0, scale=1\)N) r rr/r"r0r1rQrr#r2s& r.test_signature_inspection_2args7TestSubclassingNoShapes.test_signature_inspection_2args$sX"w/ [((!,[''0jj=$,,.CA r0c2\\\RR7R#)r r'N)rr,r rXs&r.0test_signature_inspection_2args_incorrect_shapesHTestSubclassingNoShapes.test_signature_inspection_2args_incorrect_shapes$si7;r0cn!RR\P4p\\V3/\ RR7BR#)c*a]tRtRtoRRltRtVtR#)>TestSubclassingNoShapes.test_defaults_raise.._dist_geni%$c^*#r r6r s&&&r.r CTestSubclassingNoShapes.test_defaults_raise.._dist_gen._pdf&$r0r6Nr r r s@r.r r=%$s  r0r r r'Nr4rCrr,r%rYr s& r.test_defaults_raise+TestSubclassingNoShapes.test_defaults_raise#$* ++  iAd.@Ar0cn!RR\P4p\\V3/\ RR7BR#)c&a]tRtRtoRtRtVtR#)>TestSubclassingNoShapes.test_starargs_raise.._dist_geni,$c^*#r r6)rYr^r,r4s&&&*r.r CTestSubclassingNoShapes.test_starargs_raise.._dist_gen._pdf-$r@r0r6Nr r s@r.r rH,$  r0r r r'NrArBs& r.test_starargs_raise+TestSubclassingNoShapes.test_starargs_raise*$rEr0cn!RR\P4p\\V3/\ RR7BR#)c&a]tRtRtoRtRtVtR#)._dist_geni3$c^*#r r6)rYr^r,rs&&&,r.r ATestSubclassingNoShapes.test_kwargs_raise.._dist_gen._pdf4$r@r0r6Nr r s@r.r rP3$rKr0r r r'NrArBs& r.test_kwargs_raise)TestSubclassingNoShapes.test_kwargs_raise1$rEr0r6N)rrrrr(r+rrrDOCSTRINGS_STRIPPEDr3r6r9rCrLrSrrr s@r.r$r$#s28  [[+4HIJ [[+4HIJ<BBBBr0r$r-c2.ROp\PF~p\\V4p\V\P\P ,4'gKJVF.p\ \P!W2P4RJ4K0 K R#),\s*,N)rWz\(\s*,z^\s*:) r4rGr2rArBrCrr0searchrQ)badonesdistnamerJregexs r.rl rl 9$sc-GMMuh' dU..1D1DD F F  %6$>?!"r0c\\PP\P ^ ^ 4^4\\P P\P ^R4^4R#)rbrN)rr4rrrTr+rr r r6r0r.test_infinite_inputr]C$sB ((R8!< 37;r0c\PP\PPR^4^4p\ VR4R#r{)r4lomaxrrnrrE s r.test_lomax_accuracyr`H$s/  2A6AAvr0c\PP\PPR^4^4p\ VR4R#r{)r4rrrnrrE s r.test_truncexpon_accuracyrbN$s7 U--11&!>,,Q2A6A3+r0c\P!RR7;_uu_4p\P!R4\PP R^4\PP R^4\PPR^4\PP\P)R4\V4p\V^4RRR4R# +'giR#;i)zregression test for gh-6219T)recordalwaysrzrN) r{r|r}r4rrnrfrrrTr+r#r)rnumber_of_warnings_throwns r. test_genextreme_give_no_warningsriZ$s   - -h' R# R# R# -$'F!.2 . - - -s B=C'' C8 c~\P!\PPR\P^R.RR74p\P p\P ^,^, p^ \P!^4,\P!^4,\P ^,, pRpWW4.p\P.^,;rg.ROp\VR ,VRR7\VR ,V4\VR ,VRR7\VR ,V4R#) rrrr rNrr)r9rr r)rr9)rr)rr )rr!) rTrr4rr* euler_gammargrrzetarr) rrrrrref_0ref_1ref_3ref_2s r.test_moments_gh22400rqh$s **U%%++S"&&!T,BF+S TC >>D %%(Q,C  ?W\\!_ ,ruuax 7D D  #EVVHQJE ECIu51TE"CIu51TE"r0cRp\PPR4p\V^V,^,RR7\PP^4p\W^,RR7\PPR4p\ V^4\PPR^ R7p\W^,\ P !^ 4,^,RR7\PP^ 4p\VRV,^,RR7\PPR 4p\V^ V,^,RR7R#) gox?rSr rr^Nrrrgr)r4rrErrrTr)rkrss r.test_genextreme_entropyrs~$s$K   &AAq}q(u5   #AAQU3   %AA   R 0AA1}rvvbz1A5EB   $AAr+~)6   %AAr+~)6r0cRp\PPVR4p\VR4\PP VR4p\W 4Rp\PPVR4p\VR4\PP VR4p\W 4Rp\PPV^4p\VR4\PP V^4p\W 4R#)rrg'b%4gQ@gؗҜ;g.l6?Nr)r4rrrr)r^r?rs r.test_genextreme_sf_isfru$s2 A Av&AA-.     a (BB A Au%AA()     a 'BB A Aq!AA-.     a #BBr0cdRp\PPV^^4p\VR4R#)rgH&8>N)r4r rr)probrs r.test_burr12_ppf_small_argrx$s, D||a+HH45r0cRRlp\P!.R O4p\PP V^VR7wr#p\ VRRR7V^8XgQh\ VRRR7R#) a Test fitting invweibull to data. Here is a the same calculation in R: > library(evd) > library(fitdistrplus) > x = c(1, 1.25, 2, 2.5, 2.8, 3, 3.8, 4, 5, 8, 10, 12, 64, 99) > result = fitdist(x, 'frechet', control=list(reltol=1e-13), + fix.arg=list(loc=0), start=list(shape=2, scale=3)) > result Fitting of the distribution ' frechet ' by maximum likelihood Parameters: estimate Std. Error shape 1.048482 0.2261815 scale 3.099456 0.8292887 Fixed parameters: value loc 0 c  \WW#RRR7#)r)r4dispxtolftol)r!)rx0r4r{s&&&&r.r&test_invweibull_fit..optimizer$sD4UKKr0)rrg`?rrg [@N)r6r9)rrPr rrYr!rrrrrbr@c)rTr r4 invweibullrr)rr^rrcrds r.test_invweibull_fitr$s].L HIA$$((i(HMAEAxd+ !8O8E8$/r0c^\PPW4p\W2RR7R#r~)r4rrr)r^rrNrs&&& r.test_invweibull_sfr$s$ ""1(HHU3r0zp, c, expectedc^\PPW4p\W2RR7R#r~)r4rrr)rzrrNrs&&& r.test_invweibull_isfr$s$##A)HHU3r0z df1,df2,xc>^p\PPW V4p\PPW W4p\ WERR7\PP W V4p\PP W W4p\ WgRR7R#)r9r rrN)r4r$rnr rrf)df1df2r^rw expected_cdfcalculated_cdfrcalculated_pdfs&&& r.test_ncf_edge_caser%so B77;;qs+LYY]]133NLu=77;;qs+LYY]]133NLt@g&4I,)?)decimalsN)r4r rnrrTround) scipy_val check_vals r.test_ncf_cdf_spotcheckr%%s6  b!R.IIIrxx A>?r0c \P!^^R4pRV^&Rp\PP!V.VO5!pVUu.uF%p\PP!V.VO5!NK' pp\ W$4R#uupi)r9rrN)rr rr9r)rTrWr4r rr)r^parrxiq0s r.test_ncf_ppf_issue_17026r/%sl Aq#A AaD C a#A,- .Ab%))-- !S !AB .A /s+BcDa]tRtRtoRtRtRtRtRtRt Rt Vt R #) TestHistogrami9%cx\PPR4Vn\P!.RO^R7p\ P !V4Vn\ PPRRRVPR7p\P!V^2R7p\ P !V4Vn R#)l&'binsrSrrrN)rr r r!r!r!rrrrrrrrrrrrrrrrrrrZ) rTrUrVrW histogramr4rDtemplater8r norm_template)rYrrnorm_histograms& r.rZTestHistogram.setup_method:%s99((4LL"@FGI **95 zz~~#Su488~Td4"//?r0c \P!.R O4p\P!.R O4p\VPP V4V4\ VPP R4R 4\ VPP R4R 4\ VPP R4R4\ VPP R4R4\P!R^^ 4p\VPP V4\PP VRRR7RR 7R #)rrSrr!@r4rrbrrNrrzrSrrrr rr3rrT@rVr-r. rYrrr4rQ?)rrrKrK{Gz?rrr{Gz?rrwrwrrrrrrrrr ) rTr rrrrfrrWrr4r8)rYrv pdf_valuesr^s& r.rPTestHistogram.test_pdfJ%sMNZZ!IJ  ))&1:> DMM--c2H=DMM--c2H= DMM--c2H=DMM--d3X> KKAr "**..q1 qc=C Ir0c \P!.RO4p\P!.R O4p\VPP V4V4\VPP V^R 4V^R 4\P!RR^d4p\VPP VPP V44V4\P!RR^d4p\VPP VPP V44V4\P!R ^^ 4p\VPP V4\PP VRRR7RR7R#) rrSrr4rbrrNr)rrrrrKrr ףp= ?Q?{Gz?rrzr`(\?RQ?= ףp=?)\(?rk rSrSrr ) rTr rrrrnrrWrr4r8)rYrv cdf_valuesr^s& r.rTestHistogram.test_cdf_ppfd%s)MNZZ!MN  ))&1:> ))*Qr*:;VAb\J KKS# & ))$--*;*;A*>?C KKS# & ))$--*;*;A*>?C KKAr "**..q1 qc=C Ir0cbRpVPPWPR7p\\P !VR84R4\ \P !VR8*4RV,RR7\ \P !VR8*4RV,RR7\ \P !VR8*4RV,R R7\ \P !VR 8*4RV,R R7\ \P !VR 8*4RV,R R7\ \P !VR 8*4RV,R R7\ \P !VR8*4RV,RR7\ \P !VR8*4RV,RR7\ \P !VR 8*4R V,RR7\ \P !VR8*4R!V,RR7\ \P !VR8*4R"V,RR7\ \P !VR8*4R#V,RR7\ \P !VR8*4R$V,RR7\ \P !VR8*4R%V,RR7\ \P !VR8*4RV,RR7\ \P !VR8*4RV,RR7\\P !VR84R4R#)&rrrSrrrwrrr rrrr3rVrrTrjrr-r. rYrr4NrKrrrrrrrzr`rrrrrk )rrrWrrTrpr)rYrLrus& r.rGTestHistogram.test_rvsz%s  """ARVVFSL)3/v}-x!|#Fv}-x!|#Fv}-x!|#Fv}-x!|#Fv}-x!|#Fv}-x!|#Fv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4Hv}-y1}4HRVVFSL)3/r0c\^4FNp\VPPV4\P !RR4P V4RR7KP R#)rrSrrjrN)rrrr5r4r8rrs& r.r- TestHistogram.test_munp%sGqA D..44Q7!JJsC077: Gr0c\VPP4\PPRRR7RR7R#)rSrrbrjrN)rrrEr4r8rXs&r.rtTestHistogram.test_entropy%s5**224 **s#*>T Kr0)rrWrN) rrrrrZrPrrGr- rtrrr s@r.rr9%s.@ I4I,0,G KKr0rc^^..R Or\P!W3RR7p\PP VP R^.4RR.4VP 4^8XgQh\P!W3RR7p\PP VP R^.4R 4VP 4R 8XgQhRp\P!\VR7;_uu_4\P!W34pVP 4R 8XgQhRRR4\P!V.R O34pVP 4^8XgQhR# +'giLC;i) rF)r+ rzrTz(Bin widths are not constant. Assuming...rN)r9rrg'^P?gH@r) r4rDrTrIrrfrrrM r~)countsrrJr s r.test_histogram_non_uniformr%sFLD   vne \PPR4p\PP R^RVR7p\PP W1R7wrErgV^8XgQh\PP V^VR7wrErgV^8XgQhR#)rrrrrrN)rTrUrVr4rrr)rYrrWrr,r-rcrds&& r.test_fit_override TestLogUniform.test_fit_override%sii##K0""33"G ++///Cczz ++//Af/Mczzr0c\PPR4pRRr2\P!W#4pVP ^^RR7p\ VPVPV44V4VPRVR7p\ VPVPV44V4R\P!R ^4,pVPV4p\ VRR VR,, ^ 4W2, \P!V4\P!V4, , p \ VP4V 4R#) leVi PgN~hrrrTrrNr,i8r)rTrUrVr4rrrrnrrr'rfrr) rYrWr,r-rJrnrr^rfrs & r. test_overflowTestLogUniform.test_overflow%sii##$78u1%kk!QTk*#/5hhDsh3#/5 4% %hhqkCRR("-q BFF1I-. T*r0r6N) rrrrrrrrrrrrr s@r.rr%s=4" [[Xt}5 6 ++r0rca]tRtRtoRt]P PRRR.R^.R^..4R4t]P PR R R .4R 4t ]P PR .RO4R4t ]P PR.RO4R4t ]P PR.RO4R4t ]P PR.RO4R4t ]P PR.R O4R4t]P PR.R!O4R4t]P PR.R"O4R4t]P PR.R#O4R4tRtVtR#)$ TestArgusi%c\PP^2RRR7p\\P!^24P 4VP 4^R7R#)r=rErrUN)r4argusrrrrs& r.test_argus_rvs_large_chi"TestArgus.test_argus_rvs_large_chi%s= KKOOBSsO ;EKKO002AFFHaHr0zchi, random_staterrrg rc\PPVRVR7p\P!VRV34wrE\ VR84R#)rrrrjNr4rrrmr)rYrrr^rrzs&&& r.rGTestArgus.test_rvs%s> KKOOCc O E||Aw0Dr0rrrc\PPVRRR7p\P!VR4wr4\ VR84R#)rie rc>^^V^,, R,, #)rrr6rs&r.r@.TestArgus.test_rvs_small_chi..%sa!Q$h#->)>r0rjNr)rYrrrrzs&& r.test_rvs_small_chiTestArgus.test_rvs_small_chi%s; KKOOCcO ?||A>?Dr0zchi, expected_meancb\PPV^R7p\W2RR7R#)rr^rrN)r4rrr)rYr expected_meanr s&&& r.rTestArgus.test_mean&s' KK  S  *u5r0zchi, expected_var, rtolcb\PPV^R7p\WBVR7R#)rr^rN)r4rrr)rYr expected_varrrGs&&&& r.rTestArgus.test_var &s$ KKOOCqO )d3r0zchi, expected, rtolc2\\V4W#R7R#r)rr)rYrrNrs&&&&r.test_argus_phi_small_chi"TestArgus.test_argus_phi_small_chi&s  3=r0z chi, expectedc\P!.RO4p\\PP W14VRR7R#)rrrNrrzr)rTr rr4rrfrYrrNr^s&&& r.test_pdf_small_chiTestArgus.test_pdf_small_chi$&+ HH_ % /Fr0c\P!.RO4p\\PP W14VRR7R#)rr rNr)rTr rr4rrrs&&& r.test_sf_small_chiTestArgus.test_sf_small_chi3&s+ HH_ % q.uEr0zx, chi, expectedc^\PPW4p\WCRR7R#r)r4rrr)rYr^rrNrs&&&& r.test_sf_near_1TestArgus.test_sf_near_1B&s [[^^A #51r0c\P!.RO4p\\PP W14VRR7R#)rrrNr)rTr rr4rrnrs&&& r.test_cdf_small_chiTestArgus.test_cdf_small_chiM&rr0cb\PPVRR7p\WBVR7R#)rJrrN)r4rr)rYrrNrrFs&&&& r.test_stats_small_chiTestArgus.test_stats_small_chi\&s'kkT2D1r0r6N))rgщi?)rbg.憃?)rg\paP?)rg ?)rg,6?))rgLH'B?r)rbgo$?r)rg =>r)rg2W>r)rgb~P>r))rgwM?r )rzgj?r )rgHEA`!?r )rgc ف>r)r2gF=r )rg߹Q#B=r )rg0ޟ)rrg=~L]=)rrgD,G '=))rz)g+&?gk?gμm/?)rw)g?CL\f?g^N?g]AͭR?)r)ga>$?g;Lf?gaW?)r)gx?g43n?gwY?)r2)go7🤎?g +-oGn?gIQ9‹Y?)r)g T?gGn?g@̋Y?)r)gU?g{Gn?giG͋Y?)r)gDU?g{Gn?gH͋Y?))rz)g Q)?g6{?r)gB`"۹?)g?g8%?r)r)g388?gV%?r)r)g*j?gRS%?r)r2)g! |?gZ%?r)r)gy[|?gG%?r)r)g"3|?gW%?r)r)g!3|?gW%?r)rrrrrrrrrGrrrrrrrrrrrr s@r.rr%sI  [[0 #J #J #J3      [[UT4L12 [[189 6 9 6  [[6BC 4 C 4  [[2DE>E> [[ M NG NG  [[ N OF OF  [[ 67 2 7 2  [[ P QG QG  [[ C D2 D2r0rca]tRtRtoRtRtRtRtRt] PPR.RO4R4t ] PPR.RO4R 4t R t] PPR .RO4R 4t] PP!R R7] PPR.RO4] PPR.RO4] PPR.RO4R4444t] PPR.RO4] PPR.RO4R44tRtVtR#) TestNakagamiik&cN\PPR4VnR#)li{NrSrXs&r.rZTestNakagami.setup_methodl&r\r0cdRp^p\PPW!4p\VR4R#)rNgM+)r4nakagamirr)rYr r^rs& r.rTestNakagami.test_logpdfo&s- ~~$$Q+12r0cRpRp\PPW!4p\VRRR7\PP W14p\WBRR7R#)rrTg"S+:rrN)r4rrrr)rYr r~rrs& r.ruTestNakagami.test_sf_isf&sM  ^^  r &1> ^^   'U+r0cj^pRpRp\PPW4p\WCRR7R#)rrzrrNgU iּ)r4rrr)rYr^r rrs& r.rTestNakagami.test_logcdf&s/  %&&q-%0r0cjRp^ pRp\PPW4p\WCRR7R#)rjrrNgD&`U)r4rrr)rYr^r rrs& r.rTestNakagami.test_logsf&s/  %$$Q+/r0zm, refc\\\PPV4VRR7R#)g^~z=rNrr4rrErYr rs&&&r.rtTestNakagami.test_entropy&s ..q13WEr0cX\\PPV4V4R#rtr r s&&&r.test_extreme_nuTestNakagami.test_extreme_nu&s ..q137r0c\P!\PP R44'gQh\P!\PP R44'gQhR#)rr|N)rTrr4r_entropyrXs&r.test_entropy_overflow"TestNakagami.test_entropy_overflow&sJ{{5>>2259::::{{5>>226:;;;;r0znu, refc\\\PPV4VRR7R#r)rr4rr)rYr rs&&&r.rTestNakagami.test_mean&s ++B/5Ar0z+Fit of nakagami not reliable, see gh-10908.rr rcrdca a ^do \PPS WW0PR7o \PP S 4wrEp\ WARR7\ WRRR7\ WcRR7V 3RlpV V 3RlpV V 3Rlp \ V!WEV4^RR7\ V!WEV4^RR7\ V !WEV4^RR7R #) rrr rcrdrrwrc<RV,^,\P!^SV, , 4,^V,V^,, \P!SV, 4,,#r r rTrp)r rcrdrzs&&&r. dlogl_dnu(TestNakagami.test_fit..dlogl_dnu&sPUQY"&&GcM):";;d5!8mbffWs]&;;< =r0cV<S^\P!V4,\^V4, ,^\P!\P!SV, V, 44,,\P!SV, V, ^,4, #r)rTrrrpr rcrdrLrzs&&&r. dlogl_dloc)TestNakagami.test_fit..dlogl_dloc&spRVVBZ)Ar*::;rvvw}&=>??@ffw}59:; .dlogl_dscale&sF!GbL5("fuz)BFFGcMa3G,HHI Jr0r2rN)r4rrrWrr) rYr rcrdnu_estloc_est scale_estrr r#rLrzs &&&& @@r.r TestNakagami.test_fit&s ..$$!+0xx%I%*^^%7%7%@"-3/ s3 = <  J  &9=qtL 6I>M Vi@!$Or0cRp^dp\PPWCVW PR7p\PP WSR7wrgp\ P !V4p \ P!\ P!WW, ^,44p \WcRR7\WyRR7\WRR7R#)rzrrwrGrN) r4rrrWrrTrrrr) rYrcrdr ryrzr%r&r'loc_theo scale_theos &&& r. test_fit_nuTestNakagami.test_fit_nu&s  ..$$!+0xx%I%*^^%7%7%7%G"66'?WWRWWg&7A%=>? .5 D9r0rN))rg|Bd븿)rzr)rbg/ۿ))r|g}Ô%I)rgq_)gScAgB4U)gcAgu1U)rgb >%)rgWҙ\))rg2H?)rg:U?)rg>)r( rr)rrbr")rrrQ)rrrrrZrrurrrrrrtrrrr_ r r,rrr s@r.rrk&s[53",&10 [[X %& F & F [[X DE8 E8< [[Y>?B ?B [[KL [[T?3 [[UN3 [[Wk2P344MP8 [[UN3 [[Wk2:34:r0rc |a]tRtRtoRtRtRtRt]PPRRR.4]PPR RR] P,] P)^] P^] P,^d.4]PPR .RO4R 444t R tVtR #)TestWrapCauchyi'cN\PPR4VnR#)l3R"NrSrXs&r.rZTestWrapCauchy.setup_method'r\r0c\P!RR.RR..4p\P!R.R..4p\PP W!4pVP R 8XgQh\P !W!34UUu.uF$wrE\PP WE4NK& ppp\VP4VRR7R#uuppi) Q?rorzrrSr3rrNr) rTr r4 wrapcauchyrnr?nditerrrc )rYrr^rzrrg scalar_valuess& r.test_cdf_shape_broadcasting*TestWrapCauchy.test_cdf_shape_broadcasting's HHtTlS$K0 1 HHsecU^ $     &ww&   )+A6):<):Xb))--b5): < =u=?!a"))Brvvbeebdl/C,C"DRUU"JJKr0rrrrcrdc\PPWVRVPR7p\P !V^84'gQh\P !V^\P ,8*4'gQhR#)r)rrcrdrrN)r4r4rrWrTrrg)rYrrcrdr^s&&&& r.test_rvs_lie_on_circle%TestWrapCauchy.test_rvs_lie_on_circle'sd    1U&* ! Cvva1f~~~vva1ruu9n%%%%r0rNrr )rrr)rrrrrZr7r:rZrrrrTrgrBrrr s@r.r/r/'s5 >,L [[S5$-0 [[UT2bee8beeVQqwPS$TU [[W&67&8V1&r0r/c!RR\P4pV!RR7p\PP R4p\ \ RR7;_uu_4VPVR7RRR4R# +'giR#;i) c&a]tRtRtoRtRtVtR#)/test_rvs_no_size_error..rvs_no_size_geni*'c^#rr6rXs&r._rvs4test_rvs_no_size_error..rvs_no_size_gen._rvs+'sr0r6N)rrrrrHrrr s@r.rvs_no_size_genrF*'s  r0rJ rvs_no_sizer'i_Oz_rvs\(\) got (an|\d) unexpectedrrN)r4rCrTrUrVrr,r)rJrKrWs r.test_rvs_no_size_errorrL('s`%--"}5K ))   +C y(J K KS) L K K Ks A88 B zdistname, argscV\9d\P!RV R24\\V4p\ V\P 4'd\V4^8wdGVP!V!wr4\V\P4\V\P4^RreVP!.VOVNVN5!wrx\V\P4\V\P4R#VP!V!wr\V \P4\V \P4R#)z6skipping test for the support method for distribution .Nr) $skip_test_support_gh13294_regressionrrr2r4rArCr#r rrTr*) rZr4rJa0b0r r a1b1r,r-s && r.test_support_gh13294_regressionrT4's77 $$,:Q0 1 5( #D$++,, t9>\\4(FB RVV $ RVV $ "f2t2T262R R ||T"QQr0c\PP.RO.RO4wr\P!\P )\P )\P )\P .4p\P!\P \P \P \P .4p\W4\W4VPVP8XgQhVPVP8XgQh\PP..4wrE\P!.4\P!.4rv\WF4\WW4VPVP8XgQhVPVP8XgQh\PP.ROR.4wr\P!^\P .,4p \P!^\P .,4p \W4\W4VPV P8XgQhV PV P8XgQhR#)r9N)r9r9r9rr)rrrr) r4r8r rTr r+r*rr?) rPrQex_a0ex_b0rRrSex_a1ex_b1rr ex_a2ex_b2s r.,test_support_broadcasting_gh13294_regressionr\O's ZZ   m ># % + +v 55 5r0c\\PPV4V4\\PP V)4V4R#rt)rr4cosinernr)r^rNs&&r.test_cosine_cdf_sfrav's6ELL$$Q'2ELLOOQB'2r0c\\PPV4V4\\PP V4V)4R#rt)rr4r`rr)rzrNs&&r.test_cosine_ppf_isfrc~'s8 ELL$$Q'2ELL$$Q'(3r0c\PP\P)\P.4p\ VR4R#)grB@NgrB)r4r`rrTrgr )rs r.test_cosine_logpdf_endpointsre's1 <<   /Dd./r0cN\UUu0uFwr\V\4'gKVkK! ppp\UUu0uFwrVkK pppW#8XgQh\UUu0uFwrVkK ppp\ UUu0uFwrVkK pppWE8XgQhR#uuppiuuppiuuppiuuppirt)rrAr<rrr)rrdiscrete_distnamesinvdiscrete_distnamescont_distnamesinvcont_distnamess r.test_distr_params_listsrk's/;4l74'c2$l41@AgdTA  66 6*23(wtd(N3-89['$[9  .. .4A49sBB B B5 B!c\PP^RR7\PP^^R7R8XgQh\PP ^^4R8XgQhR#)r9r]r\)rr,r N)r4r_statsrr5r6r0r.test_moment_order_4rn's` NNAs+ >> qA .# 55 5 >>  1 % ,, ,r0c :a]tRtRto]P R4t]PPR.RO4R 4t ]PPR .RO4R 4t ]PP]PPR ]P!RR]PPR7]P!RR]PPR7]P!RR ]PPR7.4R44tRtVtR#)TestRelativisticBWi'c \P!\\4PR, 4p\P P VPRR7pV#)zSample data points for pdf computed with CERN's ROOT See - https://root.cern/ Uses ROOT.TMath.BreitWignerRelativistic, available in ROOT versions 6.27+ pdf calculated for Z0 Boson, W Boson, and Higgs Boson for x in `np.linspace(0, 200, 401)`. z-data/rel_breitwigner_pdf_sample_data_ROOT.npyzx,pdf,rho,gammarr r$s& r.ROOT_pdf_sample_data'TestRelativisticBW.ROOT_pdf_sample_data'sLww N ! ! ; < vv  /@ A r0zrho,gamma,rtolqVEB@gj+@‰FC@Gz@ND@_LU?c VVR,V8HVR,V8H,,pVR,VR,rv\V\PPWbVR7VR7R#)rhorr^rfr^rN)rr4rel_breitwignerrf)rYrrr{rrrr^rfs&&&&& r.test_pdf_against_ROOT(TestRelativisticBW.test_pdf_against_ROOT'sc$ !% (C /#G,5 7 cDK3 &&**1*?d r0zrho, Gamma, rtolcRp\P!RR^ 4p\PP WQVR7p\PP WaVR7pV!WaV,V4p\ WxVR7R#)c\P!V^,V^,V^,,,4p^\P!^4,V,V,V,\P\P!V^,V,4,, pW@^,V^,, ^,V^,V^,,,, #r)rTrrg)ErGammarr]s&&& r.rfFTestRelativisticBW.test_pdf_against_simple_implementation..pdf'sGGAqDAqD5!8O45ERWWQZ!#e+e3EEBGGAqD5L113AA1q(1a4q=89 9r0rjrr^rN)rTrWr4r|rrfr) rYr{rrrfrzr^rrs &&&& r.&test_pdf_against_simple_implementation9TestRelativisticBW.test_pdf_against_simple_implementation'sk : KKdB '  ! ! % %aE % :##''e'<!Y&t,r0z rho,gammar c Rp\PPV4p\PP WRVR7p\PP V^R7p\V^,V^,3W3RR7V^,^8XgQh\PP V^VR7p\V^,VRR7V^,V^,3^V38XgQhR #) z`Tests fit for cases where floc is set. `rel_breitwigner` has special handling for these cases. ls ; rl#rrcrrwrrrN)rTrUrVr4r|rrr)rYr{rrrWrrs&&& r. test_fit_floc TestRelativisticBW.test_fit_floc's&#ii##D)$$(( 4c) ##''1'5QQ(3,TB1v{{##''1U'CA$/AAAu:---r0r6N))rtrur)rvrwr )rxryrt))rtrurt)rvrwrt)rxryg& .>)rrrrrr\rrrrr}rrr-rrrrr s@r.rprp's  ^^$ [[    [[/2  -  - [[ [[ LL"F&++2B2B  LL!5 0A0A  LL & 0A0A    . .r0rpcha]tRtRto]P P RRR.4R4tRtVt R#) TestJohnsonSUi (r@c|\PP!VR,RR/p\W!R,RR7R#):NrNrr:rNNr rN)r4rr)rYr@rs&& r.test_moment_gh18071!TestJohnsonSU.test_moment_gh18071 (s/oo##T"X>v>"XE2r0r6N)rd rRrrg={?g?5mV>gn?g\}\1@)g`fo@g$ں _@r9rgvUg N#2QsV )=>r0rTFrrlrcJ\PPR4pR wrgr\P!WgWR7p V P RVR7p /p V'dWR&V'dWR&V'dWlR&V'dW|R&V'dkV'dcV'd[V'dSRp \ P!\V R 7;_uu_4\PP!V 3/V BR R R 4R #\\PV 3/V BR # +'giR #;i) lZ #LmrbrrrrrrrrN)rrrr) rTrUrVr4rrrr r(rr)rYrrrlrrWr-rrcrdrJrrr s&&&&& r.r TestTruncPareto.test_fit:(s ii##$78+c  3<xxSsx3 L "N J J yUuKG|7;;!!%%d3d3<; *%*;*;T JT J<;;s "D D" r6N) rrrrrPrrrrr rrr s@r.rr(s"?. [[Yu 6 [[[4-8 [[WtUm4 [[WtUm4K5597Kr0rc&a]tRtRtoRtRtVtR#) TestKappa3iW(c^\PPRR4, p\PPRR4p\ W!4R#)rrzgj@N)r4kappa3rnrr)rYsf0rs& r.rTestKappa3.test_sfX(s;%,,""3,,llooc3'!r0r6N)rrrrrrrr s@r.rrW(s""r0rca]tRtRto]P !^^4t]P!^4t]P!^ 4t Rt Rt Rt Rt RtRtRtR tR tR tVtR #) TestIrwinHallia(cP\VPPR4R4R#)rN)rg?r9gQ)r ih10r4rXs&r.test_stats_ih10TestIrwinHall.test_stats_ih10f(s TYY__V46JKr0c.ROp\\V44Uu.uF$q PPV^,4NK& pp\ W14VPP^24pRp\ WE4R#uupi)rN) rgUUUUU9@g0a@igx@gn۶mˠ@g:PAg8)AgSAg2⍚H)rr#rrr )rYrEryrm50 m50_exacts& r.test_moments_ih10TestIrwinHall.test_moments_ih10l(sj-38D 2BC2BQ99##AaC(2BCW+iir"H S, Ds*A<c\P!^^^d4pVPPV4pVPPV4p\ W2^ R7R#r9maxulpN)rTrWunifrfih1r )rYptspdf_unifpdf_ih1s& r.test_pdf_ih1_unifTestIrwinHall.test_pdf_ih1_unifz(sDkk!Q$99==%((,,s#Wr:r0c$\P!^4p^ep\P!^^V4p\P!^^V4p^WB^,^,R, WB^,^,R%VP V4p\ WT^ R7R#)r Nr)r4r rTrWrfr )rYih2nptsrrNpdf_ih2s& r.test_pdf_ih2_triangle#TestIrwinHall.test_pdf_ih2_triangle(svooa kk!Q%;;q!T*%&!8/2B)C%C(q!"''#,Wr:r0c\P!^^^d4pVPPV4pVPPV4p\ W2^ R7R#r)rTrWrrnrr )rYrcdf_unifcdf_ih1s& r.test_cdf_ih1_unifTestIrwinHall.test_cdf_ih1_unif(sBkk!Q$99==%((,,s#Wr:r0c\P!^^ 4p\P!V4pVP V^, 4p\P !R\ V44p\W4^ R7R#)rrNrz)rTr'r4r rnryr#r )rYryihih_cdfexacts& r.rZTestIrwinHall.test_cdf(sO IIa  __Q A #s1v&V26r0c$.ROp\VPP\P!^ 44V^ R7\VPPR4R^ R7Rp\VPPR4V^ R7R#)r9rN) r9g\xO~>g-P2?g8ʝi?gg{Ȇ?rzg&M?g׈Y?gB?gD l?rrg+|gP3NV?g*8f?gNgX?gjS?r)r rrfrTr'r+s& r.test_pdf_ih10_exact!TestIrwinHall.test_pdf_ih10_exact(sDJ $tt*,,TYY]]299R=94Kr0c 2\VPP\P!^ 44^VPP \P!^ 44, 4Rp\ VPPR4V^ R7R#)rrNrSr)rrrrTr'rnr )rYrs& r.test_sf_ih10_exact TestIrwinHall.test_sf_ih10_exact(s\ RYYr]3QryyQS}9U5UV2TYY\\$/R@r0r6N)rrrrr4rrr rrrrrrrrZrrrrrr s@r.rra(sh ==A D //! C ??2 DL - ;;;7C LAAr0rc&a]tRtRtoRtRtVtR#)TestDParetoLognormi(cRwrr4p\P!W#WE4p\PP VP V4R4\PP VP V4R4R#)rRg6Z?g@# ?N)rRr!rrr )r4dpareto_lognormrTrIrrfrn)rYr^r r?r,r-rJs& r.test_against_R!TestDParetoLognorm.test_against_R(s[, aA$$Q10 ""488A;0CD ""488A;0CDr0r6N)rrrrrrrr s@r.rr(sEEr0rr@cVwrr4p\P!R4pT;'gRpT;'gRpT;'g^pT;'gRp\\V4p\ \ 4V,pV!V!p \P !W&4p V PV PV 44p \WWER7^\P !W6^4, p V PV PV 44p \WWER7R#)rzrrNii) rTrr2r4r%rrrrr) r@rZlp1lp2rrlpmrJr dist_frozenrrs & r.test_sf_isf_overridesr(s&*"H3d ((3-C ++C **C 991D ==5D 5( #D (^H %F-K ++c C ..- .CC43 bkk#B' 'C ..- .CC43r0rt)NF)rr r rrrrr re rrrrrj waldg@gףp= 1@gRZ@gS㥛Ԋ@g&@r iiiki1i4lEI.gʎK$/gԆX1))r rg|1q?)r r g~?)rrV<3T1=)7yQCru lY<)rgvHGBg^b!*1=)rrgS?))r g?)rr)rr)g6ްPg8R#/))g?g*wӄZ>)333333?gN9?)rzrS)rg8<,O@)r g:B)gl ׶w/gMb}[P))r rgwo {?)rQrgGd-P/)rbrgof?)rrg Oul";))r!rg7Ըh?)rrgz]r>)rg"@goBڱ9)rrg$=cC;)rzrg9?)geK r@ rC rF rH rJ r\ rb rf ri rl rp rt rw ry r{ r~ r r r r r r r r rCr r r r r r$rl r]r`rbrdrirqrsrurxrrrrrrrrrrrrr/rLrTr\r^rarcrerkrnrprrrrrrr6r0r.rs   DDD*"4/338( 00;>C==ii((1,||x'MMh.>.>.@H.L )6z'B$$ 8 J@J@Z4D..>? 6 ? 699>G,G,T**>O,O,d<B<B~D/D/N##>*7D7Dt C CA2A2H44DI6I6X**. G GFJ"J"Z=I=I@`&`&FjDjDZDD,DD,Y<Yw>t 3P6P6fB=B=J&&D$$*] /] /@))   , lGlG^b,b,Jw;w;t]H]H@",",JE&E&P668="="@u0u0p...SFSFl..2J/J/Z ? ?Fn%n%b66(00D N NRBRBj2,2,jCCL>6>6BD0D0N <'0EF/G/ <*585.:#4L ,)T)TX:72 -& 9%&P-M"+ BHHo7 89 BHHo7 89 rxx|456 BHHk;/ 01 . 2  226 16 7:);<)<):; ) ; );<!<! ();< ) < ) %*$$ %%% %%%%%%u9u9p8B8Bv'0DE@F@<   , 3#,7.)X 60D)FG 4 G 4 )<<>?4?4  !" "#_ =  = *@^K^KB.3+3+lA2A2HR:R:j%&%&P *)?[+HI J 4#26;467373 784 84 0 /-(W.W.t 3 3 :K:Kz""]A]A@EE*"GH 4 H 4r0