--- nipy.orig/nipy/algorithms/statistics/rft.py
+++ nipy/nipy/algorithms/statistics/rft.py
@@ -20,8 +20,7 @@
 from numpy.linalg import pinv
 
 from scipy import stats
-from scipy.misc import factorial
-from scipy.special import gamma, gammaln, beta, hermitenorm
+from scipy.special import gamma, gammaln, beta, hermitenorm, factorial
 
 # Legacy repr printing from numpy.
 from nipy.testing import legacy_printing as setup_module  # noqa
--- nipy.orig/nipy/algorithms/statistics/tests/test_rft.py
+++ nipy/nipy/algorithms/statistics/tests/test_rft.py
@@ -4,9 +4,8 @@
 
 import numpy as np
 
-from scipy.special import gammaln, hermitenorm
+from scipy.special import gammaln, hermitenorm, factorial
 import scipy.stats
-from scipy.misc import factorial
 
 from .. import rft
 
--- nipy.orig/nipy/labs/group/permutation_test.py
+++ nipy/nipy/labs/group/permutation_test.py
@@ -7,6 +7,7 @@
 # Third-party imports
 import numpy as np
 import scipy.misc as sm
+from scipy.special import comb
 import warnings
 
 # Our own imports
@@ -374,7 +375,7 @@
         elif self.nsamples == 2:
             n1,p = self.data1.shape[self.axis], self.data1.shape[1-self.axis]
             n2 = self.data2.shape[self.axis]
-            max_nperms = sm.comb(n1+n2,n1,exact=1)
+            max_nperms = comb(n1+n2,n1,exact=1)
             data = np.concatenate((self.data1,self.data2), self.axis)
             if self.vardata1 is not None:
                 vardata = np.concatenate((self.vardata1,self.vardata2), self.axis)
