pysal.explore.giddy.rank.Theta

class pysal.explore.giddy.rank.Theta(y, regime, permutations=999)[源代码]

制度流动性措施。 [Rey04]

对于时间序列,theta测量在n个位置上测量的变量在n个位置的互斥和穷尽分区(区域)内等级变化的程度。

θ的定义是,在所有绝对等级变化的总和之上,政权内等级变化的绝对和的总和。

参数:
y : 数组

(n,k)当k大于等于2时,y的连续列是时间上较迟的时刻(年、月等)。

regime : 数组

(n,)值,与每个观测所属的状态相对应。

permutations : 利息

为基于计算的推理生成的随机空间排列数。

实例

>>> import pysal.lib as ps
>>> from pysal.explore.giddy.rank import Theta
>>> import numpy as np
>>> f=ps.io.open(ps.examples.get_path("mexico.csv"))
>>> vnames=["pcgdp%d"%dec for dec in range(1940,2010,10)]
>>> y=np.transpose(np.array([f.by_col[v] for v in vnames]))
>>> regime=np.array(f.by_col['esquivel99'])
>>> np.random.seed(10)
>>> t=Theta(y,regime,999)
>>> t.theta
array([[0.41538462, 0.28070175, 0.61363636, 0.62222222, 0.33333333,
        0.47222222]])
>>> t.pvalue_left
array([0.307, 0.077, 0.823, 0.552, 0.045, 0.735])
>>> t.total
array([130., 114.,  88.,  90.,  90.,  72.])
>>> t.max_total
512
属性:
ranks : 数组

原始Y数组的列组(按列)。

regimes : 数组

原始状态数组。

total : 数组

(k-1,)k个期间中每个期间的排名变化总数。

max_total : 利息

n个观测的理论最大等级变化数。

theta : 数组

(k-1,)k-1区间的θ统计。

permutations : 利息

排列的数目。

pvalue_left : 浮动

在完全空间随机性条件下,观测到的θ值明显低于预期值。

pvalue_right : 浮动

在完全空间随机性下,观测到的θ的检验P值明显大于其预期值。

__init__(y, regime, permutations=999)[源代码]

初始化自身。请参阅帮助(键入(self))以获得准确的签名。

方法

__init__ \(Y,状态[,排列]) 初始化自身。