policyengine-us 1.413.0__py3-none-any.whl → 1.439.0__py3-none-any.whl

policyengine_us/variables/gov/fcc/lifeline/is_lifeline_eligible.py

@@ -21,10 +21,9 @@ class is_lifeline_eligible(Variable):
         )
         categorically_eligible = np.where(
             is_on_tribal_land,
-            np.any(tribal_lifeline_programs),
-            np.any(non_tribal_lifeline_programs),
+            tribal_lifeline_programs > 0,
+            non_tribal_lifeline_programs > 0,
         )
-        fpg_ratio = spm_unit("fcc_fpg_ratio", period)
-        fpg_limit = p.fpg_limit
-        fpg_eligible = fpg_ratio <= fpg_limit
-        return categorically_eligible | fpg_eligible
+        # Use the new unified income eligibility variable
+        income_eligible = spm_unit("is_lifeline_income_eligible", period)
+        return categorically_eligible | income_eligible

policyengine_us/variables/gov/fcc/lifeline/is_lifeline_income_eligible.py

@@ -0,0 +1,35 @@
+from policyengine_us.model_api import *
+
+
+class is_lifeline_income_eligible(Variable):
+    value_type = bool
+    entity = SPMUnit
+    label = "Income-eligible for Lifeline"
+    documentation = (
+        "Meets income requirements for Lifeline (federal or state-expanded)"
+    )
+    definition_period = YEAR
+    reference = (
+        "https://www.law.cornell.edu/cfr/text/47/54.409",
+        "https://statutes.capitol.texas.gov/Docs/UT/htm/UT.55.htm#55.015",
+    )
+
+    def formula(spm_unit, period, parameters):
+        p = parameters(period).gov.fcc.lifeline
+        fpg_ratio = spm_unit("fcc_fpg_ratio", period)
+
+        # Some states expand eligibility beyond federal limit
+        state = spm_unit.household("state_code_str", period)
+        is_tx = state == "TX"
+
+        # Texas expands to 150% FPG per Texas Utilities Code § 55.015(d-1)
+        tx_p = parameters(period).gov.states.tx.uct.lifeline
+
+        # Use state-specific limit where applicable, otherwise federal (135%)
+        # Convert to float32 to match fpg_ratio dtype for comparison
+        # Without this, float32(1.35) > float64(1.35) due to precision differences
+        fpg_limit = where(is_tx, tx_p.fpg_limit, p.fpg_limit).astype(
+            np.float32
+        )
+
+        return fpg_ratio <= fpg_limit

policyengine_us/variables/gov/fcc/lifeline/lifeline.py

@@ -33,11 +33,14 @@ class lifeline(Variable):
                     base_amount,
                 )
 
-        max_monthly_amount = (
+        federal_monthly_amount = (
             base_amount
             + is_rural_tribal * p.fcc.lifeline.amount.rural_tribal_supplement
         )
-        max_amount = max_monthly_amount * MONTHS_IN_YEAR
+        # Get TX supplement for the specific month (now properly defined as MONTH period)
+        tx_supplement = spm_unit("tx_lifeline_supplement", period.first_month)
+        total_monthly_amount = federal_monthly_amount + tx_supplement
+        max_amount = total_monthly_amount * MONTHS_IN_YEAR
         phone_broadband_cost = add(
             spm_unit, period, ["phone_cost", "broadband_cost"]
         )

policyengine_us/variables/gov/hhs/head_start/early_head_start.py

@@ -4,7 +4,8 @@ from policyengine_us.model_api import *
 class early_head_start(Variable):
     value_type = float
     entity = Person
-    label = "Amount of Early Head Start benefit"
+    label = "Early Head Start"
+    unit = USD
     definition_period = YEAR
     defined_for = "is_early_head_start_eligible"
     reference = "https://headstart.gov/program-data/article/head-start-program-facts-fiscal-year-2022"

policyengine_us/variables/gov/hhs/head_start/head_start.py

@@ -4,7 +4,8 @@ from policyengine_us.model_api import *
 class head_start(Variable):
     value_type = float
     entity = Person
-    label = "Amount of Head Start benefit"
+    label = "Head Start"
+    unit = USD
     definition_period = YEAR
     defined_for = "is_head_start_eligible"
     reference = "https://headstart.gov/program-data/article/head-start-program-facts-fiscal-year-2022"

policyengine_us/variables/gov/hhs/head_start/is_head_start_categorically_eligible.py

@@ -12,7 +12,5 @@ class is_head_start_categorically_eligible(Variable):
     )
 
     def formula(person, period, parameters):
-        tax_unit = person.tax_unit
         p = parameters(period).gov.hhs.head_start
-        programs = add(tax_unit, period, p.categorical_eligibility)
-        return np.any(programs)
+        return add(person.tax_unit, period, p.categorical_eligibility) > 0

policyengine_us/variables/gov/hhs/tanf/cash/income/tanf_gross_earned_income.py

@@ -0,0 +1,13 @@
+from policyengine_us.model_api import *
+
+
+class tanf_gross_earned_income(Variable):
+    value_type = float
+    entity = Person
+    label = (
+        "Temporary Assistance for Needy Families (TANF) gross earned income"
+    )
+    unit = USD
+    definition_period = MONTH
+
+    adds = "gov.hhs.tanf.cash.income.sources.earned"

policyengine_us/variables/gov/hhs/tanf/cash/income/tanf_gross_unearned_income.py

@@ -0,0 +1,13 @@
+from policyengine_us.model_api import *
+
+
+class tanf_gross_unearned_income(Variable):
+    value_type = float
+    entity = Person
+    label = (
+        "Temporary Assistance for Needy Families (TANF) gross unearned income"
+    )
+    unit = USD
+    definition_period = MONTH
+
+    adds = "gov.hhs.tanf.cash.income.sources.unearned"

policyengine_us/variables/gov/hhs/tanf/tanf_fpg.py

@@ -0,0 +1,28 @@
+from policyengine_us.model_api import *
+
+
+class tanf_fpg(Variable):
+    value_type = float
+    entity = SPMUnit
+    label = "TANF federal poverty guideline"
+    unit = USD
+    documentation = (
+        "The federal poverty guideline used to determine TANF eligibility."
+    )
+    definition_period = MONTH
+    reference = "https://www.law.cornell.edu/uscode/text/42/9902"
+
+    def formula(spm_unit, period, parameters):
+        # We can use tanf_unit_size for more accurate model in the future
+        n = spm_unit("spm_unit_size", period.this_year)
+        state_group = spm_unit.household("state_group_str", period.this_year)
+        year = period.start.year
+        month = period.start.month
+        if month >= 10:
+            instant_str = f"{year}-10-01"
+        else:
+            instant_str = f"{year - 1}-10-01"
+        p_fpg = parameters(instant_str).gov.hhs.fpg
+        p1 = p_fpg.first_person[state_group] / MONTHS_IN_YEAR
+        pn = p_fpg.additional_person[state_group] / MONTHS_IN_YEAR
+        return p1 + pn * (n - 1)

policyengine_us/variables/gov/irs/tax/federal_income/alternative_minimum_tax/amt_tax_including_cg.py

@@ -34,7 +34,7 @@ class amt_tax_including_cg(Variable):
         income_taxes_at_amt_rates = p.income.amt.brackets.calc(excess_income)
         # Line 19 - First CG tax bracket threshold
         filing_status = tax_unit("filing_status", period)
-        cg_bracket = p.capital_gains.brackets.thresholds["1"][filing_status]
+        cg_bracket = p.capital_gains.thresholds["1"][filing_status]
         # Line 20 - Schedule D Line 14
         lt_capital_loss_carryover = tax_unit("dwks14", period)
         # Line 21 - Line 20 minus Line 19
@@ -42,7 +42,7 @@ class amt_tax_including_cg(Variable):
         # Line 22 - smaller of Line 12 or Line 13
         smaller_of_income_or_cg = min_(reduced_income, cg_distributions)
         # Line 23 - smaller of Line 22 or Line 21 (amount is taxed at 0%)
-        cg_first_rate = p.capital_gains.brackets.rates["1"]
+        cg_first_rate = p.capital_gains.rates["1"]
         disregarded_gains = (
             min_(smaller_of_income_or_cg, reduced_cg_bracket) * cg_first_rate
         )
