ó áp7]c@sBdZddlmZmZddlZdefd„ƒYZdS(sgBase classes for statistical test results Created on Mon Apr 22 14:03:21 2013 Author: Josef Perktold i˙˙˙˙(tlziptzipNtAllPairsResultscBsGeZdZdddd„Zdd„Zd„Zd„Zd„ZRS(sŸResults 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 : string 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 : None or list of strings 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). thscCsž||_||_|dkr7tj|ƒd|_n ||_||_||_|dkr„g|D]}d|f^qe|_n6g|D]&}d||d||df^q‹|_dS(Nis%rs%s-%si( t pvals_rawt all_pairstNonetnptmaxtn_levelstmultitest_methodtlevelstall_pairs_names(tselfRRR R R tpairs((s5lib/python2.7/site-packages/statsmodels/stats/base.pyt__init__'s       &cCsDddljj}|dkr*|j}n|j|jd|ƒdS(sşp-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. i˙˙˙˙Ntmethodi(tstatsmodels.stats.multitesttstatst multitestRR t multipletestsR(R Rtsmt((s5lib/python2.7/site-packages/statsmodels/stats/base.pytpval_corrected:s  cCs |jƒS(N(tsummary(R ((s5lib/python2.7/site-packages/statsmodels/stats/base.pyt__str__GscCs;|j}tj||fƒ}|jƒ|t|jŒ<|S(screate a (n_levels, n_levels) array with corrected p_values this needs to improve, similar to R pairwise output (R RtzerosRRR(R tkt pvals_mat((s5lib/python2.7/site-packages/statsmodels/stats/base.pyt pval_tableJs cCsddljj}td„|jDƒƒ}d|j|j}|dd|ddd 7}|d jd „t|j|j ƒƒDƒƒ7}|S( s•returns text summarizing the results uses the default pvalue correction of the instance stored in ``self.multitest_method`` i˙˙˙˙Ncss|]}t|ƒVqdS(N(tlen(t.0tss((s5lib/python2.7/site-packages/statsmodels/stats/base.pys \ss0Corrected p-values using %s p-value correction tPairst iis p-values s css%|]\}}d||fVqdS(s %s %6.4gN((RRtpv((s5lib/python2.7/site-packages/statsmodels/stats/base.pys as( RRRRR tmultitest_methods_namesR tjoinRR(R Rtmaxlevelttext((s5lib/python2.7/site-packages/statsmodels/stats/base.pyRUs N( t__name__t __module__t__doc__RRRRRR(((s5lib/python2.7/site-packages/statsmodels/stats/base.pyR s  (R)tstatsmodels.compat.pythonRRtnumpyRtobjectR(((s5lib/python2.7/site-packages/statsmodels/stats/base.pyts