pysal.explore.giddy.markov.kullback¶

pysal.explore.giddy.markov.kullback(F)[源代码]¶

马尔可夫同质性的基于Kullback信息的检验。

参数:	F : 数组 （s，r，r），值是s地层，r初始状态，r终端状态的转换（不是概率）。
返回:	Results : 词典 （键值） 条件同质性（浮动）测试跨地层过渡概率同质性的统计。 测试统计的条件同质性pValue-（float）p值。 条件同质性dof-（int）自由度=r（s-1）（r-1）。

笔记

基于 [KKK62] . 下例摘自表9.2。

实例

>>> import numpy as np
>>> from pysal.explore.giddy.markov import kullback
>>> s1 = np.array([
...         [ 22, 11, 24,  2,  2,  7],
...         [ 5, 23, 15,  3, 42,  6],
...         [ 4, 21, 190, 25, 20, 34],
...         [0, 2, 14, 56, 14, 28],
...         [32, 15, 20, 10, 56, 14],
...         [5, 22, 31, 18, 13, 134]
...     ])
>>> s2 = np.array([
...     [3, 6, 9, 3, 0, 8],
...     [1, 9, 3, 12, 27, 5],
...     [2, 9, 208, 32, 5, 18],
...     [0, 14, 32, 108, 40, 40],
...     [22, 14, 9, 26, 224, 14],
...     [1, 5, 13, 53, 13, 116]
...     ])
>>>
>>> F = np.array([s1, s2])
>>> res = kullback(F)
>>> "%8.3f"%res['Conditional homogeneity']
' 160.961'
>>> "%d"%res['Conditional homogeneity dof']
'30'
>>> "%3.1f"%res['Conditional homogeneity pvalue']
'0.0'