Preserve household weight fuel uprating

Author: MaxGhenis

changelog.d/fuel-uprating-household-weight.md

@@ -0,0 +1 @@
+Preserve household-weight uprating when calibrating fuel litre proxies.

policyengine_uk_data/storage/uprating_factors.csv

@@ -12,7 +12,7 @@ childcare_expenses,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,
 clothing_and_footwear_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 communication_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 corporate_wealth,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384
-diesel_spending,1.0,1.138,1.56,1.531,1.579,1.537,1.509,1.467,1.422,1.363,1.286,1.28,1.275,1.27,1.265
+diesel_spending,1.0,1.138,1.57,1.539,1.589,1.54,1.507,1.46,1.415,1.355,1.284,1.284,1.284,1.284,1.284
 dividend_income,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384
 domestic_energy_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 education_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
@@ -30,7 +30,7 @@ free_school_milk,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.
 gross_financial_wealth,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384
 health_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 household_furnishings_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
-household_weight,1.0,1.0,1.009,1.023,1.033,1.041,1.045,1.049,1.053,1.058,1.062,1.067,1.071,1.076,1.08
+household_weight,1.0,1.0,1.003,1.017,1.027,1.039,1.046,1.054,1.058,1.064,1.064,1.064,1.064,1.064,1.064
 housing_benefit_reported,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 housing_service_charges,1.0,1.0,1.0,1.064,1.137,1.191,1.235,1.262,1.289,1.318,1.318,1.318,1.318,1.318,1.318
 housing_water_and_electricity_consumption,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
@@ -56,7 +56,7 @@ owned_land,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.3
 pension_credit_reported,1.0,1.04,1.144,1.209,1.237,1.277,1.301,1.327,1.353,1.38,1.38,1.38,1.38,1.38,1.38
 pension_income,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384
 personal_pension_contributions,1.0,1.059,1.127,1.205,1.261,1.308,1.337,1.365,1.396,1.431,1.431,1.431,1.431,1.431,1.431
-petrol_spending,1.0,1.104,1.624,1.786,1.521,1.48,1.454,1.413,1.37,1.313,1.239,1.233,1.228,1.223,1.218
+petrol_spending,1.0,1.104,1.635,1.796,1.531,1.483,1.452,1.406,1.363,1.305,1.237,1.237,1.237,1.237,1.237
 private_pension_income,1.0,1.003,1.053,1.106,1.161,1.216,1.261,1.288,1.315,1.346,1.346,1.346,1.346,1.346,1.346
 private_transfer_income,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384
 property_income,1.0,1.0,1.092,1.147,1.19,1.223,1.258,1.297,1.34,1.384,1.384,1.384,1.384,1.384,1.384

policyengine_uk_data/storage/uprating_growth_factors.csv

@@ -12,7 +12,7 @@ childcare_expenses,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0
 clothing_and_footwear_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 communication_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 corporate_wealth,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0
-diesel_spending,0,0.138,0.371,-0.019,0.031,-0.027,-0.018,-0.028,-0.031,-0.041,-0.056,-0.005,-0.004,-0.004,-0.004
+diesel_spending,0,0.138,0.379,-0.019,0.033,-0.031,-0.021,-0.032,-0.031,-0.043,-0.052,0.0,0.0,0.0,0.0
 dividend_income,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0
 domestic_energy_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 education_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
@@ -30,7 +30,7 @@ free_school_milk,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0
 gross_financial_wealth,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0
 health_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 household_furnishings_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
-household_weight,0,0.0,0.009,0.014,0.01,0.008,0.004,0.004,0.004,0.005,0.004,0.005,0.004,0.005,0.004
+household_weight,0,0.0,0.003,0.014,0.01,0.012,0.007,0.008,0.004,0.006,0.0,0.0,0.0,0.0,0.0
 housing_benefit_reported,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 housing_service_charges,0,0.0,0.0,0.064,0.069,0.047,0.037,0.022,0.021,0.022,0.0,0.0,0.0,0.0,0.0
 housing_water_and_electricity_consumption,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
@@ -56,7 +56,7 @@ owned_land,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,
 pension_credit_reported,0,0.04,0.1,0.057,0.023,0.032,0.019,0.02,0.02,0.02,0.0,0.0,0.0,0.0,0.0
 pension_income,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0
 personal_pension_contributions,0,0.059,0.064,0.069,0.046,0.037,0.022,0.021,0.023,0.025,0.0,0.0,0.0,0.0,0.0
