python - 求马尔可夫链状态转移概率矩阵

Question

一个时间序列如下： {代码...} 这个时间序列一共有4032个值，这些值的取值空间有115个取法a=[-0.49，-0.47,-0.45······1.75,1.77,1.79]，现在想求出来这个115*115的状态转移概率矩阵，写了一个matlab程序但是计算...

漂亮男人 · Answer

不太清楚马尔可夫链如何计算，猜测是状态转移，看看下面的有没有用。

a = [1,2,3,0,0,0,3,3,2,3,3,2,3,2,1,2,3]


l = 4
N = [[0] * l for i in range(l)]

amount = 0
for i in data:
    print i
    amount += 1
    N[i[0]][i[1]] += 1

print amount
for i in range(l):
    for j in range(l):
        N[i][j] /= float(amount)
print N