@@ -51,9 +51,7 @@ class amt_tax_including_cg(Variable):
             smaller_of_income_or_cg - disregarded_gains, 0
         )
         # Line 25 - Second CG tax bracket threshold
-        second_cg_bracket = p.capital_gains.brackets.thresholds["2"][
-            filing_status
-        ]
+        second_cg_bracket = p.capital_gains.thresholds["2"][filing_status]
         # Line 26 - same as line 21
         # Line 27 - Schedule D Line 21
         loss_limited_net_capital_gains = tax_unit(
@@ -73,16 +71,14 @@ class amt_tax_including_cg(Variable):
         )
         # Line 31 - multiply Line 30by second CG tax rate
         cg_second_bracket_tax = (
-            capped_income_including_cg * p.capital_gains.brackets.rates["2"]
+            capped_income_including_cg * p.capital_gains.rates["2"]
         )
         # Line 32 - Line 23 plus Line 30
         taxed_gains = disregarded_gains + capped_income_including_cg
         # Line 33 - Line 22 minus Line 32
         excess_taxed_gains = max_(0, smaller_of_income_or_cg - taxed_gains)
         # Line 34 - multiply Line 33 by third CG tax rate
-        cg_third_bracket_tax = (
-            excess_taxed_gains * p.capital_gains.brackets.rates["3"]
-        )
+        cg_third_bracket_tax = excess_taxed_gains * p.capital_gains.rates["3"]
         # Line 35 - sum of Line 17, Line 32, Line 33
         final_taxed_income = excess_income + taxed_gains + excess_taxed_gains
         # Line 36 Line 12 minus Line 35

policyengine_us/variables/gov/irs/tax/federal_income/alternative_minimum_tax/regular_tax_before_credits.py

@@ -10,29 +10,28 @@ class regular_tax_before_credits(Variable):
     unit = USD
 
     def formula(tax_unit, period, parameters):
+        p = parameters(period).gov.irs
         filing_status = tax_unit("filing_status", period)
         dwks1 = tax_unit("taxable_income", period)
 
-        capital_gains = parameters(period).gov.irs.capital_gains.brackets
-
-        dwks16 = min_(capital_gains.thresholds["1"][filing_status], dwks1)
+        dwks16 = min_(p.capital_gains.thresholds["1"][filing_status], dwks1)
         dwks17 = min_(tax_unit("dwks14", period), dwks16)
         dwks20 = dwks16 - dwks17
-        lowest_rate_tax = capital_gains.rates["1"] * dwks20
+        lowest_rate_tax = p.capital_gains.rates["1"] * dwks20
         # Break in worksheet lines
         dwks13 = tax_unit("dwks13", period)
         dwks21 = min_(dwks1, dwks13)
         dwks22 = dwks20
         dwks23 = max_(0, dwks21 - dwks22)
-        dwks25 = min_(capital_gains.thresholds["2"][filing_status], dwks1)
+        dwks25 = min_(p.capital_gains.thresholds["2"][filing_status], dwks1)
         dwks19 = tax_unit("dwks19", period)
         dwks26 = min_(dwks19, dwks20)
         dwks27 = max_(0, dwks25 - dwks26)
         dwks28 = min_(dwks23, dwks27)
-        dwks29 = capital_gains.rates["2"] * dwks28
+        dwks29 = p.capital_gains.rates["2"] * dwks28
         dwks30 = dwks22 + dwks28
         dwks31 = dwks21 - dwks30
