import bnpy
import MoCap6
Data = MoCap6.get_data()
import os
import numpy as np
import glob
import PlotUtil as PU
import scipy.io
StateSeqUtil = bnpy.util.StateSeqUtil
def PrintContigBlocks(jobpath='', taskid='',
zBySeq=None, seqIDs=[1,2,3,4,5,6],
zSuffix='Aligned', minBlockLen=5):
if zBySeq is None:
taskpath = PU.MakePath(os.path.join(jobpath, str(taskid)))
print taskpath
zpaths = glob.glob(os.path.join(taskpath, '*MAPStateSeqs' + zSuffix + '.mat'))
zpaths.sort()
zpath = zpaths[-1]
zBySeq = scipy.io.loadmat(zpath)
zBySeq = zBySeq['zHatBySeq' + zSuffix]
zBySeq = bnpy.util.StateSeqUtil.convertStateSeq_MAT2list(zBySeq)
seqIDs = np.asarray(seqIDs) -1
for seqID in seqIDs:
print 'Seq #%d' % (seqID + 1)
z_n = zBySeq[seqID]
def Make_data_struct():
Data = MoCap6.get_data()
from IPython import embed; embed()
Data.save_to_mat('HMMdataset.mat')
