U mdYg@sdZddlmZddlZddlZddlZddlm Z ej Z dZ dddd d gd d d ddgdddddgdddddgdddddgdd d!d"d#gd$d%d&d'd(gd)d*d+d,d-gd.d/d0d1d2gd3d4d5d6d7gd8d9d:d;dd?d@dAgdBdCdDdEdFgdGdHdIdJdKgdLdMdNdOdPgdQdRdSdTdUgdVdWdXdYdZgd[d\d]d^d_gd`dadbdcddgdedfdgdhdigdjdkdldmdngdodpdqdrdsgdtdudvdwdxgdydzd{d|d}gd~e fddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgde fddddgdddddgddddd gd d d d dgdddddgdddddgdddddgddd d!d"gd#d$d%d&d'gd(d)d*d+d,gd-d.d/d0d1gd2d3d4d5d6gd7d8d9d:d;gd<d=d>d?d@gdAdBdCdDdEgdFdGdHdIdJgdKdLdMdNdOgdPdQdRdSdTgdUdVdWdXdYgdZd[d\d]d^gd_d`dadbdcgdddedfdgdhgdidjdkdldmgdndodpdqdrgdsdtdudvdwgdxe fdydzd{d|gd}d~dddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdÐdĐdŐdƐdgdȐdɐdʐdːdgd͐dΐdϐdАdgdҐdӐdԐdՐdgdאdؐdِdڐdgdܐdݐdސdߐdgdddddgdddddgdddddgdddddgde fddddgdddddgdddddgdddddgd d d d d gdddddgdddddgdddddgdddd d!gd"d#d$d%d&gd'd(d)d*d+gd,d-d.d/d0gd1d2d3d4d5gd6d7d8d9d:gd;d<d=d>d?gd@dAdBdCdDgdEdFdGdHdIgdJdKdLdMdNgdOdPdQdRdSgdTdUdVdWdXgdYdZd[d\d]gd^d_d`dadbgdcdddedfdggdhdidjdkdlgdmdndodpdqgdre fdsdtdudvgdwdxdydzd{gd|d}d~ddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgddÐdĐdŐdgdǐdȐdɐdʐdgd̐d͐dΐdϐdgdѐdҐdӐdԐdgd֐dאdؐdِdgdېdܐdݐdސdgdddddgdddddgdddddgde fddddgdddddgdddddgdddddgdddddgdd d d d gd ddddgdddddgdddddgddddd gd!d"d#d$d%gd&d'd(d)d*gd+d,d-d.d/gd0d1d2d3d4gd5d6d7d8d9gd:d;d<d=d>gd?d@dAdBdCgdDdEdFdGdHgdIdJdKdLdMgdNdOdPdQdRgdSdTdUdVdWgdXdYdZd[d\gd]d^d_d`dagdbdcdddedfgdgdhdidjdkgdldmdndodpgdqe fdrdsdtdugdvdwdxdydzgd{d|d}d~dgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddÐdĐdgdƐdǐdȐdɐdgdːd̐d͐dΐdgdАdѐdҐdӐdgdՐd֐dאdؐdgdڐdېdܐdݐdgdߐddddgdddddgdddddgdddddgde fddddgdddddgdddddgdddddgddd d d gd d dddgdddddgdddddgdddddgd d!d"d#d$gd%d&d'd(d)gd*d+d,d-d.gd/d0d1d2d3gd4d5d6d7d8gd9d:d;d<d=gd>d?d@dAdBgdCdDdEdFdGgdHdIdJdKdLgdMdNdOdPdQgdRdSdTdUdVgdWdXdYdZd[gd\d]d^d_d`gdadbdcdddegdfdgdhdidjgdkdldmdndogdpdqdrdsdtgdue fdvdwdxdygdzd{d|d}d~gdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgddddÐdgdŐdƐdǐdȐdgdʐdːd̐d͐dgdϐdАdѐdҐdgdԐdՐd֐dאdgdِdڐdېdܐdgdސdߐdddgdddddgdddddgdddddgdddddgde fddddgdddddgdddddgdddd d gd d d ddgdddddgdddddgdddddgdd d!d"d#gd$d%d&d'd(gd)d*d+d,d-gd.d/d0d1d2gd3d4d5d6d7gd8d9d:d;dd?d@dAgdBdCdDdEdFgdGdHdIdJdKgdLdMdNdOdPgdQdRdSdTdUgdVdWdXdYdZgd[d\d]d^d_gd`dadbdcddgdedfdgdhdigdjdkdldmdngdodpdqdrdsgdtdudvdwdxgdye fdzd{d|d}gd~ddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdddddgdĐdŐdƐdǐdgdɐdʐdːd̐dgdΐdϐdАdѐdgdӐdԐdՐd֐dgdؐdِdڐdېdgdݐdސdߐddgdddddgdddddgdddddgdddddgdddddgde fddddgi2a d~ddxddrddqdduddydg a eddddddde gaddZd d Zd d Zd dZddZddZddZddZddZeeZde_ddZddZddZeeZd e_d!d"ZdS(#a# Implementation of Gleason's (1999) non-iterative upper quantile studentized range approximation. According to Gleason this method should be more accurate than the AS190 FORTRAN algorithm of Lund and Lund (1983) and works from .5 <= p <= .999 (The AS190 only works from .9 <= p <= .99). It is more efficient then the Copenhaver & Holland (1988) algorithm (used by the _qtukey_ R function) although it requires storing the A table in memory. (q distribution) approximations in Python. see: Gleason, J. R. (1999). An accurate, non-iterative approximation for studentized range quantiles. Computational Statistics & Data Analysis, (31), 147-158. Gleason, J. R. (1998). A table of quantile points of the Studentized range distribution. http://www.stata.com/stb/stb46/dm64/sturng.pdf )lrangeN) fminboundz0.2.3)皙?@gw g8g"qH?g$/)r@gfOg2Vvg,/F?g{ =گ)r@g9oYgW'g8.I?g;y )r@gs4<)g襹Ngޡf*X?g@gүM=g!&g[{zp?g'tT)rD@gtN ga/t g棝?g)rN@gX0ÈgU=Yg&7?grMld)r^@geb.g7gF]Z\?g@)rg%Tg%#!ҎgV>9`?g]cI )?rgu3AgtTgR5?g hg)rrg^?gOW,Ag 5?gfOio3f)rrg&SEgSA%\gG^?g[Aze)rrg'tgrc渿gB?g/ :e)rr g7{yg%"g}?gUVe)rr g(|H^ggݸgSu*?gFcXe)rr gVogЈgwe)rr gg(rf g?gid)rrg9ȷgvH?gSr(e)皙?rg80gY0xڸgqqϣũ?g&/a)rrg1,g X۴g`b2?gkcu\)rrg")g8Hgla=?gZ%IX)rrg^0ŵgrg8Qf?ge S4VV)rr glvgPjg|[Y?g.@XIU)rr gu^g̬gD+ k?g=n>ST)rr glgqGQg?G[ސ?g齥S)rr ggKgVa?gn`Ѫ(S)rr gZg, &g͵;?g R)rrgbASg%Ug_P/g?g7dQR)rrgPbcgLg? p?gP\R)rrgz(g[6g`Hŏ?grg4#̡gq:?g~woSrga,mg+Nug8>uE?g5`BT)?rg8C=kgj gw۸4?g2=ȱZ`)rrg9gK>.g'oʕ?ga˓Y)rrg8(.gD챎gİ?gʑ*9T)rrgVDg D}#gm(n?g˟R)rr g:g#NG4g֨?gcǴQDtP)rr g6Wg~sg; ?gvq,ۡN)rr go*69gCIJ2gi?glM)rr ghgT~ǥg d?g-\K)rr g>g͋\ğgtR̀?g& k K)rrgYIܙg,^gV %͇?g o_tJ)rrgfg"eܢgן??g %J)rrgRڏgف-g%͆?g[I)rrg%gNg \Cg)q?g'<~I)rrgkS!gDyZgv;&?gSXI)rrgARͅgZ9gۼ{?gX~?I)rrgQ >gMS Dg+3k?gRġa1I)rrgg.0-ڟg "-?gFKX+I)rrgIkb~gll?g}^{k?g$A-+I)rrgc}glv`xNghO?gs6/I)rrg7wgIN߶ g g?gۆ_I)rrgpg [ךzgEbn?g0 *I)rrg `ig.Cgý ?g_GqjJ)rrgmbg+^gW%xM?gBP9iNK)rrgPӗZg^:gm?gwTaLrg.Qg-sg:?gBX8JfM)皙?rg`4gc֨g|o?gtZ_)r rg63[HgܳgICP?gh)5UW)r rgבgAS g;K?g R)r rg$%@DIg~Eݬg C?gvO)r r gPtg|gwkv?gHK)r r gx߿g@gܒ,|?g{nzI)r r g*߿gGXgօچ?gvtG)r r gc[߿gУge8?glIXF)r r g'H߿g]&gTʭ?gmvE)r rg,k߿g򇉡g ?gE)r rgNj]o߿g ."g‹חQ?geԑƍD)r rgL+a߿gM_2៿giV(oׂ?gf=:D)r rg TCU߿gg럞gu?gK?C)r rg1<꣟K߿gcWgLeV%?gbo7C)r rgx|B߿g5d{g)t-?gyC)r rgp:߿gWSLƛg嚮?g[#YC)r rgd̨3߿gÏ7 gEBN?gW >C)r rg+y-߿gMqfg,s]?g.C)r rg)i"(߿g̺,әgQCXu??g-C)r rgHG߿g/L gER7?g[yEEC)r rg0t^߿gm:3YgPŀ?gp~g@7@?gFEK)r!r gGܿg0!gl0ԇ?g@G)r!r g d̂ܿg|1g(ba}?g%E)r!r gU`Gܿgz8 gqWNƒ?gyvaTC)r!r guܿgL^g?u?s?g B)r!r g羀ۿgzj g1+_yq?g:4 A)r!rgۿgE쟿g 5T?g"wx@)r!rgJۿg:ugy݀ ?gɢb@)r!rg jۿgb)eg+uQ?gqV?)r!rgg[߱ۿg>gmI~?g!B>)r!rgۿg(gJ]}?gLTT>)r!rgVDۿgzoJg7S &}?g}[>)r!rg?{ۿg)@g j}>|?g>)r!rgtzbnqۿg'W}ghsQ|?g|$i1>)r!rgd'hۿg.{gЖg⴮|?gFû0>)r!rg_`ۿg6kD$6g +{?gk9>)r!rgel CGۿg>kXgJ%{?g6t>)r!rgLE-ۿg1zggTz?g zY?)r!rgI$kۿg6䐿g֥ibz?gYksE@)r!rgDzڿg*x{0gdh z?ggjA)r!rg]gXڿgmV⋿gAf{i{?guAr!gy_mdڿgo.h,?gz* |?g醁_B)??gg<8g`?g]ng)r"rgy giDʧԸg#Q?g .dT])r"rgwvܿgsy³g8l?gà`T)r"rgO=/*ڿg/)r"r g^w+qؿg- ig ?gg =)r"r go@ؿgR gP0S}?g;1D;)r"rg&ؿg)xYxgن]j{?g2 8)r"rgxN# ؿgzg͈`Fz?gqA7)r"rgr9|Ѵ׿gy gSsLy?g a7)r"rgBdF׿gNቘg7V3Fx?g%$C,6)r"rgO׿gogT=gzߵÕw?g쎁A6)r"rgd7Y׿ge2$–gVwSiw?g036)r"rg^VD׿g&égR2v?gtqw6)r"rg<ե[э׿gHjӔg`Y"cv?giy5)r"rg׿gu\xՓg 1u?gu.`5)r"rguxs׿g̷ΓgPnNv?gG6)r"rgcCT׿gng,ou?g2,f6)r"rgHʖ3׿gXOtg"Ut?gajP7)r"rg:, ׿gO۴gLy[WEu?gFm9)r"rg.ֿgH " g.ȳ+t?g C?9)r"rg#Z&wֿgCW-g'b\u?g16GƚgN'z?g񟥚u5)r$r gJmԿgSK)՝gw?gG2)r$rg쨳ӿg}W~gTu?g?ެ=$1)r$rgMS'Lӿg CA酪gP3Ct?g[~0e0)r$rg6ӿgLog5?kRt?gAȤ0)r$rgNh-|ӿgvgæ=r?g 8-)r$rg:dZB^ӿg]OCэgxFnܙ`r?g˖vd.)r$rgdSԔFӿg\4⤔g< ^r?g~;T/)r$rg$6z6ӿg).jg-7֖0,q?g;1-)r$rg@ӿg[1gAnq?gK'$/)r$rgӿg) '(.gE1np?g( #-)r$rgbߝҿgi2]÷gwhep?g ]W/)r$rgJ]iҿg[@fgb\j*p?gG8v/)r$rg!ҿgln,FZg"p?gRZ1)r$rg\ rҿg=\g7UKo?gY>1)r$rg.7Aҿg ƾg@o?g M2)r$rgo'ҿge Rg۩o?g:i 3r$gʣѿggeVq?g[ؼ6)333333?r#gߐzg>#gu3H?gZH[g)r%rgݿgIҸg#"?g35\)r%rgCb}ؿgS>gґ?gS)r%rg^տgb[Kgi<?g ;arhK)r%rg?{sVԿgY۹wg=?gfPAA)r%r g0oQ^9PӿgxIgN_ _?gQR<)r%r gpbURҿgc LWg <0t~?gf5|8)r%r g9ҿg񮉣̡gdy?g'N:3)r%r g9bFIѿgV7ݢgķrv?gs()0)r%r gK"npѿgK' gVqt?g͊:,)r%rg_.ѿg3ƚgeIbs?gdYBP*)r%rg Dпgi gIB r?g&)bv*)r%rg#пgo:.Ug\ҕp?gƤ()r%rg (пgRXЕgf# p?gXNz]')r%rgEuFпg\4Yeg?.|o?g g]()r%rg*0=Hpпg4J|WgZfm m?gG%)r%rgRпg8ؼ͒gP+5m?g^x&)r%rg`tdA5пg;Զlgju]Dm?g~m()r%rgH&пg [gsk?gL%\&)r%rgA пg3qͿg݁C}{g@Rj?gF31)Gz?r#gv&lgKbh6g7ʠex?gu#|g)r&rgXݿg{Dg@@6?g>Xw\)r&rgȐL׿g#eIg,=(D?gS)r&rg/YkԿg\xZCg@KĈ?g( F)r&rgÁҿgtM g뻁?g÷{wB)r&r gvѿg`Eåg==}?g.Ō4)r&r gpۙ~пg/X**g\O|?g8nv|5)r&r gϿg$5_Rg:Lv?gRY-)r&r gP* Ͽgf-2g3/r?g]*k")r&r gŨGοgP-g>a! q?g@")r&rg5 Ϳgc!fgom;= o?gb^* )r&rgeͿg 9gؑo?gbPCa#)r&rgo\\̿g{{5gOrk?gD&M )r&rgC榷i̿gbg~ag i5i?gWW8)r&rg ̿g,~ʓg@!i?gd޶7 )r&rgN$˿gɜrU-gMQGh?g&)r&rg=w&˿gsݑgO90g?gO#}Tv)r&rgLe?gx!)r&rg{6ɿg>ig4,c?d?g*bE#)r&rgZqbɿgggfQXIe?gvqm &)r&rgC 5ȿg/tg*-e?ggd$Jd?gPo)r'rgDOS@ȿgCJcgMB[jc?gx&k )r'rgh_mϛǿgYi3g`"d?g. )r'rg9n ǿgc;gfb?gp'wn )r'rg*wzƿg> D̈́gPNb?gG!)r'rg` JſgbzgjAo2c?gϺWI%)r'rg^ 2ZOſg11؆~g3!c?g'Uo&r'gXĿgUX!)xg~Ȅb?g8m%&)+?r#gg&ggU^`bX?g[/KRg)r(rgܿg.*)g~t?gfeK\)r(rgp[TֿgN=g?gdAwQ)r(rgK @PҿgwC}g#?gw̃G)r(rg=vпgc_gFFi?gd9A)r(r gnMοg@ng-$y?gvBD,)r(r g˿g5-gbv?g K*)r(r gԑWʿg7gwgs?g P[#)r(r gusɿgbBݝg_ޥp?gl )r(r g;\Yg+#ȿgӉU9g{'Y8hm?gSP)V)r(rgY 6Nǿg;`+Șg(,!Zi?g E)r(rgh0ƿgVBg+xc?g[} )r(rgzb r<ƿgWz* g6-Wb?gO#/D)r(rgi0ſg400$gnre?gb")r(rg,,Kſgz]Ugk$Db?gk(= )r(rg!8FĿg=qg^)`?g%')r(rgĿgyI=bgE]M``?g:GJB )r(rg>SĿgf6O5gKeca?g3")r(rguĿgA[gێb?g<¿g<@gE4^?gǣRF)r(rgg僿gřX̴]?gg~I՞f)r(rg,R.gd~ӀggrL]?g7)r(rg ́g\\|gsC^?g{ r(gN!LgWNwg2Z]?guyb!(<xcCs z t|WnYdSXdS)zC returns True if x is a float, returns False otherwise FT)float)xr2^/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/statsmodels/stats/libqsturng/qsturng_.py_isfloatxs  r4c Cs|dks|dkrtd|d}d}d}d}d}d|}||krtd t|}|d||d||d ||d ||d ||d  |d||d||d ||d |dS||krhtd td|}|d||d||d ||d ||d ||d |d||d||d ||d |dS|d}||}|d||d||d ||d ||d ||d  ||d||d||d ||d ||d |dS)a Modified from the author's original perl code (original comments follow below) by dfield@yahoo-inc.com. May 3, 2004. Lower tail quantile for standard normal distribution function. This function returns an approximation of the inverse cumulative standard normal distribution function. I.e., given P, it returns an approximation to the X satisfying P = Pr{Z <= X} where Z is a random variable from the standard normal distribution. The algorithm uses a minimax approximation by rational functions and the result has a relative error whose absolute value is less than 1.15e-9. Author: Peter John Acklam Time-stamp: 2000-07-19 18:26:14 E-mail: pjacklam@online.no WWW URL: http://home.online.no/~pjacklam rz5Argument to ltqnorm %f must be in open interval (0,1))g%1Cg4pFk@g;->qg@ rKa@g͋40>gyTW @)g0rEcCsdddtd|dS)z+function for p-value abcissa transformationgr#g?r)rEr>r2r2r3 _ptransformsrGcCs|dt|d|dt|dd|dt|dd|dt|dd}|dkr|dddt|d7}|d krt|s|nd }|d dd |7}n"t|s|nd }|dd |7}| S)z| calculates f-hat for the coefficients in a, probability p, sample mean difference r, and degrees of freedom v. rr#r5r)r7r8gMb`rg~jt@籡*Gg9ư_?gs@gg@)r;r=rEnpisinf)r?r>rDvfr2r2r3_funcsrMcCst|dkr dS|dkrdS|dkr$dS|dkr0dS|d krrDrKp0p1p2Zy0y1y2Zy_log0Zy_log1Zy_log2Zp_tZp0_tZp1_tZp2_td2d1d0Zy_logyq0q1rCr2r2r3_interpolate_psf   " $ " $ (rbcCs|dkrddtfS|dkrdS|dkr*dS|dkr6d S|d krBd S|d krNd S|dkrd|dkrpdSn |dkrpdStt|}|d||dfS)z-returns the points to use for interpolating vrr.r/r)r-r.r/r)r,r-r.r)r+r,r-r)r+r,g3@)rcr+r"g@)r5r)r7g @)r)r7r8r5)infintround)rKr>vir2r2r3 _select_vsWs&  ricCsJt||\}}}tt||f|||dd}tt||f|||dd}tt||f|||dd}|dkrvd}d|d|d|d|f\} } } } d||| | ||| | | | } | | | | kr||| | d| | | }n ||| | d| | | }|}t| d| | d|| | |}|S)z` interpolates v based on the values in the A table for the scalar value of r and th r#rrHr)rirMrOr;r<)r>rDrKv0v1v2Zy0_sqZy1_sqZy2_sqv_v0_v1_v2_r\r]r^r_r2r2r3_interpolate_vss$ $" *rqcCsB|dks|dkrtd|dkr2|dkrBtdn|dkrBtdt|}t|tjr^|}||ftkrtt||f|||d }n|tkr|t gdgf|dkkrt ||\}}}t |\}}} t |||d} t |||d} t |||d} d |d |d |d |f\} }}}d | | ||| | ||||}||||kr| | ||d |||}n | | ||d |||}| }t |d | |d || ||}n@|t gdgf|dkkrt|||}n|tkrt |||}t|d }t d| tjjd |d |S) zscalar version of qsturngrr(zp must be between .1 and .999r"r)zv must be > 2 when p < .9r5zv must be > 1 when p >= .9r#rrrH)r:r0 isinstancerIZndarrayitemrOrMp_keysv_keysrirNrbr;r<rqrRrSrTrUrV)r>rDrKr_rjrkrlrWrXrYZr0_sqZr1_sqZr2_sqrmrnrorpr\r]r^r2r2r3_qsturngsF   &$" ,   rvzvector version of qsturngcCs,ttt|||gr t|||St|||S)aApproximates the quantile p for a studentized range distribution having v degrees of freedom and r samples for probability p. Parameters ---------- p : (scalar, array_like) The cumulative probability value p >= .1 and p <=.999 (values under .5 are not recommended) r : (scalar, array_like) The number of samples r >= 2 and r <= 200 (values over 200 are permitted but not recommended) v : (scalar, array_like) The sample degrees of freedom if p >= .9: v >=1 and v >= inf else: v >=2 and v >= inf Returns ------- q : (scalar, array_like) approximation of the Studentized Range )allmapr4rv _vqsturngr>rDrKr2r2r3qsturngs r{csdkrtdfdd}|dkrntd|dkr8dStd|dkrLd Sd t|dd||fd }t|Std||krdStd||krd Sd t|dd||fd }t|Sd S) zscalar version of psturnggzq should be >= 0csttt|||SN)rIsqueezeabsrvrzrCr2r3opt_func8sz_psturng..opt_funcr5r"rr(gMbP?r#)argsN)r:rvrrIZ atleast_1d)rCrDrKrZsolnr2rr3_psturng3s   rcCstt|||Sr|)rIr}rrCrDrKr2r2r3_psturng_scalarJsrzvector version of psturngcCs,ttt|||gr t|||St|||S)aJEvaluates the probability from 0 to q for a studentized range having v degrees of freedom and r samples. Parameters ---------- q : (scalar, array_like) quantile value of Studentized Range q >= 0. r : (scalar, array_like) The number of samples r >= 2 and r <= 200 (values over 200 are permitted but not recommended) v : (scalar, array_like) The sample degrees of freedom if p >= .9: v >=1 and v >= inf else: v >=2 and v >= inf Returns ------- p : (scalar, array_like) 1. - area from zero to q under the Studentized Range distribution. When v == 1, p is bound between .001 and .1, when v > 1, p is bound between .001 and .9. Values between .5 and .9 are 1st order appoximations. )rwrxr4r _vpsturngrr2r2r3psturngPs r)__doc__Zstatsmodels.compat.pythonrr;Z scipy.statsrSnumpyrIZscipy.optimizerre __version__rOrtrur4rErGrMrNrbrirqrvZ vectorizeryr{rrrrr2r2r2r3s                                                    60(  @    T  * Z  9