def run(self, xs, params = None, *args, **kwargs): """ Run with default args """ if params == None: # Randomly initialize with means near uniform. M = dirichlet(ones(self.d) * 1., self.k).T w = ones(self.k)/self.k params = M, w return EMAlgorithm.run(self, xs, params, *args, **kwargs)
def __init__(self, k, d): self.k, self.d = k, d EMAlgorithm.__init__(self)