code: def predict(self, Xnew, full_cov=False, Y_metadata=None, kern=None, Yindex=0):
          
          b = self.bgplvms[Yindex]
          self.posterior = b.posterior
          self.kern = b.kern
          self.likelihood = b.likelihood
          return super(MRD, self).predict(Xnew, full_cov, Y_metadata, kern)
docs: Prediction for data set Yindex[default=0].
          This predicts the output mean and variance for the dataset given in Ylist[Yindex]