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])))