-        dwks32 = capital_gains.rates["3"] * dwks31
+        dwks32 = p.capital_gains.rates["3"] * dwks31
         # Break in worksheet lines
         dwks33 = min_(
             tax_unit("dwks09", period),
@@ -43,18 +42,18 @@ class regular_tax_before_credits(Variable):
         dwks36 = max_(0, dwks34 - dwks1)
         dwks37 = max_(0, dwks33 - dwks36)
 
-        p = parameters(period).gov.irs.income
-
-        dwks38 = p.amt.capital_gains.capital_gain_excess_tax_rate * dwks37
+        dwks38 = (
+            p.income.amt.capital_gains.capital_gain_excess_tax_rate * dwks37
+        )
         # Break in worksheet lines
         dwks39 = dwks19 + dwks20 + dwks28 + dwks31 + dwks37
         dwks40 = dwks1 - dwks39
-        dwks41 = p.amt.brackets.rates[-1] * dwks40
+        dwks41 = p.income.amt.brackets.rates[-1] * dwks40
 
         # Compute regular tax using bracket rates and thresholds
         reg_taxinc = max_(0, dwks19)
-        bracket_tops = p.bracket.thresholds
-        bracket_rates = p.bracket.rates
+        bracket_tops = p.income.bracket.thresholds
+        bracket_rates = p.income.bracket.rates
         reg_tax = 0
         bracket_bottom = 0
         for i in range(1, len(list(bracket_rates.__iter__())) + 1):

policyengine_us/variables/gov/irs/tax/federal_income/capital_gains/capital_gains_excluded_from_taxable_income.py

@@ -14,15 +14,15 @@ class capital_gains_excluded_from_taxable_income(Variable):
     )
 
     def formula(tax_unit, period, parameters):
+        p = parameters(period).gov.irs.capital_gains
         net_capital_gain = tax_unit("net_capital_gain", period)
         adjusted_net_capital_gain = tax_unit(
             "adjusted_net_capital_gain", period
         )
         taxable_income = tax_unit("taxable_income", period)
         filing_status = tax_unit("filing_status", period)
