A layered perturbation method for the protection of tabular outputs
Section 4. Empirical investigation
We applied the LPM and EZS methods to personal data from a taxation file. Two variables were used: income and (to increase skewness). Cells of between 15 and 148 units were generated by combining age groups within postal code, sex and marital status. Different levels of noise from a split triangular distribution were tried. Results presented are those with Following the Risk-Utility Framework (Duncan, Keller-McNulty and Stokes 2001) the impacts of the methods on data accuracy and on risk were examined.
Table 4.1 shows the impact of the LPM on the quality of cell totals by cell size range. The LPM was applied 500 times in each cell. For each cell size range the table gives the number of cells, their average coefficient of variation (CV) after perturbation, and the percentage of times that the perturbed total was within 2%, 5%, 8% and 12% of the original cell total. For this study we assumed that cells that failed a percent sensitivity rule with would be suppressed so they were not included in the results. There were more such cells with variable (which may resemble business data more). As expected, the impact of the perturbation was higher for smaller cells, and for variable All cells perturbed by more than 8% were near-sensitive and would have been suppressed with
| Cell size |
Variable = Income (x) | Variable = Income2 (y) | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Num. cells |
Avg. CV |
% times relative distance ≤ | Num. cells |
Avg. CV |
% times relative distance ≤ | |||||||
| 2% | 5% | 8% | 12% | 2% | 5% | 8% | 12% | |||||
| 15 – 18 | 1,822 | 2.37 | 58.5 | 95.1 | 99.5 | 100.0 | 1,777 | 4.09 | 34.5 | 72.0 | 92.4 | 99.6 |
| 19 – 25 | 2,230 | 2.03 | 66.2 | 97.2 | 99.7 | 100.0 | 2,185 | 3.71 | 38.1 | 77.1 | 94.4 | 99.7 |
| 26 – 40 | 1,920 | 1.57 | 78.2 | 99.1 | 99.9 | 100.0 | 1,899 | 3.24 | 44.2 | 82.8 | 96.0 | 99.8 |
| 41 – 148 | 1,312 | 1.05 | 92.1 | 99.5 | 99.9 | 100.0 | 1,301 | 2.53 | 57.1 | 90.0 | 97.7 | 99.9 |
| All | 7,284 | 1.82 | 72.1 | 97.6 | 99.7 | 100.0 | 7,162 | 3.47 | 42.3 | 79.7 | 94.9 | 99.7 |
| Note: values of 100.0 represent values above 99.95 that were rounded to 100. | ||||||||||||
Table 4.2 gives the impact of the EZS multiplicative noise, for the same on the cell totals. Results for income are fairly similar, while results for are noticeably better. Similar results were obtained when a value for near 0.014 was used (LPM was slightly better with EZS slightly better with
| Cell size |
Variable = Income (x) | Variable = Income2 (y) | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Num. cells |
Avg. CV |
% times relative distance ≤ | Num. cells |
Avg. CV |
% times relative distance ≤ | |||||||
| 2% | 5% | 8% | 12% | 2% | 5% | 8% | 12% | |||||
| 15 – 18 | 1,822 | 2.33 | 58.7 | 97.1 | 100.0 | 100.0 | 1,777 | 3.19 | 41.2 | 86.4 | 99.8 | 100.0 |
| 19 – 25 | 2,230 | 2.08 | 64.5 | 98.5 | 100.0 | 100.0 | 2,185 | 2.93 | 45.2 | 90.0 | 99.9 | 100.0 |
| 26 – 40 | 1,920 | 1.74 | 73.9 | 99.6 | 100.0 | 100.0 | 1,899 | 2.59 | 51.4 | 93.8 | 99.9 | 100.0 |
| 41 – 148 | 1,312 | 1.30 | 86.9 | 99.9 | 99.9 | 100.0 | 1,301 | 2.09 | 63.4 | 97.1 | 100.0 | 100.0 |
| All | 7,824 | 1.91 | 69.6 | 98.7 | 99.9 | 100.0 | 7,162 | 2.76 | 49.2 | 91.4 | 99.9 | 100.0 |
| Note: values of 100.0 in the 8% columns represent values above 99.95 that were rounded to 100. | ||||||||||||
We next examined the amount of protection offered to the largest units in each cell. For each cell, an estimate for unit was obtained by differencing perturbed cell totals with and without the unit. Relative differences were calculated and incorporated in a score equal to where if if and otherwise. Table 4.3 shows the quartiles of and the scores for variables and for the largest twelve units in each cell with LPM, and for the largest unit with EZS (EZS offers the same level of protection to all units).
With the LPM the largest three units tend to be protected the most, as expected. Patterns for variables and are different. If one looks at the quartiles of for variable the level of protection gradually declines until unit 10 and increases afterwards. Since the are the same for results should keep improving after the 10th largest unit. The scores give a similar story. For variable the descent is not as regular, with unit 5 being protected the least (unit 10 if one looks at Q1 only). The weaker protection around units 5 and 10 is predicted by the formulas for whose basic form changes around those two units. Unit 10 is vulnerable to repeated targeted attacks the most, where an attack consists of obtaining an estimate from totals for units 1 to 10, and for units 1 to 9, with some set of smaller units (e.g., obtain from totals excluding unit and excluding units and 10, for Averaging the if there are enough of them, may give good estimates of Such attacks require carefully set up tabulation requests, which a semi-controlled custom tabulation environment could discourage.
| Cells | Q1 | Med | Q3 | Score (%) | |
|---|---|---|---|---|---|
| Variable = Income (x) | |||||
| Unit 1 | 7,962 | 7.9 | 15.7 | 26.6 | 3,196 (40) |
| Unit 2 | 7,962 | 8.6 | 17.5 | 29.3 | 2,895 (36) |
| Unit 3 | 7,962 | 8.1 | 16.9 | 28.7 | 3,021 (38) |
| Unit 4 | 7,962 | 7.2 | 15.5 | 26.2 | 3,314 (42) |
| Unit 5 | 7,962 | 6.4 | 13.9 | 23.8 | 3,647 (46) |
| Unit 6 | 7,962 | 6.4 | 13.9 | 23.3 | 3,614 (45) |
| Unit 7 | 7,962 | 6.2 | 13.3 | 22.4 | 3,765 (47) |
| Unit 8 | 7,962 | 6.3 | 13.4 | 22.3 | 3,731 (47) |
| Unit 9 | 7,962 | 5.1 | 11.5 | 19.9 | 4,267 (54) |
| Unit 10 | 7,962 | 3.3 | 10.7 | 20.9 | 4,373 (55) |
| Unit 11 | 7,962 | 3.8 | 11.8 | 22.4 | 4,121 (52) |
| Unit 12 | 7,962 | 3.8 | 12.2 | 24.7 | 4,031 (51) |
| U1/EZS | 7,962 | 6.7 | 7.5 | 8.4 | 7,941 (100) |
| Variable = Income2 (y) | |||||
| Unit 1 | 7,823 | 7.6 | 14.4 | 23.2 | 3,365 (43) |
| Unit 2 | 7,782 | 7.2 | 15.0 | 25.2 | 3,311 (43) |
| Unit 3 | 7,782 | 6.6 | 14.1 | 24.2 | 3,522 (45) |
| Unit 4 | 7,799 | 6.1 | 13.3 | 22.5 | 3,726 (48) |
| Unit 5 | 7,808 | 5.5 | 11.9 | 20.5 | 4,052 (52) |
| Unit 6 | 7,811 | 6.0 | 12.6 | 21.6 | 3,885 (50) |
| Unit 7 | 7,814 | 6.0 | 12.6 | 22.2 | 3,868 (50) |
| Unit 8 | 7,818 | 6.5 | 13.8 | 23.7 | 3,581 (46) |
| Unit 9 | 7,818 | 5.7 | 13.0 | 24.2 | 3,750 (48) |
| Unit 10 | 7,818 | 4.4 | 13.5 | 27.4 | 3,704 (47) |
| Unit 11 | 7,818 | 4.8 | 15.7 | 32.1 | 3,422 (44) |
| Unit 12 | 7,820 | 5.8 | 17.9 | 37.9 | 3,110 (40) |
| U1/EZS | 7,823 | 6.7 | 7.5 | 8.5 | 7,803 (100) |
In contrast, results for EZS show that the level of protection offered to unit 1 (or for any unit for that matter) is fairly constant, and it is generally much poorer than that with the LPM. The score for EZS is almost 100%, a very poor outcome. But EZS was designed to offer protection for totals, not to protect from differencing. If protection from differencing is required then the level of noise would have to be set much higher to protect values at levels comparable to the LPM. But with EZS units around unit 10 would not be more vulnerable to repeated targeted attacks.
To investigate the roles of and we generated random values from a uniform distribution, but created an outlier in each cell by setting the value of as the highest value that would not make the cell sensitive, i.e., for set The LPM was used with set to 1, and and either calculated as suggested above or set to 1. For our generated data the calculated value of never left 1. Table 4.4 shows that factor is useful because when it is set to 1 the level of protection for the outlier is not high enough when
| Q1 | Med. | Q3 | Score | |
|---|---|---|---|---|
| Standard LPM (K ≥ 1) | 11.1 | 12.6 | 14.2 | 472 |
| LPM with K = L = M = 1 | 6.7 | 7.5 | 8.6 | 996 |
- Date modified: