03 yellow green eda

%22%22%22%0ANYC%20Taxi%20Marimo%20-%20Notebook%2003%3A%20Yellow%20%26%20Green%20Taxi%20EDA%0A%0AExploratory%20Data%20Analysis%20of%20NYC%20Yellow%20and%20Green%20taxi%20data.%0AFocus%20on%20understanding%20trip%20patterns%20and%20congestion%20pricing%20impact.%0A%0ALearning%20Objectives%3A%0A%20%20%20%201.%20Perform%20temporal%20analysis%20(hourly%2C%20daily%2C%20weekly%20patterns)%0A%20%20%20%202.%20Analyze%20fare%20structures%20and%20pricing%20components%0A%20%20%20%203.%20Compare%20Yellow%20vs%20Green%20taxi%20behaviour%0A%20%20%20%204.%20Investigate%20congestion%20pricing%20impact%20(pre%2Fpost%20Jan%205%2C%202025)%0A%20%20%20%205.%20Build%20interactive%20visualisations%20with%20Plotly%0A%0APrerequisites%3A%0A%20%20%20%20-%20Completed%2002_data_loading.py%0A%20%20%20%20-%20Understanding%20of%20DuckDB%20queries%0A%0ACongestion%20Pricing%20Context%3A%0A%20%20%20%20NYC's%20congestion%20pricing%20started%20January%205%2C%202025.%20The%20cbd_congestion_fee%0A%20%20%20%20field%20captures%20per-trip%20charges%20for%20the%20MTA%20Congestion%20Relief%20Zone%0A%20%20%20%20(Manhattan%20below%2060th%20Street).%20This%20is%20one%20of%20the%20first%20large-scale%0A%20%20%20%20urban%20congestion%20pricing%20implementations%20in%20the%20US.%0A%0ATime%20estimate%3A%2090-120%20minutes%0A%22%22%22%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.19.5%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20duckdb%0A%20%20%20%20import%20polars%20as%20pl%0A%20%20%20%20import%20plotly.express%20as%20px%0A%20%20%20%20import%20plotly.graph_objects%20as%20go%0A%20%20%20%20from%20datetime%20import%20date%2C%20datetime%2C%20timedelta%0A%20%20%20%20return%20datetime%2C%20duckdb%2C%20mo%2C%20pl%2C%20px%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20from%20utils.data_urls%20import%20BASE_URL%2C%20TRIP_TYPES%2C%20build_trip_url%2C%20CONGESTION_PRICING_START%0A%20%20%20%20from%20utils.styling%20import%20(%0A%20%20%20%20%20%20%20%20TAXI_COLORS%2C%0A%20%20%20%20%20%20%20%20BOROUGH_COLORS%2C%0A%20%20%20%20%20%20%20%20CONGESTION_COLORS%2C%0A%20%20%20%20%20%20%20%20format_currency%2C%0A%20%20%20%20%20%20%20%20format_large_number%2C%0A%20%20%20%20)%0A%20%20%20%20from%20utils.zone_lookup%20import%20MANHATTAN_CBD_ZONES%2C%20AIRPORT_ZONES%2C%20is_cbd_zone%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20BASE_URL%2C%0A%20%20%20%20%20%20%20%20CONGESTION_COLORS%2C%0A%20%20%20%20%20%20%20%20MANHATTAN_CBD_ZONES%2C%0A%20%20%20%20%20%20%20%20TAXI_COLORS%2C%0A%20%20%20%20%20%20%20%20format_currency%2C%0A%20%20%20%20%20%20%20%20format_large_number%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell%0Adef%20_(duckdb)%3A%0A%20%20%20%20conn%20%3D%20duckdb.connect()%0A%20%20%20%20return%20(conn%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20analysis_options%20%3D%20mo.ui.dictionary(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22include_dec_2024%22%3A%20mo.ui.checkbox(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20label%3D%22Include%20Dec%202024%20(pre-congestion%20baseline)%22%2C%20value%3DTrue%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22sample_pct%22%3A%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start%3D1%2C%20stop%3D100%2C%20step%3D1%2C%20value%3D10%2C%20label%3D%22Sample%20%25%20(for%20large%20queries)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20)%0A%0A%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20analysis_options%5B%22include_dec_2024%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20analysis_options%5B%22sample_pct%22%5D%2C%0A%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20gap%3D2%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(analysis_options%2C)%0A%0A%0A%40app.cell%0Adef%20_(BASE_URL%2C%20analysis_options%2C%20conn%2C%20datetime%2C%20mo%2C%20pl)%3A%0A%20%20%20%20%23%20Build%20URL%20list%0A%20%20%20%20_urls%20%3D%20%5Bf%22%7BBASE_URL%7D%2Fyellow_tripdata_2025-01.parquet%22%5D%0A%20%20%20%20if%20analysis_options.value%5B%22include_dec_2024%22%5D%3A%0A%20%20%20%20%20%20%20%20_urls.append(f%22%7BBASE_URL%7D%2Fyellow_tripdata_2024-12.parquet%22)%0A%0A%20%20%20%20_sample_pct%20%3D%20analysis_options.value%5B%22sample_pct%22%5D%0A%20%20%20%20_sample_clause%20%3D%20f%22USING%20SAMPLE%20%7B_sample_pct%7D%20PERCENT%22%20if%20_sample_pct%20%3C%20100%20else%20%22%22%0A%20%20%20%20_query%20%3D%20f%22%22%22%0A%20%20%20%20%20%20%20%20SELECT%20%0A%20%20%20%20%20%20%20%20%20%20%20%20tpep_pickup_datetime%20as%20pickup_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tpep_dropoff_datetime%20as%20dropoff_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20PULocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20DOLocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20passenger_count%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20trip_distance%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tip_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tolls_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mta_tax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20improvement_surcharge%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20congestion_surcharge%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20cbd_congestion_fee%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20total_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20payment_type%0A%20%20%20%20%20%20%20%20FROM%20read_parquet(%7B_urls%7D%2C%20union_by_name%3Dtrue)%0A%20%20%20%20%20%20%20%20WHERE%20tpep_pickup_datetime%20%3E%3D%20'2024-12-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20tpep_pickup_datetime%20%3C%20'2025-02-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20fare_amount%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20AND%20trip_distance%20%3E%200%0A%20%20%20%20%20%20%20%20%7B_sample_clause%7D%0A%20%20%20%20%22%22%22%0A%20%20%20%20yellow_df%20%3D%20conn.execute(_query).pl()%0A%0A%20%20%20%20%23%20Add%20derived%20columns%0A%20%20%20%20yellow_df%20%3D%20yellow_df.with_columns(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.hour().alias(%22hour%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.weekday().alias(%22weekday%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.date().alias(%22trip_date%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22)%20%3E%3D%20pl.lit(datetime(2025%2C%201%2C%205))).alias(%22post_congestion%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.lit(%22Yellow%20Taxi%22).alias(%22taxi_type%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%0A%20%20%20%20mo.md(f%22**Yellow%20Taxi**%3A%20%7Blen(yellow_df)%3A%2C%7D%20trips%20loaded%20(%7B_sample_pct%7D%25%20sample)%22)%0A%20%20%20%20return%20(yellow_df%2C)%0A%0A%0A%40app.cell%0Adef%20_(BASE_URL%2C%20analysis_options%2C%20conn%2C%20datetime%2C%20mo%2C%20pl)%3A%0A%20%20%20%20_urls%20%3D%20%5Bf%22%7BBASE_URL%7D%2Fgreen_tripdata_2025-01.parquet%22%5D%0A%20%20%20%20if%20analysis_options.value%5B%22include_dec_2024%22%5D%3A%0A%20%20%20%20%20%20%20%20_urls.append(f%22%7BBASE_URL%7D%2Fgreen_tripdata_2024-12.parquet%22)%0A%20%20%20%20_sample_pct%20%3D%20analysis_options.value%5B%22sample_pct%22%5D%0A%20%20%20%20_sample_clause%20%3D%20f%22USING%20SAMPLE%20%7B_sample_pct%7D%20PERCENT%22%20if%20_sample_pct%20%3C%20100%20else%20%22%22%0A%20%20%20%20_query%20%3D%20f%22%22%22%0A%20%20%20%20%20%20%20%20SELECT%20%0A%20%20%20%20%20%20%20%20%20%20%20%20lpep_pickup_datetime%20as%20pickup_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20lpep_dropoff_datetime%20as%20dropoff_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20PULocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20DOLocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20passenger_count%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20trip_distance%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tip_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tolls_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mta_tax%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20improvement_surcharge%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20congestion_surcharge%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20cbd_congestion_fee%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20total_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20payment_type%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20trip_type%0A%20%20%20%20%20%20%20%20FROM%20read_parquet(%7B_urls%7D%2C%20union_by_name%3Dtrue)%0A%20%20%20%20%20%20%20%20WHERE%20lpep_pickup_datetime%20%3E%3D%20'2024-12-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20lpep_pickup_datetime%20%3C%20'2025-02-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20fare_amount%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20AND%20trip_distance%20%3E%200%0A%20%20%20%20%20%20%20%20%7B_sample_clause%7D%0A%20%20%20%20%22%22%22%0A%20%20%20%20green_df%20%3D%20conn.execute(_query).pl()%0A%20%20%20%20green_df%20%3D%20green_df.with_columns(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.hour().alias(%22hour%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.weekday().alias(%22weekday%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.date().alias(%22trip_date%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22)%20%3E%3D%20pl.lit(datetime(2025%2C%201%2C%205))).alias(%22post_congestion%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.lit(%22Green%20Taxi%22).alias(%22taxi_type%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20mo.md(f%22**Green%20Taxi**%3A%20%7Blen(green_df)%3A%2C%7D%20trips%20loaded%20(%7B_sample_pct%7D%25%20sample)%22)%0A%20%20%20%20return%20(green_df%2C)%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20hourly%20%3D%20yellow_df.group_by(%22hour%22).agg(pl.len().alias(%22trips%22)).sort(%22hour%22)%0A%0A%20%20%20%20fig_hourly%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20hourly.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22hour%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Yellow%20Taxi%20Trips%20by%20Hour%20of%20Day%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22hour%22%3A%20%22Hour%22%2C%20%22trips%22%3A%20%22Trip%20Count%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Add%20rush%20hour%20annotations%0A%20%20%20%20fig_hourly.add_vrect(x0%3D7%2C%20x1%3D9%2C%20fillcolor%3D%22red%22%2C%20opacity%3D0.1%2C%20annotation_text%3D%22AM%20Rush%22)%0A%20%20%20%20fig_hourly.add_vrect(x0%3D17%2C%20x1%3D19%2C%20fillcolor%3D%22red%22%2C%20opacity%3D0.1%2C%20annotation_text%3D%22PM%20Rush%22)%0A%0A%20%20%20%20fig_hourly.update_layout(xaxis%3Ddict(tickmode%3D%22linear%22%2C%20dtick%3D2))%0A%20%20%20%20fig_hourly%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20green_df%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20yellow_hourly%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.group_by(%22hour%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.with_columns(pl.lit(%22Yellow%20Taxi%22).alias(%22taxi_type%22))%0A%20%20%20%20)%0A%0A%20%20%20%20green_hourly%20%3D%20(%0A%20%20%20%20%20%20%20%20green_df.group_by(%22hour%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.with_columns(pl.lit(%22Green%20Taxi%22).alias(%22taxi_type%22))%0A%20%20%20%20)%0A%0A%20%20%20%20combined_hourly%20%3D%20pl.concat(%5Byellow_hourly%2C%20green_hourly%5D).sort(%22hour%22)%0A%0A%20%20%20%20fig_compare%20%3D%20px.line(%0A%20%20%20%20%20%20%20%20combined_hourly.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22hour%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20color%3D%22taxi_type%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Hourly%20Pattern%3A%20Yellow%20vs%20Green%20Taxi%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22hour%22%3A%20%22Hour%20of%20Day%22%2C%20%22trips%22%3A%20%22Trip%20Count%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_map%3DTAXI_COLORS%2C%0A%20%20%20%20%20%20%20%20markers%3DTrue%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_compare.update_layout(xaxis%3Ddict(tickmode%3D%22linear%22%2C%20dtick%3D2))%0A%20%20%20%20fig_compare%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20DOW_NAMES%20%3D%20%5B%22Mon%22%2C%20%22Tue%22%2C%20%22Wed%22%2C%20%22Thu%22%2C%20%22Fri%22%2C%20%22Sat%22%2C%20%22Sun%22%5D%0A%0A%20%20%20%20daily%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.group_by(%22weekday%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.sort(%22weekday%22)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22weekday%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20.map_elements(lambda%20x%3A%20DOW_NAMES%5Bx%20-%201%5D%2C%20return_dtype%3Dpl.Utf8)%20%20%23%20x%20-%201%20to%20convert%201-7%20to%200-6%0A%20%20%20%20%20%20%20%20%20%20%20%20.alias(%22day_name%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20fig_dow%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20daily.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22day_name%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Yellow%20Taxi%20Trips%20by%20Day%20of%20Week%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22day_name%22%3A%20%22Day%22%2C%20%22trips%22%3A%20%22Trip%20Count%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%20%20%20%20fig_dow%0A%20%20%20%20return%20(DOW_NAMES%2C)%0A%0A%0A%40app.cell%0Adef%20_(DOW_NAMES%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20heatmap_data%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.group_by(%5B%22weekday%22%2C%20%22hour%22%5D)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.sort(%5B%22weekday%22%2C%20%22hour%22%5D)%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Pivot%20for%20heatmap%0A%20%20%20%20pivot%20%3D%20heatmap_data.pivot(%0A%20%20%20%20%20%20%20%20values%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20index%3D%22weekday%22%2C%0A%20%20%20%20%20%20%20%20on%3D%22hour%22%2C%0A%20%20%20%20).fill_null(0)%0A%0A%20%20%20%20fig_heatmap%20%3D%20px.imshow(%0A%20%20%20%20%20%20%20%20pivot.select(pl.exclude(%22weekday%22)).to_numpy()%2C%0A%20%20%20%20%20%20%20%20x%3Dlist(range(24))%2C%0A%20%20%20%20%20%20%20%20y%3DDOW_NAMES%2C%0A%20%20%20%20%20%20%20%20labels%3Ddict(x%3D%22Hour%20of%20Day%22%2C%20y%3D%22Day%20of%20Week%22%2C%20color%3D%22Trips%22)%2C%0A%20%20%20%20%20%20%20%20title%3D%22Trip%20Volume%20Heatmap%3A%20Hour%20%C3%97%20Day%20of%20Week%22%2C%0A%20%20%20%20%20%20%20%20color_continuous_scale%3D%22YlOrRd%22%2C%0A%20%20%20%20%20%20%20%20aspect%3D%22auto%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_heatmap%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20hour_filter%20%3D%20mo.ui.range_slider(%0A%20%20%20%20%20%20%20%20start%3D0%2C%0A%20%20%20%20%20%20%20%20stop%3D23%2C%0A%20%20%20%20%20%20%20%20step%3D1%2C%0A%20%20%20%20%20%20%20%20value%3D%5B0%2C%2023%5D%2C%0A%20%20%20%20%20%20%20%20label%3D%22Filter%20by%20hour%20range%22%2C%0A%20%20%20%20)%0A%20%20%20%20hour_filter%0A%20%20%20%20return%20(hour_filter%2C)%0A%0A%0A%40app.cell%0Adef%20_(hour_filter%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_min_hour%2C%20_max_hour%20%3D%20hour_filter.value%0A%0A%20%20%20%20filtered_by_hour%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20(pl.col(%22hour%22)%20%3E%3D%20_min_hour)%20%26%20(pl.col(%22hour%22)%20%3C%3D%20_max_hour)%0A%20%20%20%20)%0A%0A%20%20%20%20_avg_fare%20%3D%20filtered_by_hour%5B%22fare_amount%22%5D.mean()%0A%20%20%20%20_avg_dist%20%3D%20filtered_by_hour%5B%22trip_distance%22%5D.mean()%0A%20%20%20%20_count%20%3D%20len(filtered_by_hour)%0A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20**Hours%20%7B_min_hour%7D%3A00%20-%20%7B_max_hour%7D%3A00**%0A%0A%20%20%20%20-%20Trips%3A%20%7B_count%3A%2C%7D%0A%20%20%20%20-%20Avg%20Fare%3A%20%24%7B_avg_fare%3A.2f%7D%0A%20%20%20%20-%20Avg%20Distance%3A%20%7B_avg_dist%3A.1f%7D%20miles%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Fare%20Analysis%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20%23%20Filter%20reasonable%20fares%0A%20%20%20%20fare_data%20%3D%20yellow_df.filter((pl.col(%22fare_amount%22)%20%3E%3D%202.50)%20%26%20(pl.col(%22fare_amount%22)%20%3C%3D%20100))%0A%0A%20%20%20%20_median%20%3D%20fare_data%5B%22fare_amount%22%5D.median()%0A%20%20%20%20_mean%20%3D%20fare_data%5B%22fare_amount%22%5D.mean()%0A%0A%20%20%20%20fig_fare_hist%20%3D%20px.histogram(%0A%20%20%20%20%20%20%20%20fare_data.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22fare_amount%22%2C%0A%20%20%20%20%20%20%20%20nbins%3D50%2C%0A%20%20%20%20%20%20%20%20title%3Df%22Fare%20Distribution%20(Median%3A%20%24%7B_median%3A.2f%7D%2C%20Mean%3A%20%24%7B_mean%3A.2f%7D)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22fare_amount%22%3A%20%22Fare%20(%24)%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_fare_hist.add_vline(%0A%20%20%20%20%20%20%20%20x%3D_median%2C%20line_dash%3D%22dash%22%2C%20line_color%3D%22red%22%2C%20annotation_text%3Df%22Median%3A%20%24%7B_median%3A.2f%7D%22%0A%20%20%20%20)%0A%20%20%20%20fig_fare_hist.add_vline(%0A%20%20%20%20%20%20%20%20x%3D_mean%2C%20line_dash%3D%22dot%22%2C%20line_color%3D%22blue%22%2C%20annotation_text%3Df%22Mean%3A%20%24%7B_mean%3A.2f%7D%22%0A%20%20%20%20)%0A%0A%20%20%20%20fig_fare_hist%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20%23%20Calculate%20average%20of%20each%20component%0A%20%20%20%20components%20%3D%20yellow_df.select(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22fare_amount%22).mean().alias(%22Base%20Fare%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22tip_amount%22).mean().alias(%22Tip%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22tolls_amount%22).mean().alias(%22Tolls%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22mta_tax%22).mean().alias(%22MTA%20Tax%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22congestion_surcharge%22).fill_null(0).mean().alias(%22Congestion%20Surcharge%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22cbd_congestion_fee%22).fill_null(0).mean().alias(%22CBD%20Fee%20(New)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22improvement_surcharge%22).mean().alias(%22Improvement%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Reshape%20for%20plotting%0A%20%20%20%20component_data%20%3D%20pl.DataFrame(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22component%22%3A%20components.columns%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22amount%22%3A%20components.row(0)%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20)%0A%0A%20%20%20%20fig_components%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20component_data.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22component%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22amount%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Average%20Fare%20Components%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22component%22%3A%20%22Component%22%2C%20%22amount%22%3A%20%22Average%20Amount%20(%24)%22%7D%2C%0A%20%20%20%20%20%20%20%20color%3D%22component%22%2C%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_components%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**Total%20Average%20Fare**%3A%20%24%7Bsum(components.row(0))%3A.2f%7D%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20PAYMENT_TYPES%20%3D%20%7B%0A%20%20%20%20%20%20%20%200%3A%20%22Flex%20Fare%22%2C%0A%20%20%20%20%20%20%20%201%3A%20%22Credit%20Card%22%2C%0A%20%20%20%20%20%20%20%202%3A%20%22Cash%22%2C%0A%20%20%20%20%20%20%20%203%3A%20%22No%20Charge%22%2C%0A%20%20%20%20%20%20%20%204%3A%20%22Dispute%22%2C%0A%20%20%20%20%20%20%20%205%3A%20%22Unknown%22%2C%0A%20%20%20%20%20%20%20%206%3A%20%22Voided%22%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20tip_by_payment%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(pl.col(%22fare_amount%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22tip_amount%22)%20%2F%20pl.col(%22fare_amount%22)%20*%20100).alias(%22tip_pct%22))%0A%20%20%20%20%20%20%20%20.group_by(%22payment_type%22)%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22tip_pct%22).mean().alias(%22avg_tip_pct%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22tip_amount%22).mean().alias(%22avg_tip%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.sort(%22trips%22%2C%20descending%3DTrue)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22payment_type%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20.map_elements(lambda%20x%3A%20PAYMENT_TYPES.get(x%2C%20%22Unknown%22)%2C%20return_dtype%3Dpl.Utf8)%0A%20%20%20%20%20%20%20%20%20%20%20%20.alias(%22payment_name%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20fig_tip%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20tip_by_payment.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22payment_name%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22avg_tip_pct%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Average%20Tip%20%25%20by%20Payment%20Type%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22payment_name%22%3A%20%22Payment%20Type%22%2C%20%22avg_tip_pct%22%3A%20%22Avg%20Tip%20%25%22%7D%2C%0A%20%20%20%20%20%20%20%20text%3D%22avg_tip_pct%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_tip.update_traces(texttemplate%3D%22%25%7Btext%3A.1f%7D%25%22%2C%20textposition%3D%22outside%22)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_tip%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(mo.md(%22**Note**%3A%20Cash%20tips%20are%20not%20recorded%20in%20the%20data%22)%2C%20kind%3D%22info%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(CONGESTION_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20%23%20Sample%20for%20performance%0A%20%20%20%20scatter_data%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20(pl.col(%22trip_distance%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22trip_distance%22)%20%3C%2030)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22fare_amount%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22fare_amount%22)%20%3C%20100)%0A%20%20%20%20).sample(n%3Dmin(5000%2C%20len(yellow_df)))%0A%0A%20%20%20%20fig_scatter%20%3D%20px.scatter(%0A%20%20%20%20%20%20%20%20scatter_data.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22trip_distance%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22fare_amount%22%2C%0A%20%20%20%20%20%20%20%20color%3D%22post_congestion%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Fare%20vs%20Distance%20(sampled)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22trip_distance%22%3A%20%22Distance%20(miles)%22%2C%20%22fare_amount%22%3A%20%22Fare%20(%24)%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_map%3D%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20True%3A%20CONGESTION_COLORS%5B%22post_congestion%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20False%3A%20CONGESTION_COLORS%5B%22pre_congestion%22%5D%2C%0A%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20opacity%3D0.5%2C%0A%20%20%20%20%20%20%20%20trendline%3D%22ols%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_scatter%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20fare_filters%20%3D%20mo.ui.dictionary(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22min_fare%22%3A%20mo.ui.number(start%3D0%2C%20stop%3D100%2C%20value%3D0%2C%20label%3D%22Min%20fare%20(%24)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22max_fare%22%3A%20mo.ui.number(start%3D0%2C%20stop%3D500%2C%20value%3D100%2C%20label%3D%22Max%20fare%20(%24)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22min_distance%22%3A%20mo.ui.number(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start%3D0%2C%20stop%3D50%2C%20value%3D0%2C%20step%3D0.5%2C%20label%3D%22Min%20distance%20(mi)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22max_distance%22%3A%20mo.ui.number(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start%3D0%2C%20stop%3D50%2C%20value%3D30%2C%20step%3D0.5%2C%20label%3D%22Max%20distance%20(mi)%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20)%0A%0A%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_filters%5B%22min_fare%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_filters%5B%22max_fare%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_filters%5B%22min_distance%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_filters%5B%22max_distance%22%5D%2C%0A%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20gap%3D1%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(fare_filters%2C)%0A%0A%0A%40app.cell%0Adef%20_(fare_filters%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_filters%20%3D%20fare_filters.value%0A%0A%20%20%20%20filtered_trips%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20(pl.col(%22fare_amount%22)%20%3E%3D%20_filters%5B%22min_fare%22%5D)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22fare_amount%22)%20%3C%3D%20_filters%5B%22max_fare%22%5D)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22trip_distance%22)%20%3E%3D%20_filters%5B%22min_distance%22%5D)%0A%20%20%20%20%20%20%20%20%26%20(pl.col(%22trip_distance%22)%20%3C%3D%20_filters%5B%22max_distance%22%5D)%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**%7Blen(filtered_trips)%3A%2C%7D%20trips**%20matching%20filters%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20filtered_trips.select(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pickup_datetime%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trip_distance%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22fare_amount%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22tip_amount%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22cbd_congestion_fee%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22total_amount%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20).head(20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(BASE_URL%2C%20conn%2C%20mo%2C%20pl)%3A%0A%20%20%20%20jan_2024_df%20%3D%20conn.execute(f%22%22%22%0A%20%20%20%20%20%20%20%20SELECT%20%0A%20%20%20%20%20%20%20%20%20%20%20%20tpep_pickup_datetime%20as%20pickup_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tpep_dropoff_datetime%20as%20dropoff_datetime%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20PULocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20DOLocationID%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20trip_distance%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fare_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20total_amount%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20congestion_surcharge%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20passenger_count%0A%20%20%20%20%20%20%20%20FROM%20read_parquet('%7BBASE_URL%7D%2Fyellow_tripdata_2024-01.parquet')%0A%20%20%20%20%20%20%20%20WHERE%20tpep_pickup_datetime%20%3E%3D%20'2024-01-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20tpep_pickup_datetime%20%3C%20'2024-02-01'%0A%20%20%20%20%20%20%20%20%20%20AND%20fare_amount%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20AND%20trip_distance%20%3E%200%0A%20%20%20%20%20%20%20%20USING%20SAMPLE%2010%20PERCENT%0A%20%20%20%20%22%22%22).pl()%0A%0A%20%20%20%20%23%20Add%20derived%20columns%0A%20%20%20%20jan_2024_df%20%3D%20jan_2024_df.with_columns(%5B%0A%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.day().alias(%22day%22)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.weekday().alias(%22weekday%22)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.date().alias(%22trip_date%22)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22pickup_datetime%22).dt.hour().alias(%22hour%22)%2C%20%20%23%20Add%20hour%20too%20for%20consistency%0A%20%20%20%20%20%20%20%20pl.lit(2024).alias(%22year%22)%2C%0A%20%20%20%20%5D)%0A%0A%20%20%20%20mo.md(f%22**Jan%202024%20baseline**%3A%20%7Blen(jan_2024_df)%3A%2C%7D%20trips%20loaded%20(10%25%20sample)%22)%0A%20%20%20%20return%20(jan_2024_df%2C)%0A%0A%0A%40app.cell%0Adef%20_(jan_2024_df%2C%20mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20%23%20Jan%202025%20daily%20counts%20(filter%20to%20Jan%20only)%0A%20%20%20%20_jan_2025_daily%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(pl.col(%22pickup_datetime%22).dt.year()%20%3D%3D%202025)%0A%20%20%20%20%20%20%20%20.group_by(pl.col(%22pickup_datetime%22).dt.day().alias(%22day%22))%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips_2025%22))%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Jan%202024%20daily%20counts%0A%20%20%20%20_jan_2024_daily%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df%0A%20%20%20%20%20%20%20%20.group_by(%22day%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips_2024%22))%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Join%20for%20comparison%0A%20%20%20%20yoy_daily%20%3D%20(%0A%20%20%20%20%20%20%20%20_jan_2025_daily%0A%20%20%20%20%20%20%20%20.join(_jan_2024_daily%2C%20on%3D%22day%22%2C%20how%3D%22inner%22)%0A%20%20%20%20%20%20%20%20.sort(%22day%22)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20((pl.col(%22trips_2025%22).cast(pl.Float64)%20-%20pl.col(%22trips_2024%22).cast(pl.Float64))%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%2F%20pl.col(%22trips_2024%22).cast(pl.Float64)%20*%20100.0)%0A%20%20%20%20%20%20%20%20%20%20%20%20.alias(%22pct_change%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20fig_yoy%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20yoy_daily.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22day%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22pct_change%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Daily%20Trip%20Volume%20Change%3A%20Jan%202025%20vs%20Jan%202024%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22day%22%3A%20%22Day%20of%20January%22%2C%20%22pct_change%22%3A%20%22%25%20Change%20from%202024%22%7D%2C%0A%20%20%20%20%20%20%20%20color%3D%22pct_change%22%2C%0A%20%20%20%20%20%20%20%20color_continuous_scale%3D%22RdYlGn%22%2C%0A%20%20%20%20%20%20%20%20color_continuous_midpoint%3D0%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_yoy.add_vline(x%3D5%2C%20line_dash%3D%22dash%22%2C%20line_color%3D%22black%22%2C%20annotation_text%3D%22Congestion%20Pricing%20Starts%22)%0A%20%20%20%20fig_yoy.update_layout(xaxis%3Ddict(tickmode%3D%22linear%22%2C%20dtick%3D1))%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20fig_yoy%2C%0A%20%20%20%20%20%20%20%20mo.md(%22*Negative%20%3D%20fewer%20trips%20in%202025%2C%20Positive%20%3D%20more%20trips.%20Days%201-3%20affected%20by%20New%20Year's%20holiday%20timing.*%22)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%20(yoy_daily%2C)%0A%0A%0A%40app.cell%0Adef%20_(datetime%2C%20jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_post_week_2025%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22)%20%3E%3D%20datetime(2025%2C%201%2C%206))%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22pickup_datetime%22)%20%3C%20datetime(2025%2C%201%2C%2013))%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(pl.col(%22weekday%22))%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22total_amount%22).mean().alias(%22avg_fare_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22trip_distance%22).mean().alias(%22avg_dist_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Same%20weekdays%20in%202024%20(Jan%208-14%2C%202024%20%3D%20Mon-Sun)%0A%20%20%20%20_baseline_week_2024%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22)%20%3E%3D%20datetime(2024%2C%201%2C%208))%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22pickup_datetime%22)%20%3C%20datetime(2024%2C%201%2C%2015))%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22weekday%22)%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22total_amount%22).mean().alias(%22avg_fare_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22trip_distance%22).mean().alias(%22avg_dist_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20DOW_NAME_VALUES%20%3D%20%7B1%3A%20%22Mon%22%2C%202%3A%20%22Tue%22%2C%203%3A%20%22Wed%22%2C%204%3A%20%22Thu%22%2C%205%3A%20%22Fri%22%2C%206%3A%20%22Sat%22%2C%207%3A%20%22Sun%22%7D%0A%0A%20%20%20%20weekday_yoy%20%3D%20(%0A%20%20%20%20%20%20%20%20_post_week_2025.join(_baseline_week_2024%2C%20on%3D%22weekday%22%2C%20how%3D%22full%22)%0A%20%20%20%20%20%20%20%20.sort(%22weekday%22)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20((pl.col(%22trips_2025%22)%20-%20pl.col(%22trips_2024%22))%20%2F%20pl.col(%22trips_2024%22)%20*%20100).alias(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trip_change_pct%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22avg_fare_2025%22)%20-%20pl.col(%22avg_fare_2024%22)).alias(%22fare_change%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22weekday%22).map_elements(lambda%20x%3A%20DOW_NAME_VALUES.get(x%2C%20%22%3F%22)%2C%20return_dtype%3Dpl.Utf8).alias(%22day_name%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Same-Weekday%20Comparison%20(First%20Full%20Week%20Post-Congestion)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*Jan%206-12%2C%202025%20vs%20Jan%208-14%2C%202024%20(both%20Mon-Sun)*%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20weekday_yoy.select(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22day_name%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trips_2024%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trips_2025%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trip_change_pct%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22avg_fare_2024%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22avg_fare_2025%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22fare_change%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%20(weekday_yoy%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20px%2C%20weekday_yoy)%3A%0A%20%20%20%20fig_weekday%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20weekday_yoy.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22day_name%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22trip_change_pct%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Trip%20Volume%20Change%20by%20Day%20of%20Week%20(2025%20vs%202024)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22day_name%22%3A%20%22Day%22%2C%20%22trip_change_pct%22%3A%20%22%25%20Change%22%7D%2C%0A%20%20%20%20%20%20%20%20color%3D%22trip_change_pct%22%2C%0A%20%20%20%20%20%20%20%20color_continuous_scale%3D%22RdYlGn%22%2C%0A%20%20%20%20%20%20%20%20color_continuous_midpoint%3D0%2C%0A%20%20%20%20%20%20%20%20text%3D%22trip_change_pct%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_weekday.update_traces(texttemplate%3D%22%25%7Btext%3A.1f%7D%25%22%2C%20textposition%3D%22outside%22)%0A%0A%20%20%20%20_avg_change%20%3D%20weekday_yoy%5B%22trip_change_pct%22%5D.mean()%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_weekday%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**Average%20change%20across%20the%20week%3A%20%7B_avg_change%3A.1f%7D%25**%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22success%22%20if%20_avg_change%20%3E%200%20else%20%22warn%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(CONGESTION_COLORS%2C%20mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20_post_df%20%3D%20yellow_df.filter(pl.col(%22post_congestion%22))%0A%0A%20%20%20%20fee_dist%20%3D%20(%0A%20%20%20%20%20%20%20%20_post_df.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.when(pl.col(%22cbd_congestion_fee%22).is_null()%20%7C%20(pl.col(%22cbd_congestion_fee%22)%20%3D%3D%200))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.then(pl.lit(%22%240%20(No%20fee)%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.otherwise(pl.lit(%22%240.75%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.alias(%22fee_category%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22fee_category%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22trips%22)%20%2F%20pl.col(%22trips%22).sum()%20*%20100).alias(%22pct%22))%0A%20%20%20%20)%0A%0A%20%20%20%20fig_fee%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20fee_dist.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22fee_category%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22CBD%20Congestion%20Fee%20Distribution%20(Post%20Jan%205%2C%202025)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22fee_category%22%3A%20%22Fee%22%2C%20%22trips%22%3A%20%22Trip%20Count%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BCONGESTION_COLORS%5B%22post_congestion%22%5D%5D%2C%0A%20%20%20%20%20%20%20%20text%3D%22pct%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_fee.update_traces(texttemplate%3D%22%25%7Btext%3A.1f%7D%25%22%2C%20textposition%3D%22outside%22)%0A%0A%20%20%20%20_pct_with_fee%20%3D%20_post_df.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200).height%20%2F%20len(_post_df)%20*%20100%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20fig_fee%2C%0A%20%20%20%20%20%20%20%20mo.md(f%22**%7B_pct_with_fee%3A.1f%7D%25**%20of%20post-Jan%205%20trips%20have%20a%20CBD%20fee%20(%240.75)%22)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(MANHATTAN_CBD_ZONES%2C%20jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_cbd_zones%20%3D%20set(MANHATTAN_CBD_ZONES)%0A%0A%20%20%20%20%23%202025%20CBD%20trips%20(post-congestion%20only)%0A%20%20%20%20_jan_2025_cbd%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22PULocationID%22).is_in(_cbd_zones)%20%7C%20pl.col(%22DOLocationID%22).is_in(_cbd_zones)%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22cbd_trip%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22cbd_trip%22)%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22total_amount%22).mean().alias(%22avg_fare_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20%23%202024%20CBD%20trips%20(same%20date%20range%3A%20Jan%205-31)%0A%20%20%20%20_jan_2024_cbd%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df.filter(pl.col(%22day%22)%20%3E%3D%205)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22PULocationID%22).is_in(_cbd_zones)%20%7C%20pl.col(%22DOLocationID%22).is_in(_cbd_zones)%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22cbd_trip%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22cbd_trip%22)%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22total_amount%22).mean().alias(%22avg_fare_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20cbd_yoy%20%3D%20(%0A%20%20%20%20%20%20%20%20_jan_2025_cbd.join(_jan_2024_cbd%2C%20on%3D%22cbd_trip%22%2C%20how%3D%22full%22)%0A%20%20%20%20%20%20%20%20.with_columns(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20((pl.col(%22trips_2025%22).cast(pl.Float64)%20-%20pl.col(%22trips_2024%22).cast(pl.Float64))%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%2F%20pl.col(%22trips_2024%22).cast(pl.Float64)%20*%20100.0).alias(%22trip_change_pct%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22avg_fare_2025%22)%20-%20pl.col(%22avg_fare_2024%22)).alias(%22fare_change%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.when(pl.col(%22cbd_trip%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.then(pl.lit(%22CBD%20Trip%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.otherwise(pl.lit(%22Non-CBD%20Trip%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20.alias(%22trip_type%22)%2C%0A%20%20%20%20%20%20%20%20%5D)%0A%20%20%20%20)%0A%0A%20%20%20%20_cbd_change%20%3D%20cbd_yoy.filter(pl.col(%22trip_type%22)%20%3D%3D%20%22CBD%20Trip%22)%5B%22trip_change_pct%22%5D%5B0%5D%0A%20%20%20%20_non_cbd_change%20%3D%20cbd_yoy.filter(pl.col(%22trip_type%22)%20%3D%3D%20%22Non-CBD%20Trip%22)%5B%22trip_change_pct%22%5D%5B0%5D%0A%20%20%20%20_diff%20%3D%20_cbd_change%20-%20_non_cbd_change%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20CBD%20vs%20Non-CBD%20Trips%3A%20Year-over-Year%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22*Using%20%7Blen(_cbd_zones)%7D%20zones%20classified%20as%20Manhattan%20CBD*%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20cbd_yoy.select(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22trip_type%22%2C%20%22trips_2024%22%2C%20%22trips_2025%22%2C%20%22trip_change_pct%22%2C%20%22fare_change%22%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22Compare%20CBD%20trip%20change%20vs%20non-CBD%20to%20isolate%20congestion%20pricing%20effect%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22info%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20**Difference-in-Differences%3A%20%7B_diff%3A.1f%7D%20percentage%20points**%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20CBD%20trips%20fell%20while%20non-CBD%20trips%20rose%20%E2%80%94%20consistent%20with%20congestion%20pricing%20%0A%20%20%20%20%20%20%20%20%20%20%20%20discouraging%20travel%20into%20Manhattan%20below%2060th%20St.%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22success%22%20if%20_diff%20%3C%20-5%20else%20%22info%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20green_df%2C%20mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20def%20_calc_stats(df%2C%20name)%3A%0A%20%20%20%20%20%20%20%20_post%20%3D%20df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20%20%20%20%20_with_fee%20%3D%20_post.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20return%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22taxi_type%22%3A%20name%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22total_trips%22%3A%20len(_post)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22trips_with_fee%22%3A%20_with_fee.height%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22pct_with_fee%22%3A%20_with_fee.height%20%2F%20len(_post)%20*%20100%20if%20len(_post)%20%3E%200%20else%200%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22avg_fee%22%3A%20_with_fee%5B%22cbd_congestion_fee%22%5D.mean()%20if%20_with_fee.height%20%3E%200%20else%200%2C%0A%20%20%20%20%20%20%20%20%7D%0A%0A%20%20%20%20taxi_comparison%20%3D%20pl.DataFrame(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20_calc_stats(yellow_df%2C%20%22Yellow%20Taxi%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_calc_stats(green_df%2C%20%22Green%20Taxi%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%0A%20%20%20%20fig_taxi%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20taxi_comparison.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22taxi_type%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22pct_with_fee%22%2C%0A%20%20%20%20%20%20%20%20color%3D%22taxi_type%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22%25%20of%20Trips%20with%20CBD%20Fee%20(Post%20Jan%205)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22taxi_type%22%3A%20%22Taxi%20Type%22%2C%20%22pct_with_fee%22%3A%20%22%25%20with%20Fee%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_map%3DTAXI_COLORS%2C%0A%20%20%20%20%20%20%20%20text%3D%22pct_with_fee%22%2C%0A%20%20%20%20)%0A%20%20%20%20fig_taxi.update_traces(texttemplate%3D%22%25%7Btext%3A.1f%7D%25%22%2C%20textposition%3D%22outside%22)%0A%0A%20%20%20%20fig_taxi.update_layout(yaxis%3Ddict(range%3D%5B0%2C%2085%5D))%20%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_taxi%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*Green%20taxis%20serve%20outer%20boroughs%20-%20fewer%20trips%20enter%20the%20CBD*%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20taxi_comparison%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20%23%20Jan%205%2C%202025%20(Sun)%20vs%20Jan%207%2C%202024%20(Sun)%20-%20same%20weekday%0A%20%20%20%20_day_2025_df%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22).dt.year()%20%3D%3D%202025)%20%26%20%0A%20%20%20%20%20%20%20%20(pl.col(%22pickup_datetime%22).dt.day()%20%3D%3D%205)%0A%20%20%20%20)%0A%20%20%20%20_day_2024_df%20%3D%20jan_2024_df.filter(pl.col(%22day%22)%20%3D%3D%207)%20%20%23%20First%20Sunday%20in%20Jan%202024%0A%0A%20%20%20%20_trips_2025%20%3D%20len(_day_2025_df)%0A%20%20%20%20_trips_2024%20%3D%20len(_day_2024_df)%0A%20%20%20%20_change%20%3D%20(_trips_2025%20-%20_trips_2024)%20%2F%20_trips_2024%20*%20100%0A%0A%20%20%20%20_fare_2025%20%3D%20_day_2025_df%5B%22total_amount%22%5D.mean()%0A%20%20%20%20_fare_2024%20%3D%20_day_2024_df%5B%22total_amount%22%5D.mean()%0A%0A%20%20%20%20_pct_cbd%20%3D%20(%0A%20%20%20%20%20%20%20%20_day_2025_df.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200).height%20%2F%20_trips_2025%20*%20100%0A%20%20%20%20%20%20%20%20if%20_trips_2025%20%3E%200%20else%200%0A%20%20%20%20)%0A%0A%20%20%20%20_change_icon%20%3D%20%22%F0%9F%93%88%22%20if%20_change%20%3E%200%20else%20%22%F0%9F%93%89%22%0A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20%23%23%23%20Congestion%20Pricing%20Launch%20Day%20(Sunday%20Comparison)%0A%0A%20%20%20%20*Jan%205%2C%202025%20vs%20Jan%207%2C%202024%20%E2%80%94%20both%20Sundays*%0A%0A%20%20%20%20%7C%20Metric%20%7C%202024%20%7C%202025%20%7C%20Change%20%7C%0A%20%20%20%20%7C--------%7C------%7C------%7C--------%7C%0A%20%20%20%20%7C%20Trips%20%7C%20%7B_trips_2024%3A%2C%7D%20%7C%20%7B_trips_2025%3A%2C%7D%20%7C%20%7B_change_icon%7D%20%7B_change%3A%2B.1f%7D%25%20%7C%0A%20%20%20%20%7C%20Avg%20Fare%20%7C%20%24%7B_fare_2024%3A.2f%7D%20%7C%20%24%7B_fare_2025%3A.2f%7D%20%7C%20%24%7B_fare_2025%20-%20_fare_2024%3A%2B.2f%7D%20%7C%0A%20%20%20%20%7C%20%25%20with%20CBD%20Fee%20%7C%20N%2FA%20%7C%20%7B_pct_cbd%3A.1f%7D%25%20%7C%20%E2%80%94%20%7C%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20pl%2C%20yoy_daily)%3A%0A%20%20%20%20_pre_avg%20%3D%20yoy_daily.filter(pl.col(%22day%22)%20%3C%205)%5B%22pct_change%22%5D.mean()%0A%20%20%20%20_post_avg%20%3D%20yoy_daily.filter(pl.col(%22day%22)%20%3E%3D%205)%5B%22pct_change%22%5D.mean()%0A%20%20%20%20_diff%20%3D%20_post_avg%20-%20_pre_avg%0A%0A%20%20%20%20_interpretation%20%3D%20(%0A%20%20%20%20%20%20%20%20%22Trip%20volume%20increased%20post-congestion%20vs%202024%20%E2%80%94%20other%20factors%20(weather%2C%20economy%2C%20seasonal%20variation)%20%22%0A%20%20%20%20%20%20%20%20%22may%20outweigh%20the%20%240.75%20congestion%20fee%20impact.%22%0A%20%20%20%20%20%20%20%20if%20_diff%20%3E%200%20else%0A%20%20%20%20%20%20%20%20%22Trip%20volume%20decreased%20post-congestion%20vs%202024%20%E2%80%94%20consistent%20with%20congestion%20pricing%20reducing%20demand.%22%0A%20%20%20%20)%0A%0A%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20**Year-over-Year%20Summary**%0A%0A%20%20%20%20-%20**Pre-congestion%20(Jan%201-4)**%3A%20Avg%20%7B_pre_avg%3A%2B.1f%7D%25%20vs%202024%0A%20%20%20%20-%20**Post-congestion%20(Jan%205-31)**%3A%20Avg%20%7B_post_avg%3A%2B.1f%7D%25%20vs%202024%0A%20%20%20%20-%20**Difference**%3A%20%7B_diff%3A%2B.1f%7D%20percentage%20points%0A%0A%20%20%20%20*%7B_interpretation%7D*%0A%20%20%20%20%20%20%20%20%22%22%22)%2C%0A%20%20%20%20%20%20%20%20kind%3D%22success%22%20if%20_diff%20%3E%200%20else%20%22warn%22%2C%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Trip%20Distance%20and%20Duration%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20_distance_data%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20(pl.col(%22trip_distance%22)%20%3E%3D%200.1)%20%26%20(pl.col(%22trip_distance%22)%20%3C%3D%2030)%0A%20%20%20%20)%0A%0A%20%20%20%20_p25%20%3D%20_distance_data%5B%22trip_distance%22%5D.quantile(0.25)%0A%20%20%20%20_p50%20%3D%20_distance_data%5B%22trip_distance%22%5D.quantile(0.50)%0A%20%20%20%20_p75%20%3D%20_distance_data%5B%22trip_distance%22%5D.quantile(0.75)%0A%20%20%20%20_p90%20%3D%20_distance_data%5B%22trip_distance%22%5D.quantile(0.90)%0A%0A%20%20%20%20fig_dist%20%3D%20px.histogram(%0A%20%20%20%20%20%20%20%20_distance_data.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22trip_distance%22%2C%0A%20%20%20%20%20%20%20%20nbins%3D60%2C%0A%20%20%20%20%20%20%20%20title%3Df%22Trip%20Distance%20Distribution%20(Median%3A%20%7B_p50%3A.1f%7D%20mi)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22trip_distance%22%3A%20%22Distance%20(miles)%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_dist.add_vline(x%3D_p25%2C%20line_dash%3D%22dot%22%2C%20annotation_text%3D%2225th%20%25ile%22)%0A%20%20%20%20fig_dist.add_vline(x%3D_p50%2C%20line_dash%3D%22dash%22%2C%20line_color%3D%22red%22%2C%20annotation_text%3D%22Median%22)%0A%20%20%20%20fig_dist.add_vline(x%3D_p75%2C%20line_dash%3D%22dot%22%2C%20annotation_text%3D%2275th%20%25ile%22)%0A%20%20%20%20fig_dist.add_vline(x%3D_p90%2C%20line_dash%3D%22dot%22%2C%20annotation_text%3D%2290th%20%25ile%22)%0A%0A%20%20%20%20fig_dist%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20%23%202025%20stats%20(post-congestion)%0A%20%20%20%20_dist_2025%20%3D%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20pl.col(%22post_congestion%22)%20%26%20(pl.col(%22pickup_datetime%22).dt.day()%20%3E%3D%206)%0A%20%20%20%20)%5B%22trip_distance%22%5D%0A%0A%20%20%20%20%23%202024%20stats%20(same%20period%3A%20Jan%205-31)%0A%20%20%20%20_dist_2024%20%3D%20jan_2024_df.filter(pl.col(%22day%22)%20%3E%3D%208)%5B%22trip_distance%22%5D%0A%0A%20%20%20%20distance_yoy%20%3D%20pl.DataFrame(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22metric%22%3A%20%5B%22Mean%22%2C%20%22Median%22%2C%20%2225th%20%25ile%22%2C%20%2275th%20%25ile%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2024%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2024.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2024.median()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2024.quantile(0.25)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2024.quantile(0.75)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2025%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2025.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2025.median()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2025.quantile(0.25)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_dist_2025.quantile(0.75)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20).with_columns((pl.col(%22jan_2025%22)%20-%20pl.col(%22jan_2024%22)).alias(%22change%22))%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Trip%20Distance%3A%20Year-over-Year%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*Comparing%20same%20calendar%20dates%3B%20weekday%20mix%20differs%20slightly%20between%20years*%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20distance_yoy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*Positive%20change%20%3D%20longer%20trips%20in%202025*%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20duration_df%20%3D%20yellow_df.with_columns(%0A%20%20%20%20%20%20%20%20((pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%2060).alias(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22duration_mins%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20).filter((pl.col(%22duration_mins%22)%20%3E%201)%20%26%20(pl.col(%22duration_mins%22)%20%3C%20120))%0A%0A%20%20%20%20_median_dur%20%3D%20duration_df%5B%22duration_mins%22%5D.median()%0A%0A%20%20%20%20fig_duration%20%3D%20px.histogram(%0A%20%20%20%20%20%20%20%20duration_df.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22duration_mins%22%2C%0A%20%20%20%20%20%20%20%20nbins%3D60%2C%0A%20%20%20%20%20%20%20%20title%3Df%22Trip%20Duration%20Distribution%20(Median%3A%20%7B_median_dur%3A.0f%7D%20min)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22duration_mins%22%3A%20%22Duration%20(minutes)%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_duration%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**Median%20trip**%3A%20%7B_median_dur%3A.0f%7D%20minutes%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_dur_2025%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%2060%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_mins%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_mins%22)%20%3E%201)%20%26%20(pl.col(%22duration_mins%22)%20%3C%20120))%0A%20%20%20%20)%5B%22duration_mins%22%5D%0A%0A%20%20%20%20%23%202024%20durations%20(same%20period)%0A%20%20%20%20_dur_2024%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df.filter(pl.col(%22day%22)%20%3E%3D%205)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%2060%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_mins%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_mins%22)%20%3E%201)%20%26%20(pl.col(%22duration_mins%22)%20%3C%20120))%0A%20%20%20%20)%5B%22duration_mins%22%5D%0A%0A%20%20%20%20duration_yoy%20%3D%20pl.DataFrame(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22metric%22%3A%20%5B%22Mean%20(mins)%22%2C%20%22Median%20(mins)%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2024%22%3A%20%5B_dur_2024.mean()%2C%20_dur_2024.median()%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2025%22%3A%20%5B_dur_2025.mean()%2C%20_dur_2025.median()%5D%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20).with_columns(%5B%0A%20%20%20%20%20%20%20%20(pl.col(%22jan_2025%22)%20-%20pl.col(%22jan_2024%22)).alias(%22change%22)%2C%0A%20%20%20%20%20%20%20%20((pl.col(%22jan_2025%22)%20-%20pl.col(%22jan_2024%22))%20%2F%20pl.col(%22jan_2024%22)%20*%20100).alias(%22change_pct%22)%2C%0A%20%20%20%20%5D).with_columns(%5B%0A%20%20%20%20%20%20%20%20pl.col(%22jan_2024%22).round(1)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22jan_2025%22).round(1)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22change%22).round(2)%2C%0A%20%20%20%20%20%20%20%20pl.col(%22change_pct%22).round(1)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Trip%20Duration%3A%20Year-over-Year%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20duration_yoy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22Shorter%20durations%20in%202025%20could%20indicate%20less%20congestion%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22info%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20speed_df%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%203600%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_hours%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_hours%22)%20%3E%200.01)%20%26%20(pl.col(%22duration_hours%22)%20%3C%203))%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22trip_distance%22)%20%2F%20pl.col(%22duration_hours%22)).alias(%22speed_mph%22))%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22speed_mph%22)%20%3E%201)%20%26%20(pl.col(%22speed_mph%22)%20%3C%2060))%0A%20%20%20%20)%0A%0A%20%20%20%20_speed_by_hour%20%3D%20(%0A%20%20%20%20%20%20%20%20speed_df.group_by(%22hour%22).agg(pl.col(%22speed_mph%22).mean().alias(%22avg_speed%22)).sort(%22hour%22)%0A%20%20%20%20)%0A%0A%20%20%20%20fig_speed%20%3D%20px.line(%0A%20%20%20%20%20%20%20%20_speed_by_hour.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22hour%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22avg_speed%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Average%20Trip%20Speed%20by%20Hour%20of%20Day%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22hour%22%3A%20%22Hour%22%2C%20%22avg_speed%22%3A%20%22Avg%20Speed%20(mph)%22%7D%2C%0A%20%20%20%20%20%20%20%20markers%3DTrue%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_speed.add_vrect(x0%3D7%2C%20x1%3D9%2C%20fillcolor%3D%22red%22%2C%20opacity%3D0.1%2C%20annotation_text%3D%22AM%20Rush%22)%0A%20%20%20%20fig_speed.add_vrect(x0%3D17%2C%20x1%3D19%2C%20fillcolor%3D%22red%22%2C%20opacity%3D0.1%2C%20annotation_text%3D%22PM%20Rush%22)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_speed%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(mo.md(%22Rush%20hour%20slowdowns%20clearly%20visible!%22)%2C%20kind%3D%22info%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(CONGESTION_COLORS%2C%20jan_2024_df%2C%20mo%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20_speed_2025%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%203600%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_hours%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_hours%22)%20%3E%200.01)%20%26%20(pl.col(%22duration_hours%22)%20%3C%203))%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22trip_distance%22)%20%2F%20pl.col(%22duration_hours%22)).alias(%22speed_mph%22))%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22speed_mph%22)%20%3E%201)%20%26%20(pl.col(%22speed_mph%22)%20%3C%2060))%0A%20%20%20%20%20%20%20%20.group_by(%22hour%22)%0A%20%20%20%20%20%20%20%20.agg(pl.col(%22speed_mph%22).mean().alias(%22speed_2025%22))%0A%20%20%20%20)%0A%0A%20%20%20%20%23%202024%20speed%20by%20hour%20(same%20period)%0A%20%20%20%20_speed_2024%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df.filter(pl.col(%22day%22)%20%3E%3D%205)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%20%2F%203600%0A%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_hours%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_hours%22)%20%3E%200.01)%20%26%20(pl.col(%22duration_hours%22)%20%3C%203))%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22trip_distance%22)%20%2F%20pl.col(%22duration_hours%22)).alias(%22speed_mph%22))%0A%20%20%20%20%20%20%20%20.filter((pl.col(%22speed_mph%22)%20%3E%201)%20%26%20(pl.col(%22speed_mph%22)%20%3C%2060))%0A%20%20%20%20%20%20%20%20.group_by(pl.col(%22pickup_datetime%22).dt.hour().alias(%22hour%22))%0A%20%20%20%20%20%20%20%20.agg(pl.col(%22speed_mph%22).mean().alias(%22speed_2024%22))%0A%20%20%20%20)%0A%0A%20%20%20%20speed_yoy_hourly%20%3D%20(%0A%20%20%20%20%20%20%20%20_speed_2025.join(_speed_2024%2C%20on%3D%22hour%22%2C%20how%3D%22full%22)%0A%20%20%20%20%20%20%20%20.sort(%22hour%22)%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22speed_2025%22)%20-%20pl.col(%22speed_2024%22)).alias(%22speed_change%22))%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Melt%20for%20plotting%0A%20%20%20%20_plot_data%20%3D%20pl.concat(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20speed_yoy_hourly.select(%5B%22hour%22%2C%20%22speed_2024%22%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20.rename(%7B%22speed_2024%22%3A%20%22speed%22%7D)%0A%20%20%20%20%20%20%20%20%20%20%20%20.with_columns(pl.lit(%222024%22).alias(%22year%22))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20speed_yoy_hourly.select(%5B%22hour%22%2C%20%22speed_2025%22%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20.rename(%7B%22speed_2025%22%3A%20%22speed%22%7D)%0A%20%20%20%20%20%20%20%20%20%20%20%20.with_columns(pl.lit(%222025%22).alias(%22year%22))%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%0A%20%20%20%20fig_speed_yoy%20%3D%20px.line(%0A%20%20%20%20%20%20%20%20_plot_data.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22hour%22%2C%0A%20%20%20%20%20%20%20%20y%3D%22speed%22%2C%0A%20%20%20%20%20%20%20%20color%3D%22year%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Average%20Speed%20by%20Hour%3A%202025%20vs%202024%20(Jan%205-31)%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22hour%22%3A%20%22Hour%22%2C%20%22speed%22%3A%20%22Avg%20Speed%20(mph)%22%2C%20%22year%22%3A%20%22Year%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_map%3D%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%222024%22%3A%20CONGESTION_COLORS%5B%22pre_congestion%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%222025%22%3A%20CONGESTION_COLORS%5B%22post_congestion%22%5D%2C%0A%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20markers%3DTrue%2C%0A%20%20%20%20)%0A%0A%20%20%20%20fig_speed_yoy.add_vrect(x0%3D7%2C%20x1%3D9%2C%20fillcolor%3D%22gray%22%2C%20opacity%3D0.1)%0A%20%20%20%20fig_speed_yoy.add_vrect(x0%3D17%2C%20x1%3D19%2C%20fillcolor%3D%22gray%22%2C%20opacity%3D0.1)%0A%0A%20%20%20%20_avg_improvement%20%3D%20speed_yoy_hourly%5B%22speed_change%22%5D.mean()%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig_speed_yoy%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22**Average%20speed%20change%3A%20%7B_avg_improvement%3A%2B.1f%7D%20mph**%20across%20all%20hours%20*(aggregated%20over%2027%20days)*%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(MANHATTAN_CBD_ZONES%2C%20jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_cbd_zones%20%3D%20set(MANHATTAN_CBD_ZONES)%0A%0A%20%20%20%20def%20_calc_cbd_speed(df%2C%20day_filter%3DNone)%3A%0A%20%20%20%20%20%20%20%20_filtered%20%3D%20df%0A%20%20%20%20%20%20%20%20if%20day_filter%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_filtered%20%3D%20df.filter(day_filter)%0A%20%20%20%20%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20_filtered.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22PULocationID%22).is_in(_cbd_zones)%20%7C%20pl.col(%22DOLocationID%22).is_in(_cbd_zones)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22dropoff_datetime%22)%20-%20pl.col(%22pickup_datetime%22)).dt.total_seconds()%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2F%203600%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20).alias(%22duration_hours%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20.filter((pl.col(%22duration_hours%22)%20%3E%200.01)%20%26%20(pl.col(%22duration_hours%22)%20%3C%203))%0A%20%20%20%20%20%20%20%20%20%20%20%20.with_columns((pl.col(%22trip_distance%22)%20%2F%20pl.col(%22duration_hours%22)).alias(%22speed_mph%22))%0A%20%20%20%20%20%20%20%20%20%20%20%20.filter((pl.col(%22speed_mph%22)%20%3E%201)%20%26%20(pl.col(%22speed_mph%22)%20%3C%2060))%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%23%202025%20CBD%20speeds%20(post-congestion)%0A%20%20%20%20_cbd_2025%20%3D%20_calc_cbd_speed(yellow_df%2C%20pl.col(%22post_congestion%22))%0A%0A%20%20%20%20%23%202024%20CBD%20speeds%20(same%20period)%0A%20%20%20%20_cbd_2024%20%3D%20_calc_cbd_speed(jan_2024_df%2C%20pl.col(%22day%22)%20%3E%3D%205)%0A%0A%20%20%20%20%23%20Cast%20everything%20to%20float%20to%20avoid%20type%20mismatch%0A%20%20%20%20cbd_speed_yoy%20%3D%20pl.DataFrame(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22metric%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Trips%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Mean%20Speed%20(mph)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Median%20Speed%20(mph)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Rush%20Hour%20Mean%20(7-9%2C%2017-19)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2024%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20float(len(_cbd_2024))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2024%5B%22speed_mph%22%5D.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2024%5B%22speed_mph%22%5D.median()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2024.filter(pl.col(%22pickup_datetime%22).dt.hour().is_in(%5B7%2C%208%2C%2017%2C%2018%5D))%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22speed_mph%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22jan_2025%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20float(len(_cbd_2025))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2025%5B%22speed_mph%22%5D.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2025%5B%22speed_mph%22%5D.median()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_cbd_2025.filter(pl.col(%22hour%22).is_in(%5B7%2C%208%2C%2017%2C%2018%5D))%5B%22speed_mph%22%5D.mean()%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20).with_columns(%5B%0A%20%20%20%20%20%20%20%20(pl.col(%22jan_2025%22)%20-%20pl.col(%22jan_2024%22)).alias(%22change%22)%2C%0A%20%20%20%20%5D).with_columns(%5B%0A%20%20%20%20%20%20%20%20%23%20Round%20for%20display%20(trips%20as%20int-like%2C%20speeds%20to%201%20decimal)%0A%20%20%20%20%20%20%20%20pl.when(pl.col(%22metric%22)%20%3D%3D%20%22Trips%22)%0A%20%20%20%20%20%20%20%20%20%20.then(pl.col(%22jan_2024%22).round(0))%0A%20%20%20%20%20%20%20%20%20%20.otherwise(pl.col(%22jan_2024%22).round(1))%0A%20%20%20%20%20%20%20%20%20%20.alias(%22jan_2024%22)%2C%0A%20%20%20%20%20%20%20%20pl.when(pl.col(%22metric%22)%20%3D%3D%20%22Trips%22)%0A%20%20%20%20%20%20%20%20%20%20.then(pl.col(%22jan_2025%22).round(0))%0A%20%20%20%20%20%20%20%20%20%20.otherwise(pl.col(%22jan_2025%22).round(1))%0A%20%20%20%20%20%20%20%20%20%20.alias(%22jan_2025%22)%2C%0A%20%20%20%20%20%20%20%20pl.when(pl.col(%22metric%22)%20%3D%3D%20%22Trips%22)%0A%20%20%20%20%20%20%20%20%20%20.then(pl.col(%22change%22).round(0))%0A%20%20%20%20%20%20%20%20%20%20.otherwise(pl.col(%22change%22).round(2))%0A%20%20%20%20%20%20%20%20%20%20.alias(%22change%22)%2C%0A%20%20%20%20%5D)%0A%0A%20%20%20%20_speed_improvement%20%3D%20cbd_speed_yoy.filter(pl.col(%22metric%22)%20%3D%3D%20%22Mean%20Speed%20(mph)%22)%5B%22change%22%5D%5B0%5D%0A%20%20%20%20_rush_improvement%20%3D%20cbd_speed_yoy.filter(pl.col(%22metric%22)%20%3D%3D%20%22Rush%20Hour%20Mean%20(7-9%2C%2017-19)%22)%5B%22change%22%5D%5B0%5D%0A%0A%20%20%20%20%23%20Determine%20significance%20(less%20than%200.5%20mph%20is%20essentially%20no%20change)%0A%20%20%20%20_is_significant%20%3D%20abs(_speed_improvement)%20%3E%3D%200.5%0A%0A%20%20%20%20if%20_is_significant%3A%0A%20%20%20%20%20%20%20%20_interpretation%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Positive%20change%20suggests%20reduced%20congestion!%22%20if%20_speed_improvement%20%3E%200%20%0A%20%20%20%20%20%20%20%20%20%20%20%20else%20%22Negative%20change%20-%20congestion%20may%20have%20worsened%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_kind%20%3D%20%22success%22%20if%20_speed_improvement%20%3E%200%20else%20%22warn%22%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_interpretation%20%3D%20%22No%20significant%20speed%20change%20detected%20in%20the%20first%20month%20of%20congestion%20pricing.%22%0A%20%20%20%20%20%20%20%20_kind%20%3D%20%22info%22%0A%0A%20%20%20%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20CBD%20Zone%20Speed%3A%20Year-over-Year%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20mo.md(%22*This%20is%20a%20key%20metric%20for%20congestion%20pricing%20effectiveness*%22)%2C%0A%20%20%20%20%20%20%20%20cbd_speed_yoy%2C%0A%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20**Key%20Finding%3A**%0A%0A%20%20%20%20-%20Overall%20CBD%20speed%20change%3A%20**%7B_speed_improvement%3A%2B.2f%7D%20mph**%0A%20%20%20%20-%20Rush%20hour%20speed%20change%3A%20**%7B_rush_improvement%3A%2B.2f%7D%20mph**%0A%0A%20%20%20%20%7B_interpretation%7D%0A%0A%20%20%20%20*Note%3A%20Speed%20improvements%20may%20take%20months%20to%20materialise%20as%20driver%20behaviour%20adjusts.*%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20kind%3D_kind%2C%0A%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%5D)%0A%20%20%20%20return%20(cbd_speed_yoy%2C)%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Passenger%20Analysis%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(TAXI_COLORS%2C%20pl%2C%20px%2C%20yellow_df)%3A%0A%20%20%20%20_passenger_dist%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22passenger_count%22).is_not_null()%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3C%3D%206)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips%22))%0A%20%20%20%20%20%20%20%20.sort(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22passenger_count%22).cast(pl.Utf8).alias(%22passengers%22)%20%20%23%20Convert%20to%20string%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20)%0A%0A%20%20%20%20fig_pax%20%3D%20px.bar(%0A%20%20%20%20%20%20%20%20_passenger_dist.to_pandas()%2C%0A%20%20%20%20%20%20%20%20x%3D%22passengers%22%2C%20%20%23%20Use%20string%20column%0A%20%20%20%20%20%20%20%20y%3D%22trips%22%2C%0A%20%20%20%20%20%20%20%20title%3D%22Trips%20by%20Passenger%20Count%22%2C%0A%20%20%20%20%20%20%20%20labels%3D%7B%22passengers%22%3A%20%22Passengers%22%2C%20%22trips%22%3A%20%22Trip%20Count%22%7D%2C%0A%20%20%20%20%20%20%20%20color_discrete_sequence%3D%5BTAXI_COLORS%5B%22Yellow%20Taxi%22%5D%5D%2C%0A%20%20%20%20)%0A%20%20%20%20fig_pax%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(jan_2024_df%2C%20mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_pax_2025%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22post_congestion%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20pl.col(%22passenger_count%22).is_not_null()%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3C%3D%206)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips_2025%22))%0A%20%20%20%20)%0A%0A%20%20%20%20%23%202024%20passenger%20distribution%20(same%20period)%0A%20%20%20%20_pax_2024%20%3D%20(%0A%20%20%20%20%20%20%20%20jan_2024_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22day%22)%20%3E%3D%205)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20pl.col(%22passenger_count%22).is_not_null()%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3C%3D%206)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.group_by(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.agg(pl.len().alias(%22trips_2024%22))%0A%20%20%20%20)%0A%0A%20%20%20%20pax_yoy%20%3D%20(%0A%20%20%20%20%20%20%20%20_pax_2025.join(_pax_2024%2C%20on%3D%22passenger_count%22%2C%20how%3D%22full%22)%0A%20%20%20%20%20%20%20%20.sort(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.with_columns(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22trips_2025%22)%20%2F%20pl.col(%22trips_2025%22).sum()%20*%20100).alias(%22pct_2025%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22trips_2024%22)%20%2F%20pl.col(%22trips_2024%22).sum()%20*%20100).alias(%22pct_2024%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22pct_2025%22)%20-%20pl.col(%22pct_2024%22)).alias(%22pct_change%22))%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Passenger%20Count%3A%20Year-over-Year%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20pax_yoy.select(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22passenger_count%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trips_2024%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trips_2025%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct_2024%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct_2025%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct_change%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22*Shows%20%25%20of%20trips%20by%20passenger%20count%20and%20year-over-year%20shift*%22)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20pl%2C%20yellow_df)%3A%0A%20%20%20%20_pax_fare%20%3D%20(%0A%20%20%20%20%20%20%20%20yellow_df.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22passenger_count%22).is_not_null()%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3E%200)%0A%20%20%20%20%20%20%20%20%20%20%20%20%26%20(pl.col(%22passenger_count%22)%20%3C%3D%206)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.with_columns((pl.col(%22fare_amount%22)%20%2F%20pl.col(%22passenger_count%22)).alias(%22fare_per_pax%22))%0A%20%20%20%20%20%20%20%20.group_by(%22passenger_count%22)%0A%20%20%20%20%20%20%20%20.agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.len().alias(%22trips%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22fare_amount%22).mean().alias(%22avg_fare%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20pl.col(%22fare_per_pax%22).mean().alias(%22avg_fare_per_pax%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.sort(%22passenger_count%22)%0A%20%20%20%20)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Fare%20by%20Passenger%20Count%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_pax_fare%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.callout(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22Taxi%20fares%20don't%20scale%20with%20passenger%20count%20-%20unlike%20rideshare%20pooling%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20kind%3D%22info%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Summary%20Dashboard%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(%0A%20%20%20%20MANHATTAN_CBD_ZONES%2C%0A%20%20%20%20format_currency%2C%0A%20%20%20%20format_large_number%2C%0A%20%20%20%20green_df%2C%0A%20%20%20%20jan_2024_df%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20pl%2C%0A%20%20%20%20yellow_df%2C%0A)%3A%0A%20%20%20%20_cbd_zones%20%3D%20set(MANHATTAN_CBD_ZONES)%0A%0A%20%20%20%20%23%202025%20stats%20(post-congestion)%0A%20%20%20%20_y25%20%3D%20yellow_df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20_y25_trips%20%3D%20len(_y25)%0A%20%20%20%20_y25_fare%20%3D%20_y25%5B%22total_amount%22%5D.mean()%0A%20%20%20%20_y25_dist%20%3D%20_y25%5B%22trip_distance%22%5D.mean()%0A%20%20%20%20_y25_cbd_pct%20%3D%20_y25.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200).height%20%2F%20_y25_trips%20*%20100%0A%20%20%20%20_y25_cbd_fee%20%3D%20(%0A%20%20%20%20%20%20%20%20_y25.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200)%5B%22cbd_congestion_fee%22%5D.mean()%0A%20%20%20%20%20%20%20%20if%20_y25.filter(pl.col(%22cbd_congestion_fee%22)%20%3E%200).height%20%3E%200%0A%20%20%20%20%20%20%20%20else%200%0A%20%20%20%20)%0A%0A%20%20%20%20%23%202024%20stats%20(same%20period)%0A%20%20%20%20_y24%20%3D%20jan_2024_df.filter(pl.col(%22day%22)%20%3E%3D%205)%0A%20%20%20%20_y24_trips%20%3D%20len(_y24)%0A%20%20%20%20_y24_fare%20%3D%20_y24%5B%22total_amount%22%5D.mean()%0A%0A%20%20%20%20%23%20Calculate%20changes%0A%20%20%20%20_trip_change%20%3D%20(_y25_trips%20-%20_y24_trips)%20%2F%20_y24_trips%20*%20100%0A%20%20%20%20_fare_change%20%3D%20_y25_fare%20-%20_y24_fare%0A%0A%20%20%20%20%23%20Green%20stats%0A%20%20%20%20_g25%20%3D%20green_df.filter(pl.col(%22post_congestion%22))%0A%20%20%20%20_g_trips%20%3D%20len(_g25)%0A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%22%23%23%23%20Key%20Metrics%20(Jan%205-31)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Yellow%20Trips%202025%22%2C%20value%3Dformat_large_number(_y25_trips))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22vs%202024%22%2C%20value%3Df%22%7B_trip_change%3A%2B.1f%7D%25%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Avg%20Fare%202025%22%2C%20value%3Dformat_currency(_y25_fare))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Fare%20Change%22%2C%20value%3Df%22%7B_fare_change%3A%2B.2f%7D%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22%25%20with%20CBD%20Fee%22%2C%20value%3Df%22%7B_y25_cbd_pct%3A.1f%7D%25%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20mo.stat(label%3D%22Avg%20CBD%20Fee%22%2C%20value%3Dformat_currency(_y25_cbd_fee))%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20justify%3D%22space-around%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(cbd_speed_yoy%2C%20mo%2C%20pl%2C%20yoy_daily)%3A%0A%20%20%20%20_pre_trip_change%20%3D%20yoy_daily.filter(pl.col(%22day%22)%20%3C%205)%5B%22pct_change%22%5D.mean()%0A%20%20%20%20_post_trip_change%20%3D%20yoy_daily.filter(pl.col(%22day%22)%20%3E%3D%205)%5B%22pct_change%22%5D.mean()%0A%0A%20%20%20%20_speed_change%20%3D%20cbd_speed_yoy.filter(pl.col(%22metric%22)%20%3D%3D%20%22Mean%20Speed%20(mph)%22)%5B%22change%22%5D%5B0%5D%0A%20%20%20%20_rush_speed_change%20%3D%20cbd_speed_yoy.filter(pl.col(%22metric%22)%20%3D%3D%20%22Rush%20Hour%20Mean%20(7-9%2C%2017-19)%22)%5B%0A%20%20%20%20%20%20%20%20%22change%22%0A%20%20%20%20%5D%5B0%5D%0A%0A%20%20%20%20mo.md(f%22%22%22%0A%20%20%20%20%23%23%20Key%20Findings%3A%20Congestion%20Pricing%20Impact%0A%0A%20%20%20%20%23%23%23%20Trip%20Volume%20(Year-over-Year)%0A%20%20%20%20%7C%20Period%20%7C%20Avg%20Change%20vs%202024%20%7C%0A%20%20%20%20%7C--------%7C-------------------%7C%0A%20%20%20%20%7C%20Pre-congestion%20(Jan%201-4)%20%7C%20%7B_pre_trip_change%3A%2B.1f%7D%25%20%7C%0A%20%20%20%20%7C%20Post-congestion%20(Jan%205-31)%20%7C%20%7B_post_trip_change%3A%2B.1f%7D%25%20%7C%0A%20%20%20%20%7C%20**Difference**%20%7C%20**%7B_post_trip_change%20-%20_pre_trip_change%3A%2B.1f%7D%20pp**%20%7C%0A%0A%20%20%20%20%23%23%23%20CBD%20Zone%20Speed%0A%20%20%20%20%7C%20Metric%20%7C%20Change%20vs%202024%20%7C%0A%20%20%20%20%7C--------%7C---------------%7C%0A%20%20%20%20%7C%20Overall%20CBD%20speed%20%7C%20%7B_speed_change%3A%2B.1f%7D%20mph%20%7C%0A%20%20%20%20%7C%20Rush%20hour%20CBD%20speed%20%7C%20%7B_rush_speed_change%3A%2B.1f%7D%20mph%20%7C%0A%0A%20%20%20%20%23%23%23%20Interpretation%0A%20%20%20%20%7B%22**Trip%20volume%20increased**%20despite%20congestion%20pricing%20-%20taxis%20remain%20competitive%22%20if%20_post_trip_change%20%3E%200%20else%20%22**Trip%20volume%20decreased**%20-%20possible%20demand%20reduction%20from%20pricing%22%7D%0A%0A%20%20%20%20%7B%22**Speeds%20improved**%20-%20congestion%20pricing%20may%20be%20reducing%20traffic%22%20if%20_speed_change%20%3E%200%20else%20%22**Speeds%20unchanged%2Fdecreased**%20-%20early%20days%2C%20more%20data%20needed%22%7D%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%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