-        cg = parameters(period).gov.irs.capital_gains.brackets
         income_taxed_below_first_rate = clip(
-            taxable_income, 0, cg.thresholds["1"][filing_status]
+            taxable_income, 0, p.thresholds["1"][filing_status]
         )
         reduced_taxable_income = max_(
             taxable_income - net_capital_gain,

policyengine_us/variables/gov/irs/tax/federal_income/capital_gains/capital_gains_tax.py

@@ -27,8 +27,8 @@ class capital_gains_tax(Variable):
 
         filing_status = tax_unit("filing_status", period)
 
-        first_threshold = cg.brackets.thresholds["1"][filing_status]
-        second_threshold = cg.brackets.thresholds["2"][filing_status]
+        first_threshold = cg.thresholds["1"][filing_status]
+        second_threshold = cg.thresholds["2"][filing_status]
 
         income_ordinarily_under_second_rate = clip(
             taxable_income, 0, first_threshold
@@ -55,9 +55,9 @@ class capital_gains_tax(Variable):
         )
 
         main_cg_tax = (
-            cg_in_first_bracket * cg.brackets.rates["1"]
-            + cg_in_second_bracket * cg.brackets.rates["2"]
-            + cg_in_third_bracket * cg.brackets.rates["3"]
+            cg_in_first_bracket * cg.rates["1"]
+            + cg_in_second_bracket * cg.rates["2"]
+            + cg_in_third_bracket * cg.rates["3"]
         )
 
         unrecaptured_s_1250_gain = tax_unit(

policyengine_us/variables/gov/irs/tax/federal_income/capital_gains/dwks19.py

@@ -9,11 +9,11 @@ class dwks19(Variable):
     unit = USD
 
     def formula(tax_unit, period, parameters):
+        p = parameters(period).gov.irs.capital_gains
         dwks14 = tax_unit("dwks14", period)
-        capital_gains = parameters(period).gov.irs.capital_gains.brackets
         filing_status = tax_unit("filing_status", period)
         dwks01 = tax_unit("taxable_income", period)
-        dwks16 = min_(capital_gains.thresholds["1"][filing_status], dwks01)
+        dwks16 = min_(p.thresholds["1"][filing_status], dwks01)
         dwks17 = min_(dwks14, dwks16)
         dwks10 = tax_unit("dwks10", period)
         dwks18 = max_(0, dwks01 - dwks10)

policyengine_us/variables/gov/local/tx/harris/in_harris_county_tx.py

@@ -0,0 +1,12 @@
+from policyengine_us.model_api import *
+
+
+class in_harris_county_tx(Variable):
+    value_type = bool
+    entity = Household
+    definition_period = YEAR
+    label = "Is in Harris County, Texas"
+
+    def formula(household, period, parameters):
+        county = household("county_str", period)
+        return county == "HARRIS_COUNTY_TX"

policyengine_us/variables/gov/local/tx/harris/rides/tx_harris_rides_eligible.py

@@ -0,0 +1,22 @@
+from policyengine_us.model_api import *
+
+
+class tx_harris_rides_eligible(Variable):
+    value_type = bool
+    entity = Person
+    definition_period = YEAR
+    label = "Harris County RIDES eligibility"
+    reference = "https://rides.harriscountytx.gov/Registration/Registration-Information"
+    defined_for = "in_harris_county_tx"
+
+    def formula(person, period, parameters):
+        p = parameters(period).gov.local.tx.harris.rides
+
+        # Age eligibility
+        age = person("age", period)
+        age_eligible = age >= p.age_threshold
+
+        # Disability eligibility
+        has_disability = person("is_disabled", period)
+
+        return age_eligible | has_disability

policyengine_us/variables/gov/local/tx/harris/rides/tx_harris_rides_subsidy.py

@@ -0,0 +1,26 @@
+from policyengine_us.model_api import *
+
+
+class tx_harris_rides_subsidy(Variable):
+    value_type = float
+    entity = Person
+    definition_period = YEAR
+    label = "Harris County RIDES transportation subsidy"
+    reference = "https://rides.harriscountytx.gov/"
+    unit = USD
+    defined_for = "tx_harris_rides_eligible"
+
+    def formula(person, period, parameters):
+        p = parameters(period).gov.local.tx.harris.rides
+
+        # Get estimated transportation costs
+        # In a real implementation, this would be based on actual trip data
+        # For now, we'll use a simplified estimate
+        transportation_costs = person(
+            "pre_subsidy_transportation_expense", period
+        )
+
+        # RIDES subsidizes 60% of trip costs (customer pays 40%)
+        subsidy_rate = 1 - p.customer_payment_rate
+
+        return transportation_costs * subsidy_rate

policyengine_us/variables/gov/ssa/ssi/eligibility/income/deemed/from_ineligible_spouse/is_ssi_spousal_deeming_applies.py

@@ -0,0 +1,41 @@
+from policyengine_us.model_api import *
+
+
+class is_ssi_spousal_deeming_applies(Variable):
+    value_type = bool
+    entity = Person
+    label = "SSI spousal deeming applies"
+    definition_period = YEAR
+    reference = "https://www.law.cornell.edu/cfr/text/20/416.1163"
+    defined_for = "is_ssi_eligible_individual"
+    documentation = """
+    Returns True when spousal deeming applies according to 20 CFR §416.1163.
+
+    Spousal deeming applies when:
+    1. Person is an eligible individual (not part of an eligible couple)
+    2. Spouse is ineligible (not aged/blind/disabled)
+    3. Spouse's leftover income (after child allocations) exceeds the difference
+       between couple FBR and individual FBR
+
+    When deeming applies, the benefit calculation uses the couple FBR instead
+    of the individual FBR, recognizing that two people have higher expenses.
+    """
+
+    def formula(person, period, parameters):
+        # Get spouse's deemed income (after child allocations)
+        spouse_deemed_income = add(
+            person,
+            period,
+            [
+                "ssi_earned_income_deemed_from_ineligible_spouse",
+                "ssi_unearned_income_deemed_from_ineligible_spouse",
+            ],
+        )
+
+        # Compare to FBR differential (couple rate - individual rate)
+        p = parameters(period).gov.ssa.ssi.amount
+        diff = (p.couple - p.individual) * MONTHS_IN_YEAR
+
+        # Deeming applies when spouse's deemed income exceeds the differential
+        # Note: regulation says "not more than" means ≤, so we use > for deeming
+        return spouse_deemed_income > diff

policyengine_us/variables/gov/ssa/ssi/eligibility/income/deemed/from_ineligible_spouse/ssi_income_deemed_from_ineligible_spouse.py

@@ -22,37 +22,37 @@ class ssi_income_deemed_from_ineligible_spouse(Variable):
     unit = USD
 
     def formula(person, period, parameters):
-        # 1. Sum leftover spouse earned/unearned (post-child allocations).
+        # 1. Check if spousal deeming applies (avoids duplicating threshold logic)
+        deeming_applies = person("is_ssi_spousal_deeming_applies", period)
+
+        # 2. If deeming applies, calculate deemed amount using "difference" approach:
+        #    (couple combined countable) - (individual alone countable)
+
+        # Get spouse's leftover income (post-child allocations)
         spouse_earned = person(
             "ssi_earned_income_deemed_from_ineligible_spouse", period
         )
         spouse_unearned = person(
             "ssi_unearned_income_deemed_from_ineligible_spouse", period
         )
-        leftover_spouse = spouse_earned + spouse_unearned
-
-        # 2. Compare leftover to FBR difference
-        p = parameters(period).gov.ssa.ssi.amount
-        diff = (p.couple - p.individual) * MONTHS_IN_YEAR
-        if_leftover_exceeds = leftover_spouse > diff
-
-        # 3. If leftover <= diff => no deeming
-        #    If leftover > diff => difference of incomes approach
-        #    (couple combined) - (individual alone).
-        #   (a) individual's own income
-        indiv_earned = person("ssi_earned_income", period)
-        indiv_unearned = person("ssi_unearned_income", period)
+
+        # Get individual's own income
+        individual_earned = person("ssi_earned_income", period)
+        individual_unearned = person("ssi_unearned_income", period)
+
+        # Calculate individual's countable income (alone)
         alone_countable = _apply_ssi_exclusions(
-            indiv_earned, indiv_unearned, parameters, period
+            individual_earned, individual_unearned, parameters, period
         )
 
-        #   (b) couple combined income
+        # Calculate couple's combined countable income
         couple_countable = _apply_ssi_exclusions(
-            indiv_earned + spouse_earned,
-            indiv_unearned + spouse_unearned,
+            individual_earned + spouse_earned,
+            individual_unearned + spouse_unearned,
             parameters,
             period,
         )
 
-        deemed_if_exceeds = max_(0, couple_countable - alone_countable)
-        return if_leftover_exceeds * deemed_if_exceeds
+        # Deemed amount is the difference (only when deeming applies)
+        deemed_amount = max_(0, couple_countable - alone_countable)
+        return deeming_applies * deemed_amount

policyengine_us/variables/gov/ssa/ssi/ssi.py

@@ -13,4 +13,25 @@ class ssi(Variable):
     def formula(person, period, parameters):
         if parameters(period).gov.ssa.ssi.abolish_ssi:
             return 0
-        return max_(0, person("uncapped_ssi", period))
+
+        uncapped = person("uncapped_ssi", period)
+
+        # Apply floor: SSI cannot be negative
+        benefit = max_(0, uncapped)
+
+        # Apply cap when spousal deeming: cannot exceed individual FBR
+        # per 20 CFR §416.1163
+        # This cap is necessary when spouse's gross income is just above $483:
+        # - Deeming applies (uses couple FBR)
+        # - After exclusions, countable may be low
+        # - Benefit could exceed individual FBR without this cap
+        deeming_applies = person("is_ssi_spousal_deeming_applies", period)
+        p = parameters(period).gov.ssa.ssi.amount
+        individual_max = p.individual * MONTHS_IN_YEAR
+        capped_benefit = min_(benefit, individual_max)
+
+        return where(
+            deeming_applies,
+            capped_benefit,
+            benefit,
+        )

policyengine_us/variables/gov/ssa/ssi/ssi_amount_if_eligible.py

@@ -12,12 +12,34 @@ class ssi_amount_if_eligible(Variable):
     def formula(person, period, parameters):
         p = parameters(period).gov.ssa.ssi.amount
         is_dependent = person("is_tax_unit_dependent", period)
+
+        # Three scenarios for adults:
+        # 1. Both spouses eligible (joint claim) → couple rate / 2
+        # 2. One eligible, no deeming → individual rate
+        # 3. One eligible, deeming applies → couple rate (capped in ssi)
+        #
+        # Note: In scenario 3, deeming applies when spouse's GROSS income > $483.
+        # Income exclusions are applied to the COMBINED income afterwards.
+        # This means countable income can be much lower than $483, and the
+        # benefit can exceed individual FBR, requiring the cap in ssi.
+
+        is_joint_claim = person("ssi_claim_is_joint", period)
+        deeming_applies = person("is_ssi_spousal_deeming_applies", period)
+
+        # Determine FBR to use based on scenario
+        individual_or_deeming_amount = where(
+            deeming_applies,
+            p.couple,  # Scenario 3: Deeming applies - use couple rate!
+            p.individual,  # Scenario 2: No deeming
+        )
+
         head_or_spouse_amount = where(
-            person("ssi_claim_is_joint", period),
-            p.couple / 2,
-            p.individual,
+            is_joint_claim,
+            p.couple / 2,  # Scenario 1: Both eligible
+            individual_or_deeming_amount,
         )
-        # Adults amount is based on whether it is a joint claim
+
+        # Adults amount is based on scenario (see above)
         # Dependents always use individual amount.
         ssi_per_month = where(
             is_dependent, p.individual, head_or_spouse_amount

policyengine_us/variables/gov/states/ca/cpuc/care/eligibility/ca_care_categorically_eligible.py

@@ -17,6 +17,6 @@ class ca_care_categorically_eligible(Variable):
         tribal_lifeline_programs = add(household, period, p.tribal_categorical)
         return np.where(
             is_on_tribal_land,
-            np.any(tribal_lifeline_programs),
-            np.any(non_tribal_lifeline_programs),
+            tribal_lifeline_programs > 0,
+            non_tribal_lifeline_programs > 0,
         )

policyengine_us/variables/gov/states/co/tax/income/credits/care_worker/co_care_worker_credit.py

@@ -0,0 +1,26 @@
+from policyengine_us.model_api import *
+
+
+class co_care_worker_credit(Variable):
+    value_type = float
+    entity = TaxUnit
+    label = "Colorado Care Worker Tax Credit"
+    defined_for = "co_care_worker_credit_eligible"
+    unit = USD
+    definition_period = YEAR
+    reference = "https://law.justia.com/codes/colorado/title-39/specific-taxes/income-tax/article-22/part-5/section-39-22-566/"
+
+    def formula(tax_unit, period, parameters):
+        p = parameters(period).gov.states.co.tax.income.credits.care_worker
+        joint = tax_unit("tax_unit_is_joint", period)
+        eligible_care_worker_count = tax_unit(
+            "co_care_worker_credit_eligible_care_worker_count", period
+        )
+        joint_with_two_eligible_care_workers = joint & (
+            eligible_care_worker_count == 2
+        )
+        return where(
+            joint_with_two_eligible_care_workers,
+            p.amount.joint,
+            p.amount.single,
+        )

policyengine_us/variables/gov/states/co/tax/income/credits/care_worker/co_care_worker_credit_eligible.py

@@ -0,0 +1,27 @@
+from policyengine_us.model_api import *
+
+
+class co_care_worker_credit_eligible(Variable):
+    value_type = float
+    entity = TaxUnit
+    label = "Colorado Care Worker Tax Credit eligible"
+    defined_for = StateCode.CO
+    unit = USD
+    definition_period = YEAR
+    reference = "https://law.justia.com/codes/colorado/title-39/specific-taxes/income-tax/article-22/part-5/section-39-22-566/"
+
+    def formula(tax_unit, period, parameters):
+        p = parameters(period).gov.states.co.tax.income.credits.care_worker
+        joint = tax_unit("tax_unit_is_joint", period)
+        agi = tax_unit("adjusted_gross_income", period)
+        income_limit = where(
+            joint, p.income_limit.joint, p.income_limit.single
+        )
+        income_eligible = agi <= income_limit
+        eligible_care_worker_present = (
+            tax_unit(
+                "co_care_worker_credit_eligible_care_worker_count", period
+            )
+            > 0
+        )
+        return income_eligible & eligible_care_worker_present

policyengine_us/variables/gov/states/co/tax/income/credits/care_worker/co_care_worker_credit_eligible_care_worker.py

@@ -0,0 +1,10 @@
+from policyengine_us.model_api import *
+
+
+class co_care_worker_credit_eligible_care_worker(Variable):
+    value_type = bool
+    entity = Person
+    label = "Eligible Care Worker for the Colorado Care Worker Tax Credit"
+    defined_for = StateCode.CO
+    definition_period = YEAR
+    reference = "https://law.justia.com/codes/colorado/title-39/specific-taxes/income-tax/article-22/part-5/section-39-22-566/"