ÿØÿà JFIF    ÿÛ „ !.%+&8&+/1555$;@;4?.451 4,$,44444444444414444444444444444444444444444444444444ÿÀ  á á" ÿÄ     ÿÄ ?    !1AQaq"2‘¡±ÁðBRbrÑá#‚’¢²3S CñÿÄ   ÿÄ !    !1QAa‘2ÿÚ   ? 5˜Z¯V¦cø)›t/? z¨±>Õ5€¶‹Á¤·¼z¼Ü¬+ñ®v¤¨_ˆR­BFn©—˜ý®ç̝P8gýt·ÉSTŦˆìät?þé¼íìN/Þa)ì–í6ô… Ï¿øÃj´¿KÇü]ÿ ªô¹-eKànëÕHTx}ýSÜ›ÿ ”7Ø×&µ<¦  ¥ÑO¶[Ù¯ä¨ÞÃÿ PZ-¬;#õ|•oaÿ ©CìÞz3˜öː/¤­ñTûIØ}š^ mÓ%ªxˆ¥ÉŸu=Z+ISe¿45™¼u;ú&WØ÷€æßQ™®{|íx*TC“#ZŠìZ§²‹ 6pv…³¿¡äª*áZÐ%ÒOáˆo"x«OHk w±æ+¬V(kMúŸ5Vö«$ ÁrÏbàb57/luR ¸ÑÛj Òµì`Мq­û žICÀÊ•©4€Âcà¨Ï€O´<èÐ:›ù(Ë^L8þ‘ÍÌ#¸Ð_Ì©ÙK(Öz 4¬û+¸;ü’V’84‘¬ÃŽ:[â‡ÔÌáõp¢~§ªlæ£ö{®G>J¼"°‡7¯ÆÉèßû ‹É‹§ÁòÃýâßî ^ƾÙõ‹×óH#«LP½ïX=xÑÍ$|W?•~• îëÔ©ª‹ {ÝT…Kÿ ”hûâá)J*ö˜–ÔU;iÇ€/ ÆþjóZ\ýwØ=Ìm ºèËL9 ýèÆð/¨’¥öo=nË.%Îì ŽÕ¯È|{Oj²ƒE6e/ßdÄõ²Ìâ1O®ò×TsəԸhOMýíMˆ¿¼H˜l²,7Â¥#MF/Úf°Ö½± ¸–dr‹NýÊ íjqx{œÉ ä-È ¦ øÄër¨q°ð †nцýÑÄÆ’mä…n<0È™;ÁÝá¯ÁZƒ7FÀmì­ É&9ˆîéi¶ùN§Y• ÃZãAâ?•‡©‰ , ó¾IŸŠc1 4â&y­&pŠ­6;M À 0¹qç»p.á …ŸÅáK@%6·y6ƒ‰3?”úºŽ‰éX5ªPT §µ!=Mž«Ú½‹ÅgÂSâÉaþÓoö–¯ÁÔìR>5éÿ üs¶ÆUcÌ kÇR ]ÿ ù¬¼«VŽ;Â|‡~¢¦”ÏŰæ {L™Õ°Óv¹ò¸írޡעCÃ!íVÕ {¶»sŒNPg/ "uÕbkm²“$ďå¿é¹§°½æz¯6 †s¿!s–wÚÝ“™Œ °.ûj>·+™Òa…©Œ&rÝÎtÛë긪Ît’LAVp%c Úý[ÄzJ¾ÇàXXç@˜ó<êL]·T˜¾¥1Ó©V‡g´æ½¦Ý@¹óø!_@´ÞâSÁ —S3™•& ]@JHÚý©ZŽ €×æÔr»Áf!‡yÞ4Mv*èÓã_{‘åóUuљØ«Oïé*®EvÑ Œ÷‡U \"㪒ÍK+À 4“M¡ï:0¥5í!'<@î´”>Ç»&Z–ïCCV˜Ì5Šo&îhè.žû |ÓK©h$s6KìŒëã)¹hI¦GïOåóI;ììü#É$Š0…Ææ¥TØ.5­¾gn´ “ÂÖ\:hœ89G)J@„}œ:’Ò{/Š"¦_Æ×7Æ3VÇŠÊa]ÚŒÙ€Ä–=®uÁßâACZƒ§§£ Qnâ:«,×{tyø¬iÛcœÜÄ€H½ÄÍCk´÷šß .W'b¤Íåh]÷€=,Žv×cÚEÚHXJX¶îo¨FÒtèöŸ>ªª6[J®Fµ£sGÁeqõfe\íjÒÐïÄÐGˆe1Ø‹.Ø”‘Ëuø Y­ˆÜ ŽG|zùªüMpDnQWÄ”%JŠ™)â*p@Örš«ÕT2Ð%ˆG#ª„ ·¤!°ŸOTÂT¸aÚ%4&h™LµšØüÐ.F¿²ÐÞ_Ç‚¾ÅÃaÜ÷09Æ q€öy˜v‡85õN÷]¬äѼóS{°_MެúÔ#°Ç¸0åÞè2ëôPcvÆw9®ií1Ä8F™˜à‰´+‰Ik1òÝ7“Ñ×ÒsÝ\x‚h`ÞÑ`ó"|µEcý£n˜h`}GÞ !±ù²Ápü²ß6 0ïi󜵩SÈÇ7˜-ÕURO˜¦´f$ªž-Í6(œ}<„ éc øs]ŽŽ„*—¾ ìdŽ„)méª\¿êÎIg¾ØÞ~I#C/¼¼´EÁÈŽi8“©õådô·>euä ƒ'Ê×लR1ÉJE1ÐAát`t;ÇР%Ý<‡¥„ÍÆ`×Oyó)õiI€ñQaŸ4Ûù\áàaÃÔ¹HÃu¹*k€¦<„e S‡&õÏ B!ŽhüÞ`yj}mªf×\¿ Ç~æ­9‡û\՞Ǖg²1Žû5V7 !àöšm° c`ܬøÇìµÒ'P"?…´Ö,"§^•õލsÔ)6˜sæéÍR¼ ò|Sl”‹7 nPW Gòú÷½§O¯‡„l¡kSÞŒr½PÊ@æ¢pŽ-mÿ #Ÿ˜Àº¶Áä¦;ïÔæ$1££`“Õ>„—·ž)ßð³ñ#Ï Ô$¶œ‰ÊE‹À;÷º ¯«P:Ñ”8–IÊtpÞ3ª“>ê“þës4ò2OÏÕ­±zô†Õ§‰.÷ä¸;¿˜“'œ›žª}«Œ{ª±Ì 9ÔóÞÕ‡0 $íWV3Üì¬ —@kÝ4@¿r¼±½¬™›?øØæ´'Áé®CË3-g$˜ö‡×auÚi´Žp/êÛ æF›Ú2v‹ã¿¿,nB1̨ƃqÞa5͝@&Æû“él÷ \C²½UÍc ¯k×¢U ÖéQå™—-r wô ÞÏ<Ò=&=ÿ Ôê Òêˈt,i—;LîÜ á¸*ÚÃ1$êL•LÍ <É)ýÐà’ ;F™{ƒ™˜€&'}‚ãÄK`¡ÞT@I;®žZóè‚s’7®°›+§O­Åq©é»²9<Ô J ¼9O’HL»Ùïì¸rk¼Ž_ý‘TŸu[²ßÚŒ·ü÷B%¯E ŸÔX5êO´ Ç•€’I0 ÉJX` ñ¹õ%;µŸD‘«´€àwÒ™U ûئžÖö\×®×´8 ½‡ºÐÆÓ§?Àkmœ=;d5*@-ì0F Rªýš[Ü6âö̃ڸr*KA9· u*µæ£?U¸Âêí†8@¦X4 e-ò„0s{ HâUpU?¼mñRa°®a%Ð'tÉ×’\¾ÊÉ]t›h>·(Ë@R¼¡Ãt h}’O÷au<+nT…Ö…MӐ??Óe95 q>í/;&JSû °¯ÊéÞ øƒ*Ã2½Ài&:nôUl=¾¿5eˆ3”ñc|Ú2V”>„»&eE;«ÚäC p¢Û úy 9š[ŒÌx¼擼A&DåÒ¯ˆ¤ÀÌ;"˜ ÏQä¸åhÊ}Ûq«Û0WžÒ|»€ø®öCm5•\ÇÀ§Pe3£]0ÃàLDÉ‰1øªxjgwT‚÷¿LΨK‹›ùs—xˆÜ±µ kæ¸f‰‰ÜGk/LÛØ6d9ò¶ùA{ƒA3š/¬D¬khÓk‰`˜"㯒r¿±Óã jx‡°e}<Ñø\3y:'À•/h½Í€Ç4~g ?Û(¼]v‘ªlKÎâ~?O‚W%{Ì:“'©úNq¾›úo(X’¥¯ˆ nFê{Ç€ü?º'ë ø‹ì Þ09ŒÌç9Æ —ËC`j@ÓÄ(+a‹un¸#ÂꟋ{K`‘ÑÍÍ'à´»/Û,KW;Þ4²þð ï Nm|~fGÏ(…³Ã)«1ö­Õ ¥‡¨©ƒÃ™ü-s=à=U66Ï«Ýc蓦W¹íž®›nÔ%êÇìŒ<#Ü×84ån®Ð ÒåOC` ñânÑs‡¢ç 1õ%Îhì½Ã½® e:ݼUZo™`  ÅZŸŒÊ«ê1ÏÄo$q¹Þ€©ˆhÐÉä¯ñ[!…Ú˜àJ:x2$Íß&PåT£6ç— ‡Í*4Ýšçjÿ ‰É nófÐ ó(L5C•åÆ\rMÒ@ò }y-W}™üýVù—ú¢=Ù”c®‘< M ž ´Phr ¦©TD ‘ù.$´÷O‡‘V2Æò.=IUŒ=ž‡â¬i™aþÓåÙ?òUø'ØÖ•.~* šTŒ!•-×áºTâ®ä#õü'´ eýlYÅÓeÕKÂrT"CÚ@u!Óxƒ{š3€}1¿(r}%«nËamjÑ%ÑNEò v ˜à  σöK³,*º.àzù¨™Ó ÚçâU¦*¿ 9{%Ö¹ njûdaXöb) kÛÆ±ûÓ\°M7ˆÂ=û›ç¿Ã‚­V»Cg–8ÙêE- j)k$º`Ã-ùEýeBÆÇ]c¡°ñty&Òd0nõ'¡W+ƒ*|–øµFa\GQªEAÔp5\Ǽ·¼Ç8·õ -â§Ú[ ‡ uZeÖ 3}×d'+¹:ð+K†Û®s!Ï$úe€<Û”x)1»a­¡LC]¸µík…ÚàA»AYº{†ªS[¦5HÒ7ù --,ísòDØ€èk ÞÀîÜ ò@â( ËNˆë›4ô½•/¦o‡€Û7 ê•ÆêòðÜy'Án½µ á˜ݦ ndeo…[ì¶Ê,¥R³Ä=À±—–ß;£™´ñSâ*g§”ïaið‘Jå~™ÓÞ ß³Õ¢»8x埒²52>AÊb&-÷\7´éÄù€T˜,w;3{ï˜k…à¹ÄqÀ«œ{€\ ˆ¾[´¨јr &Úé„Ívˆ±8†¿]|¬ņ4I×pÞS1ÈÖz‰#Ìv‡G!YNògñ:màTz¢Ý1ô©^O=~ë|5Bã™ç•¼µõ•bÆ@úÕS¬ÈŒ#¬zünrŸ û” Z²•èðV"ÁHÚý©wÝ €7¼Ìu1hÑa3Éä û f$o¿É ™Ú›ÝçnpÒ3äÌ3†Í§,Äï]$‰/pê †«À¼¸e9­Æê_C]žƒ·ý·frÁN«, E=›Çq -‰öŒ:aÏ¿±í&£Í:-} 84‘ÿ eƒQÑeëSsuiA ³g㟥ú£?ÿ ʼn*”“÷aühe:ÊWa@ÒÞk±eØ] F Ô—r.åä˜ @ö¥ªZoÐýYL·¥S²G/‡ñ <~*ZÆ´è>JlòàÛÆ½ÿ 窘ìGN¢:I®KšJp/`íIÁÀõ#Ä-€ö­šµŒoF4|ÆQØÆ@Ì|£Ô…¢À{9˜è½Üó›€ôYÒÎYsið;ís¤€à²ˆ‚4qÉVŒI$ ‰"° æµ8cXGjœˏ¡Aâý•ËÜ¢ûï e·çLx']á"oÅÎê3¯Ç—¹”ó0nå‚âg{Œñ> S´˜îè°g238‚ãköÝfÚd´6Ò€;ò÷±¢™¼›º ¢Æ'¥Ðx'e¬ç ]bÈÆV¢ó‹kýBO ðÊâ$Ÿ!×T 3Mýמ žìٍàÌü‘8÷€àæØ8æ©6‰©L´«…oãpð„~Çk‰!ñ;‹”ÛžÍ àž±z Ÿôû øŸÝužÏ;ÿ #|u6™Þ¬ÚˆÐõA4¶â|ôl|Ê2ŽÇ¤ÝÅÇY.<#Aí.k§hóF‚”Y; M½Ö4hŸ4&›­¿tès´%FìL¥£Ãk‰ÇT¤haÁ¤ÚxfÉ`ÑìË›>i 3t‚:,–+^÷´–{Û–Nxi"x‘Ûg î¨>¥Õ܁ùZH,2Û“:8xÊ¢Çí9.É-Ìâã-=çjwµS˜dütžçwýGòú®®ûº_ˆýx$–¡ãøO EÚÛÏ÷R„×w+3£Á£öUMyR²¹âŒ°š›¸Ñãò9§Ó_Dl+Ùßc›úšGÅÌc†Ž!Ko=¶.‘Îÿ c²(2®V mª.ÿ ¹B›¹å ù„öŸSV>™ü¯$y:G¢Z×àøúdî¹û­·ýÇ´:•c LÍõi_‹ö+ÎæGÊè>OŠ•äž´§Þ{X}¨1ÚTc›»Qþ•êô°t¿OP?eæ~É{5]•ÙR£r5†nZ\ã@ &îJõ ¾àC°þV>fé¥/ü5ñÊIº_é5 ;e­h<@ Ä&æÃëE%;X,ÒãÆÞ`Oò¦kŸm#˜!ÀyÄ¢| óLšò¥Ä` ¶R=|ÈCâh5ò3DˆïF†ðÒ#ÅìÛœ?¸yhBãœí ZxßÎÄhºRK„`Þödvײ™ÀÈÑÒgŒuY w³%†ƒÓzõ ÖÏp‚dH®¦A´ù§»ÓÇMæ~)ˆð‡û:ù&Ä •vGD´À n ݇¼Ö8Fö óáà£~Ë¥x`oK|Ä?fxiØü%pìR>éò+Û±éÎ>núlFŤ'tq8LZÏvÃ?„¡ß±È⽆¯³íü@x|PöUäèØã¡ð‚ŒAìÏ"vÍwóŸÍ{ ý0.z È•Ö{,N¡£¡ŸKÕÙž>Ýœþ ÍÀ°<×EA!Å‚D™IúOÍ¡>ôG}Â` ÍßkÜL™Ž Þð™ {IøF²¹òQ3&!ÃÂÞz.d&Ï-sH¸,Ôõ˜ŽP€ 77ˆÝ¼ÊëÜw =cÕ Ú,ØÐ5ÎYÐ)ì´öœgŒ[¤ßv㙑8心>h]§µháYš£²ºÑ.{Ï7Sð•?´~×SÃKýJÛ˜ ™Íäiúu<µX¶1õ^kâçIÑ£sZ4h>j*ÔšD:4­¿_ ÷¸ Õxæÿ ¸?Mù _•­ÊÐ ä ÷ý ÑwL œ­ïnTkÛUÍN©ë:¦fV ¶ÜÔÜMªÅâA½–¿R×TXš-%iTÊT•‡Ù‚JôϐZxWÑè‰f‰òG º ×Õû2aZ7OU3[“×AT–ÞŒ…-‘¤”Ì ì&(ˆ¿­•ƒkï’:ðY¦W‘ Å)“†‘˜³Åtcø˜ñTÂwÚÇ4|üLÇªí–v- qˆèU qPE.†â‘˜µ Æ,ÐÅs]8¾„oúÑ i>ÜxxÈó)ƒ ´æÁâØ$À‰vžŸf$Ž |ãw;ÀÁIJ»b` {¦Ó¤Ú$©YÀ‘n@Óïž«9J¼êG m¤ ܯ¹ÌW4€ÐÒÅÛ‡#褕Ÿn-?í|с¥÷Ú¹¬'´ÞÜ9ÓK `hê£SÄSà?7—Wí_´…óB›»:=Ãïq`<8ñÓŒÑlú2d¬ê³£hÖ[l|$vÝro~'R®‰§°ñmY ͧäP |PUª¹·:3Œ[Û{Xÿ ºâ@‚W–Äé u‚ ¯´*=íή.pûÒdt @G‰¬ s¸ ëÉücr ÞæÑ¨Ê@>¤¢Ö±. Þ'¯°ÌME[YéïĵÂCå½ Ué©Áû'Ê9%eÔðNU”ë‘ÌsD3/®+UI˜9h.WC”빓$#:pz:YÓ ¿xž* ³$Í +$kñAŠ‹†¢ Uê>¸)_š¬÷©ßAÂÔb9ÇU ¯¾á•9¯ÏÏ÷O÷¼¼Fähal1‰3Ì[Ïr•´UCksNÐ] R‘¸¥H+§Šé†c©vÖÞ0iÓ76s†î!§=ß ¼~Ô'°Ãmäoäš³ªøi1úÉ)³yV8 CLÄØÁ‘WYïi€H6ÖÑiámø^ÈY´°Ñ7¥Û*—Ñ©L«Qƒï—Ùrÿ ›£Ð*š¸ˆL©ˆ$ˆ ÷¾D§9È®«qbqC)–ˆïv´çñsÑVT­Ø, <àïºÀO«Jý·õ àfPìð .wFšir´þ’2_Y *Æ€x\« ì€9š@ Ž|F⇥ˆkZ@hÖÄ0t¿-<“‹qµ¾*ZL¤Ú)&BJpÓF5=$„at*Zš$’ÑtdûÝRI1 2މ$€$I$#‰SÞ’Hë¬ï;Á$¡t$’`<(ñÇt)$‡Ð.Êf¢X’Kt=Éé$‚ˆªè¢oÝëòI%Rgcª÷ŠyI%¡‰ÿ !ñ)´õ $¤ Ô’IIGÿÙa 3g@sgdZddlZddlZddlZddlmZddlmZddlm Z ddl m Z m Z ddlm Z mZmZmZmZmZmZmZddlmZdd lmZGd d d eZdFd d ZddZddZddZddZddZ ddZ!dGddZ"ddZ#dd Z$d!d"Z%d#d$Z&d%d&Z'd'd(Z(d)d*Z)dHd,d-Z*d.d/Z+d0d1Z,d2d3d4d5d6Z-dId7d8Z.dJd9d:Z/dKd;d<Z0dLd=d>Z1dMd?d@Z2dAdBZ3zddCl4m3Z3Wne5yYn0GdDdEdEZ6dS)N) NormalDistStatisticsErrorfmeangeometric_mean harmonic_meanmeanmedianmedian_grouped median_high median_lowmode multimodepstdev pvariance quantilesstdevvarianceNFraction)Decimal)groupby) bisect_left bisect_right)hypotsqrtfabsexperftaulogfsum) itemgetter)Counterc@s eZdZdS)rN)__name__ __module__ __qualname__r&r&)/opt/python39/lib/python3.9/statistics.pyrusrc Csd}t|\}}||i}|j}ttt|}t|tD]@\}} t||}tt| D]"\}}|d7}||d|||<qRq6d|vr|d} ntddt| D} || |fS)Nrcss|]\}}t||VqdSNr).0dnr&r&r' z_sum..) _exact_ratioget_coerceinttypermapsumsorteditems) datastartcountr,r+ZpartialsZ partials_getTtypvaluestotalr&r&r'_sum{s$  r?cCs,z |WSty&t|YS0dSr))Z is_finiteAttributeErrormathisfinite)xr&r&r' _isfinites  rDcCs||ur |S|tus|tur |S|tur,|St||r:|St||rH|St|trV|St|trd|St|tr|t|tr||St|trt|tr|Sd}t||j|jfdS)Nz"don't know how to coerce %s and %s)r2bool issubclassrfloat TypeErrorr#)r;Smsgr&r&r'r1s  r1c Csznt|tust|tur$|WSz|j|jfWWStyjz|WYWStydYn0Yn0Wnttfy|dfYS0d}t | t|j dS)Nz0can't convert type '{}' to numerator/denominator) r3rGras_integer_ratio numerator denominatorr@ OverflowError ValueErrorrHformatr#)rCrJr&r&r'r/s   r/cCsnt||ur|St|tr(|jdkr(t}z ||WStyht|trb||j||jYSYn0dS)Nr()r3rFr2rMrGrHrrL)valuer;r&r&r'_converts    rRcCs.t||}|t|kr&|||kr&|StdSr))rlenrO)arCir&r&r' _find_lteq s rVcCs>t|||d}|t|dkr6||d|kr6|dStdS)N)lor()rrSrO)rTlrCrUr&r&r' _find_rteqs rYnegative valueccs$|D]}|dkrt||VqdS)Nr)r)r=errmsgrCr&r&r' _fail_negsr\cCsHt||urt|}t|}|dkr,tdt|\}}}t|||S)Nr(z%mean requires at least one data point)iterlistrSrr?rR)r8r,r;r>r:r&r&r'r's rcspz t|Wn.ty:dfdd}t||}Yn 0t|}z |WStyjtddYn0dS)Nrc3s t|ddD]\}|Vq dS)Nr()r9) enumerate)iterablerCr,r&r'r:Oszfmean..countz&fmean requires at least one data point)rSrHr ZeroDivisionErrorr)r8r:r>r&rar'rAs      rcCs6ztttt|WSty0tddYn0dS)NzGgeometric mean requires a non-empty dataset containing positive numbers)rrr4rrOr)r8r&r&r'r\s  rcCst||urt|}d}t|}|dkr2tdn<|dkrn|d}t|tjtfrf|dkrbt||Stdz"t ddt ||D\}}}Wnt yYdS0t |||S)Nz.harmonic mean does not support negative valuesr(z.harmonic_mean requires at least one data pointrzunsupported typecss|]}d|VqdS)r(Nr&r*rCr&r&r'r-r.z harmonic_mean..) r]r^rSr isinstancenumbersZRealrrHr?r\rbrR)r8r[r,rCr;r>r:r&r&r'ros$  " rcCs\t|}t|}|dkr td|ddkr8||dS|d}||d||dSdSNrno median for empty datar(r6rSr)r8r,rUr&r&r'rs   rcCsLt|}t|}|dkr td|ddkr8||dS||ddSdSrfrir8r,r&r&r'r s   r cCs,t|}t|}|dkr td||dS)Nrrgrhrirjr&r&r'r s r r(c Cst|}t|}|dkr"tdn|dkr2|dS||d}||fD]}t|ttfrFtd|qFz||d}Wn&tyt|t|d}Yn0t||}t |||}|}||d} |||d|| S)Nrrgr(rhzexpected number but got %r) r6rSrrdstrbytesrHrGrVrY) r8intervalr,rCobjLl1l2cffr&r&r'rs&       rcCsBtt|d}z|ddWSty<tddYn0dS)Nr(rzno mode for empty data)r"r] most_common IndexErrorr)r8pairsr&r&r'r s  r cCs@tt|}tt|tdddgf\}}tttd|S)Nr()keyr)r"r]rtnextrr!r^r4)r8countsZmaxcountZ mode_itemsr&r&r'r 4s r  exclusive)r,methodc CsB|dkrtdt|}t|}|dkr0td|dkr|d}g}td|D]D}t|||\}}||||||d||} || qN|S|dkr0|d}g}td|D]r}|||}|dkrdn||dkr|dn|}||||}||d||||||} || q|Std|dS)Nr(zn must be at least 1rhz"must have at least two data pointsZ inclusiver{zUnknown method: )rr6rSrangedivmodappendrO) r8r,r|ZldmresultrUjZdeltaZ interpolatedr&r&r'rks2$   $$ rcsdur,tfdd|D\}}}||fSt|tfdd|D\}}}tfdd|D\}}}||dt|8}||fS)Nc3s|]}|dVqdSrhNr&rccr&r'r-r.z_ss..c3s|]}|dVqdSrr&rcrr&r'r-r.c3s|]}|VqdSr)r&rcrr&r'r-r.rh)r?rrS)r8rr;r>r:UZtotal2Zcount2r&rr'_sssrcCsLt||urt|}t|}|dkr,tdt||\}}t||d|S)Nrhz*variance requires at least two data pointsr(r]r^rSrrrR)r8xbarr,r;ssr&r&r'rs& rcCsHt||urt|}t|}|dkr,tdt||\}}t|||S)Nr(z*pvariance requires at least one data pointr)r8mur,r;rr&r&r'rs# rcCs6t||}z |WSty0t|YS0dSr))rrr@rA)r8rvarr&r&r'rs   rcCs6t||}z |WSty0t|YS0dSr))rrr@rA)r8rrr&r&r'r $s   r cCs|d}t|dkrd||}d|d|d|d|d|d |d |d |}d |d |d|d|d|d|d|d}||}|||S|dkr|nd|}tt| }|dkr^|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}||}|dkr| }|||S)5N?g333333?gQ?g^}o)@gE.kR@g Ul@g*u>l@gN@g"]Ξ@gnC`@gu @giK~j@gv|E@gd|1@gfRr@gu.2@g~y@gn8(E@?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?gEF?gPn@g&>@gigtcI,\>gŝI?g*F2v?gC4?gO1?)rrr)prsigmaqrnumZdenrCr&r&r'_normal_dist_inv_cdf7sd      r)rc@seZdZdddZd9ddZedd Zd d d d ZddZddZ ddZ d:ddZ ddZ ddZ eddZeddZedd Zed!d"Zed#d$Zd%d&Zd'd(Zd)d*Zd+d,Zd-d.Zd/d0ZeZd1d2ZeZd3d4Zd5d6Zd7d8Zd S);rz(Arithmetic mean of a normal distributionz+Standard deviation of a normal distribution)_mu_sigmarrcCs(|dkrtdt||_t||_dS)Nrzsigma must be non-negative)rrGrr)selfrrr&r&r'__init__s zNormalDist.__init__cCs.t|ttfst|}t|}||t||Sr))rdr^tuplerr)clsr8rr&r&r' from_samplesszNormalDist.from_samplesN)seedcsB|durtjn t|j|j|jfddt|DS)Ncsg|]}qSr&r&r*rUgaussrrr&r' r.z&NormalDist.samples..)randomrRandomrrr})rr,rr&rr'samplesszNormalDist.samplescCs<|jd}|stdt||jdd|tt|S)N@z$pdf() not defined when sigma is zerog)rrrrrr)rrCrr&r&r'pdfs zNormalDist.pdfcCs2|jstdddt||j|jtdS)Nz$cdf() not defined when sigma is zerorrr)rrrrrrrCr&r&r'cdfszNormalDist.cdfcCs:|dks|dkrtd|jdkr*tdt||j|jS)Nrrz$p must be in the range 0.0 < p < 1.0z-cdf() not defined when sigma at or below zero)rrrr)rrr&r&r'inv_cdfs  zNormalDist.inv_cdfrzcsfddtdDS)Ncsg|]}|qSr&)rrr,rr&r'rr.z(NormalDist.quantiles..r()r})rr,r&rr'rs zNormalDist.quantilesc Cst|tstd||}}|j|jf|j|jfkr>||}}|j|j}}|rT|s\td||}t|j|j}|sdt|d|jt dS|j||j|}|j|jt |d|t ||} || |} || |} dt| | | | t| | | | S)Nz$Expected another NormalDist instancez(overlap() not defined when sigma is zerorr) rdrrHrrrrrrrrr) rotherXYZX_varZY_varZdvZdmrTbx1x2r&r&r'overlaps"   (  zNormalDist.overlapcCs|jstd||j|jS)Nz'zscore() not defined when sigma is zero)rrrrr&r&r'zscoreszNormalDist.zscorecCs|jSr)rrr&r&r'rszNormalDist.meancCs|jSr)rrr&r&r'rszNormalDist.mediancCs|jSr)rrr&r&r'r szNormalDist.modecCs|jSr)rrr&r&r'rszNormalDist.stdevcCs |jdS)Nrrrr&r&r'rszNormalDist.variancecCs8t|tr&t|j|jt|j|jSt|j||jSr)rdrrrrrrr&r&r'__add__s zNormalDist.__add__cCs8t|tr&t|j|jt|j|jSt|j||jSr)rrr&r&r'__sub__'s zNormalDist.__sub__cCst|j||jt|Sr)rrrrrr&r&r'__mul__5szNormalDist.__mul__cCst|j||jt|Sr)rrr&r&r' __truediv__=szNormalDist.__truediv__cCst|j|jSr)rrrrr&r&r'__pos__EszNormalDist.__pos__cCst|j |jSr)rrr&r&r'__neg__IszNormalDist.__neg__cCs || Sr)r&rr&r&r'__rsub__OszNormalDist.__rsub__cCs&t|tstS|j|jko$|j|jkSr))rdrNotImplementedrrrr&r&r'__eq__Us zNormalDist.__eq__cCst|j|jfSr))hashrrrr&r&r'__hash__[szNormalDist.__hash__cCs t|jd|jd|jdS)Nz(mu=z, sigma=))r3r#rrrr&r&r'__repr___szNormalDist.__repr__)rr)rz) r#r$r% __slots__r classmethodrrrrrrrrpropertyrrr rrrrrrrr__radd__r__rmul__rrrr&r&r&r'rsF   "      r)r)rZ)r()N)N)N)N)N)7__all__rArerZ fractionsrZdecimalr itertoolsrbisectrrrrrrrrrr operatorr! collectionsr"rOrr?rDr1r/rRrVrYr\rrrrrr r rr r rrrrrr rZ _statistics ImportErrorrr&r&r&r'SsR   (   :  / 778  / ,  K