HMM前向-后向算法中的欠流



我正在学习HMM,并希望在python中实现前向-后向算法(Baum-Welch算法(。我使用来自NCBI人类基因组资源的人类基因组作为我的输入数据。但它太大了,我在正向算法中的alpha和反向算法中的beta将下溢(数字太小,无法保存(。我发现我也有同样的问题,但它的资源链接断了,我无法理解答案。我也找到了这个pdf,但我不能完全理解如何计算C_n和alpha^(zn(,以及如何在python中实现它。有人能帮我解决这个问题吗?非常感谢。

试着看看这个,在这里很难写出C_t的方程。您还应该考虑使用帮助您计算Baum-Welch算法的库,如石榴或torchegranate

最新更新