E = S[2]
if perturbation:
for obs in listObservables:
for state in listState:
if obs in E and not state=="":
E[obs][state]+=(random.random()-0.5)*coef
###########################################################
hmm = HMM_BW(listObservables, listState, 1)
if hmm_determine:
hmm.Pi = I
hmm.Pi[""] = 0.0
T[""] = hmm.Pi
for state in hmm.Pi:
T[state][""] = 0.01
hmm.transitions = T
hmm.transitions[""] = {}
E[""] = hmm.Pi
for obs in hmm.listObservables:
E[obs][""] = 0.0
E[""][""] = 1.0
hmm.emissions = E
hmm.emissions[""] = {}
for key in hmm.Pi:
hmm.transitions[""][key] = hmm.Pi[key]
hmm.emissions[""][key] = hmm.Pi[key]
for state in hmm.emissions[""]:
if state == "":
hmm.emissions[""][state] = 1.0
