mgkit.counts.glm module

New in version 0.3.3.

GLM models with metagenomes and metatranscriptomes. Experimental

mgkit.counts.glm.fit_lowess_interpolate(endog, exog, frac=0.2, it=3, kind='slinear')
mgkit.counts.glm.lowess_ci_bootstrap(endog, exog, num=100, frac=0.2, it=3, alpha=0.05, delta=0.0, min_alpha=0.001, kind='slinear')

Performance increase with the value of delta.

mgkit.counts.glm.variance_to_alpha(mu, func, min_alpha=0.001)

Based on the variance defined in the Negative Binomial in statsmodels

var = mu + alpha * (mu ** 2)