# Imports
#
import matplotlib.pyplot as plt
from libpysal.weights.contiguity import Queen
from libpysal import examples
import geopandas as gpd
from esda.moran import Moran
from splot.esda import plot_moran_simulation
#
# Load data and calculate weights
#
guerry = examples.load_example('Guerry')
link_to_data = guerry.get_path('guerry.shp')
gdf = gpd.read_file(link_to_data)
y = gdf['Donatns'].values
w = Queen.from_dataframe(gdf)
w.transform = 'r'
#
# Calculate Global Moran
#
moran = Moran(y, w)
#
# plot
#
plot_moran_simulation(moran)
plt.show()
#
# customize plot
#
plot_moran_simulation(moran, fitline_kwds=dict(color='#4393c3'))
plt.show()