factorial, comb and logsumexp have been moved: trying to import them
from the old location is an error in scipy 1.3+.
https://ci.debian.net/data/autopkgtest/unstable/amd64/s/statsmodels/
2949940/log.gz
Author: Rebecca N. Palmer <rebecca_palmer@zoho.com>
Forwarded: accepted https://github.com/statsmodels/statsmodels/pull/6164
Gbp-Pq: Name scipy13_compat.patch
stats.binom(5, 1. / 6).pmf(2)
-from scipy.misc import comb
+from scipy.special import comb
comb(5, 2) * (1 / 6.)**2 * (5 / 6.)**3
from statsmodels.formula.api import glm
>>> import matplotlib.pyplot as plt
>>> from scipy import stats
- >>> from scipy.misc import factorial
+ >>> from scipy.special import factorial
>>> df = 12
>>> chi2_c = [2**(j-1) * factorial(j-1) * df for j in range(1, 5)]
>>> edgw_chi2 = ExpandedNormal(chi2_c, name='edgw_chi2', momtype=0)
from scipy import stats
import numpy as np
from matplotlib import pyplot as plt
-from scipy.misc import logsumexp as sp_logsumexp
+from statsmodels.compat.scipy import logsumexp as sp_logsumexp
#TODO: change these to use maxentutils so that over/underflow is handled
#with the logsumexp.