import pysal
import numpy as np


f = pysal.open("soja.dbf")
y = np.array(f.by_col["DryYield ("])
w = pysal.knnW_from_shapefile("soja.shp", k=6)
lm = pysal.Moran_Local(y,w)
print lm