def test0(): random.set_state(state) (x, k, w) = initialize() (mu, S, priors, gamma, err) = kcluster0.gmm(x, k, weights=w, nreplicates=10)
def test0(): random.set_state(state) (x,k,w) = initialize() (mu,S,priors,gamma,err) = kcluster0.gmm(x,k,weights=w,nreplicates=10)
(mu, S, priors, gamma, err) = kcluster0.gmm(x, k, weights=w, nreplicates=10) #print "mu = " + str(mu) #print "idx = " + str(idx) if __name__ == "__main__": for trial in range(10): print "trial = " + str(trial) (x, k, w) = initialize() state = random.get_state() if True: (mu0, S0, priors0, gamma0, err0) = kcluster0.gmm(x, k, weights=w, nreplicates=10) random.set_state(state) (mu, S, priors, gamma, err) = kcluster.gmm(x, k, weights=w, nreplicates=10) D = (num.tile(mu[:,0].reshape((k,1)),(1,k)) - mu0[:,0])**2 + \ (num.tile(mu[:,1].reshape((k,1)),(1,k)) - mu0[:,1])**2 (a1, a2) = hungarian.hungarian(D) if True or \ num.max(num.abs(mu - mu0[a1,:])) > .001 or \ num.max(num.abs(S-S0[:,:,a1])) > .001 or \ num.max(num.abs(priors - priors0[a1])) > .001 or \ num.max(num.abs(gamma-gamma0[:,a1])) > .001 or \ num.abs(err - err0) > .001:
# D = kcluster.clusterdistfun1(x,c) def test0(): random.set_state(state) (x,k,w) = initialize() (mu,S,priors,gamma,err) = kcluster0.gmm(x,k,weights=w,nreplicates=10) #print "mu = " + str(mu) #print "idx = " + str(idx) if __name__ == "__main__": for trial in range(10): print "trial = " + str(trial) (x,k,w) = initialize() state = random.get_state() if True: (mu0,S0,priors0,gamma0,err0) = kcluster0.gmm(x,k,weights=w,nreplicates=10) random.set_state(state) (mu,S,priors,gamma,err) = kcluster.gmm(x,k,weights=w,nreplicates=10) D = (num.tile(mu[:,0].reshape((k,1)),(1,k)) - mu0[:,0])**2 + \ (num.tile(mu[:,1].reshape((k,1)),(1,k)) - mu0[:,1])**2 (a1,a2) = hungarian.hungarian(D) if True or \ num.max(num.abs(mu - mu0[a1,:])) > .001 or \ num.max(num.abs(S-S0[:,:,a1])) > .001 or \ num.max(num.abs(priors - priors0[a1])) > .001 or \ num.max(num.abs(gamma-gamma0[:,a1])) > .001 or \ num.abs(err - err0) > .001: print "cluster approx cluster0["+str(a1)+"]" print "max(|mu - mu0|) = " + str(num.max(num.abs(mu - mu0[a1,:]))) print "max(|S - S0|) = " + str(num.max(num.abs(S-S0[:,:,a1]))) print "max(|priors - priors0|) = " + str(num.max(num.abs(priors - priors0[a1])))