U e5d ï¿½ï¿½@sjdddddddddd d ddd dddgZddlZddlZddlZddlmZddlmZddlm Z ddl mZmZddlm Z mZmZmZmZmZmZmZddlmZddlmZGddï¿½deï¿½Zdcddï¿½Zddï¿½Zdd ï¿½Zd!d"ï¿½Zd#d$ï¿½Zd%d&ï¿½Z d'd(ï¿½Z!ddd*d+ï¿½Z"d,dï¿½Z#d-dï¿½Z$d.dï¿½Z%d/dï¿½Z&d0dï¿½Z'd1d ï¿½Z(d2dï¿½Z)ded4dï¿½Z*d5d ï¿½Z+d6dï¿½Z,d7d8d9ï¿½d:dï¿½Z-dfd;d<ï¿½Z.dgd=dï¿½Z/dhd>d ï¿½Z0did?dï¿½Z1djd@dï¿½Z2dAdBï¿½Z3GdCdï¿½dï¿½Z4zddDl5m3Z3Wne6k ï¿½rï¿½YnXe7dEkï¿½rfddFlm8Z8ddGlm9Z9m:Z:m;Z;mZ>e4dIdJï¿½Z?e4dKdLï¿½Z@dMZAe?ï¿½BeAï¿½ZCe@ï¿½BeAï¿½ZDe9e:fD]ï¿½Pï¿½ï¿½dS)kï¿½ NormalDistï¿½StatisticsErrorï¿½fmeanï¿½geometric_meanï¿½ harmonic_meanï¿½meanï¿½medianï¿½median_groupedï¿½median_highï¿½ median_lowï¿½modeï¿½ multimodeï¿½pstdevï¿½ pvarianceï¿½ quantilesï¿½stdevï¿½varianceï¿½Nï¿½ï¿½Fraction)ï¿½Decimal)ï¿½groupby)ï¿½bisect_leftï¿½bisect_right)ï¿½hypotï¿½sqrtï¿½fabsï¿½expï¿½erfï¿½tauï¿½logï¿½fsum)ï¿½ itemgetter)ï¿½Counterc@seZdZdS)rN)ï¿½__name__ï¿½ __module__ï¿½__qualname__ï¿½r&r&ï¿½"/usr/lib64/python3.8/statistics.pyruscCsï¿½d}t|ï¿½\}}||i}|j}ttt|ï¿½ï¿½}t|tï¿½D]@\}} t||ï¿½}tt| ï¿½D]"\}}|d7}||dï¿½|||<qRq6d|krï¿½|d} ntddï¿½t|ï¿½ ï¿½ï¿½Dï¿½ï¿½} || |fS)Nrï¿½css|]\}}t||ï¿½VqdSï¿½Nr)ï¿½.0ï¿½dï¿½nr&r&r'ï¿½ ï¿½sz_sum..) ï¿½_exact_ratioï¿½getï¿½_coerceï¿½intï¿½typerï¿½mapï¿½sumï¿½sortedï¿½items)ï¿½dataï¿½startï¿½countr,r+ZpartialsZpartials_getï¿½Tï¿½typï¿½valuesï¿½totalr&r&r'ï¿½_sum{s$ r>cCs.z |ï¿½ï¿½WStk r(tï¿½|ï¿½YSXdSr))Z is_finiteï¿½AttributeErrorï¿½mathZisfinite)ï¿½xr&r&r'ï¿½ _isfiniteï¿½s rBcCsï¿½||kr|S|tks|tkr |S|tkr,|St||ï¿½r:|St||ï¿½rH|St|tï¿½rV|St|tï¿½rd|St|tï¿½r|t|tï¿½r||St|tï¿½rï¿½t|tï¿½rï¿½|Sd}t||j|jfï¿½ï¿½dS)Nz"don't know how to coerce %s and %s)r1ï¿½boolï¿½ issubclassrï¿½floatï¿½ TypeErrorr#)r:ï¿½Sï¿½msgr&r&r'r0ï¿½s( r0cCsï¿½zrt|ï¿½tkst|ï¿½tkr$|ï¿½ï¿½WSz|j|jfWWStk rnz|ï¿½ï¿½WYWStk rhYnXYnXWn ttfk rï¿½|dfYSXd}t |ï¿½ t|ï¿½jï¿½ï¿½ï¿½dS)Nz0can't convert type '{}' to numerator/denominator)r2rErï¿½as_integer_ratioï¿½ numeratorï¿½denominatorr?ï¿½ OverflowErrorï¿½ ValueErrorrFï¿½formatr#)rArHr&r&r'r.ï¿½s r.cCspt|ï¿½|kr|St|tï¿½r(|jdkr(t}z ||ï¿½WStk rjt|tï¿½rd||jï¿½||jï¿½YSï¿½YnXdS)Nr()r2rDr1rKrErFrrJ)ï¿½valuer:r&r&r'ï¿½_convertï¿½s rPcCs.t||ï¿½}|t|ï¿½kr&|||kr&|Stï¿½dSr))rï¿½lenrM)ï¿½arAï¿½ir&r&r'ï¿½ _find_lteq s rTcCs>t|||dï¿½}|t|ï¿½dkr6||d|kr6|dStï¿½dS)N)ï¿½lor()rrQrM)rRï¿½lrArSr&r&r'ï¿½ _find_rteqs rWï¿½negative valueccs$|D]}|dkrt|ï¿½ï¿½|VqdS)Nr)r)r<ï¿½errmsgrAr&r&r'ï¿½ _fail_negsrZcCsHt|ï¿½|krt|ï¿½}t|ï¿½}|dkr,tdï¿½ï¿½t|ï¿½\}}}t|||ï¿½S)Nr(z%mean requires at least one data point)ï¿½iterï¿½listrQrr>rP)r7r,r:r=r9r&r&r'r'scstzt|ï¿½ï¿½Wn0tk r<dï¿½ï¿½fddï¿½}t||ï¿½ï¿½}Yn Xt|ï¿½}z |ï¿½WStk rntdï¿½dï¿½YnXdS)Nrc3s t|ddï¿½D]\ï¿½}|VqdS)Nr()r8)ï¿½ enumerate)ï¿½iterablerAï¿½r,r&r'r9Oszfmean..countz&fmean requires at least one data point)rQrFr ï¿½ZeroDivisionErrorr)r7r9r=r&r_r'rAs cCs8ztttt|ï¿½ï¿½ï¿½WStk r2tdï¿½dï¿½YnXdS)NzHgeometric mean requires a non-empty dataset containing positive numbers)rrr3rrMr)r7r&r&r'r\sï¿½cCsï¿½t|ï¿½|krt|ï¿½}d}t|ï¿½}|dkr2tdï¿½ï¿½n<|dkrn|d}t|tjtfï¿½rf|dkrbt|ï¿½ï¿½|Stdï¿½ï¿½z"t ddï¿½t ||ï¿½Dï¿½ï¿½\}}}Wntk rï¿½YdSXt|||ï¿½S)Nz.harmonic mean does not support negative valuesr(z.harmonic_mean requires at least one data pointrzunsupported typecss|]}d|VqdS)r(Nr&ï¿½r*rAr&r&r'r-ï¿½sz harmonic_mean..) r[r\rQrï¿½ isinstanceï¿½numbersZRealrrFr>rZr`rP)r7rYr,rAr:r=r9r&r&r'ros$ "cCs\t|ï¿½}t|ï¿½}|dkr tdï¿½ï¿½|ddkr8||dS|d}||d||dSdSï¿½Nrï¿½no median for empty dataï¿½r(ï¿½r5rQr)r7r,rSr&r&r'rï¿½s cCsLt|ï¿½}t|ï¿½}|dkr tdï¿½ï¿½|ddkr8||dS||ddSdSrdrgï¿½r7r,r&r&r'r ï¿½scCs,t|ï¿½}t|ï¿½}|dkr tdï¿½ï¿½||dS)Nrrerfrgrhr&r&r'r ï¿½s r(c Csï¿½t|ï¿½}t|ï¿½}|dkr"tdï¿½ï¿½n|dkr2|dS||d}||fD]}t|ttfï¿½rFtd|ï¿½ï¿½qFz||d}Wn(tk rï¿½t|ï¿½t|ï¿½d}YnXt||ï¿½}t |||ï¿½}|}||d} |||d|| S)Nrrer(rfzexpected number but got %r) r5rQrrbï¿½strï¿½bytesrFrErTrW) r7Zintervalr,rAï¿½objï¿½Lï¿½l1ï¿½l2Zcfï¿½fr&r&r'rï¿½s& cCsHt|ï¿½}t|ï¿½ï¿½dï¿½}z|ddWStk rBtdï¿½dï¿½YnXdS)Nr(rzno mode for empty data)r[r"ï¿½most_commonï¿½ IndexErrorr)r7Zpairsr&r&r'rscCs@tt|ï¿½ï¿½ï¿½ï¿½}tt|tdï¿½dï¿½dgfï¿½\}}tttdï¿½|ï¿½ï¿½S)Nr()ï¿½keyr)r"r[rpï¿½nextrr!r\r3)r7ZcountsZmaxcountZ mode_itemsr&r&r'r5s ï¿½ï¿½ exclusive)r,ï¿½methodc CsL|dkrtdï¿½ï¿½t|ï¿½}t|ï¿½}|dkr0tdï¿½ï¿½|dkrï¿½|d}g}td|ï¿½D]N}|||}||||}||||||d||} |ï¿½| ï¿½qN|S|dkï¿½r:|d}g}td|ï¿½D]r}|||}|dkrï¿½dn||dkrï¿½|dn|}||||}||d||||||} |ï¿½| ï¿½qï¿½|Std|ï¿½ï¿½ï¿½ï¿½dS)Nr(zn must be at least 1rfz"must have at least two data pointsZ inclusiveruzUnknown method: )rr5rQï¿½rangeï¿½appendrM) r7r,rvZldï¿½mï¿½resultrSï¿½jZdeltaZinterpolatedr&r&r'rls4$ $$csï¿½ï¿½dk r,tï¿½fddï¿½|Dï¿½ï¿½\}}}||fSt|ï¿½ï¿½tï¿½fddï¿½|Dï¿½ï¿½\}}}tï¿½fddï¿½|Dï¿½ï¿½\}}}||dt|ï¿½8}||fS)Nc3s|]}|ï¿½dVqdSï¿½rfNr&raï¿½ï¿½cr&r'r-ï¿½sz_ss..c3s|]}|ï¿½dVqdSr|r&rar}r&r'r-ï¿½sc3s|]}|ï¿½VqdSr)r&rar}r&r'r-ï¿½srf)r>rrQ)r7r~r:r=r9ï¿½UZtotal2Zcount2r&r}r'ï¿½_ssï¿½srï¿½cCsLt|ï¿½|krt|ï¿½}t|ï¿½}|dkr,tdï¿½ï¿½t||ï¿½\}}t||d|ï¿½S)Nrfz*variance requires at least two data pointsr(ï¿½r[r\rQrrï¿½rP)r7ï¿½xbarr,r:ï¿½ssr&r&r'rï¿½s&cCsHt|ï¿½|krt|ï¿½}t|ï¿½}|dkr,tdï¿½ï¿½t||ï¿½\}}t|||ï¿½S)Nr(z*pvariance requires at least one data pointrï¿½)r7ï¿½mur,r:rï¿½r&r&r'rï¿½s#cCs8t||ï¿½}z |ï¿½ï¿½WStk r2tï¿½|ï¿½YSXdSr))rrr?r@)r7rï¿½ï¿½varr&r&r'rs cCs8t||ï¿½}z |ï¿½ï¿½WStk r2tï¿½|ï¿½YSXdSr))rrr?r@)r7rï¿½rï¿½r&r&r'r &s cCs|d}t|ï¿½dkrï¿½d||}d|d|d|d|d|d |d |d|}d|d |d|d|d|d|d|d}||}|||S|dkrï¿½|nd|}tt|ï¿½ï¿½}|dkï¿½r^|d}d|d|d|d|d|d|d|d}d|d |d!|d"|d#|d$|d%|d}nï¿½|d}d&|d'|d(|d)|d*|d+|d,|d-}d.|d/|d0|d1|d2|d3|d4|d}||}|dkï¿½rï¿½|}|||S)5Nï¿½ï¿½?g333333ï¿½?gï¿½ï¿½Qï¿½ï¿½?g^ï¿½}o)ï¿½ï¿½@gï¿½E.kï¿½Rï¿½@g ï¿½ï¿½Ulï¿½@g*uï¿½ï¿½>lï¿½@gï¿½Nï¿½ï¿½ï¿½ï¿½@gï¿½ï¿½"]Îž@gnCï¿½ï¿½ï¿½`@guï¿½ï¿½@giKï¿½ï¿½~jï¿½@gvï¿½ï¿½|Eï¿½@gï¿½ï¿½dï¿½|1ï¿½@gfRï¿½ï¿½rï¿½ï¿½@gï¿½ï¿½u.2ï¿½@gï¿½ï¿½ï¿½ï¿½~yï¿½@gï¿½n8(E@ï¿½ï¿½?ï¿½g@gï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½?gé¬·ï¿½ZaI?ggï¿½Elï¿½Dï¿½?g7\ï¿½ï¿½ï¿½ï¿½ï¿½?gï¿½uSï¿½Sï¿½?gï¿½=ï¿½. @gj%bï¿½@gï¿½ï¿½ï¿½Hwï¿½@gjRï¿½ï¿½eï¿½ï¿½?gï¿½9dh? >g('ß¿ï¿½ï¿½A?gï¿½ï¿½~z ï¿½?g@ï¿½3ï¿½ï¿½?gÉ…3ï¿½ï¿½?g3fRï¿½xï¿½ï¿½?gIï¿½Fï¿½ï¿½l@gï¿½ï¿½ï¿½ï¿½tï¿½ï¿½>g*ï¿½Yï¿½ï¿½nï¿½>gESB\T?gï¿½N;A+ï¿½?gï¿½UR1ï¿½ï¿½ï¿½?gEï¿½Fï¿½ï¿½ï¿½?gPï¿½nï¿½ï¿½@g&ï¿½>ï¿½ï¿½ï¿½@gï¿½ï¿½ï¿½ï¿½iï¿½gï¿½tcI,\ï¿½>gï¿½Åï¿½ï¿½ï¿½I?g*F2ï¿½vï¿½?gï¿½C4ï¿½ï¿½?gï¿½ï¿½Oï¿½1ï¿½?)rrr)ï¿½prï¿½ï¿½sigmaï¿½qï¿½rZnumZdenrAr&r&r'ï¿½_normal_dist_inv_cdf9sdï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ rï¿½c@sï¿½eZdZdddï¿½Zd7ddï¿½Zedd ï¿½ï¿½Zd dï¿½dd ï¿½Zddï¿½Zddï¿½Z ddï¿½Z d8ddï¿½Zddï¿½Ze ddï¿½ï¿½Ze ddï¿½ï¿½Ze ddï¿½ï¿½Ze dd ï¿½ï¿½Ze d!d"ï¿½ï¿½Zd#d$ï¿½Zd%d&ï¿½Zd'd(ï¿½Zd)d*ï¿½Zd+d,ï¿½Zd-d.ï¿½ZeZd/d0ï¿½ZeZd1d2ï¿½Zd3d4ï¿½Zd5d6ï¿½Zd S)9rz(Arithmetic mean of a normal distributionz+Standard deviation of a normal distribution)ï¿½_muï¿½_sigmarï¿½rï¿½cCs(|dkrtdï¿½ï¿½t|ï¿½|_t|ï¿½|_dS)Nrï¿½zsigma must be non-negative)rrErï¿½rï¿½)ï¿½selfrï¿½rï¿½r&r&r'ï¿½__init__ï¿½s zNormalDist.__init__cCs.t|ttfï¿½st|ï¿½}t|ï¿½}||t||ï¿½ï¿½Sr))rbr\ï¿½tuplerr)ï¿½clsr7rï¿½r&r&r'ï¿½from_samplesï¿½szNormalDist.from_samplesN)ï¿½seedcsB|dkrtjn tï¿½|ï¿½jï¿½|j|jï¿½ï¿½ï¿½ï¿½ï¿½fddï¿½t|ï¿½Dï¿½S)Ncsg|]}ï¿½ï¿½ï¿½ï¿½ï¿½qSr&r&ï¿½r*rSï¿½ï¿½gaussrï¿½rï¿½r&r'ï¿½ ï¿½sz&NormalDist.samples..)ï¿½randomrï¿½ZRandomrï¿½rï¿½rw)rï¿½r,rï¿½r&rï¿½r'ï¿½samplesï¿½szNormalDist.samplescCs<|jd}|stdï¿½ï¿½t||jdd|ï¿½tt|ï¿½S)Nï¿½@z$pdf() not defined when sigma is zerogï¿½)rï¿½rrrï¿½rr)rï¿½rArr&r&r'ï¿½pdfï¿½s zNormalDist.pdfcCs2|jstdï¿½ï¿½ddt||j|jtdï¿½ï¿½S)Nz$cdf() not defined when sigma is zerorï¿½rï¿½rï¿½)rï¿½rrrï¿½r)rï¿½rAr&r&r'ï¿½cdfï¿½szNormalDist.cdfcCs:|dks|dkrtdï¿½ï¿½|jdkr*tdï¿½ï¿½t||j|jï¿½S)Nrï¿½rï¿½z$p must be in the range 0.0 < p < 1.0z-cdf() not defined when sigma at or below zero)rrï¿½rï¿½rï¿½)rï¿½rï¿½r&r&r'ï¿½inv_cdfï¿½s zNormalDist.inv_cdfrtcsï¿½ï¿½fddï¿½tdï¿½ï¿½Dï¿½S)Ncsg|]}ï¿½ï¿½|ï¿½ï¿½ï¿½qSr&)rï¿½rï¿½ï¿½r,rï¿½r&r'rï¿½ï¿½sz(NormalDist.quantiles..r()rw)rï¿½r,r&rï¿½r'rï¿½s zNormalDist.quantilescCst|tï¿½stdï¿½ï¿½||}}|j|jf|j|jfkr>||}}|j|j}}|rT|s\tdï¿½ï¿½||}t|j|jï¿½}|sï¿½dt|d|jt dï¿½ï¿½S|j||j|}|j|jt |d|t ||ï¿½ï¿½} || |} || |}dt|ï¿½| ï¿½|ï¿½| ï¿½ï¿½t|ï¿½|ï¿½|ï¿½|ï¿½ï¿½S)Nz$Expected another NormalDist instancez(overlap() not defined when sigma is zerorï¿½rï¿½)rbrrFrï¿½rï¿½rrrrrrrï¿½)rï¿½ï¿½otherï¿½Xï¿½YZX_varZY_varZdvZdmrRï¿½bï¿½x1ï¿½x2r&r&r'ï¿½overlapï¿½s" (zNormalDist.overlapcCs|jSr)ï¿½rï¿½ï¿½rï¿½r&r&r'rï¿½szNormalDist.meancCs|jSr)rï¿½rï¿½r&r&r'rï¿½szNormalDist.mediancCs|jSr)rï¿½rï¿½r&r&r'rï¿½szNormalDist.modecCs|jSr)ï¿½rï¿½rï¿½r&r&r'rï¿½szNormalDist.stdevcCs |jdS)Nrï¿½rï¿½rï¿½r&r&r'rszNormalDist.variancecCs8t|tï¿½r&t|j|jt|j|jï¿½ï¿½St|j||jï¿½Sr)ï¿½rbrrï¿½rrï¿½ï¿½rï¿½rï¿½r&r&r'ï¿½__add__ s zNormalDist.__add__cCs8t|tï¿½r&t|j|jt|j|jï¿½ï¿½St|j||jï¿½Sr)rï¿½rï¿½r&r&r'ï¿½__sub__s zNormalDist.__sub__cCst|j||jt|ï¿½ï¿½Sr)ï¿½rrï¿½rï¿½rrï¿½r&r&r'ï¿½__mul__%szNormalDist.__mul__cCst|j||jt|ï¿½ï¿½Sr)rï¿½rï¿½r&r&r'ï¿½__truediv__-szNormalDist.__truediv__cCst|j|jï¿½Sr)ï¿½rrï¿½rï¿½ï¿½rï¿½r&r&r'ï¿½__pos__5szNormalDist.__pos__cCst|j|jï¿½Sr)rï¿½rï¿½r&r&r'ï¿½__neg__9szNormalDist.__neg__cCs ||Sr)r&rï¿½r&r&r'ï¿½__rsub__?szNormalDist.__rsub__cCs&t|tï¿½stS|j|jko$|j|jkSr))rbrï¿½NotImplementedrï¿½rï¿½rï¿½r&r&r'ï¿½__eq__Es zNormalDist.__eq__cCst|j|jfï¿½Sr))ï¿½hashrï¿½rï¿½rï¿½r&r&r'ï¿½__hash__KszNormalDist.__hash__cCs t|ï¿½jï¿½d|jï¿½d|jï¿½dï¿½S)Nz(mu=z, sigma=ï¿½))r2r#rï¿½rï¿½rï¿½r&r&r'ï¿½__repr__OszNormalDist.__repr__)rï¿½rï¿½)rt)r#r$r%ï¿½ __slots__rï¿½ï¿½classmethodrï¿½rï¿½rï¿½rï¿½rï¿½rrï¿½ï¿½propertyrrrrrrï¿½rï¿½rï¿½rï¿½rï¿½rï¿½ï¿½__radd__rï¿½ï¿½__rmul__rï¿½rï¿½rï¿½r&r&r&r'rï¿½sDï¿½ " )rï¿½ï¿½__main__)ï¿½isclose)ï¿½addï¿½subï¿½mulï¿½truediv)ï¿½repeatï¿½ ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½ï¿½iï¿½ï¿½z Test z with another NormalDist:ï¿½z with a constant:ï¿½z Test constant with ï¿½:cCsdSr)r&)ï¿½G1ï¿½G2r&r&r'ï¿½assert_closesrï¿½iï¿½ï¿½ï¿½ï¿½ï¿½Iï¿½ï¿½/g`@@cCsg|]}|tï¿½qSr&ï¿½ï¿½srar&r&r'rï¿½ï¿½srï¿½cCsg|]}|tï¿½qSr&rï¿½rar&r&r'rï¿½ï¿½scCsg|]}|tï¿½qSr&rï¿½rar&r&r'rï¿½ï¿½scCsg|]}|tï¿½qSr&rï¿½rar&r&r'rï¿½ï¿½scCsg|]\}}||ï¿½qSr&r&ï¿½r*rAï¿½yr&r&r'rï¿½ï¿½scCsg|]\}}||ï¿½qSr&r&rï¿½r&r&r'rï¿½ï¿½s)r)rX)r()N)N)N)N)N)Qï¿½__all__r@rcrï¿½Z fractionsrZdecimalrï¿½ itertoolsrZbisectrrrrrrrrrr ï¿½operatorr!ï¿½collectionsr"rMrr>rBr0r.rPrTrWrZrrrrrr r rrrrrï¿½rrrr rï¿½rZ_statisticsï¿½ImportErrorr#rï¿½rï¿½rï¿½rï¿½rï¿½rï¿½ZdoctestZg1Zg2r,rï¿½rï¿½rï¿½ï¿½funcï¿½printrï¿½r3Zconstrï¿½rï¿½rï¿½rï¿½rGï¿½zipZtestmodr&r&r&r'ï¿½Tsï¿½ï¿½( : / 779 / , JQ ï¿½ ï¿½