yeah, i have most of the distribution ones done as we speak. again i am
worried about lightweight wrappers being important. for instance we
currently depend on chi_high from skbio.

chi_high(14, 5) = 0.01560941610026691

scipy.stats.chi2.cdf(14, 5) = 1 - 0.01560941610026691 = 0.98439058389973311

so my wrapper is just:

def chi_high(x, df):
   return 1 - chi2.cdf(14,5)

Are we thinking this is an okay way to do this, or ultimately that this is
the wrong way to go?


On Thu, May 15, 2014 at 5:55 PM, Greg Caporaso <notifications@github.com>wrote:

> In qiime/stats.py:
>
> > @@ -21,6 +22,12 @@
> >  create new statistical method implementations.
> >  """
> >
> > +from scipy.stats import (spearmanr, kruskal, mannwhitneyu, kendalltau,
> > +                         power_divergence)
> > +
> > +from skbio.math.stats.special import MACHEP, ndtri
>
> Would next Sunday 5/24 be a realistic goal to make this change and the
> changes for skbio.math.stats.distribution that I noted below?
>
> —
> Reply to this email directly or view it on GitHub<https://github.com/biocore/qiime/pull/1543/files#r12719863>
> .
>