Abowd 2022 - Privacy-Loss Budget Allocation and Tuning

Summary

TDA was tuned on 2010 data (never on the 2020 data it protects) through demonstration data products. The final pre-production release used ( persons, housing units); production used ( persons, housing units). Production shifted person budget away from blocks (519/1024 down to 165/4099) onto custom block groups and tracts, enabled by spine optimization. Budgets were found by level-by-level binary search against the -percentage-point redistricting criterion.

Overview

Section 8 explains how the global budget and its allocation across levels and queries were set (iteratively, because post-processing errors are data-dependent and have no closed form), and records the numerical allocations for the April 2021 demonstration data and the production run.

Main Content

Why tune on 2010 data

It is not possible to tune a formally private system by directly using the confidential data it protects: the tuning parameters (budget allocations, query strategies) would leak information. So TDA was not tuned on the 2020 Census Edited File. Starting October 2019 the Bureau released demonstration data products from 2010 data (privacy-protected microdata files, PPMF, produced from the Microdata Detail File, MDF), each representing the then-current DAS. Quality assurance of the 2020 MDF is “a difficult conceptual task”; it received human review for unusual or unacceptable errors compared to the CEF, which does not strictly adhere to the tenets of formal privacy (no such errors were flagged for production), highlighting the need for formally private QA.

Demonstration product history

Releases and budgets

ReleaseMechanismGlobal budget
October 2019, May 2020discrete geometric, pure DP () (-allocation: 4 persons, 2 housing units)
September 2020, November 2020discrete Gaussian, zCDP
April 2021 (final pre-production)discrete Gaussian, zCDP (1.05 persons, 0.045 housing units; ) and (0.185 persons, 0.0035 housing units; )
Production (June 2021)discrete Gaussian, zCDP (2.56 persons, 0.07 housing units)

When is re-expressed as the conversion is evaluated at ; a more accurate zCDP conversion sums person and housing-unit first, then converts.

The tuning procedure (April 2021)

Experiments from December 2020 to early April 2021 estimated the minimum global budget, and its allocation to queries, needed for fitness-for-use in redistricting (Wright and Irimata 2021). Criterion: for at least 95% of geographic units, the largest of the selected CENRACE x HISPANIC populations (the major groups in OMB Statistical Policy Directive 15) as a percentage of total population must not change by more than percentage points relative to its enumerated percentage. Because of the cost of many independent runs, geographic variability was used to evaluate the 95% criterion (run-to-run consistency was observed to be high).

TDA’s error distributions are data-dependent in a way that cannot be written in closed form and depends on confidential values, so fitting was iterative rather than via ex ante zCDP properties. Six candidate query strategies were built from two dimensions: (a) measurements closely following the publication specification for redistricting data versus collapsing/excluding marginals to exploit dependence among queries; (b) per-query budgets set (i) proportional to the square of query cardinality, (ii) inversely proportional to it, or (iii) uniform. Passes were ordered heuristically with large-budget queries isolated in their own passes; strategies with almost all budget on DETAILED (proportional to squared cardinality) optimized all queries in one pass. Inverse-square strategies were quickly ruled out as noncompetitive. For the other four, the search went level by level from the top:

  1. By binary search, for each level find the approximate on that level’s queries that satisfies the criterion with all other levels at infinite budget.
  2. Set US and state to their single-level fitted , others infinite; scale US and state up proportionally until the criterion holds.
  3. Fix the county budget at its single-level fit, and scale US, state and county up together until the criterion holds.
  4. Continue down the hierarchy until every level has a .

This led to the initial identification of total . Because the criterion involves off-spine entities, every experimental run processed data down to the block level, which is why tuning took several months. The April 2021 product used the L2 and L1 multi-pass components: the US level used one pass for both phases with all queries in Table 4; state, county, tract and block group levels used two passes each (first pass: the TOTAL query as both the noisy estimate and constraint; second: the remaining queries); the block level used one pass with all queries.

Table 3: per-level proportions, April 2021

Geographic levelPersonsHousing units
US51/10241/1024
State153/10241/1024
County78/102418/1024
Tract51/102475/1024
Custom block group172/1024906/1024
Block519/102423/1024

For persons, block-level queries received about half the budget; for housing units, block-group queries received the predominant share (906/1024).

Table 4: per-query proportions within level, persons, April 2021