-petrol_spending,0,0.104,0.471,0.1,-0.148,-0.027,-0.018,-0.028,-0.03,-0.042,-0.056,-0.005,-0.004,-0.004,-0.004
+petrol_spending,0,0.104,0.481,0.099,-0.147,-0.031,-0.021,-0.032,-0.031,-0.043,-0.052,0.0,0.0,0.0,0.0
 private_pension_income,0,0.003,0.05,0.05,0.05,0.047,0.037,0.021,0.021,0.024,0.0,0.0,0.0,0.0,0.0
 private_transfer_income,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0
 property_income,0,0.0,0.092,0.05,0.037,0.028,0.029,0.031,0.033,0.033,0.0,0.0,0.0,0.0,0.0

policyengine_uk_data/tests/test_road_fuel_volume_uprating.py

@@ -78,21 +78,29 @@ def test__given_pre_pandemic_base__then_obr_forecast_volume_declines():
 def test__given_uprating_table__then_only_fuel_rows_are_overridden():
     # Given
     df = pd.DataFrame(
-        {year: [1.0, 1.0, 1.0] for year in range(START_YEAR, END_YEAR + 1)},
-        index=["petrol_spending", "diesel_spending", "some_other_variable"],
+        {year: [1.0, 1.0, 1.0, 1.0] for year in range(START_YEAR, END_YEAR + 1)},
+        index=[
+            "petrol_spending",
+            "diesel_spending",
+            "household_weight",
+            "some_other_variable",
+        ],
     )
+    df.loc["household_weight", 2024] = 1.2
 
     # When
     out = _apply_road_fuel_litre_proxy_override(df.copy())
 
     # Then
     for year in range(START_YEAR, END_YEAR + 1):
         assert out.loc["some_other_variable", year] == 1.0
+        assert out.loc["household_weight", year] == df.loc["household_weight", year]
     for variable in VOLUME_OVERRIDDEN_VARIABLES:
         expected = fuel_spending_litre_proxy_index(
             variable=variable,
             base_year=START_YEAR,
             end_year=END_YEAR,
+            household_weight_index=df.loc["household_weight"],
         )
         for year in range(START_YEAR, END_YEAR + 1):
             assert out.loc[variable, year] == round(expected[year], 3)
@@ -108,6 +116,7 @@ def test__given_storage_csv__then_fuel_rows_reflect_litre_proxy_index():
             variable=variable,
             base_year=START_YEAR,
             end_year=END_YEAR,
+            household_weight_index=df.loc["household_weight"],
         )
         assert variable in df.index
         for year in range(START_YEAR, END_YEAR + 1):

policyengine_uk_data/utils/uprating.py

@@ -88,6 +88,7 @@ def fuel_spending_litre_proxy_index(
     base_year: int = START_YEAR,
     end_year: int = END_YEAR,
     parameters=None,
+    household_weight_index: dict[int, float] | pd.Series | None = None,
 ) -> dict[int, float]:
     """Return the spending-proxy index that preserves fuel litres.
 
@@ -112,16 +113,38 @@ def fuel_spending_litre_proxy_index(
         parameters.household.consumption.fuel.prices,
         FUEL_PRICE_PARAMETER_NAME[variable],
     )
-    population = parameters.gov.economic_assumptions.indices.ons.population
     base_price = price_parameter(base_year)
-    base_population = population(base_year)
+
+    if household_weight_index is None:
+        population = parameters.gov.economic_assumptions.indices.ons.population
+        base_household_weight_index = population(base_year)
+
+        def household_weight_relative(year: int) -> float:
+            return population(year) / base_household_weight_index
+
+    else:
+
+        def household_weight_value(year: int) -> float:
+            if hasattr(household_weight_index, "index"):
+                if year in household_weight_index.index:
+                    return float(household_weight_index.loc[year])
+                return float(household_weight_index.loc[str(year)])
+            try:
+                return float(household_weight_index[year])
+            except KeyError:
+                return float(household_weight_index[str(year)])
+
+        base_household_weight_index = household_weight_value(base_year)
+
+        def household_weight_relative(year: int) -> float:
+            return household_weight_value(year) / base_household_weight_index
+
     return {
         year: (
             volume_index[year]
             * price_parameter(year)
             / base_price
-            * base_population
-            / population(year)
+            / household_weight_relative(year)
         )
         for year in volume_index
     }
@@ -141,10 +164,14 @@ def _apply_road_fuel_litre_proxy_override(df: pd.DataFrame) -> pd.DataFrame:
     for variable in VOLUME_OVERRIDDEN_VARIABLES:
         if variable not in df.index:
             continue
+        household_weight_index = (
+            df.loc["household_weight"] if "household_weight" in df.index else None
+        )
         index = fuel_spending_litre_proxy_index(
             variable=variable,
             base_year=START_YEAR,
             end_year=END_YEAR,
+            household_weight_index=household_weight_index,
         )
         missing_years = [
             year for year in range(START_YEAR, END_YEAR + 1) if year not in index