Warning: Cannot modify header information - headers already sent by (output started at /var/www/jonas-eriksen.dk/pic/private/index.php:1) in /var/www/jonas-eriksen.dk/pic/private/index.php on line 215
!A?hY r U g d Z ddlZddlZddlZddlZddlmZ ddlmZ ddl m
Z
m Z ddl m
Z
mZ ddlmZmZmZmZmZmZmZmZ ddlmZ dd lmZ dd
lm Z m Z m Z ed Z G d d
e Z!d Z"d>dZ#d Z$d Z%d Z&d Z'd?dZ(de)de)de)fdZ*dej+ j, z dz Z-e)e.d < de)de)de/fd Z0de)de)defd Z1d Z2d>d Z3d! Z4d>d"Z5d# Z6d$ Z7d% Z8d@d'Z9d( Z:d) Z;d*d+d,d-Zd.Z=d>d/Z>d>d0Z?d>d1Z@d2 ZAd3 ZBd4 ZC e d5d6 ZDd7d8d9ZEd: ZF dd;lGmFZF n # eH$ r Y nw xY w G d< d= ZIdS )A)
NormalDistStatisticsError correlation
covariancefmeangeometric_mean
harmonic_meanlinear_regressionmeanmedianmedian_grouped median_high
median_lowmode multimodepstdev pvariance quantilesstdevvariance NFraction)Decimal)groupbyrepeat) bisect_left bisect_right)hypotsqrtfabsexperftaulogfsum)reduce)mul)Counter
namedtuple defaultdict @c e Zd ZdS )r N)__name__
__module__ __qualname__ #/usr/lib64/python3.11/statistics.pyr r s Dr1 r c d}t }|j }i }|j }t | t D ]B\ }} || t t | D ] \ }} |dz
} || d |z || < Cd |v r |d }
n+t d | D }
t t |t } | |
|fS )Nr c 3 < K | ]\ }}t || V d S Nr .0dns r2 z_sum.. s. @@tq!HQNN@@@@@@r1 ) setaddgetr typemap _exact_ratiosumitemsr& _coerceint) datacounttypes types_addpartials partials_gettypvaluesr: r9 totalTs r2 _sumrP s >
E EEE IH. s0 < . s. @@DAq!Q @@@@@@r1 c 3 B K | ]\ }}t |||z V d S r6 r r7 s r2 r; z_ss.. s4 DDtq!(1ac""DDDDDDr1 )rP r< r= r* rE r r? r@ rA r rB rC r& rD )rF rT rO ssdrG rH rI sx_partials sxx_partialsrL rM r: sxsxxr9 s ` @r2 _ssr\ s } < < < < c t | |u r| S t |t r| j dk rt } || S # t
$ r: t |t r# || j || j z cY S w xY w)Nr4 )r? re rE ro rf rg r rn )valuerO s r2 _convertrr M s 0 E{{a !S e /1 4 4 qxx a ! ! 1U_%%%*;(<(<< < < < s
A A B
B negative valuec # H K | D ]}|dk rt | |V d S )Nr )r )rM errmsgrS s r2 _fail_negrv _ sD A
q55!&)) ) r1 r: mreturnc P t j | |z }|||z |z | k z S r6 )r` isqrt)r: rw as r2 _integer_sqrt_of_frac_rtor| g s/ P
16A !A
r1 _sqrt_bit_widthc | | z
t z
dz }|dk r t | |d|z z |z }d}n t | d|z z | }d| z }||z S )Nr} r r4 )
bit_lengthr r| )r: rw qrn ro s r2 _float_sqrt_of_fracr s s ;
!,,.. (? :q@AAvv-aa!e<<A -a26k1== A2g { ""r1 c
| dk r| st d S | | }} t | t | z }| \ }}| }| \ }}d| z ||z dz z |||z ||z z dz z k r|S | }| \ } }
d| z ||
z dz z ||| z |
|z z dz z k r|S |S )Nr z0.0 r} )r r rk next_plus
next_minus) r: rw rootnrdrplusnpdpminusnmdms r2 _decimal_sqrt_of_fracr s# = Avv "5>> !rA21 AJJ # ) ) + +D
"
"
$
$FB >> D
"
"
$
$FB1u2zABB 2 222 OO E
#
#
%
%FB1u2zABB 2 222 Kr1 c z t | \ }}}|dk rt d t ||z | S )Nr4 z%mean requires at least one data point)rP r rr )rF rO rN r: s r2 r
r
sF t**KAua1uu EFFF EAIq ! !!r1 c t | n"# t $ r dfd} || } Y nw xY w|%t | }st d |z S t | }n.# t $ r! t | }t | }Y nw xY wt t t | | }|k rt d t | }|st d ||z S )Nr c 3 B K t | d D ] \ }|V
d S )Nr4 )start) enumerate)iterablerS r: s r2 rG zfmean..count s< !(!444
1
r1 z&fmean requires at least one data pointz(data and weights must be the same lengthz sum of weights must be non-zero)lenrg r% r listr@ r' )rF weightsrG rN num_weightsnumdenr: s @r2 r r sC II
uT{{ T
L!"JKK Kqy#'ll # # #w--'ll # s3g&&
'
'CK HIII
w--C A ?@@@ 9s 33 A. .(BBc 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 r@ r$ rm r )rF s r2 r r se
G5S$ (())) G G G < = =BF GGs .2 A
c . t | | u rt | } d}t | }|dk rt d |dk rQ|O| d }t |t
j t f r|dk rt | |S t d |t d| }|}nmt | |u rt | }t | |k rt d t d t || D \ }}} t | | } t d t || D \ }}} n# t $ r Y dS w xY w|dk rt d t ||z | S )
Nz.harmonic mean does not support negative valuesr4 z.harmonic_mean requires at least one data pointr zunsupported typez*Number of weights does not match data sizec 3 K | ]}|V d S r6 r0 )r8 ws r2 r; z harmonic_mean.. s" G Gq G G G G G Gr1 c 3 . K | ]\ }}|r||z nd V dS )r Nr0 )r8 r rS s r2 r; z harmonic_mean.. s3 P PTQ 0q1uuq P P P P P Pr1 z Weighted sum must be positive)iterr r r
isinstancenumbersRealr rg r rP rv zipZeroDivisionErrorrr )
rF r ru r: rS sum_weights_rO rN rG s
r2 r r s , DzzTDzz
=F D A1uu NOOO
aGO G a',0 1 1 01uu%f---H .// /A,, ==G # #7mmG w<<1 !"NOO O G GIgv,F,F G G G G G;v&& P PS$=O=O P P PPP5%% qq zz =>>> K%' + ++s ";E
E,+E,c t | } t | }|dk rt d |dz dk r | |dz S |dz }| |dz
| | z dz S Nr no median for empty datar} r4 sortedr r )rF r: is r2 r r % sw $< >Bc t t | d } |d d S # t $ r t d d w xY w)Nr4 r zno mode for empty data)r( r most_common
IndexErrorr )rF pairss r2 r r sn ,
DJJ + +A . .EBQx{ B B B 677TABs
A Ac t t | }|sg S t | fd| D S )Nc & g | ]
\ }}|k |S r0 r0 )r8 rq rG maxcounts r2
z multimode.. s' J J JleU8I8IE8I8I8Ir1 )r( r maxrM rC )rF countsr s @r2 r r sa T$ZZ
F 6==??##H J J J Jfllnn J J JJr1 r exclusive)r: methodc |dk rt d t | } t | }|dk rt d |dk rg|dz
}g }t d| D ]M}t ||z | \ }}| | ||z
z | |dz |z z |z } | | N|S |dk r||dz }g }t d| D ]b}||z |z }|dk rdn||dz
k r|dz
n|}||z ||z z
}| |dz
||z
z | | |z z |z } | | c|S t
d| )Nr4 zn must be at least 1r} z"must have at least two data points inclusiver zUnknown method: )r r r rangedivmodappendrm )
rF r: r ldrw resultr r delta interpolateds
r2 r r s 1uu 4555 $< >>r1 c t | | \ }}}}|dk rt d ||z }t |t rt |j |j S t |j |j S )Nr4 z'pstdev requires at least one data pointr )rF r rO r rT r: r s r2 r r sz dB--KAr1a1uu GHHH
q&C!W E$S]CODDD s}co > >>r1 c 6 t | \ }}}}|dk rt d ||dz
z } t | t |j |j fS # t $ r1 t | t | t | z fcY S w xY wr )r\ r rf r rn ro r_ )rF rO r r r: r s r2 _mean_stdevr s PYYNAr41uu GHHH
A,C4T{{ /
s O OOO 4 4 4T{{E$KK%)) 333334s (A 8BBc @ t | }t | |k rt d |dk rt d t | |z t | |z t fdt | | D }||dz
z S )NzDcovariance requires that both inputs have same number of data pointsr} z,covariance requires at least two data pointsc 3 4 K | ]\ }}|z
|z
z V d S r6 r0 r8 xiyir ybars r2 r; z covariance.. 4 AAVRT b4i(AAAAAAr1 r4 )r r r% r )rS yr: sxyr r s @@r2 r r s AA
1vv{{ deee1uu LMMM 77Q;D 77Q;D
AAAAAs1ayyAAA
A
AC !a%=r1 c t | }t | |k rt d |dk rt d t | |z t | |z t fdt | | D }t fd| D }t fd|D } |t ||z z S # 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 pointsc 3 4 K | ]\ }}|z
|z
z V d S r6 r0 r s r2 r; z correlation.. r r1 c 3 , K | ]}|z
xz V d S r6 r0 r8 r r9 r s r2 r; z correlation.. 0 00R$Y!#000000r1 c 3 , K | ]}|z
xz V d S r6 r0 )r8 r r9 r s r2 r; z correlation.. r r1 z&at least one of the inputs is constant)r r r% r r r ) rS r r: r r[ syyr9 r r s @@@r2 r r s+ AA
1vv{{ efff1uu MNNN 77Q;D 77Q;D
AAAAAs1ayyAAA
A
AC
00000a000
0
0C
00000a000
0
0CHT#)__$$ H H H FGGGHs C' 'DLinearRegressionslope interceptF) proportionalc r
t | }t | |k rt d |dk rt d |rAt d t | | D }t d | D }njt | |z t | |z
t
fdt | | D }t fd| D } ||z }n # t $ r t d w xY w|rd n
| z z
}t ||
S ) NzKlinear regression requires that both inputs have same number of data pointsr} z3linear regression requires at least two data pointsc 3 & K | ] \ }}||z V
d S r6 r0 )r8 r r s r2 r; z$linear_regression../ s* 33vr227333333r1 c 3 K | ] }||z V
d S r6 r0 )r8 r s r2 r; z$linear_regression..0 s& ''r27''''''r1 c 3 4 K | ]\ }}|z
|z
z V d S r6 r0 r s r2 r; z$linear_regression..4 s4 EERBI"t),EEEEEEr1 c 3 , K | ]}|z
xz V d S r6 r0 r s r2 r; z$linear_regression..5 s0 44BdNAa'444444r1 z
x is constant r )r r r% r r r ) rS r r r: r r[ r r r9 r r s @@@r2 r r sk $J AA
1vv{{ klll1uu STTT 533Q33333''Q'''''Aww{Aww{EEEEE3q!99EEEEE44444!44444/c / / / o.../#< )l@g N@g"]Ξ@gnC`@gu @giK~j@gv|E@gd|1@gfRr@gu.2@g~y@gn8(E@r 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&>@gigtcI,\>gŝI?g*F2v?gC4?gO1?)r r r$ )pr sigmar rr r rS s r2 _normal_dist_inv_cdfr A s
CA Aww% 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
#XX37A c!ffW
ACxx
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 A3ww
B
U r1 )r c ( e Zd Z dddZd#dZed Zdd d
Zd Zd Z d
Z
d$dZ 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!dS )%r z(Arithmetic mean of a normal distributionz+Standard deviation of a normal distribution_mu_sigmar r c |dk rt d t | | _ t | | _ d S )Nr zsigma must be non-negative)r rf r r )selfr r s r2 __init__zNormalDist.__init__ s; N 3;;!">?? ?99Ell r1 c ( | t | S r6 )r )clsrF s r2 from_sampleszNormalDist.from_samples s ?sK%%&&r1 N)seedc | t j nt j | j | j | j cfdt | D S )Nc ( g | ]} S r0 r0 )r8 r gaussr r s r2 r z&NormalDist.samples.. s% 333Qb% 333r1 )randomr
Randomr r r )r r: r
r
r r s @@@r2 sampleszNormalDist.samples sY G $ &-2E2E2KHdk E333333%((3333r1 c | j | j z }|st 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 rS r diffs r2 pdfzNormalDist.pdf sb >;, J!"HII I48|4$;$/233d3>6J6JJJr1 c | j st 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 rS s r2 cdfzNormalDist.cdf sI 6{ J!"HII IcCTX $+2F G H HHIIr1 c |dk s|dk rt d | j dk rt d 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)r r r r )r r s r2 inv_cdfzNormalDist.inv_cdf s[
88qCxx!"HII I ;# !"QRR R#Atx===r1 r c B fdt d D S )Nc @ g | ]} |z S r0 )r )r8 r r: r s r2 r z(NormalDist.quantiles.. s) 999 QU##999r1 r4 )r )r r: s ``r2 r zNormalDist.quantiles s0 :9999U1a[[9999r1 c t |t st d | |}}|j |j f|j |j fk r||}}|j |j }}|r|st
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 | |
| |
z
t | | | | z
z z
S )Nz$Expected another NormalDist instancez(overlap() not defined when sigma is zeror r+ )
r r rg r r r r r r" r r r$ r ) r otherXYX_varY_vardvr r{ bx1x2s r2 overlapzNormalDist.overlap s
%,, D BCC CU1
Hae !% 0 0 0aqAz1:u NE N!"LMM M
U]
!%!%-
=R3>F#: ;<<< <
EEMAEEM )
Hqx $rBwc%%-6H6H1H'H"I"I I!er\!er\d15599quuRyy011DrQUU2YY9N4O4OOPPr1 c T | j st d || j z
| j z S )Nz'zscore() not defined when sigma is zero)r r r r s r2 zscorezNormalDist.zscore s6 { M!"KLL LDH ++r1 c | j S r6 r r s r2 r
zNormalDist.mean s 5xr1 c | j S r6 r* r+ s r2 r zNormalDist.median s 6xr1 c | j S r6 r* r+ s r2 r zNormalDist.mode s
xr1 c | j S r6 r r+ s r2 r zNormalDist.stdev s
8{r1 c " | j | j z S r6 r/ r+ s r2 r zNormalDist.variance s +{T[((r1 c t |t r5t | j |j z t | j |j S t | j |z | j S r6 r r r r r r$ r% s r2 __add__zNormalDist.__add__! Z b* % % L bfrvouRY /J/JKK K"&2+ry111r1 c t |t r5t | j |j z
t | j |j S t | j |z
| j S r6 r2 r3 s r2 __sub__zNormalDist.__sub__/ r5 r1 c ^ t | j |z | j t | z S r6 r r r r r3 s r2 __mul__zNormalDist.__mul__= ,
"&2+ry488';<< > > : : : : Q Q QD , , , X X X X ) ) X) 2 2 2 2 2 2 = = == = =- - -. . . H H; ; ; - - -P P P% % %& & & & &r1 r r6 )rs )r )J__all__r` r r sys fractionsr decimalr itertoolsr r bisectr r r r r r! r" r# r$ r% functoolsr& operatorr' collectionsr( r) r* r rm r rP r\ rb rD rA rr rv rE r|
float_infomant_digr __annotations__rf 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 r0 r1 r2 rg sv hT .
% % % % % % % % , , , , , , , , < < < < < < < < < < < < < < < < < < < < 8 8 8 8 8 8 8 8 8 8
c j 3 3 3 l& & & & R 4 4 4>+ + + \ $ 3>22Q6 6 6 6
#3
#3
#5
#
#
#
#S S W <" " ",# # # #LG G G&5, 5, 5, 5,p+ + +0 , &E+ E+ E+ E+PB B B<K K Kr ; (4 (4 (4 (4 (4b)% )% )% )%X& & & & R? ? ? ?$? ? ? ?$
4
4
4( 8 H H HB : 02HII 05 8> 8> 8> 8> 8>|G G G V 0000000 D \& \& \& \& \& \& \& \& \& \&s D D'&D'