diff --git a/GPy/kern/src/ODE_UYC.py b/GPy/kern/src/ODE_UYC.py
index d02eb1d9..57c41767 100644
--- a/GPy/kern/src/ODE_UYC.py
+++ b/GPy/kern/src/ODE_UYC.py
@@ -18,7 +18,7 @@ class ODE_UYC(Kern):
         self.lengthscale_U = Param('lengthscale_U', lengthscale_U, Logexp())
         self.ubias = Param('ubias', ubias, Logexp())
 
-        self.add_parameters(self.variance_Y, self.variance_U, self.lengthscale_Y, self.lengthscale_U, self.ubias)
+        self.link_parameters(self.variance_Y, self.variance_U, self.lengthscale_Y, self.lengthscale_U, self.ubias)
 
     def K(self, X, X2=None):
         # model :   a * dy/dt + b * y = U
diff --git a/GPy/kern/src/ODE_st.py b/GPy/kern/src/ODE_st.py
index f9d4e684..0b4fecae 100644
--- a/GPy/kern/src/ODE_st.py
+++ b/GPy/kern/src/ODE_st.py
@@ -38,7 +38,7 @@ class ODE_st(Kern):
         self.b = Param('b', b, Logexp())
         self.c = Param('c', c, Logexp())
 
-        self.add_parameters(self.a, self.b, self.c, self.variance_Yt, self.variance_Yx, self.lengthscale_Yt,self.lengthscale_Yx)
+        self.link_parameters(self.a, self.b, self.c, self.variance_Yt, self.variance_Yx, self.lengthscale_Yt,self.lengthscale_Yx)
 
 
     def K(self, X, X2=None):        
diff --git a/GPy/kern/src/ODE_t.py b/GPy/kern/src/ODE_t.py
index ffd349ec..d5dae665 100644
--- a/GPy/kern/src/ODE_t.py
+++ b/GPy/kern/src/ODE_t.py
@@ -17,7 +17,7 @@ class ODE_t(Kern):
                 self.a= Param('a', a, Logexp())
                 self.c = Param('c', c, Logexp())
                 self.ubias = Param('ubias', ubias, Logexp())
-                self.add_parameters(self.a, self.c, self.variance_Yt, self.lengthscale_Yt,self.ubias)
+                self.link_parameters(self.a, self.c, self.variance_Yt, self.lengthscale_Yt,self.ubias)
 
         def K(self, X, X2=None):
                 """Compute the covariance matrix between X and X2."""