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)
