pysal.model.spreg.ML_Lag_Regimes

class pysal.model.spreg.ML_Lag_Regimes(y, x, regimes, w=None, constant_regi='many', cols2regi='all', method='full', epsilon=1e-07, regime_lag_sep=False, regime_err_sep=False, cores=False, spat_diag=False, vm=False, name_y=None, name_x=None, name_w=None, name_ds=None, name_regimes=None)[源代码]

带状态的空间滞后模型的ML估计(注:未添加一致性检查、诊断或常量);Anselin(1988年)【Anselin1988年】_

参数:
y : 数组

因变量nx1数组

x : 数组

二维数组,n行,每个独立(外生)变量一列,不包括常量

regimes : 列表

n值列表,并将每个观测值映射到一个区域。假定与“x”对齐。

常量_regi:[一','多']

切换器控制恒量项设置。它可以采用以下值:

  • “one”:1的向量附加到x并保持
    各政权间的常数
  • “many”:将一个向量附加到x并考虑
    每个制度不同(默认)
cols2regi : 列出“全部”

表明x的每一列是否应视为不同的制度或保持不变的论点(错误)。如果一个列表,k布尔值表示每个变量的选项(如果每个区域一个,则为真,如果保持不变,则为假)。如果“全部”(默认),则所有变量都会随状态而变化。

w : 稀疏矩阵

空间权重稀疏矩阵

method : 一串

如果“满”,强力计算(全矩阵表达式)如果“ORD”,ORD特征值方法如果“LU”,LU稀疏矩阵分解

ε : 浮动

拟素标量函数与逆积的公差准则

状态_lag_sep:boolean

如果为真,空间滞后的空间参数也会根据不同的状态进行计算。如果为“假”(默认),则空间参数为固定交叉区域。

cores : 布尔

指定是否使用多处理默认值:无多处理,cores=false注意:多处理可能无法在所有平台上工作。

spat_diag : 布尔

如果为真,则包括空间诊断(尚未实现)

vm : 布尔

如果为真,则在汇总结果中包含方差-协方差矩阵

name_y : 一串

输出中使用的从属变量的名称

name_x : 字符串表

输出中使用的自变量的名称

name_w : 一串

输出中使用的权重矩阵的名称

name_ds : 一串

用于输出的数据集名称

name_regimes : 一串

输出中使用的状态变量的名称

属性:
summary : 一串

回归结果和诊断摘要(注:与print命令一起使用)

betas : 数组

(k+1)x1估计系数数组(rho-first)

rho : 浮动

空间自回归系数的估计仅在多重回归时字典“multi”中可用(详见下文“multi”)。

u : 数组

nx1残差数组

predy : 数组

nx1预测Y值数组

n : 整数

观测次数

k : 整数

仅当进行多次回归时,字典“multi”中才提供估算系数的变量数(包括常数,不包括rho)(详情请参见下面的“multi”)。

y : 数组

因变量nx1数组

x : 数组

二维数组,有n行,每个独立(外生)变量有一列,其中常量仅在多重回归时字典“multi”中可用(有关详细信息,请参见下面的“multi”)。

method : 一串

对数雅可比方法if‘full’:蛮力(全矩阵计算)if‘ord’,ord特征值方法if‘lu’,lu稀疏矩阵分解

ε : 浮动

极小标量函数与逆积的公差准则

mean_y : 浮动

因变量均值

std_y : 浮动

因变量标准差

vm : 数组

方差协方差矩阵(k+1 x k+1),所有系数

vm1 : 数组

方差协方差矩阵(k+2 x k+2),包括sig2,仅在多重回归时字典“multi”中可用(详情请参见下面的“multi”)。

sig2 : 浮动

sigma平方在计算中使用,只有在多重回归时字典“multi”中才可用(有关详细信息,请参阅下面的“multi”)。

logll : 浮动

最大化的日志可能性(包括常量项)仅在多重回归时字典“multi”中可用(有关详细信息,请参阅下面的“multi”)。

aic : 浮动

Akaike信息标准仅在多重回归时在字典“multi”中可用(有关详细信息,请参阅下面的“multi”)。

schwarz : 浮动

Schwarz标准仅在多重回归时在字典“multi”中可用(有关详细信息,请参阅下面的“multi”)。

predy_e : 数组

还原形式的预测值

e_pred : 数组

用简化预测值预测误差

pr2 : 浮动

伪r平方(y和ypred之间的平方相关性)仅在多重回归时字典“multi”中可用(有关详细信息,请参阅下面的“multi”)。

pr2_e : 浮动

伪r平方(y和ypred_e之间的平方相关性(使用约简形式))仅在多重回归时字典“multi”中可用(有关详细信息,请参阅下面的“multi”)。

std_err : 数组

1XK的betas标准错误数组,仅在多重回归时字典'multi'中可用(有关详细信息,请参阅下面的'multi')

z_stat : 元组列表

z statistic;每个元组都包含一对(statistic,p-value),其中每个元组都是一个浮点,只有在进行多次回归时才能在字典“multi”中使用(有关详细信息,请参见下面的“multi”)。

name_y : 一串

输出中使用的从属变量的名称

name_x : 字符串表

输出中使用的自变量的名称

name_w : 一串

输出中使用的权重矩阵的名称

name_ds : 一串

用于输出的数据集名称

name_regimes : 一串

输出中使用的状态变量的名称

title : 一串

当进行多次回归时,仅在字典“multi”中可用的回归方法的名称(有关详细信息,请参阅下面的“multi”)。

regimes : 列表

n值列表,并将每个观测值映射到一个区域。假定与“x”对齐。

常量_regi:[一','多']

如果状态为假,则忽略。制度不变的选择。切换器控制恒量项设置。它可以采用以下值:

  • “one”:1的向量附加到x并保持
    各政权间的常数
  • “many”:将一个向量附加到x并考虑
    不同体制
cols2regi : 列出“全部”

如果状态为假,则忽略。表明x的每一列是否应视为不同的制度或保持不变的论点(错误)。如果一个列表,k布尔值表示每个变量的选项(如果每个区域一个,则为真,如果保持不变,则为假)。如果“全部”,则所有变量都随状态而变化。

regime_lag_sep : 布尔

如果为真,空间滞后的空间参数也会根据不同的状态进行计算。如果为“假”(默认),则空间参数为固定交叉区域。

regime_err_sep : 布尔

始终设置为“假”,以便与其他体制模型兼容

kr : 利息

变量/列的数量将被“区域化”或按制度变化。这将导致对每个变量(即每个变量的nr参数)按状态进行一次参数估计。

kf : 利息

变量/列的数量被认为是固定的或全局的,因此只能得到一个参数估计。

nr : 利息

“制度”列表中不同制度的数量

multi : 词典

只有在估计多个回归时才可用,即当状态“err_sep=真”且没有变量跨状态固定时。包含每个回归的所有属性

实例
________
使用pysal打开data baltim.dbf并创建变量矩阵和权重矩阵。
>>>将numpy导入为np
>>>导入pysal.lib
>>>来自pysal.lib导入示例
>>>db=pysal.lib.io.open(examples.get ou path(“baltim.dbf”),'r”)。
>>>ds_name=“baltim.dbf”,
>>>y_name=“price”,
>>>Y=np.array(db.by_col(y_name)).t
>>>Y.shape=(len(y),1)
>>>X_name=[“nroom”,“age”,“sqft”]
>>>x=np.array([db.by_col(var)for var in x_name]).t
>>>ww=ps.open(ps.examples.get_path(“baltim_q.gal”))
>>>w=w w.read()。
>>>ww.close()。
>>>W_name=“Baltim_q.gal”,
>>>w.transform='r',
因为在这个例子中,我们有兴趣检查结果是否变化
根据制度,我们使用citcou来定义地点是在城市还是
城外(县内):。
>>>状态=db.by_col(“citcou”)。
现在我们可以使用所有参数运行回归:
>>>ml lag=ml_lag_状态(y,x,状态,w=w,name_y=y_名称,name_x=x_名称,name_w=w_名称,name_ds=ds_名称,name_状态=“citcou”)。
>>>np.around(mllag.betas,小数=4)
数组([-15.0059],

[4.496]、[-0.0318]、[0.35]、[-4.5404]、[3.9219]、[-0.1702]、[0.8194]、[0.5385])

>>>“0:.6f”.格式(mllag.rho)
'0.538503'
>>>“0:.6f”.格式(mllag.mean_y)
'44.307180'
>>>“0:.6f”.格式(mllag.std_y)
'23.606077'
>>>np.around(np.diag(mllag.vm1),小数=4)
数组([47.42,2.3953,0.0051,0.0648,69.6765,3.2066,]

0.0116、0.0486、0.004、390.7274])

>>>np.around(np.diag(mllag.vm),小数=4)
数组([47.42,2.3953,0.0051,0.0648,69.6765,3.2066,]

0.0116,0.0486,0.004])

>>>“0:.6f”.格式(mllag.sig2)
'200.044334'
>>>“0:.6f”.格式(mllag.logll)
“-864.985056”
>>>“0:.6f”.格式(mllag.aic)
'1747.970112'
>>>“0:.6f”.格式(mllag.schwarz)
'1778.136835'
>>>MLLAG.标题
“最大似然空间滞后-状态(方法=完全)”

方法

ML_Lag_Regimes_Multi  
__init__(y, x, regimes, w=None, constant_regi='many', cols2regi='all', method='full', epsilon=1e-07, regime_lag_sep=False, regime_err_sep=False, cores=False, spat_diag=False, vm=False, name_y=None, name_x=None, name_w=None, name_ds=None, name_regimes=None)[源代码]

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

方法

ML_Lag_Regimes_Multi \(Y,X,Wu i,W,regiu id,…)
__init__ \(Y,X,状态[,W,常量区域,…]) 初始化自身。

属性

mean_y 
sig2n 
sig2n_k 
std_y 
utu 
vm