balance-of-london-explorer

import%20marimo%0A%0A__generated_with%20%3D%20%220.14.17%22%0Aapp%20%3D%20marimo.App()%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%23%20Data%20Explorer%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20pandas%20as%20pd%0A%20%20%20%20from%20matplotlib%20import%20pyplot%20as%20plt%0A%20%20%20%20import%20seaborn%20as%20sns%0A%20%20%20%20from%20scipy%20import%20stats%0A%20%20%20%20import%20altair%20as%20alt%0A%20%20%20%20import%20pyarrow%0A%20%20%20%20return%20alt%2C%20mo%2C%20np%2C%20pd%2C%20plt%2C%20sns%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23data_url%20%3D%20'https%3A%2F%2Fbalance-of-london.wandering-axis.xyz%2Fdata'%0A%20%20%20%20data_url%20%3D%20'..%2Fwebmap%2Fdata'%0A%20%20%20%20return%20(data_url%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%23%23%23%20Accessibility%20Percentiles%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(data_url%2C%20pd)%3A%0A%20%20%20%20node_results%20%3D%20pd.read_parquet(f'%7Bdata_url%7D%2Fnode_place_mean_h3.parquet')%0A%20%20%20%20node_results.head()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%23%23%23%20Spatial%20Join%20to%20Census%20Population%20Data%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(data_url%2C%20pd)%3A%0A%20%20%20%20poi_census_sj%20%3D%20pd.read_parquet(f'%7Bdata_url%7D%2Fpoi_census.parquet')%0A%20%20%20%20poi_census_sj.head(1)%0A%20%20%20%20return%20(poi_census_sj%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%23%23%23%20Population%20Accessibility%20Distribution%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(plt%2C%20poi_census_sj%2C%20sns)%3A%0A%20%20%20%20clean_data%20%3D%20poi_census_sj.dropna(subset%3D%5B'poi_places_count_15_min_walk'%2C%20'totpop'%2C%20'poi_places_count_45_min_mm'%5D)%0A%20%20%20%20plt.figure(figsize%3D(15%2C%206))%0A%20%20%20%20%23%20sns.histplot(data%3Dclean_data%2C%20x%3D'poi_places_count_15_min_mm'%2C%20weights%3D'totpop'%2C%20stat%3D'count'%2C%20bins%3D50%2C%20alpha%3D0.6%2C%20color%3D'red'%2C%20kde%3DTrue%2C%20line_kws%3D%7B'color'%3A%20'blue'%2C%20'linewidth'%3A%202%7D)%0A%20%20%20%20sns.histplot(data%3Dclean_data%2C%20x%3D'poi_places_count_45_min_mm'%2C%20weights%3D'totpop'%2C%20stat%3D'count'%2C%20bins%3D50%2C%20alpha%3D0.6%2C%20color%3D'blue'%2C%20kde%3DTrue%2C%20line_kws%3D%7B'color'%3A%20'red'%2C%20'linewidth'%3A%202%7D)%0A%20%20%20%20plt.grid(True%2C%20alpha%3D0.3)%0A%20%20%20%20plt.xlabel('Places%20in%2045%20mins%20PT'%2C%20fontsize%3D12)%0A%20%20%20%20plt.ylabel('Residents'%2C%20fontsize%3D12)%0A%20%20%20%20plt.title('Accessibility%20Distribution'%2C%20fontsize%3D14)%0A%20%20%20%20plt.gca().xaxis.set_major_formatter(plt.FuncFormatter(lambda%20x%2C%20p%3A%20f'%7Bint(x%20%2F%201000)%7Dk'))%0A%20%20%20%20plt.gca().yaxis.set_major_formatter(plt.FuncFormatter(lambda%20y%2C%20p%3A%20f'%7Bint(y%20%2F%201000)%7Dk'))%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.show()%0A%20%20%20%20return%20(clean_data%2C)%0A%0A%0A%40app.cell%0Adef%20_(clean_data%2C%20np)%3A%0A%20%20%20%20mu%20%3D%20np.average(clean_data%5B'poi_places_count_45_min_mm'%5D%2C%20weights%3Dclean_data%5B'totpop'%5D)%0A%20%20%20%20sigma%20%3D%20np.sqrt(np.average((clean_data%5B'poi_places_count_45_min_mm'%5D%20-%20mu)%20**%202%2C%20weights%3Dclean_data%5B'totpop'%5D))%0A%20%20%20%20print(f'%5CnDistribution%20parameters%3A')%0A%20%20%20%20print(f'Mean%20(%CE%BC)%3A%20%7Bint(mu%20%2F%201000)%7Dk')%0A%20%20%20%20print(f'Standard%20deviation%20(%CF%83)%3A%20%7Bint(sigma%20%2F%201000)%7Dk')%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%23%23%23%20Node%20Place%20Balance%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(data_url%2C%20pd)%3A%0A%20%20%20%20points_gdf_h3_mean%20%3D%20pd.read_parquet(f'%7Bdata_url%7D%2Fnode_place_mean_h3.parquet')%0A%20%20%20%20return%20(points_gdf_h3_mean%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20points_gdf_h3_mean)%3A%0A%20%20%20%20dropdown_y%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dpoints_gdf_h3_mean.columns%2C%0A%20%20%20%20%20%20%20%20value%3D'poi_places_count_15_min_walk'%2C%0A%20%20%20%20%20%20%20%20label%3D%22Choose%20Local%20Measure%22)%0A%0A%20%20%20%20dropdown_x%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dpoints_gdf_h3_mean.columns%2C%0A%20%20%20%20%20%20%20%20value%3D'poi_places_count_45_min_mm'%2C%0A%20%20%20%20%20%20%20%20label%3D%22Choose%20City%20Wide%20Measure%22)%0A%0A%20%20%20%20dropdown_x%2Cdropdown_y%0A%20%20%20%20return%20dropdown_x%2C%20dropdown_y%0A%0A%0A%40app.cell%0Adef%20_(alt%2C%20dropdown_x%2C%20dropdown_y%2C%20mo%2C%20points_gdf_h3_mean)%3A%0A%20%20%20%20chart%20%3D%20mo.ui.altair_chart(%0A%20%20%20%20%20%20%20%20alt.Chart(%0A%20%20%20%20%20%20%20%20%20%20%20%20points_gdf_h3_mean%5B%5Bdropdown_x.selected_key%2Cdropdown_y.selected_key%5D%5D)%0A%20%20%20%20%20%20%20%20.mark_point()%0A%20%20%20%20%20%20%20%20.encode(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3Df'%7Bdropdown_x.selected_key%7D%3AQ'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20y%3Df'%7Bdropdown_y.selected_key%7D%3AQ'%0A%20%20%20%20%20%20%20%20)%20%2B%20%0A%20%20%20%20%20%20%20%20alt.Chart(points_gdf_h3_mean%5B%5Bdropdown_x.selected_key%2Cdropdown_y.selected_key%5D%5D)%0A%20%20%20%20%20%20%20%20.transform_regression(%0A%20%20%20%20%20%20%20%20%20%20%20%20dropdown_x.selected_key%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20dropdown_y.selected_key%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20method%3D'poly'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20order%3D3%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.mark_line(color%3D'red')%0A%20%20%20%20%20%20%20%20.encode(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3Df'%7Bdropdown_x.selected_key%7D%3AQ'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20y%3Df'%7Bdropdown_y.selected_key%7D%3AQ'%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%20%20%20%20%23%20chart.save('..%2Fwebmap%2Fchart.html')%0A%20%20%20%20chart%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A