matlab在HMM预测疾病流行状态时观测值很大的处理

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

问题同之前我发的几个帖子，主要还是根据发病数预测疾病的流行状态，发病数这个观测在流行与未流行状态下产生的概率服从poisson分布，matlab程序如下：
j=0:10                                                                  # j 为观测值，即发病数
A1=(5.^j) *exp(-1.7)
B1=factorial(j)
p1=A.1/B1                                                          #p1：未流行状态下发病数的概率                                                      
j=0:10
A2=(15.^j) *exp(-4.8)
B2=factorial(j)
p2=A2./B2                                                           #p2：流行状态下发病数的概率   
tr = [0.8874,0.1126;
0.4677,0.5323];                                                  #tr：两个状态的转移概率
e = [p;p2;];                                                            #e：两个状态下发病数的概率
seq = hmmgenerate(100,tr,e);                       
[estTR, estE] = hmmtrain(seq,tr,e);
estimatedStates = hmmviterbi(seq,tr,e);
pStates =  hmmdecode(seq, tr, e)                  #分别产生观测序列、参数估计值和状态序列

现在的问题是，j 即发病数很大，在700-3000左右，这样j的阶乘B就很大，在matlab中就不能显示，为“INF”，P 就为0了，之后的参数估计值和状态序列这些就没有意义了。但实际上在几百万人中发生千多的病例，应该也算是poisson分布，我现在就不知道应该怎么处理这种观测值 j 很大matlab给不出结果的情形呢？是数据可以在怎么转化或者分布出问题了，还是程序有其他写法呢？
本人是学生，非常希望再次得到大家的帮助，谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB atlab matla Lab Mat matlab 流行 程序

我开始本来想如果就这个程序，把观测值按照发病数取值范围划分为几个区间，用分类变量（1，2，3，4.。。）代替发病数，阶乘问题倒是解决了，但是这个新的观测区间产生的概率又不知道了，还是不行一样
困惑

Possible solutions either one or both or more.
1. large state number, not just binary (agreed with 2nd 苏攀攀).
2. observed sequence be ohter variables such as distributions, instead of frequency.
3. processing your data before hmm

请问richardma ，关于你第二点建议，我没其他的指示变量了，频数转化为分布？
第三点建议是可以咋个处理data呢？
我笨

You are smart. Would you mind talking to your professor or teaching assistants about it? Cannot tell a long story short in a blog.
By the way, if this is a course assignment please take the advice: not to overdo, just make use of what you learned.
Sorry for disappointing you.

thank you，but I really want to try it although there does have much difficulties,especially to a medical student .Fortunately,I have got much help from many teachers and friends in the web,including you.For good or ill ,I have learn much.

Are you trying to estimate prevalence from incidence?

顶起来，好好学习

是的哈，大致就是通过发病数来预测疾病流行状况