|
| 1 | +import pandas as pd |
| 2 | +from pathlib import Path |
| 3 | + |
| 4 | + |
| 5 | +def project_water_bills(): |
| 6 | + df_pre_2025 = pd.DataFrame( |
| 7 | + { |
| 8 | + "Year": [2021, 2022, 2023, 2024, 2025], |
| 9 | + "Oftwat avg bills (real)": [486, 470, 486, 492, 503], |
| 10 | + "CPIH": [113.1, 123.0, 129.9, 134.0, 139.0], |
| 11 | + } |
| 12 | + ) |
| 13 | + |
| 14 | + df_pre_2025["Oftwat avg bills (nominal)"] = ( |
| 15 | + df_pre_2025["Oftwat avg bills (real)"] * df_pre_2025["CPIH"] / 100 |
| 16 | + ) |
| 17 | + df_pre_2025["Oftwat avg bills (nominal)"] = ( |
| 18 | + df_pre_2025["Oftwat avg bills (nominal)"] |
| 19 | + / df_pre_2025["Oftwat avg bills (nominal)"].iloc[0] |
| 20 | + * 100 |
| 21 | + ).round(1) |
| 22 | + df_pre_2025["Nominal YoY change"] = ( |
| 23 | + df_pre_2025["Oftwat avg bills (nominal)"].pct_change() * 100 |
| 24 | + ).round(1) |
| 25 | + |
| 26 | + proposed_increases = pd.read_csv( |
| 27 | + Path(__file__).parent / "ofwat_increases.csv" |
| 28 | + ) |
| 29 | + avg_bills_2025_onwards = ( |
| 30 | + proposed_increases[proposed_increases.columns[1:]].mean()[1:].values |
| 31 | + ) |
| 32 | + |
| 33 | + df_post_2025 = pd.DataFrame( |
| 34 | + { |
| 35 | + "Year": [2025, 2026, 2027, 2028, 2029], |
| 36 | + "CPIH": [139.0, 142.2, 145.2, 148.2, 151.3], |
| 37 | + "Pre-inflation avg bills (nominal)": avg_bills_2025_onwards, |
| 38 | + } |
| 39 | + ) |
| 40 | + |
| 41 | + # Add CPIH to each year's change |
| 42 | + |
| 43 | + df_post_2025["Avg bills (nominal)"] = df_post_2025[ |
| 44 | + "Pre-inflation avg bills (nominal)" |
| 45 | + ].values |
| 46 | + df_post_2025["CPIH change"] = df_post_2025["CPIH"].pct_change() * 100 |
| 47 | + |
| 48 | + for year in range(2026, 2030): |
| 49 | + row = df_post_2025[df_post_2025["Year"] == year].iloc[0] |
| 50 | + cpi_change = ( |
| 51 | + row["CPIH"] |
| 52 | + / df_post_2025[df_post_2025["Year"] == year - 1]["CPIH"].values[0] |
| 53 | + - 1 |
| 54 | + ) * 100 |
| 55 | + # Increase the nominal bills by the CPIH change |
| 56 | + addition = df_post_2025.loc[ |
| 57 | + df_post_2025["Year"] == year - 1, "Avg bills (nominal)" |
| 58 | + ].values[0] * (cpi_change / 100) |
| 59 | + # Add addition to this and future years |
| 60 | + df_post_2025.loc[ |
| 61 | + df_post_2025["Year"] >= year, "Avg bills (nominal)" |
| 62 | + ] += addition |
| 63 | + |
| 64 | + df_post_2025["Relative change"] = ( |
| 65 | + df_post_2025["Avg bills (nominal)"].pct_change() * 100 |
| 66 | + ).round(1) |
| 67 | + |
| 68 | + df_post_2025 |
| 69 | + |
| 70 | + combined_water_forecast = pd.DataFrame( |
| 71 | + { |
| 72 | + "Year": list(range(2022, 2030)), |
| 73 | + "Average nominal bills YoY change": df_pre_2025[ |
| 74 | + "Nominal YoY change" |
| 75 | + ].tolist()[1:] |
| 76 | + + df_post_2025["Relative change"].tolist()[1:], |
| 77 | + } |
| 78 | + ) |
| 79 | + |
| 80 | + print(combined_water_forecast.to_markdown(index=False)) |
0 commit comments