python (3.12.0)
ˑe3 U g d Z ddlZddlZddlZddlZddlmZ ddlmZ ddl m
Z
m Z m Z ddl
mZmZ ddlmZmZmZmZmZmZmZmZmZ ddlmZ dd lm Z dd
l m Z m Z m Z ed Z! G d d
e" Z#d Z$dBdZ%d Z&d Z'd Z(d Z)dCdZ*dddddde+e, fdZ-de.d e.de.fd Z/d ej` jb z d z Z2e.e3d < de.d e.de,fd!Z4de.d e.defd"Z5d# Z6dBd$Z7d% Z8dBd&Z9d' Z:d( Z;d) Z<dDd*Z=d+ Z>d, Z?d-d.d/d0Z@dBd1ZAdBd2ZBdBd3ZCdBd4ZDd5 ZEd6 ZFd7d8d9ZG e d:d; ZHdd<d=ZId> ZJ dd?lKmJZJ G d@ dA ZMy# eL$ r Y w xY w)E)
NormalDistStatisticsError correlation
covariancefmeangeometric_mean
harmonic_meanlinear_regressionmeanmedianmedian_grouped median_high
median_lowmode multimodepstdev pvariance quantilesstdevvariance NFraction)Decimal)countgroupbyrepeat) bisect_left bisect_right) hypotsqrtfabsexperftaulogfsumsumprod)reduce)
itemgetter)Counter
namedtuple defaultdict @c e Zd Zy)r N)__name__
__module__ __qualname__ >/BuggyBox/python/3.12.0/bootstrap/lib/python3.12/statistics.pyr r s r3 r c p d}t }|j }i }|j }t | t D ]9 \ }} || t t | D ] \ }} |dz
} || d |z || < ; d |v r|d }
n t d |j D }
t t |t } | |
|fS )Nr c 3 : K | ] \ }}t || y wNr .0dns r4 <genexpr>z_sum.<locals>.<genexpr> s @tq!HQN@ ) setaddgetr typemap _exact_ratiosumitemsr( _coerceint) datar types types_addpartials partials_gettypvaluesr< r; totalTs r4 _sumrR s >
E EE IH<<L tT* 1 V# f- 1DAq QJE&q!,q0HQK 11
x @x~~/?@@ws#A
ue r3 c ~ t fd| D \ }}}|||fS d}t }|j }t t }t t }t | t D ]G \ } }
|| t t |
D ]' \ } |dz
}|xx | z
cc< |xx | | z z
cc< ) I |st d x}n_d |v r|d x}nSt d |j D } t d |j D }
||
z | | z z
|z }| |z t t |t }|||fS )Nc 3 2 K | ] }|z
xz y wr8 r2 )r: xcr; s r4 r= z_ss.<locals>.<genexpr> s <!1q5jaA - <s r r6 c 3 : K | ] \ }}t || y wr8 r r9 s r4 r= z_ss.<locals>.<genexpr> s @DAq!Q @r> c 3 @ K | ] \ }}t |||z y wr8 r r9 s r4 r= z_ss.<locals>.<genexpr> s Dtq!(1ac"Ds )rR r? r@ r, rH r rB rC rD r rE rF r( rG )rI rV rQ ssdr rJ rK sx_partials sxx_partialsrN rO r< sxsxxr; s ` @r4 _ssr^ s` } <t <<
335!!
E EE I c"K s#L tT* % V# f- %DAq QJE Na N Oq1u $O %% 1+ a
d##a @K,=,=,?@
@D|/A/A/CDD s{R"W$- Jws#A
sAu r3 c l | j S # t $ r t j | cY S w xY wr8 ) is_finiteAttributeErrormathisfinite)rU s r4 _isfiniterd s1 {{} }}Q s 33c | |u r| S |t u s|t u r| S | t u r|S t || r|S t | | r| S t | t r|S t |t r| S t | t rt |t r|S t | t rt |t r| S d}t || j |j fz )Nz"don't know how to coerce %s and %s)rH bool
issubclassr float TypeErrorr/ )rQ Smsgs r4 rG rG s AvqCx19axCx(!Q(!Q(!S1H!S1H!X :a#7!U
1h 7
.C
C1::qzz22
33r3 c | j S # t $ r Y nt t f$ r | d fcY S w xY w | j | j
fS # t $ r% dt
| j d}t | w xY w)Nzcan't convert type 'z' to numerator/denominator) as_integer_ratiora
OverflowError
ValueError numerator denominatorrB r/ ri )rU rk s r4 rD rD s :!!##
: & 4y Q]]++ $T!W%5%5$66PQn s 333A .A=c t | |u r| S t |t r| j dk7 rt } || S # t
$ r9 t |t r' || j || j z cY S w xY wNr6 )rB rg rH rq rh ri r rp )valuerQ s r4 _convertru M s} 0 E{a !Se / /1 4 x a !U__%%*;*;(<< < s ? >B?Bc # D K | D ] }|dk r t | | y w)Nr )r )rO errmsgrU s r4 _fail_negrx _ s/ A
q5!&) )s Faverager6 )keyreversetiesstartreturnc V |dk7 rt d| | t || } t t | t | }|dz
}dgt | z }t
|t d D ]: \ }} t | }
t |
} || dz dz z } |
D ]
\ }
}| ||< || z
}< |S )Nry z Unknown tie resolution method: )r{ r6 r )rz ) ro rC sortedzipr lenr r) list)rI rz r{ r| r} val_posiresult_ggroupsizerankrt orig_poss r4 _rankr g s H y:4(CDD
3~Suw'9G
AS3w<
FZ]3 1Q5zD1H>! $ $OE8 #F8 $ T
Mr3 r< mc P t j | |z }|||z |z | k7 z S r8 )rb isqrt)r< r as r4 _integer_sqrt_of_frac_rtor s. P
16A !A
r3 r _sqrt_bit_widthc | j |j z
t z
dz }|dk\ r t | |d|z z |z }d}||z S t | d|z z | }d| z }||z S )Nr r r6 )
bit_lengthr r )r< r qrp rq s r4 _float_sqrt_of_fracr s ;
!,,. (? :q@AAv-aa!e<A { "" .a26k1= A2g { ""r3 c | dk r| s t d S | | }} t | t | z j }|j \ }}|j }|j \ }}d| z ||z dz z |||z ||z z dz z kD r|S |j }|j \ } }
d| z ||
z dz z ||| z |
|z z dz z k r|S |S )Nr z0.0 r )r r rm next_plus
next_minus) r< r rootnrdrplusnpdpminusnmdms r4 _decimal_sqrt_of_fracr s = Av5> !rA21 AJ # ) ) +D
"
"
$FB >> D
"
"
$FB1u2zABB 2 22 OO E
#
#
%FB1u2zABB 2 22 Kr3 c ` t | \ }}}|dk r t d t ||z | S )Nr6 z%mean requires at least one data point)rR r ru )rI rQ rP r< s r4 r
r
s< t*KAua1u EFF EAIq !!r3 c ^ |) t | t | }s t d |z S t |t
t f s t | } t | | }t | }|s t d ||z S # t $ r dfd} || } Y w xY w# t $ r t d w xY w)Nr c 3 @ K t | d D ] \ }| y w)Nr6 r} ) enumerate)iterablerU r< s r4 r zfmean.<locals>.count s& %ha 8 DAqG s z&fmean requires at least one data pointz(data and weights must be the same lengthz sum of weights must be non-zero) r ri r& r
isinstancer tupler' ro )rI weightsr rP numdenr< s @r4 r r s D A T
!"JK Kqy ge} -w-JdG$ w-C ?@@ 9+ A
;D J HIIJs A9 B 9BBB,c | t t t t | S # t $ r
t d d w xY w)NzGgeometric mean requires a non-empty dataset containing positive numbers)r" r rC r% ro r )rI s r4 r r sJ
G5S$ ()) G < =BF GGs !% ;c z t | | u r t | } d}t | }|dk r t d |dk( rD|B| d }t |t
j t f r|dk r t | |S t d |t d| }|}nQt | |u r t | }t | |k7 r t d t d t || D \ }}} t | | } t d t || D \ }}} |dk r t d t ||z | S # t $ r Y yw xY w)
Nz.harmonic mean does not support negative valuesr6 z.harmonic_mean requires at least one data pointr zunsupported typez*Number of weights does not match data sizec 3 K | ] }| y wr8 r2 )r: ws r4 r= z harmonic_mean.<locals>.<genexpr>N s Gq Gs c 3 4 K | ] \ }}|r||z nd yw)r Nr2 )r: r rU s r4 r= z harmonic_mean.<locals>.<genexpr>Q s PTQq1uq 0 Ps z Weighted sum must be positive)iterr r r r numbersRealr ri r rR rx r ZeroDivisionErrorru )
rI r rw r< rU sum_weightsr rQ rP r s
r4 r r ! sR , DzTDz
=F D A1u NOO
aGO G a',,0 11u%f--H ./ /A, =G #7mG w<1 !"NO O GIgv,F G G;v& PS$=O PP5%
z =>> K%' ++ s #,D. . D:9D:c t | } t | }|dk( r t d |dz dk( r| |dz S |dz }| |dz
| | z dz S Nr no median for empty datar r6 r r r )rI r< r s r4 r r Y sl $<D D AAv 8991uzAF|
FQU d1g%**r3 c t | } t | }|dk( r t d |dz dk( r| |dz S | |dz dz
S r r rI r< s r4 r r q sZ
$<D D AAv 8991uzAF|AFQJ r3 c ` t | } t | }|dk( r t d | |dz S )Nr r r r r s r4 r
r
s<
$<D D AAv 899 Q<r3 c , t | } t | }|s t d | |dz }t | | }t | || } t | }t | }||dz z
}|}||z
}|||dz |z
z |z z S # t $ r t d w xY w)Nr r )loz$Value cannot be converted to a floatr- )r r r r r rh ro ri ) rI intervalr< rU r jLcffs r4 r r s )T $<D D A 899
Q!V A D! AT1#AA? !H
HsNA
B AA x1q52:&* ** A>@@As A> >Bc t t | j d } |d d S # t $ r
t d d w xY w)Nr6 r zno mode for empty data)r* r most_common
IndexErrorr )rI pairss r4 r r sU ,
DJ + +A .EBQx{ B 67TABs . Ac t t | }|sg S t |j }|j D cg c]
\ }}||k( s | c}}S c c}}w r8 )r* r maxrO rF )rI countsmaxcountrt r s r4 r r sT T$Z
F 6==?#H&,lln JleU8IE JJ Js
A A r exclusive)r< methodc * |dk r t d t | } t | }|dk r t d |dk( rW|dz
}g }t d| D ]? }t ||z | \ }}| | ||z
z | |dz |z z |z } |j | A |S |dk( rn|dz }g }t d| D ]V }||z |z }|dk rdn||dz
kD r|dz
n|}||z ||z z
}| |dz
||z
z | | |z z |z } |j | X |S t
d| )Nr6 zn must be at least 1r z"must have at least two data points inclusiver Unknown method: )r r r rangedivmodappendro )
rI r< r ldr r r r delta interpolateds
r4 r r 9 sp 1u 455 $<D TB Av BCC
Fq! (A a!eQ'HAu Gq5y 1DQK%4G G1LL MM, ' (
Fq! (AA
AUBqD1aAaC!A#IE QK1u9 5Q% G1LL MM, ' (
'z2
33r3 c j t | | \ }}}}|dk r t d t ||dz
z | S )Nr z*variance requires at least two data pointsr6 r^ r ru )rI xbarrQ ssrV r< s r4 r r j sE $J dD/KAr1a1u JKK B!a%L! $$r3 c d t | | \ }}}}|dk r t d t ||z | S )Nr6 z*pvariance requires at least one data pointr )rI murQ r rV r< s r4 r r sA !D dB-KAr1a1u JKK BFA r3 c t | | \ }}}}|dk r t d ||dz
z }t |t r t |j
|j S t |j
|j S Nr z'stdev requires at least two data pointsr6 r^ r rg r r rp rq r )rI r rQ r rV r< msss r4 r r sp dD/KAr1a1u GHH
A,C!W $S]]COODD s}}coo >>r3 c t | | \ }}}}|dk r t d ||z }t |t r t |j
|j S t |j
|j S )Nr6 z'pstdev requires at least one data pointr )rI r rQ r rV r< r s r4 r r sl dB-KAr1a1u GHH
q&C!W $S]]COODD s}}coo >>r3 c t | \ }}}}|dk r t d ||dz
z } t | t |j |j
fS # t $ r% t | t | t | z fcY S w xY wr )r^ r rh r rp rq ra )rI rQ r r r< r s r4 _mean_stdevr s PYNAr41u GHH
A,C4T{ /
s OOO 4T{E$K%) 3334s *A +BBc t | }t | |k7 r t d |dk r t d t | |z t | |z t fd| D fd|D }||dz
z S )NzDcovariance requires that both inputs have same number of data pointsr z,covariance requires at least two data pointsc 3 ( K | ] }|z
y wr8 r2 )r: xir s r4 r= z covariance.<locals>.<genexpr> s )29) c 3 ( K | ] }|z
y wr8 r2 r: yiybars r4 r= z covariance.<locals>.<genexpr> s +B"BI+Br r6 )r r r& r' )rU yr< sxyr r s @@r4 r r s} AA
1v{ dee1u LMM 7Q;D 7Q;D
)q)+B+B
CC !a%=r3 linear)r c t | }t | |k7 r t d |dk r t d |dvrt d| |dk( r#|dz
dz }t | | } t || }n@t | |z }t | |z }| D cg c] }||z
} }|D cg c] }||z
}}t | | } t | | }
t || } | t
|
| z z S c c}w c c}w # t $ r t d
w xY w) NzEcorrelation requires that both inputs have same number of data pointsr z-correlation requires at least two data points> r rankedr r r6 r r z&at least one of the inputs is constant)r r ro r r& r' r r ) rU r r r< r} r r r r r r] syys r4 r r s. , AA
1v{ eff1u MNN
))+F:677
Q" !5 ! !5 !Aw{Aw{!" #2R$Y # #!" #2R$Y # #
!Q-C
!Q-C
!Q-CHT#)_$$
$ # H FGGHs C&" C+C0 0DLinearRegressionslope intercept) proportionalc
t | }t | |k7 r t d |dk r t d |s9t | |z }t | |z
| D cg c] }||z
} }
fd|D }t | | dz }t | | } ||z }|rdn
|z z
} t || S c c}w # t $ r t d w xY w)NzKlinear regression requires that both inputs have same number of data pointsr z3linear regression requires at least two data pointsc 3 ( K | ] }|z
y wr8 r2 r s r4 r= z$linear_regression.<locals>.<genexpr>u s #2R$Y #r z
x is constantr )r r r& r' r r ) rU r r r< r r r r] r r r s @r4 r r F s $J AA
1v{ kll1u STT Aw{Aw{!" #2R$Y # # # #
!Q-#
C
!Q-C/c $ )<I %9 ==
$ / o../s B,B1 1Cc | dz
}t | dk rpd||z z
}d|z dz |z dz |z dz |z dz |z d z |z d
z |z d z |z }d |z d
z |z dz |z dz |z dz |z dz |z dz |z dz }||z }|||z z S |dk r| nd| z
}t t | }|dk r^|dz
}d|z dz |z dz |z dz |z dz |z d z |z d z |z d z }d |z d z |z d!z |z d"z |z d#z |z d$z |z d%z |z dz }n]|dz
}d&|z d'z |z d(z |z d)z |z d*z |z d+z |z d,z |z d-z }d.|z d/z |z d0z |z d1z |z d2z |z d3z |z d4z |z dz }||z }|dk r| }|||z z S )5N ?g333333?g Q ?g^}o)@gE.kR@g U l@g*u>l@g N@g"]Ξ@gnC`@gu @giK~j@gv|E@gd|1@gfRr@gu.2@g~y@gn8(E@ ?r g @g?g鬷ZaI?ggElD?g7\?guSS?g=.
@gj%b@gHw@gjRe?g9dh?
>g('߿A?g~z ?g@3?gɅ3?g3fRx?gIFl @gt>g*Yn>gESB\T?gN;A+?gUR1?gE F?gPn@g&>@gi<g@F>gtcI,\>gŝI?g*F2v?gC4?gO1?)r! r r% )pr sigmar rr r rU s r4 _normal_dist_inv_cdfr s
CA Aw% q1u 0140145601456 11 566 1 1 56 6
1 1
56 6 1
1 56
6 11 566 1140145601456 11 566 1 1 56 6
1 1
56 6 1
1 56
6
#IQY
#X37A c!fW
ACx
G1A51256712567 22 677 2 2 67 7
2 2
67 7 2
2 67
7 22 2A51256712567 22 677 2 2 67 7
2 2
67 7 2
2 67
7
G1A51256712567 22 677 2 2 67 7
2 2
67 7 2
2 67
7 22 3Q61256712567 22 677 2 2 67 7
2 2
67 7 2
2 67
7 c A3w
B
U r3 )r c e Zd Z dddZd dZed ZdddZd Zd
Z d Z
d!d Z d
Z d Z
ed Zed Zed Zed Zed Zd Zd Zd Zd Zd Zd ZeZd ZeZ d Z d Z d Z d Z d Z!y)"r z(Arithmetic mean of a normal distributionz+Standard deviation of a normal distribution_mu_sigmac f |dk r t d t | | _ t | | _ y )Nr zsigma must be non-negative)r rh r r )selfr r s r4 __init__zNormalDist.__init__ s. N 3;!">? ?9El r3 c | t | S r8 )r )clsrI s r4 from_sampleszNormalDist.from_samples s ?K%&&r3 N)seedc |t j n t j | j }| j | j }}t d | D cg c] } |||
c}S c c}w r8 )randomgaussRandomr r r )r r< r r r r r s r4 sampleszNormalDist.samples sU G $ &--2E2K2KHHdkkE*0q/:Qb% :::s A,c | j | j z }|s t d || j z
}t ||z d|z z t t
|z z S )Nz$pdf() not defined when sigma is zerog )r r r r" r r$ )r rU r diffs r4 pdfzNormalDist.pdf sY >;;,!"HI I488|4$;$/23d3>6JJJr3 c | j s t d ddt || j z
| j t z z z z S )Nz$cdf() not defined when sigma is zeror r )r r r# r _SQRT2r rU s r4 cdfzNormalDist.cdf sC 6{{!"HI IcCTXX $++2F G HHIIr3 c p |dk s|dk\ r t d t || j | j S )Nr r z$p must be in the range 0.0 < p < 1.0)r r r r )r r s r4 inv_cdfzNormalDist.inv_cdf s9
8qCx!"HI I#Atxx==r3 c f t d| D cg c] }| j ||z c}S c c}w rs )r r! )r r< r s r4 r zNormalDist.quantiles s0 .31a[9 QU#999s .c t |t s t d | |}}|j |j f|j |j fk r||}}|j
|j
}}|r|s t
d ||z
}t |j |j z
}|s%dt |d|j z t z z z
S |j |z |j |z z
}|j |j z t ||z |t ||z z z z } || z |z }
|| z
|z } dt |j |
|j |
z
t |j | |j | z
z z
S )Nz$Expected another NormalDist instancez(overlap() not defined when sigma is zeror r- )
r r ri r r r r r! r# r r r% r ) r otherXYX_varY_vardvr r bx1x2s r4 overlapzNormalDist.overlap sd
%, BC CU1
HHaee !%% 0 0aqAzz1::uE!"LM M
U]
!%%!%%-
R3>F#: ;<< <
EEEMAEEEM )
HHqxx $rBwc%%-6H1H'H"I I!er\!er\d1559quuRy01DrQUU2Y9N4OOPPr3 c j | j s t d || j z
| j z S )Nz'zscore() not defined when sigma is zero)r r r r s r4 zscorezNormalDist.zscore9 s3 {{!"KL LDHH ++r3 c | j S r8 r r s r4 r
zNormalDist.meanD s 5xxr3 c | j S r8 r1 r2 s r4 r zNormalDist.medianI s 6xxr3 c | j S r8 r1 r2 s r4 r zNormalDist.modeN s
xxr3 c | j S r8 r r2 s r4 r zNormalDist.stdevW s 8{{r3 c 6 | j | j z S r8 r6 r2 s r4 r zNormalDist.variance\ s +{{T[[((r3 c t |t rAt | j |j z t | j |j S t | j |z | j S r8 r r r r r r+ r, s r4 __add__zNormalDist.__add__a T b* % bffrvvouRYY /JK K"&&2+ryy11r3 c t |t rAt | j |j z
t | j |j S t | j |z
| j S r8 r9 r: s r4 __sub__zNormalDist.__sub__o r< r3 c b t | j |z | j t | z S r8 r r r r! r: s r4 __mul__zNormalDist.__mul__} +
"&&2+ryy48';<<r3 c b t | j |z | j t | z S r8 r@ r: s r4 __truediv__zNormalDist.__truediv__ rB r3 c D t | j | j S r8 r r r r+ s r4 __pos__zNormalDist.__pos__ s ("&&")),,r3 c F t | j | j S r8 rF rG s r4 __neg__zNormalDist.__neg__ s 2266'299--r3 c | |z
S r8 r2 r: s r4 __rsub__zNormalDist.__rsub__ s Fbzr3 c t |t st S | j |j k( xr | j |j k( S r8 )r r NotImplementedr r r: s r4 __eq__zNormalDist.__eq__ s: P"j)! !vv :BII$::r3 c F t | j | j f S r8 )hashr r r2 s r4 __hash__zNormalDist.__hash__ s MTXXt{{+,,r3 c f t | j d| j d| j dS )Nz(mu=z, sigma=))rB r/ r r r2 s r4 __repr__zNormalDist.__repr__ s. t*%%&d488,ht{{oQOOr3 c 2 | j | j fS r8 r
r2 s r4 __getstate__zNormalDist.__getstate__ s xx$$r3 c " |\ | _ | _ y r8 r
)r states r4 __setstate__zNormalDist.__setstate__ s % $+r3 )r r )r )"r/ r0 r1 __slots__r classmethodr r r r r! r r- r/ propertyr
r r r r r; r> rA rD rH rJ __radd__rL __rmul__rO rR rU rW rZ r2 r3 r4 r r s .
:?I
# ' ' "& ; KJ > : QD , ) ) 2 2 ==-. H H; -P%&r3 r r8 )znegative value)r )N__all__rb r r sys fractionsr decimalr itertoolsr r r bisectr r r r r! r" r# r$ r% r& r' functoolsr( operatorr) collectionsr* r+ r, r ro r rR r^ rd rG rD ru rx r rh r rH r
float_infomant_digr __annotations__r r r
r r r r r r
r r r r r r r r r r r r r r _statistics ImportErrorr r2 r3 r4 <module>rn s hT .
, , , E E E 8 8
c j 3 l& R 4>+ \$ IQ 34PU; 3l 3>>222Q6 6
#3
#3
#5
#S S W <",!HG&5,p+0 ,&E+PB<Kr ; (4b)%X& R?$?$
4(8 $, -H` 02HI 05 7>zG V 0
Z& Z& s /E EE
Browse
Usage
Versions