用PySal计算莫兰的I？

Question

对于Uni的一个赋值，我们必须使用文件来计算莫兰的I。我已经找到了一个示例代码，并从那里开始操作，但我得到了错误消息“‘PySal’object has no attribute '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不能为你解析它。只有当PySAL将您的文件识别为DataTable时，"by_col“方法才可用。

支持的DataTable类型包括带有特殊"GeoDa“头的.csv、.dbf或.txt文件。GeoDa标题包含第一行中的数据行数和列数，以及第二行中的列名。

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

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

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