!gRdZddlmZddlZddlmZGddeZGddZy) zgBase classes for statistical test results Created on Mon Apr 22 14:03:21 2013 Author: Josef Perktold )lzipN)Holderc>eZdZdZdfd ZdZdZdZddZxZ S) HolderTuplez Holder class with indexing c tdi||tfd|D_yjjf_y)Nc36K|]}t|ywN)getattr).0attselfs S/var/www/dash_apps/app1/venv/lib/python3.12/site-packages/statsmodels/stats/base.py z'HolderTuple.__init__..sDcwtS1Ds)super__init__tuple statisticpvalue)r tuple_kwds __class__s` rrzHolderTuple.__init__s?  4  DVDDDJ..$++6DJc#8K|jEd{y7wr rr s r__iter__zHolderTuple.__iter__s::s c |j|Sr r)r idxs r __getitem__zHolderTuple.__getitem__szz#rc,t|jSr )lenrrs r__len__zHolderTuple.__len__s4::rcXtjt|j||S)N)dtypecopy)nparraylistr)r r%r&s r __array__zHolderTuple.__array__!sxxTZZ(DAArr )NT) __name__ __module__ __qualname____doc__rrr r#r* __classcell__)rs@rrr s"7Brrc6eZdZdZ ddZd dZdZdZdZy) AllPairsResultsaResults class for pairwise comparisons, based on p-values Parameters ---------- pvals_raw : array_like, 1-D p-values from a pairwise comparison test all_pairs : list of tuples list of indices, one pair for each comparison multitest_method : str method that is used by default for p-value correction. This is used as default by the methods like if the multiple-testing method is not specified as argument. levels : {list[str], None} optional names of the levels or groups n_levels : None or int If None, then the number of levels or groups is inferred from the other arguments. It can be explicitly specified, if the inferred number is incorrect. Notes ----- This class can also be used for other pairwise comparisons, for example comparing several treatments to a control (as in Dunnet's test). Nc 2||_||_|tj|dz|_n||_||_||_||Dcgc]}|c}|_y|Dcgc] }dj||d||d"c}|_ycc}wcc}w)Nz{}-{}r) pvals_raw all_pairsr'maxn_levelsmultitest_methodlevelsall_pairs_namesformat)r r4r5r8r9r7pairss rrzAllPairsResults.__init__@s""  FF9-1DM$DM 0 >*3$!&5) $D  )$58$fU1X&6$D $$s B!%Bcpddlmcm}| |j}|j |j |dS)zp-values corrected for multiple testing problem This uses the default p-value correction of the instance stored in ``self.multitest_method`` if method is None. rN)methodr3)statsmodels.stats.multiteststats multitestr8 multipletestsr4)r r>smts rpval_correctedzAllPairsResults.pval_correctedWs9 21 >**F   ?BBrc"|jSr )summaryrs r__str__zAllPairsResults.__str__ds||~rc|j}tj||f}|j|t |j <|S)zcreate a (n_levels, n_levels) array with corrected p_values this needs to improve, similar to R pairwise output )r7r'zerosrDrr5)r k pvals_mats r pval_tablezAllPairsResults.pval_tablegsA MMHHaV$ +/+>+>+@ $'(rc *ddlmcm}td|jD}d|j |j z}|dd|dz dzzzd zz }|d jd t|j|jDz }|S) zreturns text summarizing the results uses the default pvalue correction of the instance stored in ``self.multitest_method`` rNc32K|]}t|ywr )r")r sss rrz*AllPairsResults.summary..ys>2s2w>z0Corrected p-values using %s p-value correction Pairs r3z p-values  c32K|]\}}|d|dyw)z z6.4gNr)r r<pvs rrz*AllPairsResults.summary..~s'LKUBUG2bY/LrP) r?r@rAr6r:multitest_methods_namesr8joinziprD)r rCmaxleveltexts rrFzAllPairsResults.summaryrs 21>)=)=>>D--d.C.CDE 3(Q,"234|CC  Ld22D4G4G4IJLL L r)hsNNr ) r+r,r-r.rrDrGrLrFrrrr1r1%s*4?C'+. C rr1) r.statsmodels.compat.pythonrnumpyr'statsmodels.tools.testingrrr1rrrr`s/ +,B&B2[[r