'Predict label for data.
Parameters
X : array-like, shape = [n_samples, n_features]
Returns
C : array, shape = (n_samples,) component memberships'

def predict(self, X):

    (logprob, responsibilities) = self.score_samples(X)
    return responsibilities.argmax(axis=1)
