Skip to content

Instantly share code, notes, and snippets.

@ljwolf

ljwolf/induce.py

Last active April 15, 2016 23:19
Show Gist options
  • Save ljwolf/e131fb43e059ed6c54a635c6ac8c5bdc to your computer and use it in GitHub Desktop.
Save ljwolf/e131fb43e059ed6c54a635c6ac8c5bdc to your computer and use it in GitHub Desktop.
Download ZIP
inducing bug #587
Raw
induce.py
import pysal as ps
import numpy as np
df = ps.pdio.read_files(ps.examples.get_path('columbus.shp'))
W = ps.queen_from_shapefile(ps.examples.get_path('columbus.shp'))
y = df['HOVAL'].values
y.reshape(W.n,1)
X = df[['CRIME', 'INC']].values
y[5] = np.nan
X[5] = np.array([np.nan, np.nan])
ps.spreg.ML_Lag(y,X,w=W)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment