pysal.explore.giddy.ergodic.fmpt

pysal.explore.giddy.ergodic.fmpt(P)[源代码]

计算遍历转移概率矩阵的第一平均通过时间矩阵。

参数:
P : 数组

(k,k),遍历马尔可夫转移概率矩阵。

返回:
M : 数组

(k,k),元素是从状态i开始的链首先进入状态j所需的间隔数的预期值。如果i=j,则这是重复时间。

笔记

[KS67] .

实例

>>> import numpy as np
>>> from pysal.explore.giddy.ergodic import fmpt
>>> p=np.array([[.5, .25, .25],[.5,0,.5],[.25,.25,.5]])
>>> fm=fmpt(p)
>>> fm
array([[2.5       , 4.        , 3.33333333],
       [2.66666667, 5.        , 2.66666667],
       [3.33333333, 4.        , 2.5       ]])

因此,如果今天在oz下着雨,我们可以预计在4天内会有一个好天气,平均来说,3.33天会下雪。我们预计2.5天后还会下雨。如果今天在奥兹天气很好,我们将在2.67天内经历天气变化(雨或雪)。(那个邪恶的女巫只能死一次,所以我认为那是终极的吸收状态)。