diff --git a/GPy/models/gp_offset_regression.py b/GPy/models/gp_offset_regression.py
index 3699223e..2fd264f4 100644
--- a/GPy/models/gp_offset_regression.py
+++ b/GPy/models/gp_offset_regression.py
@@ -75,6 +75,10 @@ class GPOffsetRegression(GP):
             #print Gs[i]
             #print w
             dr_doffsets.append(dr_doffset)
+            
+        #lastly we need to divide by the lengthscale: So far we've found d(X_i - X_j)/dOffsets
+        #we want dr/dOffsets. (X_i - X_j)/lengthscale = r
+        dr_doffsets /= self.kern.lengthscale 
         return dr_doffsets
         
     def parameters_changed(self):