用PySal计算莫兰的I值？

Question

对于Uni的赋值，我们必须使用PySal来计算Moran的I。我找到了一个示例代码，并从中得到错误消息“'file'object has no attribute'by'u col'”

任何想法都将不胜感激。数据集是纽约的一个文本文件，“T0帴MINC”属性是包含不同地区家庭收入中值的列。脚本附在下面

import os
import pysal
import numpy as np

os.chdir = ("C:\\Users\\hsv0336\\Desktop")
wd = os.getcwd()

w = pysal.open(wd + "\\Rook.gal")
f = pysal.open(wd + "\\NewYorkData.txt")
y = np.array(f.by_col["T0_MINC"])

mi = pysal.Moran(y, w, two_tailed=False)

我使用的示例脚本位于以下链接的“Moran”下： http://pysal.readthedocs.io/en/latest/users/tutorials/autocorrelation.html#moran-s-i

Answer 1

在pysal.打开尝试根据扩展名并通过检查文件来确定文件类型。返回的类型pysal.打开不总是清晰的，应使用内置类型函数进行检查。在您的例子中，它返回一个纯python文件对象，这意味着pysal无法为您解析它。只有当PySAL将文件识别为DataTable时，“by_col”方法才可用。在

支持的数据表类型包括带有特殊“GeoDa”头的.csv、.dbf或.txt文件。GeoDa头包含第一行中的数据行数和列数，第二行中包含列的名称。在

例如STL示例的前几行_宏.txt数据集。它有78个数据行和4个列。在

78,4
"FIPSNO","HR8488","HR8893","HC8488"
17107,1.290722,1.624458,2

您可以将文件重新格式化为受支持的数据表类型之一，也可以编写自己的代码来解析文本文件。在