diff --git a/george/modeling.py b/george/modeling.py
index 3879a1be..6fe441b6 100644
--- a/george/modeling.py
+++ b/george/modeling.py
@@ -120,10 +120,10 @@ def compute_gradient(self, *args, **kwargs):
         grad = np.empty([len(vector)] + list(value0.shape), dtype=np.float64)
         for i, v in enumerate(vector):
             vector[i] = v + _EPS
-            self.set_vector(vector)
+            self.set_parameter_vector(vector)
             value = self.get_value(*args, **kwargs)
             vector[i] = v
-            self.set_vector(vector)
+            self.set_parameter_vector(vector)
             grad[i] = (value - value0) / _EPS
         return grad