policyengine-uk 2.52.1__tar.gz → 2.53.1__tar.gz

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/CHANGELOG.md

@@ -5,6 +5,18 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/), 
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [2.53.1] - 2025-09-29 12:26:51
+
+### Fixed
+
+- Temporarily disabled OBR participation responses.
+
+## [2.53.0] - 2025-09-25 09:22:41
+
+### Fixed
+
+- Bugs in the benefit cap exemption list.
+
 ## [2.52.1] - 2025-09-24 15:33:18
 
 ### Fixed
@@ -2236,6 +2248,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 
 
+[2.53.1]: https://github.com/PolicyEngine/openfisca-uk/compare/2.53.0...2.53.1
+[2.53.0]: https://github.com/PolicyEngine/openfisca-uk/compare/2.52.1...2.53.0
 [2.52.1]: https://github.com/PolicyEngine/openfisca-uk/compare/2.52.0...2.52.1
 [2.52.0]: https://github.com/PolicyEngine/openfisca-uk/compare/2.51.0...2.52.0
 [2.51.0]: https://github.com/PolicyEngine/openfisca-uk/compare/2.50.0...2.51.0

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: policyengine-uk
-Version: 2.52.1
+Version: 2.53.1
 Summary: PolicyEngine tax and benefit system for the UK.
 Project-URL: Homepage, https://github.com/PolicyEngine/policyengine-uk
 Project-URL: Repository, https://github.com/PolicyEngine/policyengine-uk

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/policyengine_uk/dynamics/labour_supply.py

@@ -58,7 +58,7 @@ def calculate_excluded_from_labour_supply_responses(
         | (adult_index == 0)
         | (adult_index >= count_adults + 1)
     )
-    return excluded
+    return excluded.values
 
 
 class FTEImpacts(BaseModel):
@@ -132,8 +132,9 @@ def apply_labour_supply_responses(
     )
     reform_income = sim.calculate(target_variable, year, map_to="person")
 
-    baseline_income = baseline_income.values
-    reform_income = reform_income.values
+    baseline_income = baseline_income
+    reform_income = reform_income
+    baseline_income = pd.Series(baseline_income, index=reform_income.index)
 
     # Calculate relative changes
     income_rel_change = np.where(
@@ -154,7 +155,9 @@ def apply_labour_supply_responses(
     )
 
     # Apply extensive margin responses (participation model)
-    participation_responses = apply_participation_responses(sim=sim, year=year)
+    participation_responses = (
+        None  # = apply_participation_responses(sim=sim, year=year)
+    )
 
     # Add FTE impacts to the response data
     fte_impacts = FTEImpacts(
@@ -336,11 +339,10 @@ def apply_progression_responses(
     response = response_df["total_response"].values
 
     # Apply the labour supply response to the simulation
-    # NOTE: Don't reset calculations as this breaks UC and other benefit calculations
-    # Instead, just update the employment income with the behavioral response
+    sim.reset_calculations()
     sim.set_input(input_variable, year, employment_income + response)
 
-    weight = sim.calculate("household_weight", year, map_to="person")
+    weight = sim.calculate("household_weight", year, map_to="person").values
 
     result = MicroDataFrame(df, weights=weight)
 

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/policyengine_uk/dynamics/progression.py

@@ -38,7 +38,7 @@ def calculate_derivative(
         Array of marginal rates clipped between 0 and 1
     """
     # Get baseline values for input variable and identify adults
-    input_variable_values = sim.calculate(input_variable, year)
+    input_variable_values = sim.calculate(input_variable, year).copy()
     adult_index = sim.calculate("adult_index")
     entity_key = sim.tax_benefit_system.variables[input_variable].entity.key
 
@@ -67,6 +67,10 @@ def calculate_derivative(
         ~pd.Series(adult_index).isin(range(1, count_adults + 1))
     ] = np.nan
 
+    # Reset simulation to original state
+    sim.reset_calculations()
+    sim.set_input(input_variable, year, input_variable_values)
+
     # Clip to ensure rates are between 0 and 1 (0% to 100% retention)
     return rel_marginal_wages.clip(0, 1)
 
@@ -96,6 +100,9 @@ def calculate_relative_income_change(
     reformed_target_values = sim.calculate(
         target_variable, year, map_to="person"
     )
+    original_target_values = pd.Series(
+        original_target_values, index=reformed_target_values.index
+    )
 
     # Calculate relative change, handling division by zero
     rel_change = (
@@ -159,6 +166,7 @@ def calculate_derivative_change(
         count_adults=count_adults,
         delta=delta,
     )
+    original_deriv = pd.Series(original_deriv, index=reformed_deriv.index)
 
     # Calculate relative and absolute changes in marginal rates
     rel_change = reformed_deriv / original_deriv - 1

policyengine_uk-2.53.1/policyengine_uk/variables/gov/dwp/is_benefit_cap_exempt.py

@@ -0,0 +1,19 @@
+from policyengine_uk.model_api import *
+
+
+class is_benefit_cap_exempt(Variable):
+    value_type = bool
+    entity = BenUnit
+    label = (
+        "Whether exempt from the benefits cap because of health or disability"
+    )
+    definition_period = YEAR
+    reference = "https://www.gov.uk/benefit-cap/when-youre-not-affected"
+
+    def formula(benunit, period, parameters):
+        exempt_health = benunit(
+            "is_benefit_cap_exempt_health_disability", period
+        )
+        exempt_other = benunit("is_benefit_cap_exempt_other", period)
+        exempt_earnings = benunit("is_benefit_cap_exempt_earnings", period)
+        return exempt_health | exempt_earnings | exempt_other

policyengine_uk-2.53.1/policyengine_uk/variables/gov/dwp/is_benefit_cap_exempt_earnings.py

@@ -0,0 +1,66 @@
+from policyengine_uk.model_api import *
+
+
+class is_benefit_cap_exempt_earnings(Variable):
+    value_type = bool
+    entity = BenUnit
+    label = "Whether exempt from the benefits cap for non-health/disability reasons"
+    definition_period = YEAR
+    reference = "https://www.gov.uk/benefit-cap/when-youre-not-affected"
+
+    def formula(benunit, period, parameters):
+        # Check if anyone in benefit unit is over state pension age
+        person = benunit.members
+        over_pension_age = person("is_SP_age", period)
+        has_pensioner = benunit.any(over_pension_age)
+
+        # UC-specific exemptions
+        # Limited capability for work and work-related activity
+        has_lcwra = benunit.any(
+            person("uc_limited_capability_for_WRA", period)
+        )
+
+        # Carer element in UC indicates caring for someone with disability
+        gets_uc_carer_element = benunit("uc_carer_element", period) > 0
+
+        # Earnings exemption for UC (£846/month = £10,152/year)
+        # Note: Only check earned income, not UC amount itself to avoid circular dependency
+        uc_earned = benunit.sum(
+            benunit.members("employment_income", period)
+            + benunit.members("self_employment_income", period)
+            - benunit.members("income_tax", period)
+            - benunit.members("national_insurance", period)
+        )
+        earnings_threshold = 10_152
+        meets_earnings_test = uc_earned >= earnings_threshold
+
+        # Disability and carer benefits that exempt from cap
+        QUAL_PERSONAL_BENEFITS = [
+            "attendance_allowance",
+            "carers_allowance",
+            "dla",  # Disability Living Allowance (includes components)
+            "pip_dl",  # PIP daily living component
+            "pip_m",  # PIP mobility component
+            "iidb",  # Industrial injuries disability benefit
+        ]
+
+        # ESA and Working Tax Credit
+        QUAL_BENUNIT_BENEFITS = [
+            "esa_income",  # Income-based ESA
+            "working_tax_credit",  # If getting WTC, likely working enough
+        ]
+
+        qualifying_personal_benefits = add(
+            benunit, period, QUAL_PERSONAL_BENEFITS
+        )
+        qualifying_benunit_benefits = add(
+            benunit, period, QUAL_BENUNIT_BENEFITS
+        )
+
+        # Check for Armed Forces Compensation Scheme payments
+        afcs = benunit("afcs", period) > 0
+
+        # ESA contribution-based with support component
+        esa_support_component = benunit("esa_contrib", period) > 0
+
+        return meets_earnings_test

policyengine_uk-2.52.1/policyengine_uk/variables/gov/dwp/is_benefit_cap_exempt.py → policyengine_uk-2.53.1/policyengine_uk/variables/gov/dwp/is_benefit_cap_exempt_health_disability.py

@@ -1,10 +1,12 @@
 from policyengine_uk.model_api import *
 
 
-class is_benefit_cap_exempt(Variable):
+class is_benefit_cap_exempt_health_disability(Variable):
     value_type = bool
     entity = BenUnit
-    label = "Whether exempt from the benefits cap"
+    label = (
+        "Whether exempt from the benefits cap because of health or disability"
+    )
     definition_period = YEAR
     reference = "https://www.gov.uk/benefit-cap/when-youre-not-affected"
 
@@ -25,7 +27,12 @@ class is_benefit_cap_exempt(Variable):
 
         # Earnings exemption for UC (£846/month = £10,152/year)
         # Note: Only check earned income, not UC amount itself to avoid circular dependency
-        uc_earned = benunit("uc_earned_income", period)
+        uc_earned = benunit.sum(
+            benunit.members("employment_income", period)
+            + benunit.members("self_employment_income", period)
+            - benunit.members("income_tax", period)
+            - benunit.members("national_insurance", period)
+        )
         earnings_threshold = 10_152
         meets_earnings_test = uc_earned >= earnings_threshold
 
@@ -59,10 +66,8 @@ class is_benefit_cap_exempt(Variable):
         esa_support_component = benunit("esa_contrib", period) > 0
 
         return (
-            has_pensioner
-            | has_lcwra
+            has_lcwra
             | gets_uc_carer_element
-            | meets_earnings_test
             | (qualifying_personal_benefits > 0)
             | (qualifying_benunit_benefits > 0)
             | afcs

policyengine_uk-2.53.1/policyengine_uk/variables/gov/dwp/is_benefit_cap_exempt_other.py

@@ -0,0 +1,66 @@
+from policyengine_uk.model_api import *
+
+
+class is_benefit_cap_exempt_other(Variable):
+    value_type = bool
+    entity = BenUnit
+    label = "Whether exempt from the benefits cap for non-health/disability reasons"
+    definition_period = YEAR
+    reference = "https://www.gov.uk/benefit-cap/when-youre-not-affected"
+
+    def formula(benunit, period, parameters):
+        # Check if anyone in benefit unit is over state pension age
+        person = benunit.members
+        over_pension_age = person("is_SP_age", period)
+        has_pensioner = benunit.any(over_pension_age)
+
+        # UC-specific exemptions
+        # Limited capability for work and work-related activity
+        has_lcwra = benunit.any(
+            person("uc_limited_capability_for_WRA", period)
+        )
+
+        # Carer element in UC indicates caring for someone with disability
+        gets_uc_carer_element = benunit("uc_carer_element", period) > 0
+
+        # Earnings exemption for UC (£846/month = £10,152/year)
+        # Note: Only check earned income, not UC amount itself to avoid circular dependency
+        uc_earned = benunit.sum(
+            benunit.members("employment_income", period)
+            + benunit.members("self_employment_income", period)
+            - benunit.members("income_tax", period)
+            - benunit.members("national_insurance", period)
+        )
+        earnings_threshold = 10_152
+        meets_earnings_test = uc_earned >= earnings_threshold
+
+        # Disability and carer benefits that exempt from cap
+        QUAL_PERSONAL_BENEFITS = [
+            "attendance_allowance",
+            "carers_allowance",
+            "dla",  # Disability Living Allowance (includes components)
+            "pip_dl",  # PIP daily living component
+            "pip_m",  # PIP mobility component
+            "iidb",  # Industrial injuries disability benefit
+        ]
+
+        # ESA and Working Tax Credit
+        QUAL_BENUNIT_BENEFITS = [
+            "esa_income",  # Income-based ESA
+            "working_tax_credit",  # If getting WTC, likely working enough
+        ]
+
+        qualifying_personal_benefits = add(
+            benunit, period, QUAL_PERSONAL_BENEFITS
+        )
+        qualifying_benunit_benefits = add(
+            benunit, period, QUAL_BENUNIT_BENEFITS
+        )
+
+        # Check for Armed Forces Compensation Scheme payments
+        afcs = benunit("afcs", period) > 0
+
+        # ESA contribution-based with support component
+        esa_support_component = benunit("esa_contrib", period) > 0
+
+        return has_pensioner | afcs | esa_support_component

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/policyengine_uk/variables/household/demographic/is_disabled_for_benefits.py

@@ -14,4 +14,16 @@ class is_disabled_for_benefits(Variable):
             "dla",
             "pip",
         ]
-        return add(person, period, QUALIFYING_BENEFITS) > 0
+
+        p_claims_lcwra_if_on_pip_dla = 0.8
+        p_claims_lcwra_if_not_on_pip_dla = 0.13
+
+        random_seed = random(person)
+
+        on_qual_benefits = add(person, period, QUALIFYING_BENEFITS) > 0
+
+        return np.where(
+            on_qual_benefits,
+            random_seed < p_claims_lcwra_if_on_pip_dla,
+            random_seed < p_claims_lcwra_if_not_on_pip_dla,
+        )

{policyengine_uk-2.52.1 → policyengine_uk-2.53.1}/pyproject.toml

@@ -1,6 +1,6 @@
 [project] 
 name = "policyengine-uk"
-version = "2.52.1"
+version = "2.53.1"
 description = "PolicyEngine tax and benefit system for the UK."
 readme = "README.md"
 authors = [