diff --git a/pybasicbayes/distributions/negativebinomial.py b/pybasicbayes/distributions/negativebinomial.py
index 26958a1ccd67652056b339e267da217930e53cda..7caed4edf798db3f57462c97a01fccc262a113c0 100644
--- a/pybasicbayes/distributions/negativebinomial.py
+++ b/pybasicbayes/distributions/negativebinomial.py
@@ -11,7 +11,7 @@ __all__ = [
 import numpy as np
 from numpy import newaxis as na
 import scipy.special as special
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 from warnings import warn
 
 from pybasicbayes.abstractions import Distribution, GibbsSampling, \
diff --git a/pybasicbayes/models/mixture.py b/pybasicbayes/models/mixture.py
index 43461c4a7c78ae06921d73ac8c2d4f188700311d..55292c644336f0dcb8c3185f9a0a03bb1d6c5730 100644
--- a/pybasicbayes/models/mixture.py
+++ b/pybasicbayes/models/mixture.py
@@ -10,7 +10,7 @@ na = np.newaxis
 import scipy.special as special
 import abc, copy
 from warnings import warn
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 
 from pybasicbayes.abstractions import ModelGibbsSampling, ModelMeanField, ModelEM
 from pybasicbayes.abstractions import Distribution, GibbsSampling, MeanField, Collapsed, \
diff --git a/pybasicbayes/util/stats.py b/pybasicbayes/util/stats.py
index 9ea9b479a449c970aaffc7cab0e6f1a124461c6d..7a133d4ff2725023f9f61e92a52aa95498ddb402 100644
--- a/pybasicbayes/util/stats.py
+++ b/pybasicbayes/util/stats.py
@@ -7,7 +7,7 @@ na = np.newaxis
 import scipy.stats as stats
 import scipy.special as special
 import scipy.linalg
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 from numpy.core.umath_tests import inner1d
 
 from .general import any_none, blockarray