L i#ddlZddlZddlmZddlmZmZddlZddl m Z ddl m Z ddl mZddlmZddlmZdd lmZgd Zd Zd Zgd ZgdZgdZgdZgdZdgZeezedZeedZdhZ gdZ!dZ"ejFjHejFjKde"ejFjKdddgdZ&d1dZ'dZ(dZ)dZ*ejFjKd e jVd!fe jXd"fgd#Z-d$Z.d%Z/d&Z0ejFjKd'e0d(Z1 d2d)Z2Gd*d+Z3gd,Z4Gd-d.Z5Gd/d0Z6y)3N)assert_allclose assert_equal)statsdifferential_evolution)distcont)FitError) distdiscrete)goodness_of_fit)i'??)dpareto_lognorm gausshypergenexpongengamma irwinhallkappa4ksonekstwoncfncx2 truncexpon tukeylambdavonmises levy_stable trapezoidtruncweibull_minstudentized_range) betaprime crystalball exponweibf geninvgauss jf_skew_tnct recipinvgaussrel_breitwigner vonmises_line) burrchichi2mielkepearson3genhalflogisticrdistparetopowerlaw powerlognorm wrapcauchyr*)4alphar"r,burr12cauchyr-r.r#dgammardweibullr% fatiguelifefisk foldcauchy genextremer genhyperbolicgennorm genpareto halfcauchyinvgamma invweibullrr' johnsonsukappa3rrlandaulevylevy_lr loglaplacelomaxr/nakagamirr(rr3r5 powernormr* skewcauchytrtriang truncparetor rr!)argusexponpowr$rrr1 halfgennormgompertz johnsonsbr kstwobignr)rrr+r6)MMMLErR)erlangr@ norminvgaussc#FKtD]\}}|tvs||fywN)r skip_fit)distnameargs `/mnt/ssd/data/python-lab/Trading/venv/lib/python3.12/site-packages/scipy/stats/tests/test_fit.pycases_test_cont_fitrcas- " # 8 #C-  s! !z distname,argmethodrZrYc 2ttjdd}ttjdd}|t|vr|sd}t j ||t |vr|sd}t j|tt|}tj|ddgg}tjtj|tztj|j d zt"gd }t$D]*} tj&j)d } tj*d 5|j,|| | d} |dk(r|j.t0vrdd i} ni} |j2| fd|i| } |dk(rLtj4| }|j2|fi| }d|j.d| d|}t7|| d||dk(r|t8vrd}tj:| | f}|d|d fxxdzcc<|d|dfxxdzcc<|j=dtj4|}|j2|fi| }d|j.d| d|}t7|| d|ddd |z }tjtj> jAtzt"g|d <tjBtjD| r tGd!tjHtj>||ks+yd"tK|d#}|d$tK d#z }|d%tKd#z }tGd&|j.d#|z#1swYxYw)'N SCIPY_XFAILF)default SCIPY_XSLOWz@Failure expected; set environment variable SCIPY_XFAIL=1 to run.z9Very slow; set environment variable SCIPY_XSLOW=1 to run.g?riignore)allsize random_staterZflocrdzCDifferent results fitting uncensored data wrapped as CensoredData: z: est=z est1=g|=)rtolerr_msgGz?r)\(?axis)intervalz2Different results fitting interval-censored data: z est2=皙?znan returned in fitz parameter:  z estimated: z diff : zfit not very good in )&intosgetenv failing_fitspytestxfail xslow_fitsskipgetattrrnphstackmaxvstackthresh_percentfullnumargs thresh_min fit_sizesrandom default_rngerrstatervsname mle_use_floc0fit CensoredDatarfail_interval_censored column_stacksortabsmeananyisnanAssertionErrorrlstr)r`rard run_xfail run_xslowmsgdistfntruearg diffthresholdfit_sizerngrkwdsestdata1est1nicrxdata2est2difftxts rb test_cont_fitrkshBIImU; &K&MNM3Lii##D) [[X & C&**S(sCC6;;-#?{&**S8848C**3/!vzz%040))/ VC5tfNcsC83I#I??C:6#q!T)!#q!T)! 1 %**H=!vzz%040!!' VC5tfFccB= C@W}FFBFF388:$6~$E$.$01 b 66"((3-  !67 7vvbffTlm34]3L`CL>, SXJb)) SYKr**4V[[MDsJKK_ C Cs D?N  N cptt|}|j|dd}t|||d|y)Nrzz poor mle fit of (loc, scale) in )atolrr)rrrr)rdatadesiredrdactuals rb_check_loc_scale_mle_fitrs;tA UU4[ FFG$>tfEGcrtjgd}td|ddgdtd|ddgdy)N)ru{Gz?rrGz?rr@uniformrurtMbP?expongp= ף?)rarrayr)rs rb"test_non_default_loc_scale_mle_fitrs5 88D EDYtTlDAWdT7OTBrcngd}tjj|d}t|ddgdy)zgh-6167)rrrrrjrjrjrjr)rprirrN)rrrr)rphats rbtest_expon_fitrs- #D ;;??4a? (DD1c(.rc|tjtjdtjdg}d}t j t |5t jt5tjj|ddddddy#1swYxYw#1swYyxYw)Nz0Optimization converged to parameters that are...match) r concatenatezerosonesrraisesr warnsRuntimeWarningrbetar)rmessages rbtest_fit_errorrs~ >>288B<5 6D@G xw / LL ( ts$B25 B&B2&B/ +B22B;z dist, params)?@) 333333?rjcZtjjd}t|dr |j}n |j }|j |d|d}||g|j }|j||}|j||}t||t||y)Nig:pdfdrm) rrrhasattrlogpdflogpmfrsumnnlf_penalized_nnlfr)distparamsrlogpxfxrefres1res2s rbtest_nnlf_and_related_methodsrs ))   *CtU&s5A ! f  ! ! # #C 99VQ D    *DD#D#rc#Khd}hd}hd}tttzD]}||vst|ts:d}t j |t jj|Q||vr:d}t j |t jj|||vr:d}t j |t jj||yw) N> rrSrrfoldnormrMr truncnormrrR poisson_binomr r!>3r%rPr-r=ricer7binomgammarLr2r8r[rGr/nbinomr3arcsinerAlognormrandintrbradfordr;rTrrVrDinvgaussr0r4skewnormr" exponnormrBr'rWrFrXrNr>r?rErKr6r#r<rU weibull_max weibull_minr\r)>rr(r,rrrQskellamzipfianrr betabinomr$ hypergeomr betanbinomr loguniform nhypergeom reciprocalr&rr5r@rr1r*nchypergeom_fishernchypergeom_walleniusztested separatelyreasonmarkstoo slow (>= 0.25s)too slow (>= 1.0s)) dictr r isinstancerrparammarkrslowxslow)skip_basic_fitslow_basic_fitxslow_basic_fitrrs rbcases_test_fit_mlers'N RN=O\H,-  > !D#)>(F,,t6;;+;+;6+;+JK K ^ #*F,,t6;;+;+;6+;+JK K _ $)F,,t6;;+<+rSrrr2rr4rrrrr&rrr\r@rr!>1rPr-r.r=waldr7rrLanglitr[rGr3rQrmaxwellrrrr;rTrrrVloggammarrr$rBrrFrXrNr>r?rErKrrr6r#r< genlogisticrRrr semicircularr)r1laplace_asymmetric>"r%rr(rr,rrrSr8r:rr/rArrrDrrMr0r"r'rWrrrrrUr5r+rr*r rrrzFails. Oh well.rrrr zignore::RuntimeWarning) r r r r rrr r rrrfilterwarnings)rrrwarns_basic_fitrrr s rbcases_test_fit_msers  NJN$ 3O!kO\H,- > !D#)>&F,,t6;;+;+;6+;+JK K ^ #*F,,t6;;+;+;6+;+JK K _ $)F,,t6;;+<+r)r!)r r itemsr r)r`shapess rbcases_test_fitstartr"UsH N002&8S)BB  sAAzdistname, shapesc&tt|}tjj d}|jd}tj dd5|j |}ddd|jddsJy#1swYxYw)NiV rrk)invaliddividerz)rrrrrr _fitstart _argcheck)r`r!rrrguesss rb test_fitstartr)]s| 5( #D ))   *C ::b>D Xh 7%t$% 4>>5": && &%%s BBct||}|||}|||} || ks$tjj|| ||yy)N)rqr)rrtestingr) rrparams1params0rqr nlff_namenlffnlff1nlff0s rbassert_nlff_less_or_closer2isM 4 #D $ E $ E EM ""5%d"F rceZdZejZdZejjeZ ejjddde Z ddgZ dddZd Zd Zeed Zd d dZdZdZdZdZd)dZej2j5dedZej2j5dedZej2j>dZ ej2j5dddZ!ej2jDdZ#dZ$dZ%ej2j>dZ&dZ'ej2j>d Z(d!Z)d"Z*d#Z+d$Z,ej2jDd%Z-d&Z.d'Z/y()*TestFitlVarrrmrrrrnpry{Gz?)rrqrrct|d|i|SNrr)selfrargsrs rboptz TestFit.opt}s%t====rcd}tjt|5tjd|j |j dddy#1swYyxYw)Nz `dist` must be an instance of...rr)rr ValueErrorrrrshape_bounds_ar?rs rb test_dist_ivzTestFit.test_dist_ivsD4 ]]:W 5 : IIb$))T%8%8 9 : : :s ,AAc8d}tjt|5tj|j gdg|j dddd}tjt|5tj|j dddtjg|j dddtjt|5tj|j dddtjg|j dddtjt|5tj|j gd|j dddy#1swY xYw#1swYxYw#1swYtxYw#1swYyxYw) Nz'`data` must be exactly one-dimensional.rrrjz.All elements of `data` must be finite numbers.rrjrI)123) rrrCrrrrDrnaninfrEs rb test_data_ivzTestFit.test_data_ivs5; ]]:W 5 C IIdii)d.A.A B CC ]]:W 5 I IIdii!Q266!2D4G4G H I ]]:W 5 I IIdii!Q266!2D4G4G H I ]]:W 5 G IIdii$2E2E F G G C C I I I I G Gs//E+2>E8>F4.F+E58FF Fcd}dddd}tjt|5tj|j |j |dddd}ddg}tjt|5tj|j |j |dddd }dd g}tjt|5tj|j |j |dddd d g}tjt|5tj|j |j |dddd }dg}tjt|5tj|j |j |dddd}gd}tjt|5tj|j |j |dddd}ddd}tjt|5tj|j |j |dddd}ddg}tjt|5tj|j |j |dddd}ddg}tjt|5tj|j |j |dddd}tjt|5tj|j |j dddtj tjfdg}tjt|5tj|j |j |dddy#1swYsxYw#1swY+xYw#1swYxYw#1swYxYw#1swYVxYw#1swYxYw#1swYxYw#1swY}xYw#1swY5xYw#1swYxYw#1swYyxYw)Nz1Bounds provided for the following unrecognized...r6r7rrr9r:rJrz6Each element of a `bounds` sequence must be a tuple...)rrrIz6Each element of `bounds` must be a tuple specifying...)rrrrrz7A `bounds` sequence must contain at least 2 elements...z;A `bounds` sequence may not contain more than 3 elements...)r6r6r6r6z.There are no values for `p` on the interval...)rrr8z.There are no values for `n` on the interval...)rrz6There are no integer values for `n` on the interval...)gffffff?g?z0The intersection of user-provided bounds for `n`) rrrrrrrrrCrrN)r?r shape_boundsboundss rbtest_bounds_ivzTestFit.test_bounds_ivs@E$6@ \\. 8 : IIdiiL 9 :K"F+ ]]:W 5 : IIdiiL 9 :K"K0 ]]:W 5 : IIdiiL 9 :1v ]]:W 5 : IIdiiL 9 :Ly ]]:W 5 : IIdiiL 9 :P5 ]]:W 5 4 IIdiiF 3 4C$62 ]]:W 5 : IIdiiL 9 :C( ]]:W 5 : IIdiiL 9 :K"F+ ]]:W 5 : IIdiiL 9 :E ]]:W 5 , IIdii + ,&&"&&)62 ]]:W 5 : IIdiiL 9 : :[ : :  : :  : : : :  : :  4 4  : :  : :  : : , , : :s,O9,O,O!,O(5,O5 ,P ,P5,P ,P)+P6 ,QO OO%(O25O?P PP&)P36P?Q cd}dddd}tjt|5tj|j |j |j|dddd}dd d }tjt|5tj|j |j |j|ddddd g}tjt|5tj|j |j |j|ddddd gg}tjt|5tj|j |j |j|dddd }dg}tjt|5tj|j |j |j|dddd}gd}tjt|5tj|j |j |j|dddd}ddd }tjt|5tj|j |j |j|dddd}gd}tjt|5tj|j |j |j|dddd}ddd }tjt|5tj|j |j |j|dddd}gd}tjt|5tj|j |j |j|dddy#1swYsxYw#1swYxYw#1swYxYw#1swYyxYw#1swY&xYw#1swYxYw#1swY}xYw#1swY)xYw#1swYxYw#1swYyxYw)Nz2Guesses provided for the following unrecognized...rrrRr)r(z+Each element of `guess` must be a scalar...hir8r%rjz-A `guess` sequence must contain at least 2...z1A `guess` sequence may not contain more than 3...)rrjrIzCGuess for parameter `n` rounded.*|Guess for parameter `p` clipped.*g@gz$Guess for parameter `loc` rounded...)r5rrz"Guess for parameter `p` clipped...r5z$Guess for parameter `loc` clipped...)r5rr) rrrrrrrshape_bounds_drrC)r?rr(s rb test_guess_ivzTestFit.test_guess_ivs]Fc, \\. 8 N IIdiiD,?,?u M N@d# ]]:W 5 N IIdiiD,?,?u M NC ]]:W 5 N IIdiiD,?,?u M NQ ]]:W 5 N IIdiiD,?,?u M NB ]]:W 5 N IIdiiD,?,?u M NF ]]:W 5 N IIdiiD,?,?u M NX% \\. 8 N IIdiiD,?,?u M N9 \\. 8 N IIdiiD,?,?u M N7d# \\. 8 N IIdiiD,?,?u M N9 \\. 8 N IIdiiD,?,?u M N NS N N  N N N N N N  N N  N N  N N  N N  N N  N Nsw8O 8O%8O'8O4%8P8P(8P 8P(+8P5 8Q OO$'O14O>P PP%(P25P>Q c2d}tttz}tjj |j }tt|}tj||}tjt|dzdftj}|dtj|zz |dddf<|dtj|zz|dddf<d|d<d |d <|j|d} |j|d } t||| | gz} t|d d r5| dd } tj | | d <|j"| ||d } |dd }t|dd r|j"| ||d } t%j&5} | j)t*dtj,| |||j.}dddddd}||}t1| j2| fi|j4d|iy#1swY;xYw)Nrrj)dtypeg$@rzrrrQ)gؗҜ))D$v&*hh0C 0 $,>? v& !$cjj# 7 7,5 7 0 0s ;H  Hrqc,|j|ddy)Nrdr5r<ryr?rqs rbtest_basic_fit_mlezTestFit.test_basic_fit_mle Iu!4rc,|j|ddy)Nrerjr<r{r|s rbtest_basic_fit_msezTestFit.test_basic_fit_mser~rc@d}tjj|j}tj }d}|j |||d}ddd}t j|||d|j}t|||j|fdd i|jy) Nr )rirrm皙?rrtrurer`r.rb) rrrrfrrrrrAr2rrpr?rrrrr!rrSrws rb test_arcsinezTestFit.test_arcsines ii##DII.}}txxac:(9= iidL$((S!$cjj& M,> MBF)) MrrdrccVd}tjj|j}tj }d}|j |||d}dddd}t j||||j|}ddd |} t|||j|fi|jd | iy) Nr )rir@rmr)r-rtru)rardrrbrcr.) rrrrfrrSrrrAr2rrp) r?rdrrrrr!rrSrwr.s rb test_arguszTestFit.test_argus-s ii##DII.{{txxac:(YO iidLDHHVT"+=>vF !$cjj& 7DII 7,5 7rcDd}tjj|j}tj }d}|j |||d}ddddd}t j|||d|j}t|||j|fdd i|jy) Nr )gdsz@g@e?rirrmr)abrtrurdr`r.r) rrrrfrrrrrAr2rrprs rb test_betazTestFit.test_beta=s ii##DII.zzBtxxac:&I(9> iidL$((S!$cjj& A,2 A6:ii Arc<d}tjj|j}tj }d}|j |||d}dddd}t j||||j}t|||j|fi|jy)Nr )g);?rrrmr)crtrura) rrrrfrrrrrAr2rrprs rb test_foldnormzTestFit.test_foldnormMs ii##DII.~~,txxac:&y9M iidLDHHE!$cjj&NDIINrcBd}tjj|j}tj }d}|j |||d}dddd}t j|||d|j}t|||j|fdd i|jy) Nr )g7h@rirrmr)nurtrurdr`r.r) rrrrfrrMrrrAr2rrprs rb test_nakagamizTestFit.test_nakagami[s ii##DII.~~-txxac:&y9M iidL$((S!$cjj& A,2 A6:ii ArcBd}tjj|j}tj }d}|j |||d}dddd}t j|||d|j}t|||j|fdd i|jy) Nr )gk?rirrmr)rrtrurer`r.rb) rrrrfrr4rrrAr2rrprs rb test_powerlawzTestFit.test_powerlawis ii##DII.~~-txxac:&y9M iidL$((S!$cjj& M,> MBF)) Mrc<d}tjj|j}tj }d}|j |||d}dgdz}t j||||j}t|||j|fi|jy)Nr )g?g333333@gffffff@gffffff@rmrrYr) rrrrfrrRrrrAr2rrprs rbtest_truncparetozTestFit.test_truncparetoxs ii##DII.  %txxac:!{1} iidLDHHE!$cjj&NDIINrc<d}tjj|j}tj }d}|j |||d}dgdz}t j||||j}t|||j|fi|jy)Nr )rrg?rrrmrr5r) rrrrfrr rrrAr2rrprs rbtest_truncweibull_minzTestFit.test_truncweibull_mins ii##DII.%%*txxac:!{1} iidLDHHE!$cjj&NDIINrcPd}tjj|j}tj }d\}}}|j |||||}dtjddgi}t j||||j} t| j|||ffi|jtj}d\}}|j ||||}t j|||j} t| j||ffi|jy) Nr r?rrtrnror9rr)gY.?r)rrrrfrrrrrrArrrp bernoulli) r?rrrrr9r:rtrrSrws rbtest_missing_shape_boundsz!TestFit.test_missing_shape_boundss ii##DII.{{ 1cxx1#ACx@RXXq"g./ iidLDHHE Q3K=499=3xxsx=iiddhh7 QH: :rcd}tjj|j}tj }d\}}|j |||}d}d|i}t j||||j} t| j||ffi|jd\}}|j |||}d } d | i}t j||||j} t| j||ffi|jtj }d \}}|j |||| }|| d }t j||||j} t| j||ffi|jy)Nr)?rr)rr5rtr)rr)rurnro)r;r5rurrrtrurnror) rrrrfrnormrrrArrrp) r?rrrrrtrur loc_boundsrTrw scale_boundss rbtest_fit_only_loc_scalezTestFit.test_fit_only_loc_scalesO ii##DII.zz UxxCacx: $iidFdhh? S%L>DII> Uxxe!#x> <(iidFdhh? S%L>DII>zz UxxCu13xG#l;iidFdhh? S%L>DII>rcd}tjj|j}tj }d\}}|j ||||}t j||}t|jdfi|j||f||fd}t j|||}t|j||ffi|jtj}d\} } }|j | | |||}dd d } t j||| |j }t|j| | |ffi|jy) Nrrrr7rrr)rr)rrr8r) rrrrfrrrrrrrprrA) r?rrrrrtrurrwrTr9r:rSs rbtest_everything_fixedzTestFit.test_everything_fixeds ii##DII.zz UxxCu13xGiid# F8dii8suen=iidF+ S%L>DII>{{ 1cxx1#ACx@$<8 iidLDHHE Q3K=499=rcbd}tjj|j}tj }d}|j |||d}|jdk(sJgd}t j|||}d}|jj|sJ|jdusJy)Nr)r5rrmr)rr7)rrz3Optimization converged to parameter values that areF) rrrrfrrrminrr startswithsuccess) r?rrrrr!rrTrwrs rb test_failurezTestFit.test_failures ii##DII.||txxac:xxzQ-iidF+G{{%%g...{{e###rcd}tjjd}tj}d}gd}|j |||d}tj ||||j}tj|j|fi|jrJtj |||||j}t|j|fi|jy)Nil^s)r r))rj)ffffff?F)g333333?xrQrmr)r(ra) rrrrrrrrAallcloserrpr)r?rrrrrrTrrws rb test_guesszTestFit.test_guesss ii##L1;txxac:iidFdhh?;;szz6?TYY???iidF&DHHM F8dii8rc"ddg}tj}ddd}tj|||d}t|jj dd tj|||d }t|jj d d y) NrjrY)rr:0yE>rrrdrTrdrIrrregRQ@)rrrrrru)r?rrrTres_mleres_mses rbtest_mse_accuracy_1zTestFit.test_mse_accuracy_1st1v{{*5))D$veD,,ad;))D$veD,,e$?rctjjd}tj}d}|ddj ||}ddd}tj |||d }tj|}||d z|d z |d z z }||d z|d z |d z z } || |z f} t|j| dy)Nl4gcvDrrIrmrQrrrerrr^r-C6?rq) rrrrrrrrrr) r?rrr9rrTrwrrrrs rbtest_mse_accuracy_2zTestFit.test_mse_accuracy_2 sii##$78}} Aqz~~13~7 :6iid6%@ GGDM qtVae^a!e $ quWqt^a!e $1f Cd3rN)r)0__name__ __module__ __qualname__rrrrfrrrrrrrDrZrrqrprArFrOrUr[ryrr  parametrizerr}rrrrrrrrrrrrrrrrrrrrrbr4r4rs ;;D D ))   %C ;;??1c#? >Dv&N"0N D D$ 'D>: G1:f-N^ 7D [[[*<*>?5@5 [[[*<*>?5@5 [[ M M [[X~6 77 7 [[ A A O A [[ M M O [[ O O;*?:>2$  [[99" @4rr4)xA=QXEY7TVrGrrrNC`BIK;rr?OLrUWrPrrrrrrr@rrMrrrrr\Zrrrrr:rrrrRrr6r9rrrJrrrrrrDrrSrrrrrrrrrrrrrrrrrrr]rrrrrr<rrrrrrrrrrrrrHcheZdZdZdZdZdZdZdZdZ dZ e jjd Zd Zd Ze jj#d d gdfdgdfdgdfgdZe jje jj#d gddZdZdZy)TestGoodnessOfFitctj}gd}d}tjt|5t tj|dddd}tjt |5t |gdgdddd}tjt |5t ||ddddd}tjt|5t ||d dddd }tjt|5t ||d dddy#1swYxYw#1swYxYw#1swYxYw#1swYbxYw#1swYyxYw)NrHz.`dist` must be a \(non-frozen\) instance of...rz2`data` must be a one-dimensional array of numbers.z`statistic` must be one of...mm) statisticz"`n_mc_samples` must be an integer.gD@) n_mc_samplesz$SeedSequence expects int or sequenceherringr<)rrrr TypeErrorr rC)r?rrrs rb test_gof_ivzTestGoodnessOfFit.test_gof_iv-s&zz C ]]9G 4 - EJJL! , -G ]]:W 5 / D9+ . /2 ]]:W 5 5 D!t 4 57 ]]9G 4 : D!& 9 :9 ]]9G 4 4 D! 3 4 4! - - / / 5 5 : : 4 4s;D+6D7+EEE+D47EE EE$ctjjd}t}tj|tj |dd}t tj||d|}tj|tjd i|jd}t|j|jt|j|jd y) NlH]Vh1rddofrks known_paramsrrexactrd{Gzt?rr)rrr examgradesrstdr rrkstestcdfrrpvaluer?rrrrwrs rbtest_against_ksz!TestGoodnessOfFit.test_against_ksEsii##$78 !wwqzBFF114EF ejj!,(,#7ll1ejj8<8< ejj!,(,#7 u- Dt4rctjjd}tjdd}ddi}t t j ||d|}t|jdt|jd d y) NlH'b4:/rrrug{U=@rrguV?rrrrrs rbtest_against_anderson_case_1z.TestGoodnessOfFit.test_against_anderson_case_1nsiii##$78 IIa !34 ejj!,(,#7 u- Cd3rctjjd}tjdd}ddi}t t j ||d|}t|jdt|jd d y) NlP^QqPrrrtgfHF@rrgZd;@皙?rrrrs rbtest_against_anderson_case_2z.TestGoodnessOfFit.test_against_anderson_case_2ysiii##$67 IIa /0 ejj!,(,#7 u- E5rc tjjd}tjj dddd|}t tj|d|}t|jd t|jd d y) NliUngBQ_ *?rrjrrrrrgS?g333333?rr) rrrrrrr rrrr)r?rrrws rbtest_against_anderson_case_3z.TestGoodnessOfFit.test_against_anderson_case_3snii##$78 NN  1q,/  1ejj!tE u- Dt4rctjjd}tjdddj d|}t tj|d| }tj|d }t|j|jd t|j|jd d z dy)Nl/JGg Q8?rrrr rmrrgumbel_r)rrrrr)rrrrr?rr r"andersonrrcritical_valuesrsignificance_level)r?rrrwrs rbtest_against_anderson_gumbel_rz0TestGoodnessOfFit.test_against_anderson_gumbel_rsii##$45   1s#& ((+C(H enna4"%'nnQZ0 s':':1'=> C$:$:1$=c$AMrchtjjd}gd}ddd}tjtj ||d|}t |jdd d |jcxkrd ksJJt |jd dt |jddy)NlfKW7/z)rrrzr5rrrrfillibenrgG=D;?rrr?gF8y?gh㈵>rgffffff?Mb`?) rrrrr rrrr)r?ryrrws rbtest_against_filliben_normz,TestGoodnessOfFit.test_against_filliben_normsii##$78 # 1- ##EJJ .8cC  wT:cjj&3&&&&&  '7dC F6rcVtjjd}|jddd}t j tj |d|}dd d }t j tj ||d| }t|j|jd y)Nl1gDF3rrr)rtrurnr*rrrrrgV瞯? %%c&;&;[_MSt,r))r5g?g;On?)rgDŽ?g -?)rg-hC/?g\ Ac?c |\}}}tjjd}|j|}t j tj |d|}t|j|dt|j|dy) NlvTV_)rnr*rrrgQ?r) rrrr0rr rayleighrrr)r?r2r9 ref_statistic ref_pvaluerrrws rbtest_against_ppccz#TestGoodnessOfFit.test_against_ppccsr(,$=*ii##$78 JJAJ ##ENNA(+- }4@ JV>>T__++1159T__++//8qii##$78u00!!.<*4,8cC;;t5577#5577dD DD;;t55#55DB BBT__++1159T__++//8  %0 ii##$78u00!!.<*4,8cC T__++--t4T__++1159T__++//8;;t55t7M7MNNNNrc8d}tjjd}tjj d|}t tj|ddd||}gd }gd }tj|j|}t||d y) Nctj||}|j|}tj||dd}tj|dz|S)Nrvrr)rwprependappendrj)rrrrr)rrrwrr-rs rb greenwoodz:TestGoodnessOfFit.test_custom_statistic..greenwoodsH4(A Aa:A66!q&t, ,rrCr5rmrrrr) r;ryrg?rg?rg333333?rg?r+r4rt) g҈?gxܙ ?g&)?gW>?gR?gEկt>r)r@kwargsrs rbraz-TestFitResult.test_plot_iv..optimizer5s)4C6CsC Crrr7rz!`plot_type` must be one of \{'...rllama) plot_typez2matplotlib must be installed to use method `plot`.)rrrrrrr matplotlibrrrCplotModuleNotFoundError ImportError)r?rrarTrwrZrrs @rb test_plot_ivzTestFitResult.test_plot_iv1sii##$78zz~~a3~? D6"ii D&IF , :Gz9 ,7+ , , ,#[1 ,KG2'B ,7+ , , , ,sH5!B>B2)B>2B;7B>;B>>,D*D=DD DDN)rrrr^rrrbrTrT0s,rrTr^)gHz>rr)7r}numpyr numpy.testingr+rmrrrscipyrscipy.optimizertest_continuous_basicr !scipy.stats._distn_infrastructurer scipy.stats._distr_paramsr scipy.statsr rrrmle_failing_fitsmle_Xslow_fitsrmm_failing_fitsmm_XXslow_fits mm_Xslow_fitsrrrr_rcr rrrrrrrrrrrrr"r)r2r4rrrTrrrbrls 7 1+62'   *#   (=  %6?OP !. 9 '   )<)>?E4=1HL2@HLVGC / ::z2 ;; 578 8  +\<~+-@-BC'D'MN(.Gk4k4^ . CJCJJ,,r