diff --git a/GPy/likelihoods/likelihood.py b/GPy/likelihoods/likelihood.py
index bd84b72c..0f17a883 100644
--- a/GPy/likelihoods/likelihood.py
+++ b/GPy/likelihoods/likelihood.py
@@ -218,7 +218,7 @@ class Likelihood(Parameterized):
         #fi_samples = np.random.randn(num_samples)*np.sqrt(var_star) + mu_star
         fi_samples = np.random.normal(mu_star, np.sqrt(var_star), size=(mu_star.shape[0], num_samples))
 
-        from scipy.misc import logsumexp
+        from scipy.special import logsumexp
         log_p_ystar = -np.log(num_samples) + logsumexp(self.logpdf(fi_samples, y_test, Y_metadata=Y_metadata), axis=1)
         log_p_ystar = np.array(log_p_ystar).reshape(*y_test.shape)
         return log_p_ystar
diff --git a/setup.py b/setup.py
index 523c4433..4d1fa0b5 100644
--- a/setup.py
+++ b/setup.py
@@ -164,7 +164,7 @@ setup(name = 'GPy',
       py_modules = ['GPy.__init__'],
       test_suite = 'GPy.testing',
       setup_requires = ['numpy>=1.7'],
-      install_requires = ['numpy>=1.7', 'scipy>=0.16', 'six', 'paramz>=0.9.0'],
+      install_requires = ['numpy>=1.7', 'scipy>=1.3.0', 'six', 'paramz>=0.9.0'],
       extras_require = {'docs':['sphinx'],
                         'optional':['mpi4py',
                                     'ipython>=4.0.0',