!g(HdZddlZddlZddlmZddlmZmZm Z ddZ ddZ y) zb Created on Sun Sep 25 21:23:38 2011 Author: Josef Perktold and Scipy developers License : BSD-3 N)stats) array_like bool_likeint_likect|dd}t|d}t|d}tj||}|j |}|r|dk(rtj tj|||}tj tj|d| |}||z |z } tjj| } net|r%|j|}|j|g|} n5td t|r|j|g|} n td tjd|dz} dg|j z} t#d| |<t%| } t#dg|j z} t#ddd | |<t%| } t'j(5t'j*d dd| | zdz |z tj,| tj.| |  zz}tj0||}ddd| z }|S#1swYxYw)a Calculate the Anderson-Darling a2 statistic. Parameters ---------- x : array_like The data to test. dist : {'norm', callable} The assumed distribution under the null of test statistic. fit : bool If True, then the distribution parameters are estimated. Currently only for 1d data x, except in case dist='norm'. params : tuple The optional distribution parameters if fit is False. axis : int If dist is 'norm' or fit is False, then data can be an n-dimensional and axis specifies the axis of a variable. Returns ------- {float, ndarray} The Anderson-Darling statistic. xN)ndimfitaxis)r norm)ddofr z!dist must be 'norm' or a Callablez+if fit is false, then dist must be callableignorez#divide by zero encountered in log1p)messageg?)rrrnpsortshape expand_dimsmeanstdrr cdfcallabler ValueErroraranger slicetuplewarningscatch_warningsfilterwarningsloglog1psum)rdistr paramsr ynobsxbarswzisl1sl2 ad_valuesa2s V/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/stats/_adnorm.pyanderson_statisticr3s0 1c%A C C D& !D A 774=D 6>>>"''!$"7>Drvvaad;TBATQA q!A d^XXa[F$V$A@A A D>$V$AJK K !TAXA &166/Cd CI *C ;-!&& CdD"%CI *C  ") C 3Z#%-RXXqvg=N1NO FF94 ( ) B I))s A-H??IcJt|dd|}|j|}|dd|z zd|dzz zz}tj|dk(r|dk\r0|d kr+dtjd d |zzd |dzzz z }||fS|d kr+dtjdd|zzd|dzzz z }||fS|dkr(tjdd|zz d|dzzz }||fS|dkr(tjdd|zz d|dzzz}||fSd}||fStj gd}d}d}d} d} d} ||| | | g} tj ||d } tjtj|z}td!D]}| |k(}| |||||<||fS)"a Anderson-Darling test for normal distribution unknown mean and variance. Parameters ---------- x : array_like The data array. axis : int The axis to perform the test along. Returns ------- ad2 : float Anderson Darling test statistic. pval : float The pvalue for hypothesis that the data comes from a normal distribution with unknown mean and variance. See Also -------- statsmodels.stats.diagnostic.anderson_statistic The Anderson-Darling a2 statistic. statsmodels.stats.diagnostic.kstest_fit Kolmogorov-Smirnov test with estimated parameters for Normal or Exponential distributions. r T)r%r r r g?g@r皙?Zd;*)\HY@(\k@(\?A`Т SeE@MbM@333333? Ac]?V-@Gz? \m?t@Έ ?)r5r6r:r>cNtjtj|zS)N)rnan ones_likead2as r2znormal_ad..}sRVVbll4&88cPdtjdd|zzd|dzzz z S)Nr r7r8r9rrexprIs r2rKznormal_ad..~s0Q ftm #ftqy&8 8"::rLcPdtjdd|zzd|dzzz z S)Nr r;r<r=rrNrIs r2rKznormal_ad..s0Q Vd] "Vdai%7 7"99rLcJtjdd|zz d|dzzz S)Nr?r@rArrNrIs r2rKznormal_ad..s&RVVFUT\$9D419.s&RVVFUT\$9FTQYrls& HH<~?rL