Query (cells)USStateCountyTractCBGBlock
TOTAL (1)n/a678/1024 (*)342/10241/1024572/10241/1024
CENRACE (63)2/10241/10241/10242/10241/10242/1024
HISPANIC (2)1/10241/10241/10241/10241/10241/1024
VOTINGAGE (2)1/10241/10241/10241/10241/10241/1024
HHINSTLEVELS (3)1/10241/10241/10241/10241/10241/1024
HHGQ (8)1/10241/10241/10241/10241/10241/1024
HISPANICxCENRACE (126)5/10242/10243/10245/10243/10245/1024
VOTINGAGExCENRACE (126)5/10242/10243/10245/10243/10245/1024
VOTINGAGExHISPANIC (4)1/10241/10241/10241/10241/10241/1024
VOTINGAGExHISPANICxCENRACE (252)17/10246/102411/102417/10248/102417/1024
HHGQxVOTINGAGExHISPANICxCENRACE = DETAILED (2,016)990/1024330/1024659/1024989/1024432/1024989/1024

(*) The TOTAL query is held invariant at the state level; the allocation shown there is for the total population of all AIAN tribal areas in the state and of the remainder of the state (36 states with AIAN areas); the US-level TOTAL cell is blank. Table 5 (housing units): Occupancy Status (2 cells) has proportion 1 at every level. The TOTAL query for counties represents of the person budget.

Production settings (8.2)

Total : persons, housing units, “considerable increases from the larger of the April 2021 settings”. The key change for persons was a shift of allocation away from the block level onto custom block group and tract level, which allowed location-confidentiality protection (controlled mainly by the block allocation) to improve while meeting even tighter accuracy targets for HISPANIC x CENRACE at tract and block-group level. This relied on the efficiency gain from spine optimization. Estimation strategy and pass orderings were the same as April 2021.

Table 6: per-level proportions, production

Geographic levelPersonsHousing units
US104/40991/205
State1440/40991/205
County447/40997/82
Tract687/4099364/1025
Custom block group1256/40991759/4100
Block165/409999/820

Table 7: per-query proportions within level, persons, production

Query (cells)USStateCountyTractCBGBlock
TOTAL (1)n/a3773/40973126/40971567/41021705/40995/4097
CENRACE (63)52/40976/409710/40974/20513/40999/4097
HISPANIC (2)26/40976/409710/40975/41023/40995/4097
VOTINGAGE (2)26/40976/409710/40975/41023/40995/4097
HHINSTLEVELS (3)26/40976/409710/40975/41023/40995/4097
HHGQ (8)26/40976/409710/40975/41023/40995/4097
HISPANICxCENRACE (126)130/409712/409728/40971933/41021055/409921/4097
VOTINGAGExCENRACE (126)130/409712/409728/409710/20519/409921/4097
VOTINGAGExHISPANIC (4)26/40976/409710/40975/41023/40995/4097
VOTINGAGExHISPANICxCENRACE (252)26/2412/241101/409767/410224/409971/4097
DETAILED (2,016)189/241230/4097754/4097241/20511288/40993945/4097

Table 8 (housing units): Occupancy Status (2 cells) has proportion 1 at every level. Normalization check (this note): the printed within-level proportions are only approximately normalized. In Table 7 the US column sums to exactly 1 while the other columns sum to between 0.9976 and 0.9993; in Table 4 the US column sums to 1024/1024 and the other columns to 1023/1024 as printed, consistent with rounding of the published fractions. The per-level proportions in Tables 3 and 6 sum to exactly 1 for both persons and housing units.

Examples

From shares to a concrete and

Setup: production persons, tract-level HISPANICxCENRACE (126 cells), group share and within-level share , persons budget 2.56. Solution: , so per cell (noise standard deviation about 2.2 people). Compare April 2021 tract HISPANICxCENRACE: , i.e. . Interpretation: the production shift of budget to HISPANICxCENRACE at the tract level cuts per-cell variance roughly 800-fold, paid for by the higher total budget and the drop in block-level share. Recomputed by this note.

person_apr = dict(zip("US state county tract CBG block".split(), [51,153,78,51,172,519]))
print(sum(person_apr.values()))   # 1024
 
# production tract HISPANICxCENRACE
rho_g = 2.56 * (687/4099) * (1933/4102)
print(rho_g, 1/rho_g, (1/rho_g) ** 0.5)       # ~0.2022, ~4.95, ~2.2
# April 2021 tract HISPANICxCENRACE
rho_a = 1.05 * (51/1024) * (5/1024)
print(rho_a, 1/rho_a)                          # ~2.55e-4, ~3.9e3

Connections

See Also