[Ghmm-list] io problem?

[Felix N. <fel...@is...>]

hi,
can anybody help me with the following listing, seems like i'm unable to 
import back after writing

thanks,
Felix

#!/usr/bin/python

from ghmm import *

sigma = IntegerRange(1,7)
A = [[0.9, 0.1], [0.3, 0.7]]
efair = [1.0 / 6] * 6
print efair
eloaded = [3.0 / 13, 3.0 / 13, 2.0 / 13, 2.0 / 13, 2.0 / 13, 1.0 / 13]
B = [efair, eloaded]
pi = [0.5] * 2
m = HMMFromMatrices(sigma, DiscreteDistribution(sigma), A, B, pi)

m.write("model.txt")
hmm=HMMOpen("model.txt")
print hmm


-- 
Felix Naef, PhD
Assistant professor
Computational systems biology
Ecole Polytechnique Federale de Lausanne (EPFL)
Swiss Insitute of Bioinformatics (SIB)