sample_motif_hhmm.m

来自「贝叶斯网络matlab源程序,可用于分类,欢迎大家下载测试」· M 代码 · 共 11 行

%bnet = mk_motif_hhmm('motif_pattern', 'acca', 'background', 't');
bnet = mk_motif_hhmm('motif_pattern', 'accaggggga', 'background', []);

chars = ['a', 'c', 'g', 't'];
Tmax = 100;

for seqi=1:5
  evidence = cell2num(sample_dbn(bnet, 'length', Tmax));
  chars(evidence(end,:))
end