python - 用 PySal 计算 Moran's I？

Question

对于 Uni 的一项作业，我们必须使用 PySal 来计算 Moran's I。我找到了一个示例代码并从那里开始，但我收到错误消息“‘文件’对象没有属性‘by_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.open 尝试根据扩展名和检查文件来确定文件类型。 pysal.open 返回的类型并不总是很清楚，应该使用内置类型函数进行检查。在您的情况下，它返回一个普通的 python 文件对象，这意味着 pysal 无法为您解析它。 “by_col”方法仅在 PySAL 将您的文件识别为 DataTable 时可用。

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

例如示例 STL_HOM.txt 数据集的前几行。其中有 78 个数据行和 4 个列。

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

您可以将您的文件重新格式化为一种受支持的 DataTable 类型，或者编写您自己的代码来解析您的文本文件。