HMM Viterbi算法_2 SAS代码_Andy的原创帖13

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

重新整理了一下Viterbi的程序，写成了一个宏，可以展示序列和每个时段的概率。（有点赶时间，就不多废话了)维特比的思路：最大可能路径


例子的参数:



观察序列:

示例结果：

我的运行结果：

我的宏(viterbi)

%macroviterbi(obs);

prociml;

/*

M=2 观察状态数为2

N=3 隐含状态数为3

*/

/*用函数创建3*3方阵-这个例子的转移矩阵不是很好*/

/*A为转移矩阵*/

A=j(3,3,0.333);

/*B为发射矩阵*/

B={0.5 0.5,

   0.75 0.25,

   0.25 0.75};

/*用函数创建初始向量*/

pi=j(1,3,0.333);

/*创建单位矩阵，用于选择选举胜出的状态*/

vec=i(3);

/*初始状态*/

/*转移矩阵-不变*/

h=pi;

/*观察序列*/

os={&obs};

/*观察状态-确定性,第一次为os[1]*/

o=b[,os[1]]`;

/*发射-产生候选比较集*/

can1=h#o;

/*获取获胜隐含状态及概率*/

ws=can1[<:>];

pr=can1`;

do i=2 to ncol(os);

h_win=vec[ws[i-1],];

h=h_win*a;

/*观察状态*/

o=b[,os]`;

/*发射-产生候选比较集*/

can1=h#o;

/*获取获胜隐含状态及概率*/

ws=ws||can1[<:>];

pr=pr||can1`;

end;

printws;

printpr;

quit;

%mend;

%viterbi(1 1 1 1 2 1 2 2 2 2)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：sas代码 ITER Andy HMM vit 程序 维特

o=b[,os]`- >o=b[,os[i]]`, 数组下标被吃掉了...