From 9f0c2ea1b7578e05db7a7642d676308f937b2771 Mon Sep 17 00:00:00 2001
From: Alan Saul <alan.daniel.saul@gmail.com>
Date: Tue, 17 Dec 2013 15:08:27 +0000
Subject: [PATCH] Changed initalise_latent to take lower of init=PCA and
 corrected import

---
 GPy/models_modules/gplvm.py        | 6 +++---
 GPy/models_modules/mrd.py          | 4 ++--
 GPy/models_modules/sparse_gplvm.py | 3 ---
 3 files changed, 5 insertions(+), 8 deletions(-)

diff --git a/GPy/models_modules/gplvm.py b/GPy/models_modules/gplvm.py
index 541b3176..da82540d 100644
--- a/GPy/models_modules/gplvm.py
+++ b/GPy/models_modules/gplvm.py
@@ -12,7 +12,7 @@ from .. import util
 
 def initialise_latent(init, input_dim, Y):
     Xr = np.random.randn(Y.shape[0], input_dim)
-    if init == 'pca':
+    if init.lower() == 'pca':
         from ..util.linalg import pca
         PC = pca(Y, input_dim)[0]
         Xr[:PC.shape[0], :PC.shape[1]] = PC
@@ -30,7 +30,7 @@ class GPLVM(GP):
     :type init: 'pca'|'random'
 
     """
-    def __init__(self, Y, input_dim, init='pca', X=None, kernel=None, normalize_Y=False):
+    def __init__(self, Y, input_dim, init='PCA', X=None, kernel=None, normalize_Y=False):
         if X is None:
             X = initialise_latent(init, input_dim, Y)
         if kernel is None:
@@ -60,7 +60,7 @@ class GPLVM(GP):
         for i in range(self.output_dim):
             target[:,:,i] = self.kern.dK_dX(np.dot(self.Ki,self.likelihood.Y[:,i])[None, :],X,self.X)
         return target
-   
+
     def magnification(self,X):
         target=np.zeros(X.shape[0])
         J = np.zeros((X.shape[0],X.shape[1],self.output_dim))
diff --git a/GPy/models_modules/mrd.py b/GPy/models_modules/mrd.py
index 2376993d..862b697b 100644
--- a/GPy/models_modules/mrd.py
+++ b/GPy/models_modules/mrd.py
@@ -42,7 +42,7 @@ class MRD(Model):
 
     """
     def __init__(self, likelihood_or_Y_list, input_dim, num_inducing=10, names=None,
-                 kernels=None, initx='pca',
+                 kernels=None, initx='PCA',
                  initz='permute', _debug=False, **kw):
         if names is None:
             self.names = ["{}".format(i) for i in range(len(likelihood_or_Y_list))]
@@ -237,7 +237,7 @@ class MRD(Model):
                                                 partial=g.partial_for_likelihood)]) \
                               for g in self.bgplvms])))
 
-    def _init_X(self, init='pca', likelihood_list=None):
+    def _init_X(self, init='PCA', likelihood_list=None):
         if likelihood_list is None:
             likelihood_list = self.likelihood_list
         Ylist = []
diff --git a/GPy/models_modules/sparse_gplvm.py b/GPy/models_modules/sparse_gplvm.py
index 44f6c7ef..04d3415a 100644
--- a/GPy/models_modules/sparse_gplvm.py
+++ b/GPy/models_modules/sparse_gplvm.py
@@ -7,9 +7,6 @@ import pylab as pb
 import sys, pdb
 from sparse_gp_regression import SparseGPRegression
 from gplvm import GPLVM, initialise_latent
-# from .. import kern
-# from ..core import model
-# from ..util.linalg import pdinv, PCA
 
 class SparseGPLVM(SparseGPRegression, GPLVM):
     """