rolling threshold marimo

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.8%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%0A%0A%0A%40app.cell%0Adef%20imports()%3A%0A%20%20%20%20import%20calendar%0A%0A%20%20%20%20import%20jinja2%20%20%23%20noqa%3A%20F401%20%E2%80%94%20required%20by%20pandas.style%20in%20Pyodide%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20pandas%20as%20pd%0A%0A%20%20%20%20COLS%20%3D%20%5B%22Jan%22%2C%20%22Feb%22%2C%20%22Mar%22%2C%20%22Apr%22%2C%20%22May%22%2C%20%22Jun%22%5D%0A%20%20%20%20return%20COLS%2C%20calendar%2C%20mo%2C%20np%2C%20pd%2C%20plt%0A%0A%0A%40app.cell%0Adef%20load_data(mo%2C%20pd)%3A%0A%20%20%20%20import%20importlib%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20with%20mo.status.spinner(subtitle%3D%22Loading%20data...%22)%3A%0A%20%20%20%20%20%20%20%20if%20sys.platform%20%3D%3D%20%22emscripten%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20import%20io%0A%20%20%20%20%20%20%20%20%20%20%20%20import%20urllib.request%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20_url%20%3D%20str(mo.notebook_location()%20%2F%20%22public%22%20%2F%20%22iri_data.csv%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20with%20urllib.request.urlopen(_url)%20as%20_resp%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20df_iri%20%3D%20pd.read_csv(io.StringIO(_resp.read().decode(%22utf-8%22)))%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_stratus%20%3D%20importlib.import_module(%22ocha_stratus%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20blob_name%20%3D%20%22ds-aa-ner-drought%2Fraw%2Firi%2Fner_maproom_export_2026-04-25_thresh35%20-%20Sheet1.csv%22%0A%20%20%20%20%20%20%20%20%20%20%20%20df_iri%20%3D%20_stratus.load_csv_from_blob(blob_name)%0A%0A%20%20%20%20df_iri.columns%20%3D%20%5B%0A%20%20%20%20%20%20%20%20%22year%22%2C%0A%20%20%20%20%20%20%20%20%22Jan%22%2C%0A%20%20%20%20%20%20%20%20%22Feb%22%2C%0A%20%20%20%20%20%20%20%20%22Mar%22%2C%0A%20%20%20%20%20%20%20%20%22Apr%22%2C%0A%20%20%20%20%20%20%20%20%22May%22%2C%0A%20%20%20%20%20%20%20%20%22Jun%22%2C%0A%20%20%20%20%20%20%20%20%22Aug%22%2C%0A%20%20%20%20%20%20%20%20%22JAS_SPI%22%2C%0A%20%20%20%20%5D%0A%20%20%20%20df_iri%5B%22year%22%5D%20%3D%20df_iri%5B%22year%22%5D.astype(int)%0A%20%20%20%20_month_cols%20%3D%20%5B%0A%20%20%20%20%20%20%20%20c%20for%20c%20in%20df_iri.columns%20if%20c%20not%20in%20(%22year%22%2C%20%22Aug%22%2C%20%22JAS_SPI%22)%0A%20%20%20%20%5D%0A%20%20%20%20df_iri%5B_month_cols%5D%20%3D%20df_iri%5B_month_cols%5D%20%2F%20100.0%0A%20%20%20%20df_iri%20%3D%20df_iri.sort_values(%22year%22).reset_index(drop%3DTrue)%0A%20%20%20%20return%20(df_iri%2C)%0A%0A%0A%40app.cell%0Adef%20params()%3A%0A%20%20%20%20ref_window%20%3D%2010%0A%20%20%20%20start_eval_year%20%3D%20(%0A%20%20%20%20%20%20%20%202001%20%20%23%20first%20year%20with%20full%2010-yr%20reference%20window%20(data%20starts%201991)%0A%20%20%20%20)%0A%20%20%20%20end_eval_year%20%3D%202025%0A%20%20%20%20pct_steps%20%3D%20list(%0A%20%20%20%20%20%20%20%20range(0%2C%20105%2C%205)%0A%20%20%20%20)%20%20%23%20fraction%20triggering%20from%20top%3A%200%25%2C%205%25%2C%20...%2C%20100%25%0A%20%20%20%20mos%20%3D%20%5B1%2C%202%2C%203%2C%204%2C%205%2C%206%5D%0A%20%20%20%20rp_target%20%3D%203.5%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20end_eval_year%2C%0A%20%20%20%20%20%20%20%20mos%2C%0A%20%20%20%20%20%20%20%20pct_steps%2C%0A%20%20%20%20%20%20%20%20ref_window%2C%0A%20%20%20%20%20%20%20%20rp_target%2C%0A%20%20%20%20%20%20%20%20start_eval_year%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell%0Adef%20note_data(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%3E%20**Data%20source%3A**%20Historical%20IRI%20forecasts%20exported%20from%20the%20Maproom%20on%20**25%20April%202026**%2C%0A%20%20%20%20%3E%20using%20the%20**OCHA%20Certification**%20model%20with%20the%20Frequency%20slider%20set%20to%20**35%25**.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_plots(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Forecast%20and%20observational%20thresholds%0A%0A%20%20%20%20Use%20the%20sliders%20to%20set%20the%20threshold%20for%20each%20trigger%20component.%20The%20**forecast%20component**%0A%20%20%20%20triggers%20if%20any%20two%20consecutive%20months%20(Jan%2BFeb%2C%20Feb%2BMar%2C%20%E2%80%A6%2C%20May%2BJun)%20both%20fall%20in%20the%0A%20%20%20%20top%20X%25%20of%20their%20rolling%2010-year%20historical%20reference%20window.%20The%20**observational%20component**%0A%20%20%20%20triggers%20if%20the%20ENACTS%20SPI%20falls%20in%20the%20bottom%20Y%25%20of%20the%20full%20historical%20record.%0A%0A%20%20%20%20The%202%C3%973%20grid%20shows%20all%20six%20forecast%20months.%20The%20observational%20chart%20shows%20ENACTS%20SPI%0A%20%20%20%20against%20the%20fixed%20historical%20threshold.%20Red%20markers%20%3D%20triggered.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20obs_ui(mo)%3A%0A%20%20%20%20obs_pct_slider%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D5%2C%0A%20%20%20%20%20%20%20%20stop%3D50%2C%0A%20%20%20%20%20%20%20%20step%3D1%2C%0A%20%20%20%20%20%20%20%20value%3D15%2C%0A%20%20%20%20%20%20%20%20label%3D%22Observational%20component%3A%20Aug%20percentile%20threshold%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(obs_pct_slider%2C)%0A%0A%0A%40app.cell%0Adef%20pct_ui(mo)%3A%0A%20%20%20%20pct_sel%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20start%3D0%2C%0A%20%20%20%20%20%20%20%20stop%3D100%2C%0A%20%20%20%20%20%20%20%20step%3D5%2C%0A%20%20%20%20%20%20%20%20value%3D35%2C%0A%20%20%20%20%20%20%20%20label%3D%22Forecast%20component%3A%20%25%20triggering%20from%20top%22%2C%0A%20%20%20%20%20%20%20%20show_value%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20return%20(pct_sel%2C)%0A%0A%0A%40app.cell%0Adef%20compute_obs_triggers(%0A%20%20%20%20df_iri%2C%0A%20%20%20%20end_eval_year%2C%0A%20%20%20%20np%2C%0A%20%20%20%20obs_pct_slider%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_eval_years%20%3D%20list(range(start_eval_year%2C%20end_eval_year%20%2B%201))%0A%20%20%20%20_all_aug%20%3D%20df_iri%5B%22Aug%22%5D.values%0A%20%20%20%20_n_hist%20%3D%20len(_all_aug)%0A%20%20%20%20if%20_obs_pct%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20_obs_thresh%20%3D%20float(%22-inf%22)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20_k_obs%20%3D%20int(np.ceil(_obs_pct%20%2F%20100%20*%20_n_hist))%0A%20%20%20%20%20%20%20%20_obs_thresh%20%3D%20float(sorted(_all_aug)%5B_k_obs%20-%201%5D)%0A%20%20%20%20_obs_rows%20%3D%20%5B%5D%0A%20%20%20%20for%20_year%20in%20_eval_years%3A%0A%20%20%20%20%20%20%20%20_actual%20%3D%20df_iri%5Bdf_iri%5B%22year%22%5D%20%3D%3D%20_year%5D.iloc%5B0%5D%0A%20%20%20%20%20%20%20%20_act_aug%20%3D%20float(_actual%5B%22Aug%22%5D)%0A%20%20%20%20%20%20%20%20_trig_obsv%20%3D%20_act_aug%20%3C%3D%20_obs_thresh%0A%20%20%20%20%20%20%20%20_obs_rows.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22year%22%3A%20_year%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22obs_threshold%22%3A%20_obs_thresh%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22actual_aug%22%3A%20_act_aug%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trig_obsv%22%3A%20_trig_obsv%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20df_obs%20%3D%20pd.DataFrame(_obs_rows).set_index(%22year%22)%0A%20%20%20%20return%20(df_obs%2C)%0A%0A%0A%40app.cell%0Adef%20compute_forecast_triggers(%0A%20%20%20%20COLS%2C%0A%20%20%20%20calendar%2C%0A%20%20%20%20df_iri%2C%0A%20%20%20%20end_eval_year%2C%0A%20%20%20%20mos%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pct_steps%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20ref_window%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20%22%22%22Forecast%20component%20only%20%E2%80%94%20does%20not%20depend%20on%20obs_pct%20so%20slider%20won't%20rerun%20this.%22%22%22%0A%20%20%20%20_eval_years%20%3D%20list(range(start_eval_year%2C%20end_eval_year%20%2B%201))%0A%20%20%20%20_consec_pairs%20%3D%20%5B(mos%5Bi%5D%2C%20mos%5Bi%20%2B%201%5D)%20for%20i%20in%20range(len(mos)%20-%201)%5D%0A%0A%20%20%20%20_result_rows%20%3D%20%5B%5D%0A%20%20%20%20_thresh_rows%20%3D%20%5B%5D%0A%0A%20%20%20%20for%20_pct%20in%20pct_steps%3A%0A%20%20%20%20%20%20%20%20for%20_year%20in%20_eval_years%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20_ref%20%3D%20df_iri%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20df_iri%5B%22year%22%5D.between(_year%20-%20ref_window%2C%20_year%20-%201)%0A%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_actual%20%3D%20df_iri%5Bdf_iri%5B%22year%22%5D%20%3D%3D%20_year%5D.iloc%5B0%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20_n_ref%20%3D%20len(_ref)%0A%20%20%20%20%20%20%20%20%20%20%20%20_k%20%3D%20int(np.ceil(_pct%20%2F%20100%20*%20_n_ref))%20if%20_pct%20%3E%200%20else%200%0A%20%20%20%20%20%20%20%20%20%20%20%20_trig_month%20%3D%20%7B%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20_col%20in%20COLS%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20k-th%20highest%20value%20in%20reference%20window%3B%20locks%20threshold%20to%20an%20actual%20historical%20value%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20_k%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_thresh%20%3D%20float(%22inf%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_thresh%20%3D%20float(sorted(_ref%5B_col%5D.values%2C%20reverse%3DTrue)%5B_k%20-%201%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_act_val%20%3D%20float(_actual%5B_col%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_trig%20%3D%20_act_val%20%3E%3D%20_thresh%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_trig_month%5B_col%5D%20%3D%20_trig%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_thresh_rows.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct%22%3A%20_pct%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22year%22%3A%20_year%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22month%22%3A%20_col%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22threshold%22%3A%20_thresh%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22actual%22%3A%20_act_val%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22triggered%22%3A%20_trig%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20_pair_trigs%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22trig_%7Bcalendar.month_abbr%5B_m1%5D%7D_%7Bcalendar.month_abbr%5B_m2%5D%7D%22%3A%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_trig_month%5Bcalendar.month_abbr%5B_m1%5D%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20and%20_trig_month%5Bcalendar.month_abbr%5B_m2%5D%5D%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%20for%20_m1%2C%20_m2%20in%20_consec_pairs%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20_trig_fcast%20%3D%20any(_pair_trigs.values())%0A%20%20%20%20%20%20%20%20%20%20%20%20_result_rows.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct%22%3A%20_pct%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22year%22%3A%20_year%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22trig_fcast%22%3A%20_trig_fcast%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20**%7Bf%22trig_%7B_col%7D%22%3A%20_trig_month%5B_col%5D%20for%20_col%20in%20COLS%7D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20**_pair_trigs%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20df_forecast%20%3D%20pd.DataFrame(_result_rows)%0A%20%20%20%20df_thresholds%20%3D%20pd.DataFrame(_thresh_rows)%0A%20%20%20%20return%20df_forecast%2C%20df_thresholds%0A%0A%0A%40app.cell%0Adef%20combine_results(df_forecast%2C%20df_obs%2C%20pct_steps%2C%20pd)%3A%0A%20%20%20%20_obs%20%3D%20df_obs%5B%5B%22trig_obsv%22%5D%5D%0A%20%20%20%20_chunks%20%3D%20%5B%5D%0A%20%20%20%20for%20_pct%20in%20pct_steps%3A%0A%20%20%20%20%20%20%20%20_fc%20%3D%20df_forecast%5Bdf_forecast%5B%22pct%22%5D%20%3D%3D%20_pct%5D.set_index(%22year%22)%0A%20%20%20%20%20%20%20%20_merged%20%3D%20_fc.join(_obs)%0A%20%20%20%20%20%20%20%20_merged%5B%22trig_either%22%5D%20%3D%20_merged%5B%22trig_fcast%22%5D%20%7C%20_merged%5B%22trig_obsv%22%5D%0A%20%20%20%20%20%20%20%20_chunks.append(_merged.reset_index())%0A%20%20%20%20df_results%20%3D%20pd.concat(_chunks%2C%20ignore_index%3DTrue)%0A%20%20%20%20return%20(df_results%2C)%0A%0A%0A%40app.cell%0Adef%20compute_summary(%0A%20%20%20%20df_results%2C%0A%20%20%20%20end_eval_year%2C%0A%20%20%20%20obs_pct_slider%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_n_years%20%3D%20end_eval_year%20-%20start_eval_year%20%2B%201%0A%20%20%20%20_rows%20%3D%20%5B%5D%0A%20%20%20%20for%20_pct%2C%20_grp%20in%20df_results.groupby(%22pct%22)%3A%0A%20%20%20%20%20%20%20%20_n_fcast%20%3D%20int(_grp%5B%22trig_fcast%22%5D.sum())%0A%20%20%20%20%20%20%20%20_n_obsv%20%3D%20int(_grp%5B%22trig_obsv%22%5D.sum())%0A%20%20%20%20%20%20%20%20_n_either%20%3D%20int(_grp%5B%22trig_either%22%5D.sum())%0A%20%20%20%20%20%20%20%20_rp_fcast%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20round((_n_years%20%2B%201)%20%2F%20_n_fcast%2C%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_n_fcast%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20else%20float(%22inf%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_rp_obsv%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20round((_n_years%20%2B%201)%20%2F%20_n_obsv%2C%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_n_obsv%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20else%20float(%22inf%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_rp_either%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20round((_n_years%20%2B%201)%20%2F%20_n_either%2C%201)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20_n_either%20%3E%200%0A%20%20%20%20%20%20%20%20%20%20%20%20else%20float(%22inf%22)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_rows.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22pct_triggering%22%3A%20_pct%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22n_fcast%22%3A%20_n_fcast%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22rp_fcast%22%3A%20_rp_fcast%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22n_obsv%20(SPI%E2%89%A4%7B_obs_pct%7D%25)%22%3A%20_n_obsv%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22rp_obsv%20(SPI%E2%89%A4%7B_obs_pct%7D%25)%22%3A%20_rp_obsv%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22n_either%22%3A%20_n_either%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22rp_either%22%3A%20_rp_either%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20df_summary%20%3D%20pd.DataFrame(_rows)%0A%20%20%20%20return%20(df_summary%2C)%0A%0A%0A%40app.cell%0Adef%20find_closest_pct(df_summary%2C%20rp_target)%3A%0A%20%20%20%20_finite%20%3D%20df_summary%5Bdf_summary%5B%22rp_either%22%5D%20%3C%20float(%22inf%22)%5D%0A%20%20%20%20closest_pct%20%3D%20int(%0A%20%20%20%20%20%20%20%20_finite.iloc%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(_finite%5B%22rp_either%22%5D%20-%20rp_target).abs().argsort().iloc%5B0%5D%0A%20%20%20%20%20%20%20%20%5D%5B%22pct_triggering%22%5D%0A%20%20%20%20)%0A%20%20%20%20return%20(closest_pct%2C)%0A%0A%0A%40app.cell%0Adef%20bad_years_data(pd)%3A%0A%20%20%20%20df_bad_years%20%3D%20pd.DataFrame(%0A%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22year%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%202021%2C%202020%2C%202019%2C%202018%2C%202017%2C%202016%2C%202015%2C%202014%2C%202013%2C%202012%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%202011%2C%202010%2C%202009%2C%202008%2C%202007%2C%202006%2C%202005%2C%202004%2C%202003%2C%202002%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%202001%2C%202000%2C%201999%2C%201998%2C%201997%2C%201996%2C%201995%2C%201994%2C%201993%2C%201992%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%201991%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%22bad_year_rank%22%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%209%2C%202%2C%2032%2C%2032%2C%207%2C%2032%2C%205%2C%2012%2C%2017%2C%2016%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%204%2C%2014%2C%201%2C%2013%2C%2032%2C%2032%2C%2011%2C%2010%2C%2015%2C%2032%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%203%2C%2032%2C%2032%2C%2032%2C%208%2C%2032%2C%2032%2C%2032%2C%2032%2C%2032%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%2032%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)%0A%20%20%20%20return%20(df_bad_years%2C)%0A%0A%0A%40app.cell%0Adef%20obs_display(df_obs%2C%20mo%2C%20obs_pct_slider%2C%20pct_sel)%3A%0A%20%20%20%20_n%20%3D%20int(df_obs%5B%22trig_obsv%22%5D.sum())%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**Forecast%20component%3A**%20%7Bpct_sel%7D%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**Observational%20component%20(ENACTS%20SPI)%3A**%20%7Bobs_pct_slider%7D%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%E2%86%92%20**%7B_n%7D%20trigger%20year%7B's'%20if%20_n%20!%3D%201%20else%20''%7D**%22%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%20rp_readout(df_summary%2C%20mo%2C%20pct_sel)%3A%0A%20%20%20%20_pct%20%3D%20pct_sel.value%0A%20%20%20%20_row%20%3D%20df_summary%5Bdf_summary%5B%22pct_triggering%22%5D%20%3D%3D%20_pct%5D.iloc%5B0%5D%0A%20%20%20%20_rp_fcast%20%3D%20_row%5B%22rp_fcast%22%5D%0A%20%20%20%20_rp_obsv_col%20%3D%20next(c%20for%20c%20in%20df_summary.columns%20if%20c.startswith(%22rp_obsv%22))%0A%20%20%20%20_rp_obsv%20%3D%20_row%5B_rp_obsv_col%5D%0A%20%20%20%20_rp_either%20%3D%20_row%5B%22rp_either%22%5D%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22Forecast%20component%20RP%3A%20**%7B_rp_fcast%7D**%20%C2%B7%20%22%0A%20%20%20%20%20%20%20%20f%22Observational%20component%20RP%3A%20**%7B_rp_obsv%7D**%20%C2%B7%20%22%0A%20%20%20%20%20%20%20%20f%22Combined%20RP%3A%20**%7B_rp_either%7D**%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20trigger_bars(%0A%20%20%20%20df_bad_years%2C%0A%20%20%20%20df_results%2C%0A%20%20%20%20df_summary%2C%0A%20%20%20%20obs_pct_slider%2C%0A%20%20%20%20pct_sel%2C%0A%20%20%20%20plt%2C%0A)%3A%0A%20%20%20%20from%20scipy.stats%20import%20spearmanr%20as%20_spearmanr%0A%0A%20%20%20%20_pct%20%3D%20pct_sel.value%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_grp%20%3D%20df_results%5Bdf_results%5B%22pct%22%5D%20%3D%3D%20_pct%5D%0A%0A%20%20%20%20_n_fcast%20%3D%20int(_grp%5B%22trig_fcast%22%5D.sum())%0A%20%20%20%20_n_obsv%20%3D%20int(_grp%5B%22trig_obsv%22%5D.sum())%0A%20%20%20%20_n_either%20%3D%20int(_grp%5B%22trig_either%22%5D.sum())%0A%0A%20%20%20%20_row%20%3D%20df_summary%5Bdf_summary%5B%22pct_triggering%22%5D%20%3D%3D%20_pct%5D.iloc%5B0%5D%0A%20%20%20%20_rp_fcast%20%3D%20_row%5B%22rp_fcast%22%5D%0A%20%20%20%20_rp_obsv_col%20%3D%20next(c%20for%20c%20in%20df_summary.columns%20if%20c.startswith(%22rp_obsv%22))%0A%20%20%20%20_rp_obsv%20%3D%20_row%5B_rp_obsv_col%5D%0A%20%20%20%20_rp_either%20%3D%20_row%5B%22rp_either%22%5D%0A%0A%20%20%20%20_merged%20%3D%20_grp.set_index(%22year%22).join(%0A%20%20%20%20%20%20%20%20df_bad_years.set_index(%22year%22)%5B%22bad_year_rank%22%5D%2C%20how%3D%22inner%22%0A%20%20%20%20)%0A%20%20%20%20_max_rank%20%3D%20int(_merged%5B%22bad_year_rank%22%5D.max())%0A%20%20%20%20_severity%20%3D%20_max_rank%20%2B%201%20-%20_merged%5B%22bad_year_rank%22%5D%0A%0A%20%20%20%20_labels%20%3D%20%5Bf%22Forecast%20(%7B_pct%7D%25)%22%2C%20f%22Observational%20(%7B_obs_pct%7D%25)%22%2C%20%22Either%22%5D%0A%20%20%20%20_colors%20%3D%20%5B%22steelblue%22%2C%20%22darkorange%22%2C%20%22crimson%22%5D%0A%20%20%20%20_counts%20%3D%20%5B_n_fcast%2C%20_n_obsv%2C%20_n_either%5D%0A%20%20%20%20_trig_cols%20%3D%20%5B%22trig_fcast%22%2C%20%22trig_obsv%22%2C%20%22trig_either%22%5D%0A%0A%20%20%20%20_rs%2C%20_ps%20%3D%20%5B%5D%2C%20%5B%5D%0A%20%20%20%20for%20_col%20in%20_trig_cols%3A%0A%20%20%20%20%20%20%20%20_r%2C%20_p%20%3D%20_spearmanr(_merged%5B_col%5D.astype(int)%2C%20_severity)%0A%20%20%20%20%20%20%20%20_rs.append(float(_r))%0A%20%20%20%20%20%20%20%20_ps.append(float(_p))%0A%0A%20%20%20%20_fig%2C%20(_ax1%2C%20_ax2)%20%3D%20plt.subplots(1%2C%202%2C%20figsize%3D(10%2C%204))%0A%0A%20%20%20%20_b1%20%3D%20_ax1.bar(_labels%2C%20_counts%2C%20color%3D_colors%2C%20alpha%3D0.8)%0A%20%20%20%20_ax1.bar_label(_b1)%0A%20%20%20%20_ax1.set_ylabel(%22Trigger%20years%22)%0A%20%20%20%20_ax1.set_title(%22Trigger%20counts%22)%0A%20%20%20%20_ax1.spines%5B%5B%22top%22%2C%20%22right%22%5D%5D.set_visible(False)%0A%0A%20%20%20%20_b2%20%3D%20_ax2.bar(_labels%2C%20_rs%2C%20color%3D_colors%2C%20alpha%3D0.8)%0A%20%20%20%20for%20_bar%2C%20_r%2C%20_p%20in%20zip(_b2%2C%20_rs%2C%20_ps)%3A%0A%20%20%20%20%20%20%20%20_stars%20%3D%20%22***%22%20if%20_p%20%3C%200.001%20else%20%22**%22%20if%20_p%20%3C%200.01%20else%20%22*%22%20if%20_p%20%3C%200.05%20else%20%22%22%0A%20%20%20%20%20%20%20%20_ax2.text(%0A%20%20%20%20%20%20%20%20%20%20%20%20_bar.get_x()%20%2B%20_bar.get_width()%20%2F%202%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_r%20%2B%20(0.03%20if%20_r%20%3E%3D%200%20else%20-0.03)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%7B_r%3A.2f%7D%7B_stars%7D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22center%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20va%3D%22bottom%22%20if%20_r%20%3E%3D%200%20else%20%22top%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20fontsize%3D9%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_ax2.axhline(0%2C%20color%3D%22black%22%2C%20lw%3D0.8)%0A%20%20%20%20_ax2.set_ylim(-1%2C%201)%0A%20%20%20%20_ax2.set_ylabel(%22Spearman%20r%22)%0A%20%20%20%20_ax2.set_title(%22Bad%20year%20correlation%20(%E2%86%91%20%3D%20triggers%20in%20worse%20years)%22)%0A%20%20%20%20_ax2.spines%5B%5B%22top%22%2C%20%22right%22%5D%5D.set_visible(False)%0A%0A%20%20%20%20_rp_line%20%3D%20(%0A%20%20%20%20%20%20%20%20f%22Forecast%20RP%3A%20%7B_rp_fcast%7D%20%20%C2%B7%20%20%22%0A%20%20%20%20%20%20%20%20f%22Observational%20RP%3A%20%7B_rp_obsv%7D%20%20%C2%B7%20%20%22%0A%20%20%20%20%20%20%20%20f%22Combined%20RP%3A%20%7B_rp_either%7D%22%0A%20%20%20%20)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20plt.subplots_adjust(bottom%3D0.18)%0A%20%20%20%20_fig.text(0.5%2C%200.03%2C%20_rp_line%2C%20ha%3D%22center%22%2C%20va%3D%22bottom%22%2C%20fontsize%3D9%2C%20color%3D%22%23555555%22)%0A%20%20%20%20_fig%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20all_months_plot(%0A%20%20%20%20COLS%2C%0A%20%20%20%20df_iri%2C%0A%20%20%20%20df_thresholds%2C%0A%20%20%20%20pct_sel%2C%0A%20%20%20%20plt%2C%0A%20%20%20%20ref_window%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20_pct%20%3D%20pct_sel.value%0A%20%20%20%20_fig%2C%20_axes%20%3D%20plt.subplots(2%2C%203%2C%20figsize%3D(14%2C%207)%2C%20sharey%3DFalse)%0A%20%20%20%20_pre%20%3D%20df_iri%5Bdf_iri%5B%22year%22%5D%20%3C%20start_eval_year%5D.sort_values(%22year%22)%0A%0A%20%20%20%20for%20_col%2C%20_ax%20in%20zip(COLS%2C%20_axes.flat)%3A%0A%20%20%20%20%20%20%20%20_df%20%3D%20df_thresholds%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(df_thresholds%5B%22month%22%5D%20%3D%3D%20_col)%20%26%20(df_thresholds%5B%22pct%22%5D%20%3D%3D%20_pct)%0A%20%20%20%20%20%20%20%20%5D.sort_values(%22year%22)%0A%20%20%20%20%20%20%20%20%23%20Pre-evaluation%20actuals%20(reference%20period)%0A%20%20%20%20%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20%20%20%20%20_pre%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_pre%5B_col%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22lightgray%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D2%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20s%3D28%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%23%20Threshold%20line%20from%20start_eval_year%0A%20%20%20%20%20%20%20%20_ax.plot(_df%5B%22year%22%5D%2C%20_df%5B%22threshold%22%5D%2C%20lw%3D1.5%2C%20color%3D%22steelblue%22)%0A%20%20%20%20%20%20%20%20_trig%20%3D%20_df%5B_df%5B%22triggered%22%5D%5D%0A%20%20%20%20%20%20%20%20_no_trig%20%3D%20_df%5B~_df%5B%22triggered%22%5D%5D%0A%20%20%20%20%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20%20%20%20%20_trig%5B%22year%22%5D%2C%20_trig%5B%22actual%22%5D%2C%20color%3D%22crimson%22%2C%20zorder%3D5%2C%20s%3D55%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20%20%20%20%20_no_trig%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20_no_trig%5B%22actual%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3D%22gray%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20zorder%3D4%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20s%3D35%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alpha%3D0.6%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ax.set_title(_col)%0A%20%20%20%20%20%20%20%20_ax.spines%5B%5B%22top%22%2C%20%22right%22%5D%5D.set_visible(False)%0A%20%20%20%20%20%20%20%20_ax.tick_params(labelsize%3D8)%0A%0A%20%20%20%20_fig.suptitle(%0A%20%20%20%20%20%20%20%20f%22Rolling%20%7Bref_window%7D-yr%20threshold%20at%20top%20%7B_pct%7D%25%20%E2%80%94%20red%20%3D%20triggered%22%2C%0A%20%20%20%20%20%20%20%20fontsize%3D12%2C%0A%20%20%20%20)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20_fig%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20aug_obs_plot(df_iri%2C%20df_obs%2C%20obs_pct_slider%2C%20plt%2C%20start_eval_year)%3A%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_obs_thresh%20%3D%20df_obs%5B%22obs_threshold%22%5D.iloc%5B0%5D%0A%20%20%20%20_all%20%3D%20df_iri.sort_values(%22year%22)%0A%20%20%20%20_pre%20%3D%20_all%5B_all%5B%22year%22%5D%20%3C%20start_eval_year%5D%0A%20%20%20%20_eval%20%3D%20_all%5B_all%5B%22year%22%5D%20%3E%3D%20start_eval_year%5D%0A%20%20%20%20_trig_eval%20%3D%20_eval%5B_eval%5B%22Aug%22%5D%20%3C%3D%20_obs_thresh%5D%0A%20%20%20%20_no_trig_eval%20%3D%20_eval%5B_eval%5B%22Aug%22%5D%20%3E%20_obs_thresh%5D%0A%0A%20%20%20%20_fig%2C%20_ax%20%3D%20plt.subplots(figsize%3D(10%2C%203.5))%0A%20%20%20%20_ax.axhline(%0A%20%20%20%20%20%20%20%20_obs_thresh%2C%0A%20%20%20%20%20%20%20%20lw%3D1.8%2C%0A%20%20%20%20%20%20%20%20color%3D%22darkorange%22%2C%0A%20%20%20%20%20%20%20%20label%3Df%22ENACTS%20SPI%20threshold%20(bottom%20%7B_obs_pct%7D%25%20of%20full%20record)%22%2C%0A%20%20%20%20)%0A%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20_pre%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20_pre%5B%22Aug%22%5D%2C%0A%20%20%20%20%20%20%20%20color%3D%22lightgray%22%2C%0A%20%20%20%20%20%20%20%20zorder%3D2%2C%0A%20%20%20%20%20%20%20%20s%3D40%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20%20%20%20%20label%3Df%22Pre-%7Bstart_eval_year%7D%20(reference%20only)%22%2C%0A%20%20%20%20)%0A%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20_trig_eval%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20_trig_eval%5B%22Aug%22%5D%2C%0A%20%20%20%20%20%20%20%20color%3D%22crimson%22%2C%0A%20%20%20%20%20%20%20%20zorder%3D5%2C%0A%20%20%20%20%20%20%20%20s%3D70%2C%0A%20%20%20%20%20%20%20%20label%3D%22Observational%20component%20triggered%22%2C%0A%20%20%20%20)%0A%20%20%20%20_ax.scatter(%0A%20%20%20%20%20%20%20%20_no_trig_eval%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20_no_trig_eval%5B%22Aug%22%5D%2C%0A%20%20%20%20%20%20%20%20color%3D%22gray%22%2C%0A%20%20%20%20%20%20%20%20zorder%3D4%2C%0A%20%20%20%20%20%20%20%20s%3D45%2C%0A%20%20%20%20%20%20%20%20alpha%3D0.7%2C%0A%20%20%20%20%20%20%20%20label%3D%22Not%20triggered%22%2C%0A%20%20%20%20)%0A%20%20%20%20_ax.set_title(%0A%20%20%20%20%20%20%20%20f%22ENACTS%20SPI%20(observational%20component)%20%E2%80%94%20full%20historical%20threshold%20at%20bottom%20%7B_obs_pct%7D%25%22%0A%20%20%20%20)%0A%20%20%20%20_ax.set_xlabel(%22Year%22)%0A%20%20%20%20_ax.set_ylabel(%22Aug%20value%22)%0A%20%20%20%20_ax.legend()%0A%20%20%20%20_ax.spines%5B%5B%22top%22%2C%20%22right%22%5D%5D.set_visible(False)%0A%20%20%20%20plt.tight_layout()%0A%20%20%20%20_fig%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_activation(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%20Activation%20record%0A%0A%20%20%20%20The%20table%20below%20shows%20the%20full%20year-by-year%20trigger%20record%20at%20the%20selected%20thresholds%2C%0A%20%20%20%20with%20return%20periods%20for%20each%20component.%20Each%20column%20shows%20whether%20that%20month%2C%0A%20%20%20%20consecutive-month%20pair%2C%20or%20component%20triggered.%20Red%20cells%20%3D%20triggered.%20Hover%20rows%20to%20highlight.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20trigger_detail_table(%0A%20%20%20%20COLS%2C%0A%20%20%20%20calendar%2C%0A%20%20%20%20df_bad_years%2C%0A%20%20%20%20df_results%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20mos%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pct_sel%2C%0A%20%20%20%20pd%2C%0A)%3A%0A%20%20%20%20_pct%20%3D%20pct_sel.value%0A%20%20%20%20_df%20%3D%20(%0A%20%20%20%20%20%20%20%20df_results%5Bdf_results%5B%22pct%22%5D%20%3D%3D%20_pct%5D%0A%20%20%20%20%20%20%20%20.copy()%0A%20%20%20%20%20%20%20%20.sort_values(%22year%22%2C%20ascending%3DFalse)%0A%20%20%20%20)%0A%20%20%20%20_df%20%3D%20_df.join(%0A%20%20%20%20%20%20%20%20df_bad_years.set_index(%22year%22)%5B%22bad_year_rank%22%5D%2C%20on%3D%22year%22%2C%20how%3D%22left%22%0A%20%20%20%20)%0A%0A%20%20%20%20_month_cols%20%3D%20%5Bf%22trig_%7Bc%7D%22%20for%20c%20in%20COLS%5D%0A%20%20%20%20_pair_cols%20%3D%20%5B%0A%20%20%20%20%20%20%20%20f%22trig_%7Bcalendar.month_abbr%5Bmos%5Bi%5D%5D%7D_%7Bcalendar.month_abbr%5Bmos%5Bi%2B1%5D%5D%7D%22%0A%20%20%20%20%20%20%20%20for%20i%20in%20range(len(mos)%20-%201)%0A%20%20%20%20%5D%0A%20%20%20%20_rename%20%3D%20%7B%0A%20%20%20%20%20%20%20%20**%7Bf%22trig_%7Bc%7D%22%3A%20c%20for%20c%20in%20COLS%7D%2C%0A%20%20%20%20%20%20%20%20**%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22trig_%7Bcalendar.month_abbr%5Bmos%5Bi%5D%5D%7D_%7Bcalendar.month_abbr%5Bmos%5Bi%2B1%5D%5D%7D%22%3A%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bcalendar.month_abbr%5Bmos%5Bi%5D%5D%7D%2B%7Bcalendar.month_abbr%5Bmos%5Bi%2B1%5D%5D%7D%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%20for%20i%20in%20range(len(mos)%20-%201)%0A%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20%22trig_fcast%22%3A%20%22Forecast%22%2C%0A%20%20%20%20%20%20%20%20%22trig_obsv%22%3A%20%22ENACTS%20SPI%22%2C%0A%20%20%20%20%20%20%20%20%22trig_either%22%3A%20%22Either%22%2C%0A%20%20%20%20%20%20%20%20%22bad_year_rank%22%3A%20%22Bad%20year%22%2C%0A%20%20%20%20%7D%0A%20%20%20%20_display%20%3D%20(%0A%20%20%20%20%20%20%20%20_df%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%22year%22%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20_month_cols%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20_pair_cols%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%5B%22trig_fcast%22%2C%20%22trig_obsv%22%2C%20%22trig_either%22%2C%20%22bad_year_rank%22%5D%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20.rename(columns%3D_rename)%0A%20%20%20%20%20%20%20%20.reset_index(drop%3DTrue)%0A%20%20%20%20)%0A%0A%20%20%20%20_rank_vals%20%3D%20df_bad_years%5B%22bad_year_rank%22%5D.values%0A%20%20%20%20_p20%20%3D%20float(np.percentile(_rank_vals%2C%2020))%0A%20%20%20%20_p33%20%3D%20float(np.percentile(_rank_vals%2C%2033))%0A%20%20%20%20_p66%20%3D%20float(np.percentile(_rank_vals%2C%2066))%0A%0A%20%20%20%20def%20_rank_style(v)%3A%0A%20%20%20%20%20%20%20%20if%20pd.isna(v)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22background-color%3A%20%23f4f4f4%3B%20color%3A%20%23cccccc%22%0A%20%20%20%20%20%20%20%20elif%20v%20%3C%3D%20_p20%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22background-color%3A%20%23ff9999%22%0A%20%20%20%20%20%20%20%20elif%20v%20%3C%3D%20_p33%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22background-color%3A%20%23ffcc88%22%0A%20%20%20%20%20%20%20%20elif%20v%20%3C%20_p66%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22%22%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20%22background-color%3A%20%23cceecc%22%0A%0A%20%20%20%20_highlight_cols%20%3D%20%5Bc%20for%20c%20in%20_display.columns%20if%20c%20not%20in%20(%22year%22%2C%20%22Bad%20year%22)%5D%0A%20%20%20%20_styled%20%3D%20(%0A%20%20%20%20%20%20%20%20_display.style%0A%20%20%20%20%20%20%20%20.map(%0A%20%20%20%20%20%20%20%20%20%20%20%20lambda%20v%3A%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22background-color%3A%20%23ffaaaa%3B%20color%3A%20%237a0000%3B%20font-weight%3A%20bold%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20v%20is%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%20%22%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%20subset%3D_highlight_cols%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.map(_rank_style%2C%20subset%3D%5B%22Bad%20year%22%5D)%0A%20%20%20%20%20%20%20%20.format(%7B%22Bad%20year%22%3A%20lambda%20v%3A%20%22%22%20if%20pd.isna(v)%20else%20str(int(v))%7D)%0A%20%20%20%20%20%20%20%20.set_uuid(%22trigger_detail%22)%0A%20%20%20%20)%0A%20%20%20%20_css%20%3D%20%22%22%22%3Cstyle%3E%0A%20%20%20%20%23T_trigger_detail%20tbody%20tr%3Ahover%20td%20%7B%0A%20%20%20%20background-color%3A%20%23dde8f8%3B%0A%20%20%20%20cursor%3A%20default%3B%0A%20%20%20%20%7D%0A%20%20%20%20%3C%2Fstyle%3E%22%22%22%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%23%23%23%20Per-year%20trigger%20detail%20(forecast%20component%3A%20top%20%7B_pct%7D%25)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.Html(_css%20%2B%20_styled.hide(axis%3D%22index%22).to_html())%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%20note_analysis(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20---%0A%0A%20%20%20%20%23%23%20Optimization%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_optimization_params(%0A%20%20%20%20end_eval_year%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20ref_window%2C%0A%20%20%20%20rp_target%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22Evaluation%20years%3A%20**%7Bstart_eval_year%7D%E2%80%93%7Bend_eval_year%7D**%20%22%0A%20%20%20%20%20%20%20%20f%22(%7Bend_eval_year%20-%20start_eval_year%20%2B%201%7D%20years%2C%20full%20%7Bref_window%7D-yr%20window%20throughout)%2C%20%22%0A%20%20%20%20%20%20%20%20f%22target%20RP%3A%20**%7Brp_target%7D**%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_obs(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%201%20%C2%B7%20Set%20the%20observational%20trigger%20threshold%0A%0A%20%20%20%20The%20**observational%20component**%20uses%20the%20ENACTS%20MON%20Jun%E2%80%93Jul%20SPI%20exported%20from%20the%0A%20%20%20%20Maproom%20as%20an%20observational%20indicator%20(not%20a%20forecast).%20Set%20the%20slider%20to%20choose%0A%20%20%20%20what%20bottom-percentile%20of%20the%20full%20historical%20record%20counts%20as%20a%20trigger.%0A%20%20%20%20This%20threshold%20is%20fixed%20across%20all%20evaluation%20years%20(not%20rolling).%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_sweep(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%202%20%C2%B7%20Forecast%20threshold%20sweep%20%E2%86%92%20return%20period%20table%0A%0A%20%20%20%20For%20every%20candidate%20percentile%20level%2C%20each%20evaluation%20year%20is%20assessed%3A%20the%0A%20%20%20%20**forecast%20component**%20triggers%20if%20any%20two%20consecutive%20months%20(Jan%2BFeb%2C%20Feb%2BMar%2C%20%E2%80%A6%2C%0A%20%20%20%20May%2BJun)%20both%20exceed%20their%20rolling%2010-year%20historical%20threshold.%20The%20table%0A%20%20%20%20below%20shows%2C%20for%20each%20percentile%2C%20how%20many%20years%20trigger%20under%20each%20component%20and%0A%20%20%20%20the%20implied%20return%20period%20when%20combined%20with%20the%20observational%20component.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20trigger_summary(df_summary%2C%20mo%2C%20obs_pct_slider%2C%20rp_target)%3A%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_near_target%20%3D%20df_summary%5B%0A%20%20%20%20%20%20%20%20df_summary%5B%22rp_either%22%5D.apply(%0A%20%20%20%20%20%20%20%20%20%20%20%20lambda%20x%3A%20isinstance(x%2C%20float)%20and%20abs(x%20-%20rp_target)%20%3C%3D%200.6%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%5D%0A%20%20%20%20_note%20%3D%20(%0A%20%20%20%20%20%20%20%20f%22Rows%20closest%20to%20target%20RP%20%7Brp_target%7D%3A%20pct%20%3D%20%7B_near_target%5B'pct_triggering'%5D.tolist()%7D%22%0A%20%20%20%20%20%20%20%20if%20len(_near_target)%0A%20%20%20%20%20%20%20%20else%20f%22No%20rows%20within%200.6%20of%20target%20RP%20%7Brp_target%7D%22%0A%20%20%20%20)%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(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%23%23%23%20Trigger%20counts%20by%20percentile%20threshold%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22Observational%20component%20fixed%20at%20bottom%20**%7B_obs_pct%7D%25**%20of%20full%20record.%20%7B_note%7D%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%20mo.ui.table(df_summary)%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%20note_auto_select(mo%2C%20rp_target)%3A%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22%22%22%0A%20%20%20%20%23%23%203%20%C2%B7%20Automatic%20threshold%20selection%0A%0A%20%20%20%20The%20forecast%20percentile%20whose%20**combined%20return%20period%20is%20closest%20to%20%7Brp_target%7D%20years**%0A%20%20%20%20is%20identified%20automatically.%20The%20years%20that%20would%20have%20triggered%20under%20each%20component%0A%20%20%20%20are%20listed%20here.%0A%20%20%20%20%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20triggered_years_detail(%0A%20%20%20%20closest_pct%2C%0A%20%20%20%20df_results%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20obs_pct_slider%2C%0A%20%20%20%20rp_target%2C%0A)%3A%0A%20%20%20%20_obs_pct%20%3D%20obs_pct_slider.value%0A%20%20%20%20_grp%20%3D%20df_results%5Bdf_results%5B%22pct%22%5D%20%3D%3D%20closest_pct%5D%0A%20%20%20%20_fcast_years%20%3D%20sorted(_grp%5B_grp%5B%22trig_fcast%22%5D%5D%5B%22year%22%5D.tolist())%0A%20%20%20%20_obsv_years%20%3D%20sorted(_grp%5B_grp%5B%22trig_obsv%22%5D%5D%5B%22year%22%5D.tolist())%0A%20%20%20%20_either_years%20%3D%20sorted(_grp%5B_grp%5B%22trig_either%22%5D%5D%5B%22year%22%5D.tolist())%0A%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20f%22At%20**pct%20%3D%20%7Bclosest_pct%7D%25**%20(closest%20to%20RP%20%7Brp_target%7D)%3A%20%20%5Cn%22%0A%20%20%20%20%20%20%20%20f%22Forecast%20component%3A%20**%7B_fcast_years%7D**%20%20%5Cn%22%0A%20%20%20%20%20%20%20%20f%22Observational%20component%20(ENACTS%20SPI%20%E2%89%A4%7B_obs_pct%7D%25)%3A%20**%7B_obsv_years%7D**%20%20%5Cn%22%0A%20%20%20%20%20%20%20%20f%22Combined%3A%20**%7B_either_years%7D**%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20note_single_month(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20%23%23%23%20Single-month%20detail%0A%0A%20%20%20%20The%20interactive%20chart%20shows%20the%20rolling%2010-year%20threshold%20(blue%20line%2C%20from%202001)%0A%20%20%20%20and%20actual%20IRI%20forecast%20values%20for%20the%20selected%20month.%20Hover%20the%20blue%20threshold%0A%20%20%20%20markers%20to%20see%20the%20reference%20years%20and%20sorted%20values%20behind%20each%20threshold.%0A%0A%20%20%20%20**Month**%20selects%20which%20forecast%20month%20is%20shown.%20**%25%20triggering%20from%20top**%20is%0A%20%20%20%20controlled%20by%20the%20slider%20at%20the%20top%20of%20the%20page.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20month_ui(COLS%2C%20mo)%3A%0A%20%20%20%20month_sel%20%3D%20mo.ui.dropdown(options%3DCOLS%2C%20value%3D%22Jan%22%2C%20label%3D%22Month%22)%0A%20%20%20%20month_sel%0A%20%20%20%20return%20(month_sel%2C)%0A%0A%0A%40app.cell%0Adef%20threshold_evolution_plot(%0A%20%20%20%20df_iri%2C%0A%20%20%20%20df_thresholds%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20month_sel%2C%0A%20%20%20%20pct_sel%2C%0A%20%20%20%20ref_window%2C%0A%20%20%20%20start_eval_year%2C%0A)%3A%0A%20%20%20%20import%20altair%20as%20alt%0A%0A%20%20%20%20_month%20%3D%20month_sel.value%0A%20%20%20%20_pct%20%3D%20pct_sel.value%0A%20%20%20%20_df%20%3D%20(%0A%20%20%20%20%20%20%20%20df_thresholds%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(df_thresholds%5B%22month%22%5D%20%3D%3D%20_month)%20%26%20(df_thresholds%5B%22pct%22%5D%20%3D%3D%20_pct)%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20.sort_values(%22year%22)%0A%20%20%20%20%20%20%20%20.copy()%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Pre-compute%20tooltip%20columns%20for%20the%20threshold%20line%0A%20%20%20%20_ref_strs%2C%20_ref_sorted%20%3D%20%5B%5D%2C%20%5B%5D%0A%20%20%20%20for%20_%2C%20_row%20in%20_df.iterrows()%3A%0A%20%20%20%20%20%20%20%20_yr%20%3D%20int(_row%5B%22year%22%5D)%0A%20%20%20%20%20%20%20%20_ref%20%3D%20df_iri%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20df_iri%5B%22year%22%5D.between(_yr%20-%20ref_window%2C%20_yr%20-%201)%0A%20%20%20%20%20%20%20%20%5D.sort_values(%22year%22)%0A%20%20%20%20%20%20%20%20_ref_strs.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bint(_ref%5B'year'%5D.min())%7D%E2%80%93%7Bint(_ref%5B'year'%5D.max())%7D%3A%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%2C%20%22.join(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bint(r%5B'year'%5D)%7D%3A%20%7Br%5B_month%5D%3A.0%25%7D%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_%2C%20r%20in%20_ref.iterrows()%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20_ref_sorted.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22sorted%20%E2%86%91%3A%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%2B%20%22%2C%20%22.join(f%22%7Bv%3A.0%25%7D%22%20for%20v%20in%20sorted(_ref%5B_month%5D.values))%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20_df%5B%22ref_years%22%5D%20%3D%20_ref_strs%0A%20%20%20%20_df%5B%22ref_sorted%22%5D%20%3D%20_ref_sorted%0A%20%20%20%20_df%5B%22thresh_fmt%22%5D%20%3D%20(_df%5B%22threshold%22%5D%20*%20100).round(1).astype(str)%20%2B%20%22%25%22%0A%20%20%20%20_df%5B%22actual_fmt%22%5D%20%3D%20(_df%5B%22actual%22%5D%20*%20100).round(1).astype(str)%20%2B%20%22%25%22%0A%20%20%20%20_df%5B%22status%22%5D%20%3D%20_df%5B%22triggered%22%5D.map(%0A%20%20%20%20%20%20%20%20%7BTrue%3A%20%22Triggered%22%2C%20False%3A%20%22Not%20triggered%22%7D%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Pre-eval%20data%0A%20%20%20%20_pre%20%3D%20(%0A%20%20%20%20%20%20%20%20df_iri%5Bdf_iri%5B%22year%22%5D%20%3C%20start_eval_year%5D%5B%5B%22year%22%2C%20_month%5D%5D%0A%20%20%20%20%20%20%20%20.rename(columns%3D%7B_month%3A%20%22actual%22%7D)%0A%20%20%20%20%20%20%20%20.copy()%0A%20%20%20%20)%0A%20%20%20%20_pre%5B%22actual_fmt%22%5D%20%3D%20(_pre%5B%22actual%22%5D%20*%20100).round(1).astype(str)%20%2B%20%22%25%22%0A%0A%20%20%20%20_pre_chart%20%3D%20(%0A%20%20%20%20%20%20%20%20alt.Chart(_pre)%0A%20%20%20%20%20%20%20%20.mark_point(color%3D%22%23cccccc%22%2C%20size%3D55%2C%20filled%3DTrue%2C%20opacity%3D0.8)%0A%20%20%20%20%20%20%20%20.encode(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3Dalt.X(%22year%3AO%22%2C%20title%3D%22Year%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20y%3Dalt.Y(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22actual%3AQ%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title%3D%22IRI%20forecast%20probability%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20axis%3Dalt.Axis(format%3D%22.0%25%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%20tooltip%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22year%3AO%22%2C%20title%3D%22Year%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22actual_fmt%3AN%22%2C%20title%3D%22Actual%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)%0A%20%20%20%20)%0A%20%20%20%20_line%20%3D%20(%0A%20%20%20%20%20%20%20%20alt.Chart(_df)%0A%20%20%20%20%20%20%20%20.mark_line(color%3D%22steelblue%22%2C%20strokeWidth%3D2)%0A%20%20%20%20%20%20%20%20.encode(x%3D%22year%3AO%22%2C%20y%3D%22threshold%3AQ%22)%0A%20%20%20%20)%0A%20%20%20%20_thresh_pts%20%3D%20(%0A%20%20%20%20%20%20%20%20alt.Chart(_df)%0A%20%20%20%20%20%20%20%20.mark_point(color%3D%22steelblue%22%2C%20size%3D65%2C%20filled%3DTrue)%0A%20%20%20%20%20%20%20%20.encode(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3D%22year%3AO%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20y%3D%22threshold%3AQ%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20tooltip%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22year%3AO%22%2C%20title%3D%22Year%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22thresh_fmt%3AN%22%2C%20title%3Df%22Threshold%20(top%20%7B_pct%7D%25)%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22ref_years%3AN%22%2C%20title%3D%22Reference%20window%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22ref_sorted%3AN%22%2C%20title%3D%22Sorted%20values%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)%0A%20%20%20%20)%0A%20%20%20%20_scatter%20%3D%20(%0A%20%20%20%20%20%20%20%20alt.Chart(_df)%0A%20%20%20%20%20%20%20%20.mark_point(size%3D90%2C%20filled%3DTrue%2C%20opacity%3D0.9)%0A%20%20%20%20%20%20%20%20.encode(%0A%20%20%20%20%20%20%20%20%20%20%20%20x%3D%22year%3AO%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20y%3D%22actual%3AQ%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20color%3Dalt.Color(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22status%3AN%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20scale%3Dalt.Scale(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20domain%3D%5B%22Triggered%22%2C%20%22Not%20triggered%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20range%3D%5B%22crimson%22%2C%20%22%23888888%22%5D%2C%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%20title%3DNone%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%20tooltip%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22year%3AO%22%2C%20title%3D%22Year%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22actual_fmt%3AN%22%2C%20title%3D%22Actual%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22thresh_fmt%3AN%22%2C%20title%3D%22Threshold%22)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip(%22status%3AN%22%2C%20title%3D%22Status%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)%0A%20%20%20%20)%0A%20%20%20%20_chart%20%3D%20(%0A%20%20%20%20%20%20%20%20(_pre_chart%20%2B%20_line%20%2B%20_thresh_pts%20%2B%20_scatter)%0A%20%20%20%20%20%20%20%20.properties(%0A%20%20%20%20%20%20%20%20%20%20%20%20title%3Df%22%7B_month%7D%20%E2%80%94%20rolling%20%7Bref_window%7D-yr%20threshold%20at%20top%20%7B_pct%7D%25%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20width%3D720%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20height%3D340%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20.configure_axis(grid%3DFalse)%0A%20%20%20%20%20%20%20%20.configure_view(strokeWidth%3D0)%0A%20%20%20%20)%0A%20%20%20%20mo.ui.altair_chart(_chart)%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

← Niger drought AA