diff --git a/pyhsmm/internals/hmm_states.py b/pyhsmm/internals/hmm_states.py
index e6ff2255bb9ec0a3391ae0ed23b542e67aabf3a6..2a5909317e89daf4de9a37d31f8490629730524c 100644
--- a/pyhsmm/internals/hmm_states.py
+++ b/pyhsmm/internals/hmm_states.py
@@ -5,7 +5,7 @@ import numpy as np
 from numpy import newaxis as na
 import abc
 import copy
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 
 from pyhsmm.util.stats import sample_discrete
 try:
diff --git a/pyhsmm/internals/hsmm_states.py b/pyhsmm/internals/hsmm_states.py
index 6a0505e8e27da8fb0ac30bfe73a4fe3228396461..a446b286d7fc468cf2a3662d9b168bb943993012 100644
--- a/pyhsmm/internals/hsmm_states.py
+++ b/pyhsmm/internals/hsmm_states.py
@@ -2,7 +2,7 @@ from __future__ import division
 from builtins import range, map
 import numpy as np
 from numpy import newaxis as na
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 
 from pyhsmm.util.stats import sample_discrete
 from pyhsmm.util.general import rle, rcumsum, cumsum
diff --git a/pyhsmm/models.py b/pyhsmm/models.py
index bf5a5ff29b8c113c6fb7c737773ea8dcb957827e..d7a246c4c09ecf7324ee417b904188a0f77ef85d 100644
--- a/pyhsmm/models.py
+++ b/pyhsmm/models.py
@@ -11,7 +11,7 @@ import matplotlib.pyplot as plt
 from matplotlib.gridspec import GridSpec, GridSpecFromSubplotSpec
 from matplotlib import cm
 from warnings import warn
-from scipy.misc import logsumexp
+from scipy.special import logsumexp
 from scipy.stats import invgamma, chi2
 
 from pyhsmm.basic.abstractions import Model, ModelGibbsSampling, \