530 lines
69 KiB
Plaintext
530 lines
69 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"# Sensitivity"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"```{admonition} Learning Objectives\n",
|
|
"After reading this chapter, you will be able to:\n",
|
|
"- Define sensitivity\n",
|
|
"- Find the sensitivity of counting queries\n",
|
|
"- Find the sensitivity of summation queries\n",
|
|
"- Decompose average queries into counting and summation queries\n",
|
|
"- Use clipping to bound the sensitivity of summation queries\n",
|
|
"```\n",
|
|
"\n",
|
|
"As we mentioned when we discussed the Laplace mechanism, the amount of noise necessary to ensure differential privacy for a given query depends on the *sensitivity* of the query. Roughly speaking, the sensitivity of a function reflects the amount the function's output will change when its input changes. Recall that the Laplace mechanism defines a mechanism $F(x)$ as follows:\n",
|
|
"\n",
|
|
"\\begin{align}\n",
|
|
"F(x) = f(x) + \\textsf{Lap}\\left(\\frac{s}{\\epsilon}\\right)\n",
|
|
"\\end{align}\n",
|
|
"\n",
|
|
"where $f(x)$ is a deterministic function (the query), $\\epsilon$ is the privacy parameter, and $s$ is the sensitivity of $f$.\n",
|
|
"\n",
|
|
"For a function $f : \\mathcal{D} \\rightarrow \\mathbb{R}$ mapping datasets ($\\mathcal{D}$) to real numbers, the *global sensitivity* of $f$ is defined as follows:\n",
|
|
"\n",
|
|
"\\begin{align}\n",
|
|
"GS(f) = \\max_{x, x': d(x,x') <= 1} |f(x) - f(x')|\n",
|
|
"\\end{align}\n",
|
|
"\n",
|
|
"Here, $d(x, x')$ represents the *distance* between two datasets $x$ and $x'$, and we say that two datasets are *neighbors* if their distance is 1 or less. How this distance is defined has a huge effect on the definition of privacy we obtain, and we'll discuss the distance metric on datasets in detail later on.\n",
|
|
"\n",
|
|
"The definition of global sensitivity says that for *any two* neighboring datasets $x$ and $x'$, the difference between $f(x)$ and $f(x')$ is at most $GS(f)$. This measure of sensitivity is called \"global\" because it is independent of the actual dataset being queried (it holds for *any* choice of neighboring $x$ and $x'$). Another measure of sensitivity, called *local sensitivity*, fixes one of the datasets to be the one being queried; we will consider this measure in a later section. For now, when we say \"sensitivity,\" we mean global sensitivity."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"## Distance\n",
|
|
"\n",
|
|
"The distance metric $d(x,x')$ described earlier can be defined in many different ways. Intuitively, the distance between two datasets should be equal to 1 (i.e. the datasets are neighbors) if they differ in the data of exactly one individual. This idea is easy to formalize in some contexts (e.g. in the US Census, each individual submits a single response containing their data) but extremely challenging in others (e.g. location trajectories, social networks, and time-series data).\n",
|
|
"\n",
|
|
"A common formal definition for datasets containing rows is to consider the number of rows which differ between the two. When each individual's data is contained in a single row, this definition often makes sense. Formally, this definition of distance is encoded as a *symmetric difference* between the two datasets:\n",
|
|
"\n",
|
|
"\\begin{align}\n",
|
|
"d(x, x') = | x - x' \\cup x' - x |\n",
|
|
"\\end{align}\n",
|
|
"\n",
|
|
"This particular definition has several interesting and important implications:\n",
|
|
"- If $x'$ is constructed from $x$ by *adding one row*, then $d(x,x') = 1$\n",
|
|
"- If $x'$ is constructed from $x$ by *removing one row*, then $d(x,x') = 1$\n",
|
|
"- If $x'$ is constructed from $x$ by *modifying one row*, then $d(x,x') = 2$\n",
|
|
"\n",
|
|
"In other words, adding or removing a row results in a neighboring dataset; *modifying* a row results in a dataset at distance *2*. \n",
|
|
"\n",
|
|
"This particular definition of distance results in what is typically called *unbounded differential privacy*. Many other definitions are possible, including one called *bounded differential privacy* in which modifying a single row in a dataset *does* result in a neighboring dataset. \n",
|
|
"\n",
|
|
"For now, we will stick to the formal definition defined above in terms of symmetric difference. We will discuss alternative definitions in a later section."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"## Calculating Sensitivity\n",
|
|
"\n",
|
|
"How do we determine the sensitivity of a particular function of interest? For some simple functions on real numbers, the answer is obvious.\n",
|
|
"\n",
|
|
"- The global sensitivity of $f(x) = x$ is 1, since changing $x$ by 1 changes $f(x)$ by 1\n",
|
|
"- The global sensitivity of $f(x) = x+x$ is 2, since changing $x$ by 1 changes $f(x)$ by 2\n",
|
|
"- The global sensitivity of $f(x) = 5*x$ is 5, since changing $x$ by 1 changes $f(x)$ by 5\n",
|
|
"- The global sensitivity of $f(x) = x*x$ is unbounded, since the change in $f(x)$ depends on the value of $x$\n",
|
|
"\n",
|
|
"For functions that map datasets to real numbers, we can perform a similar analysis. We will consider the functions which represent common aggregate database queries: counts, sums, and averages."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 2,
|
|
"metadata": {
|
|
"tags": [
|
|
"remove-cell"
|
|
]
|
|
},
|
|
"outputs": [],
|
|
"source": [
|
|
"import pandas as pd\n",
|
|
"import numpy as np\n",
|
|
"import matplotlib.pyplot as plt\n",
|
|
"plt.style.use('seaborn-whitegrid')\n",
|
|
"\n",
|
|
"adult = pd.read_csv(\"adult_with_pii.csv\")"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"### Counting Queries\n",
|
|
"\n",
|
|
"Counting queries (`COUNT` in SQL) count the number of rows in the dataset which satisfy a specific property. As a rule of thumb, **counting queries always have a sensitivity of 1**. This is because adding a row to the dataset can increase the output of the query by at most 1: either the new row has the desired property, and the count increases by 1, or it does not, and the count stays the same (the count may correspondingly decrease when a row is removed)."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"How many people are in the dataset?\"** (sensitivity: 1 - counting rows where the property = True)"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 3,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"32561"
|
|
]
|
|
},
|
|
"execution_count": 3,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult.shape[0]"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"How many people have an educational status above 10?\"** (sensitivity: 1 - counting rows with a property)"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 4,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"10516"
|
|
]
|
|
},
|
|
"execution_count": 4,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Education-Num'] > 10].shape[0]"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"How many people have an educational status equal to or below 10?\"** (sensitivity: 1 - counting rows with a property)"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 5,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"22045"
|
|
]
|
|
},
|
|
"execution_count": 5,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Education-Num'] <= 10].shape[0]"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"How many people are named Joe Near?\"** (sensitivity: 1 - counting rows with a property)"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 6,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"0"
|
|
]
|
|
},
|
|
"execution_count": 6,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Name'] == 'Joe Near'].shape[0]"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"### Summation Queries\n",
|
|
"\n",
|
|
"Summation queries (`SUM` in SQL) sum up the *attribute values* of dataset rows."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"What is the sum of the ages of people with an educational status above 10?\"**"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 7,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"422876"
|
|
]
|
|
},
|
|
"execution_count": 7,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Education-Num'] > 10]['Age'].sum()"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"Sensitivity for these queries is not as simple as it is for counting queries. Adding a new row to the dataset will increase the result of our example query by the *age of the new person*. That means the sensitivity of the query depends on the *contents* of the row we add.\n",
|
|
"\n",
|
|
"We'd like to come up with a concrete number to represent the sensitivity of the query. Unfortunately, no number really exists. We could claim, for example, that the sensitivity is 125 - but it may turn out that the row we add to the database corresponds to a person who is over 125 years old, which would violate our claim. For any number we come up with, it's possible for the added row to violate our claim.\n",
|
|
"\n",
|
|
"You might (rightly) be skeptical of this point. Say we claim the sensitivity is 1000 - it's very unlikely that we'll find a person who is 1000 years old to violate this claim. In this specific domain - ages - there's a very reasonable upper bound on how old someone can be. The oldest person ever lived to be [122 years old](https://en.wikipedia.org/wiki/List_of_the_verified_oldest_people), so an upper bound of 125 seems reasonable.\n",
|
|
"\n",
|
|
"But this is not a *proof* that nobody will ever live to be 126. And in other domains (e.g. income), it can be much harder to come up with a reasonable upper bound.\n",
|
|
"\n",
|
|
"As a rule of thumb, summation queries have **unbounded sensitivity** when no lower and upper bounds exist on the value of the attribute being summed. When lower and upper bounds do exist, the sensitivity of a summation query is equal to the **difference between them**. In the next section, we will see a technique called *clipping* for enforcing bounds when none exist, so that summation queries with unbounded sensitivity can be converted into queries with bounded sensitivity."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"### Average Queries\n",
|
|
"\n",
|
|
"Average queries (`AVG` in SQL) calculate the mean of attribute values in a particular column."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"**Example: \"What is the average age of people with an educational status above 10?\"**"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 8,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"40.21262837580829"
|
|
]
|
|
},
|
|
"execution_count": 8,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Education-Num'] > 10]['Age'].mean()"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"The easiest way to answer an average query with differential privacy is by re-phrasing it as two queries: a summation query divided by a counting query. For the above example:"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 9,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"40.21262837580829"
|
|
]
|
|
},
|
|
"execution_count": 9,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult[adult['Education-Num'] > 10]['Age'].sum() / adult[adult['Education-Num'] > 10]['Age'].shape[0]"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"The sensitivities of both queries can be calculated as described above. Noisy answers for each can be calculated (e.g. using the Laplace mechanism) and the noisy answers can be divided to obtain a differentially private mean. The total privacy cost of both queries can be calculated by sequential composition."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"## Clipping\n",
|
|
"\n",
|
|
"Queries with unbounded sensitivity cannot be directly answered with differential privacy using the Laplace mechanism. Fortunately, we can often transform such queries into equivalent queries with *bounded* sensitivity, via a process called *clipping*.\n",
|
|
"\n",
|
|
"The basic idea beind clipping is to *enforce* upper and lower bounds on attribute values. For example, ages above 125 can be \"clipped\" to exactly 125. After clipping has been performed, we are *guaranteed* that all ages will be 125 or below. As a result, the senstivity of a summation query on clipped data is equal to the difference between the upper and lower bounds used in clipping: $upper - lower$. For example, the following query has a sensitivity of 125:"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 10,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"1256257"
|
|
]
|
|
},
|
|
"execution_count": 10,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
}
|
|
],
|
|
"source": [
|
|
"adult['Age'].clip(lower=0, upper=125).sum()"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"The primary challenge in performing clipping is to determine the upper and lower bounds. For ages, this is simple - nobody can have an age less than 0, and probably nobody will be older than 125. In other domains, as mentioned earlier, it's much more difficult.\n",
|
|
"\n",
|
|
"Furthermore, there is a tradeoff between the amount of information lost in clipping and the amount of noise needed to ensure differential privacy. When the upper and lower clipping bounds are closer together, then the sensitivity is lower, and less noise is needed to ensure differential privacy. However, aggressive clipping often removes a lot of information from the data; this information loss tends to cause a *loss* of accuracy which outweighs the improvement in noise resulting from smaller sensitivity. \n",
|
|
"\n",
|
|
"As a rule of thumb, **try to set the clipping bounds to include 100% of the dataset**, or get as close as possible. This is harder in some domains (e.g. graph queries, which we will study later) than others.\n",
|
|
"\n",
|
|
"It's tempting to determine the clipping bounds by looking at the data. For example, we can look at the histogram of ages in our dataset to determine an appropriate upper bound:"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 11,
|
|
"metadata": {
|
|
"tags": [
|
|
"hide-input"
|
|
]
|
|
},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEDCAYAAADayhiNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3df1xUdaL/8dfMMPPAH5CrDHwXM3/tEqSoGeVVYdt70UTdfTSVPRTS7WHo9V7JrGsqEVleeyyIsWs/vNvdLH9tFle863K1wO1e6G5XpCvsNblL19Wyh4nBIBUFjIxwvn+UZyWiIWyGQd/Px8OHzGc4c97jMLw9nzPnHIthGAYiIiKAta8DiIhI8FApiIiISaUgIiImlYKIiJhUCiIiYlIpiIiIKaSvA1yOysrKvo4gItIv3XTTTV873q9LAbp/Yv5QU1NDXFxcwNbXW/0hZ3/ICMr5XeoPGeHqyPlN/6HW9JGIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiKmfn/wmvQPozIPfM3oe35f76ncuX5fh8iVRFsKIiJi8uuWQlFREVu3biUkJISVK1cSExPDmjVraG9vx+l0smnTJhwOB0VFRezYsQOr1cr8+fOZN28eXq+XzMxMamtrsdls5OTkMGLECH/GFRG56vltS+Hjjz9my5Yt7N69m+eff5433niDZ555hrS0NHbv3s3w4cMpLCykpaWFLVu2sH37dnbt2sXWrVv55JNP2L9/P+Hh4bzyyissXbqU/Px8f0UVEZEv+W1Loby8nKlTpzJ48GAGDx7Mhg0b+Ju/+RvWr18PQHJyMtu3b2f06NHEx8cTFhYGQEJCAlVVVZSXl+NyuQBITEwkOzvbX1GvGl8/ry8i8hd+K4UPP/wQwzB48MEHqa+vZ8WKFbS2tuJwOABwOp243W4aGhoYOnSouVxERESXcZvNhtVqpa2tzVxeRES+e37dp1BXV8dzzz1HbW0tP/vZz7BYLOZ9hmF0+vvScYvF0u34V9XU1Pgh+dfzeDwBXV9v9ZecgXC5/w795d+yP+TsDxlBOf1WCsOGDePGG28kJCSE6667jkGDBmGz2fB4PISGhlJXV0dkZCRRUVGUlZWZy9XX1zNp0iSioqJwu93Exsbi9XoxDAO73d5lPYE873n/P8+6/z8CGmwu9/Xq/6958OgPGeHqyNkn11NITEzk8OHDdHR00NjYSEtLC9OmTaOkpASAgwcPkpSUxMSJEzl27BhNTU00NzdTVVVFQkIC06dPp7i4GIDS0lKmTJnir6giIvIlv20pREVFMWvWLO69915aW1vJzs4mPj6etWvXUlBQQHR0NC6XC7vdzqpVq0hPT8disZCRkUFYWBhz5szh0KFDpKam4nA4yM3N9VfUgAvMDt+rb6tARC6fX/cpLFiwgAULFnQa27ZtW5fvS0lJISUlpdPYxWMTREQkcHREs4iImFQKIiJiUimIiIhJpSAiIiaVgoiImFQKIiJiUimIiIhJpSAiIiaVgoiImK7aazT3/lQTOn2EiFy5tKUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYvLb9RSqq6tZvnw5I0eOBCAmJoYlS5awZs0a2tvbcTqdbNq0CYfDQVFRETt27MBqtTJ//nzmzZuH1+slMzOT2tpabDYbOTk5jBgxwl9xRUQEP5ZCS0sLs2bN4tFHHzXHHnnkEdLS0pg9ezZ5eXkUFhbicrnYsmULhYWF2O12XC4XM2bMoLS0lPDwcPLz83nzzTfJz89n8+bN/oorIiL4cfqoubm5y1hFRQXJyckAJCcnU15eztGjR4mPjycsLIzQ0FASEhKoqqqivLycmTNnApCYmEhlZaW/ooqIyJf8uqVQWVnJkiVLaG1tZcWKFbS2tuJwOABwOp243W4aGhoYOnSouVxERESXcZvNhtVqpa2tzVz+opqaGn89BbkCXO7Ph8fj6Rc/Y/0hZ3/ICMrpt1KIjY0lIyOD5ORk3n//fRYvXsyFCxfM+w3D6PT3peMWi6Xb8a+Ki4vrZUJda/lq0Pufjy/U1NRc9mMEQn/I2R8ywtWR85tmXvw2fTR27Fhzqmj06NFERETQ1NSEx+MBoK6ujsjISKKiomhoaDCXq6+vx+l0EhUVhdvtBsDr9WIYBna73V9xRUQEP5ZCYWEhO3fuBMDtdnPu3DnuvPNOSkpKADh48CBJSUlMnDiRY8eO0dTURHNzM1VVVSQkJDB9+nSKi4sBKC0tZcqUKf6KKiIiX/Lb9NHMmTN5+OGHKSkpoa2tjSeeeIK4uDjWrl1LQUEB0dHRuFwu7HY7q1atIj09HYvFQkZGBmFhYcyZM4dDhw6RmpqKw+EgNzfXX1FFRORLfiuFa665hhdeeKHL+LZt27qMpaSkkJKS0mns4rEJIiISODqiWURETCoFERExqRRERMSkUhAREZNKQURETCoFERExqRRERMSkUhAREZNKQURETCoFERExqRRERMSkUhAREZPPUti3bx979uyhra2N9PR07rrrLnbv3h2IbCIiEmA+S+GVV17hjjvuoLi4mOuvv569e/ea10QQEZEri89SsFqthISEUFJSwk9/+lMAzp8/7/dgIiISeD5LYdy4ccycOROv10tcXBy7du0iOjo6ENlERCTAfF5kJzs7mxUrVnDNNdcAkJycTGpqqt+DiYhI4HVbCosWLcJisXS74MXrL4uIyJWj21JYt24dAHv27MHpdDJlyhQ6OjqoqKigqakpYAFFRCRwui2FH/7whwCcOnWKrKwsc3zSpEksWbLE/8lERCTgfO5TaGpqYteuXdx4441YrVaOHTumLQURkSuUz1J4+umn2blzJ8899xyGYTBmzBg2b94ciGwiIhJgPkvhhRdeIDs7OxBZRESkj/ksBcMwKCgoYMKECdjtdnP8Bz/4gV+DiYhI4PkshePHj3P8+HH2799vjlkslh59JNXj8TB37lwyMjKYOnUqa9asob29HafTyaZNm3A4HBQVFbFjxw6sVivz589n3rx5eL1eMjMzqa2txWazkZOTw4gRIy7vmYqIiE8+S2HXrl00NzfzwQcfYLVaGTVqFKGhoT168F/96lcMGTIEgGeeeYa0tDRmz55NXl4ehYWFuFwutmzZQmFhIXa7HZfLxYwZMygtLSU8PJz8/HzefPNN8vPztR9DRCQAfJ7moqioCJfLxebNm9m0aRO33347v//9730+8MmTJzlx4gQ//vGPAaioqCA5ORn44qjo8vJyjh49Snx8PGFhYYSGhpKQkEBVVRXl5eXMnDkTgMTERCorKy/jKYqISE/53FJ4+eWXKSoqYsCAAQA0NzeTnp5u/tLuzsaNG3nsscfYt28fAK2trTgcDgCcTidut5uGhgaGDh1qLhMREdFl3GazYbVaaWtrM5cXERH/8FkKVqvVLASAQYMGERLyzYvt27ePSZMmddoPcOkpMwzD6PT3peMWi6Xb8a9TU1Pj6ynIVexyfz48Hk+/+BnrDzn7Q0ZQTp+lMHnyZJYtW8bNN9+MYRhUVFSQkJDwjcuUlZVx+vRpysrK+Oijj3A4HAwYMACPx0NoaCh1dXVERkYSFRVFWVmZuVx9fT2TJk0iKioKt9tNbGwsXq8XwzA6ffLpUnFxcd/uGZve6+Vy0p/M3tF3r/Op3LkBW1dNTc1lvBcCoz9khKsj5zdNyfsshdWrV3PkyBGqq6sBWL58OZMnT/7GZS7dKfzss88yfPhw/vjHP1JSUsLtt9/OwYMHSUpKYuLEiWRnZ9PU1ITNZqOqqoqsrCw+//xziouLSUpKorS0lClTpvT0uYqIyGXwWQonTpzg0KFDPPDAAwD84z/+I2FhYea5kXpqxYoVrF27loKCAqKjo3G5XNjtdlatWkV6ejoWi4WMjAzCwsKYM2cOhw4dIjU1FYfDQW5ubu+enYiIfCs+S+Hxxx/noYceMm/PmzeP9evX85vf/KZHK1ixYoX59bZt27rcn5KSQkpKSqexi8cmiIhIYPn8SGp7e3unfQg33HBDlx3BIiJyZfC5pTBhwgQeeOABJk+ebF5PYcKECYHIJiIiAeazFLKysigvL+d///d/sdlsLF261Oenj0REpH/yOX104cIFGhoasFgsLF68mPDwcLxebyCyiYhIgPkshccee4x3332X4uJiAN5++23Wrl3r92AiIhJ4Pkvh7NmzrF692jwJ3sKFC6mvr/d7MBERCTyfpeD1emlqajJPM3Hy5Ena2tr8HkxERALP547mhx56iHvvvZdTp04xe/ZsAJ588km/BxMRkcDzWQoJCQn89re/5dy5c9jtdsLDwwORS0RE+kC300dut5tHH32URYsWkZOTw4ABAwgPD+f//u//SEtLC2RGEREJkG63FLKyspg7dy733XcfJSUlrFu3jqFDh1JRUcGaNWsCmVFERAKk2y2F1tZWXC4XY8eOZfny5VRWVjJy5Ej+9V//lenTpwcyo4iIBEi3WwpWa+e+uPbaa7nnnnv8HkhERPpOt6XQ0dGBx+PpdJW0S29fejU2ERG5MnRbCrW1tcydO7fTGVEv3rZYLPz7v/97QAKKiEjgdFsK//Ef/xHIHCIiEgR8HtEsIiJXD5WCiIiYui2FkpISAF5//fWAhRERkb7V7T6FX/ziF9TV1fHyyy/T2NjY5X59PFVE5MrTbSls2LCB//7v/8br9fLxxx8HMpOIiPSRbkvhlltu4ZZbbmHmzJlce+21fPDBB1gsFkaNGmVeW0FERK4sPs+SWlNTQ0ZGBqNHj6a9vZ0PP/yQhx9+mJkzZwYin4iIBJDPUti9ezdFRUXmEczNzc2kp6erFERErkA+P5JqtVo7ndJi0KBBhIT47BIREemHfP52nzx5MsuWLePmm2/GMAzefvttEhISfD5wa2srmZmZnDt3jvPnz7N8+XJiY2NZs2YN7e3tOJ1ONm3ahMPhoKioiB07dmC1Wpk/fz7z5s3D6/WSmZlJbW0tNpuNnJwcRowY8Z08aRER+Xo+S2H16tUcOXKE6upqAP7u7/6Om266yecDl5aWMn78eJYuXcqZM2e47777mDx5MmlpacyePZu8vDwKCwtxuVxs2bKFwsJC7HY7LpeLGTNmUFpaSnh4OPn5+bz55pvk5+ezefPmy3/GIiLSrR7NAyUkJPRo6+BSc+bMMb8+e/YsUVFRVFRUsH79egCSk5PZvn07o0ePJj4+nrCwMHNdVVVVlJeX43K5AEhMTCQ7O/tbrV9ERL49v+8cWLBgAR999BHPP/88ixcvxuFwAOB0OnG73TQ0NDB06FDz+yMiIrqM22w2rFYrbW1t5vIX1dTU+PspiPRKIH82PR5P0L8X+kNGUE6fpfCnP/2JG264odcrePXVV6mpqWH16tVYLBZz/NLrNFzq4qm5uxv/qri4uF4me6+Xy4n0TO9/Nr+9mpqagK6vN/pDRrg6clZWVnZ7n89PH+Xm5nLhwoVvvdLq6mrOnj0LfPHmaG9vZ8CAAXg8HgDq6uqIjIwkKiqKhoYGc7n6+nqcTidRUVG43W4AvF4vhmFgt9u/dQ4REek5n6UwcOBAbrvtNpYvX87KlSvNP74cOXKEl156CYCGhgZaWlqYNm2aeaK9gwcPkpSUxMSJEzl27BhNTU00NzdTVVVFQkIC06dPp7i4GPhip/WUKVMu53mKiEgP+Jw+uu+++3r1wAsWLODRRx8lLS0Nj8fDunXrGD9+PGvXrqWgoIDo6GhcLhd2u51Vq1aRnp6OxWIhIyODsLAw5syZw6FDh0hNTcXhcJCbm9urHCIi0nM9Ok6huLiYuro60tPTOX78OKNHj/b5wKGhoeTn53cZ37ZtW5exlJQUUlJSOo1dPDZBREQCx+f00WOPPUZNTY05lfP222+zdu1avwcTEZHA81kKZ8+eZfXq1eaZURcuXEh9fb3fg4mISOD5LAWv10tTU5P5cdCTJ0/S1tbm92AiIhJ4PvcpPPTQQ9x7772cOnWKlJQULBYLTz75ZCCyiYhIgPkshYSEBH77299y7tw5bDYbQ4YMCUQuERHpAz5LYe/evTz77LMMHjwYgJaWFv7hH/6Bn/zkJ34PJyIigeWzFHbu3Mm+ffvMLYTGxkYWL16sUhARuQL53NE8cuRIwsPDzdvf+973uO666/waSkRE+ka3WwobN240P3Hkcrm46aabsFgs/M///E+PDl4TudqNyjwQ4DV+cZLHU7lzA7xeuZJ0WwoxMTEA/PCHP+w0Hh8f/7VnKxURkf6v21K44447APj888+pqKjgs88+C1goERHpGz53NC9atIiYmJhOF8LRloKIyJXJZykMGTKEjRs3BiKLiIj0MZ+lcOedd7Jhwwbi4uIICfnLt1+8frKIiFw5fJbCCy+8QExMDCdPnjTHNH0kInJl8lkKQ4cO5amnngpEFhER6WM+S2HcuHH88pe/ZMKECZ2mj2699Va/BhMRkcDzWQqNjY0AvPHGG53GVQoiIlcen6WwYsWKQOQQEZEg0KNSuLhj2ev1cvr0acaNG8euXbv8Hk5ERAKrR6fOvpTb7ebpp5/2WyAREek7Ps+S+lVOp5N3333XH1lERKSP+dxSuOuuu8zpI8MwaGxs5K/+6q/8HkxERALPZyk888wz5tcWi4XBgwd3ur6CiIhcObothX379n3jgj05zUVeXh6VlZVcuHCBZcuWER8fz5o1a2hvb8fpdLJp0yYcDgdFRUXs2LEDq9XK/PnzmTdvHl6vl8zMTGpra7HZbOTk5DBixIhv/wxFRKTHui0FwzC6jF24cIFXX32Vuro6n6Vw+PBh/vznP1NQUMDHH3/MHXfcwdSpU0lLS2P27Nnk5eVRWFiIy+Viy5YtFBYWYrfbcblczJgxg9LSUsLDw8nPz+fNN98kPz+fzZs3X/4zFhGRbvm8nsJFr732Gjt27GDGjBncd999Ph/45ptvZsKECQBcc801tLa2UlFRwfr16wFITk5m+/btjB49mvj4eMLCwgBISEigqqqK8vJys3gSExPJzs7u3TMUEZEe87lP4fDhw2zevJlx48bx4osvMmzYsB49sM1mY+DAgQDs2bOHH/3oR7z11ls4HA7gi08xud1uGhoaOl2rISIiosu4zWbDarXS1tZmLi8iIt+9bkvh+PHj5OfnM3DgQDZt2tTr+fw33niDwsJCXnrpJWbNmmWOX5ye+uo0lWEYWCyWbse/qqample5RK5Uwfqe8Hg8QZvtUld7zm5LweVyMXbsWMaPH88//dM/dbk/JyfH54P/4Q9/4Pnnn2fr1q2EhYUxYMAAPB4PoaGh1NXVERkZSVRUFGVlZeYy9fX1TJo0iaioKNxuN7GxsXi9XgzDwG63d1lHXFxcD5/qV73Xy+VEglvv3xP+VVNTE7TZLnU15KysrOz2vm5L4fe//32vVnbRZ599Rl5eHtu3b2fIkCEATJs2jZKSEm6//XYOHjxIUlISEydOJDs7m6amJmw2G1VVVWRlZfH5559TXFxMUlISpaWlTJky5bLyiIiIb92WwvDhwy/rgV977TU+/vhjHnzwQXMsNzeX7OxsCgoKiI6OxuVyYbfbWbVqFenp6VgsFjIyMggLC2POnDkcOnSI1NRUHA4Hubm5l5VHRER887mjubfmz5/P/Pnzu4xv27aty1hKSgopKSmdxi4emyAiIoHzrc99JCIiVy6VgoiImFQKIiJiUimIiIhJpSAiIiaVgoiImFQKIiJiUimIiIhJpSAiIia/HdEsIn1jVOaBPlnvqdy5fbJe+W5pS0FEREwqBRERMakURETEpFIQERGTSkFEREwqBRERMakURETEpFIQERGTSkFEREwqBRERMakURETEpFIQERGTSkFEREwqBRERMakURETE5NdSOH78ODNmzOA3v/kNAGfPnmXRokWkpaWxcuVK2traACgqKuKuu+7i7rvvprCwEACv18uqVatITU1l4cKFnD592p9RRUQEP5ZCS0sLGzZsYOrUqebYM888Q1paGrt372b48OEUFhbS0tLCli1b2L59O7t27WLr1q188skn7N+/n/DwcF555RWWLl1Kfn6+v6KKiMiX/FYKDoeDF154gcjISHOsoqKC5ORkAJKTkykvL+fo0aPEx8cTFhZGaGgoCQkJVFVVUV5ezsyZMwFITEyksrLSX1FFRORLfiuFkJAQQkNDO421trbicDgAcDqduN1uGhoaGDp0qPk9ERERXcZtNhtWq9WcbhIREf8I6DWaLRaL+bVhGJ3+vnTcYrF0O/5VNTU1fkgqIt+Wr/eix+PpF+/Xqz1nQEthwIABeDweQkNDqaurIzIykqioKMrKyszvqa+vZ9KkSURFReF2u4mNjcXr9WIYBna7vctjxsXF9TLNe71cTkS+jq/3Yk1NzWW8XwPnasj5TdPxAf1I6rRp0ygpKQHg4MGDJCUlMXHiRI4dO0ZTUxPNzc1UVVWRkJDA9OnTKS4uBqC0tJQpU6YEMqqIyFXJb1sK1dXVbNy4kTNnzhASEkJJSQlPPfUUmZmZFBQUEB0djcvlwm63s2rVKtLT07FYLGRkZBAWFsacOXM4dOgQqampOBwOcnNz/RVVRES+5LdSGD9+PLt27eoyvm3bti5jKSkppKSkdBqz2Wzk5OT4K56IiHwNHdEsIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiIKaDXUxCRK9eozAM9+C7/XMfkVO5cvzzu1UhbCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgIiImlYKIiJhUCiIiYtIRzSLS7/XsaOpvo2dHXl+JR1JrS0FEREwqBRERMQX19NHPf/5zjh49isViISsriwkTJvR1JBGRK1rQlsLbb7/NBx98QEFBASdOnOCRRx5hz549fR1LRMT03e/L6LnX7x3jl8cN2umj8vJyZsyYAcAPfvADmpqa+Pzzz/s4lYjIlS1otxQaGhoYN26ceXvYsGG43W4GDx7c6fsqKyt79fh77/5/l5VPRKQvtbS09Pr33zcJ2lIwDKPLbYvF0mnspptuCmQkEZErXtBOH0VFRdHQ0GDerq+vJyIiog8TiYhc+YK2FKZPn05JSQkAf/rTn4iMjOwydSQiIt+toJ0+mjx5MuPGjWPBggVYLBYef/zxgGfIy8ujsrKSCxcusGzZMuLj41mzZg3t7e04nU42bdqEw+EIeK6LWltbyczM5Ny5c5w/f57ly5cTGxsbVBkv5fF4mDt3LhkZGUydOjXoclZXV7N8+XJGjhwJQExMDEuWLAm6nABFRUVs3bqVkJAQVq5cSUxMTFDl3LNnD0VFRebt6upqXnnlFZ544gkArr/+etavX99H6f6iubmZtWvX8umnn+L1esnIyMDpdAZdzo6ODh5//HH+/Oc/Y7fbeeKJJxg4cKB/XnNDvlZ5ebmxZMkSwzAMo7Gx0bj11luNzMxM47XXXjMMwzA2btxovPzyy30Z0Thw4IDx61//2jAMw/jwww+N2267LegyXuoXv/iFceeddxp79+4NypwVFRXGk08+2WksGHM2NjYat912m/HZZ58ZdXV1RnZ2dlDmvKiiosJ44oknjIULFxpHjx41DMMwHnjgAaOsrKyPkxnGrl27jKeeesowDMP46KOPjFmzZgVlzoMHDxorV640DMMwPvjgA+Nv//Zv/faaB+30UV+7+eabefrppwG45ppraG1tpaKiguTkZACSk5MpLy/vy4jMmTOHpUuXAnD27FmioqKCLuNFJ0+e5MSJE/z4xz8GCMqczc3NXcaCMWd5eTlTp05l8ODBREZGsmHDhqDMedGWLVtYunQpZ86cMQ9ADZaM3/ve9/jkk08AaGpqYsiQIUGZ89SpU2am6667jtraWr+95iqFbthsNgYOHAh8sSn8ox/9iNbWVnPzzOl04na7+zKiacGCBTz88MNkZWUFbcaNGzeSmZlp3g7GnBc/4rdkyRLuueceDh8+HJQ5P/zwQwzD4MEHHyQtLY3y8vKgzAnwzjvv8P3vfx+bzUZ4eLg5HiwZ586dS21tLTNnzmThwoWsWbMmKHPGxMTw1ltv0d7eznvvvcfp06c5c+aMX17zoN2nECzeeOMNCgsLeemll5g1a5Y5bnzlI7N96dVXX6WmpobVq1d3+thusGTct28fkyZNYsSIEeZYMOaMjY0lIyOD5ORk3n//fRYvXsyFCxfM+4MlJ0BdXR3PPfcctbW1/OxnPwvKf0+AwsJC7rjjji7jwZLxd7/7HdHR0bz44ou8++67PPDAA+Z/BiF4ct56661UVVVxzz33cP311zNmzBiOHz9u3v9d5lQpfIM//OEPPP/882zdupWwsDAGDBiAx+MhNDSUuro6IiMj+zRfdXU1w4YN4/vf/z5xcXG0t7cHXUaAsrIyTp8+TVlZGR999BEOhyMoc44dO5axY8cCMHr0aCIiIjh79mzQ5Rw2bBg33ngjISEhXHfddQwaNAibzRZ0OeGL6bfs7GwsFos5TQMETcaqqioSExOBL/5T0NLSQktLi3l/sOQEeOihh8yvZ8yYQVRUlF9ec00fdeOzzz4jLy+Pf/7nf2bIkCEATJs2zfyY7MGDB0lKSurLiBw5coSXXnoJ+OII8JaWlqDLCLB582b27t3Lv/zLv3D33XezfPnyoMxZWFjIzp07AXC73Zw7d44777wz6HImJiZy+PBhOjo6aGxsDNrXva6ujkGDBuFwOLDb7YwZM4YjR44AwZNx5MiRHD16FIAzZ84waNAgYmJigi7nu+++yyOPPALAf/7nf3LDDTf47TW3GMGyfRRkCgoKePbZZxk9erQ5lpubS3Z2NufPnyc6OpqcnBzsdnufZfR4PDz66KPm/2bvv/9+xo8fz9q1a4Mm41c9++yzDB8+nMTExKDL+emnn/Lwww/T0tJCW1sb999/P3FxcUGXE76YMjxw4ACtra38/d//PfHx8UGXs7q6ms2bN7N161YATpw4wbp16+jo6GDixInmL7m+1NzcTO3+TkIAAAJHSURBVFZWFufOnePChQusXLkSp9MZdDk7OjrIysri/fffJywsjI0bN9Le3u6X11ylICIiJk0fiYiISaUgIiImlYKIiJhUCiIiYlIpiIiISaUgchn+7d/+jXHjxtHY2NjXUUS+EyoFkcuwf/9+RowYYR5EJNLf6TQXIr30ySef8M4775CTk8OLL75Iamoqhw4d4uc//zlOp5PY2FgGDhzIihUr+OUvf8mRI0dob29n4cKF/OQnP+nr+CJfS1sKIr30+uuv89d//dckJSXx/vvvU1dXx1NPPUVeXh5bt27lj3/8I/DF6UjOnDnDyy+/zM6dO/nVr36Fx+Pp4/QiX09bCiK9tH//fjIyMrDZbKSkpPD6669z5swZbrjhBgCSkpLo6OigqqqKo0ePsmjRIuCLUxa43e5OZ40VCRYqBZFeOHv2LO+88w65ublYLBY8Hg9hYWGdvsdqtdLR0YHD4WDevHksW7asj9KK9Jymj0R6Yf/+/dxzzz0UFRXxu9/9juLiYj799FNaW1s5efIk7e3t/Nd//RcAEyZMoLS0lI6ODs6fP8+GDRv6OL1I97SlINILBw4cIC8vz7xtsVhwuVxYrVZWrFjBtddey5gxY7DZbEyePJkpU6Ywf/58DMMgLS2tD5OLfDOdJVXkO/TWW28xatQorr32WtatW8ctt9yiTxpJv6ItBZHvkGEY3H///QwaNIhhw4Zx22239XUkkW9FWwoiImLSjmYRETGpFERExKRSEBERk0pBRERMKgURETGpFERExPT/AeP175vcYe+wAAAAAElFTkSuQmCC\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"plt.hist(adult['Age'])\n",
|
|
"plt.xlabel('Age')\n",
|
|
"plt.ylabel('Number of Records');"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"It's clear from this histogram that nobody in this particular dataset is over 90, so an upper bound of 90 would suffice.\n",
|
|
"\n",
|
|
"However, it's important to note that **this approach does not satisfy differential privacy**. If we pick our clipping bounds by looking at the data, then the bounds themselves might reveal something about the data.\n",
|
|
"\n",
|
|
"Typically, clipping bounds are decided either by using a property of the dataset that can be known without looking at the data (e.g. that the dataset contains ages, which are likely to lie between 0 and 125), or by performing differentially private queries to evaluate different choices for the clipping bounds.\n",
|
|
"\n",
|
|
"To use the second approach, we typically set the lower bound to 0 and slowly increase the upper bound until the query's output stops changing (meaning we haven't included any new data by increasing the bound). For example, let's try computing the sum of ages for clipping bounds from 0 to 100, using the Laplace mechanism for each one to ensure differential privacy:"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 13,
|
|
"metadata": {
|
|
"tags": [
|
|
"hide-input"
|
|
]
|
|
},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEOCAYAAACKDawAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deUBU5f7H8fcMAwKCgGwuuKCgIrijuZVmmOVyM7Pw5tbPdi2ta7lVbpRLi1bqLTOvGVpZylVTU+9NTW/ihooiIKKCisoi+zLAMOf3hzk3riCKDDPMfF//1Dkz55zvw8iHM88553lUiqIoCCGEsHhqUxcghBCidkjgCyGElZDAF0IIKyGBL4QQVkICXwghrIQEvhBCWAmzD/yEhARCQkJYt27dHd8XHx/PiBEjGDFiBH//+99rqTohhKg7zDrwCwsLCQsLo1evXlW+97333iMsLIyNGzeSmJhIUVFRLVQohBB1h1kHvp2dHatWrcLLy8uwLjExkXHjxjF+/HgmTpxIbm4uGRkZFBYWEhgYiFqtZsmSJTg4OJiwciGEMD9mHfgajQZ7e/ty68LCwpg/fz5r166lT58+rF+/npSUFNzd3Zk3bx7PPvss33zzjWkKFkIIM6YxdQH36tSpU7z33nsAlJSU0KFDBxRFISkpic8++wx7e3tCQ0Pp3bs3bdq0MXG1QghhPupc4Ds4OPDtt9+iUqkM6y5fvoy/vz9ubm4AdOvWjcTERAl8IYT4E7Pu0qlIu3bt2L9/PwDbt28nMjKSZs2aUVBQQHZ2Nnq9nri4OFq1amXiSoUQwryozHm0zJiYGBYvXkxKSgoajQZvb2/eeOMNPvnkE9RqNfXq1eOTTz7B1dWV6OhoPv74Y4qLi3nwwQd5/fXXTV2+EEKYFbMOfCGEEDWnznXpCCGEqB4JfCGEsBJmfZdOVFSUqUsQQog6qVu3bretM+vAh4qLvhtxcXEEBATUcDXmzRrbDNbZbmtsM1hnu6vT5spOlqVLRwghrIQEvhBCWAkJfCGEsBIS+EIIYSUk8IUQwkpI4AshhJWQwBdCCCshgS+EuC96vcIH22Pp/sG/eXPDSf4dm0qxrsxoxyvWlXEsKZOLGQVGO4YpRRy/wuvfn8AYw5yZ/YNXQgjzVaLT89ZP0WyNvkrv1u7siU/jnydS8HCy47sXe9LG2/m+9p+SXURCah7Xc7RczS7iP3HXiPsuCW2pHoCQAG9e6deK4JYN71jjubQ8rmVrcbSzoX49Da6OtjRxdcDW5vZzXkVRKCotI0+rI6uwhKyCUrIKS0jJKuJKViGpucX0a+vJyG4+5bbPKijhQkY+KdlarmUX0bWFG90rqevs9Tze3x5Laq6WT0O70L5JAwA2n0hh6k/RDGrfqNycHzVFAl8IUU6etpS9Z9NxdbDlQX+PSoOnoFjHK+uiOHAugxmPt+Plh1qh0yv8JzGD6RtP8dw/jvDPSX3wbmBf4fZ3ci2niE//dY6foi6j/+NEV6WC5i62jOrenJ6t3Im9lkt4ZBIjv0ylR8uGvP1YW0PAZheWEB6ZzO7YVM5ez6OkTH/bMWzUKpq6OtCogT2FpTpyikrJLdKRX6yjTF/x2bVzPQ0NHGzZeeY6X+w7z2sD/NCWlrHj9DWOXMzkz5upVDCpvx9vhPij+eMPw438YpbtSST8UDJO9TTYadQ8+fffmf9EIA3sbZn6UzQ9fd35dFTne/6Z3Q0JfCEsxI38YlbsPc+NgmK0pWUoCjzXuyW9/TwqfH9BsY5Nx6+QnldMPY0ajY2a48lZ7EtIp0R3MyCDW7gx7bF29PAtf6aq1yu89t1xfk/M4MOnOvJM92YA2NqoeLitF/94rjuhKyP5vzVH+fGVXjjVKx81OUWl7Dh9jUYN7Ono44K7Uz1ytaUcT85if0IG6w8noyjwf318GdyhMY1d7PF0rkdiwlnDMAOPBTXilX6t+PHoZVbsO8/TX0YyoJ0XLdwd2XD0MoUlZfTwbcj/9W1JUBMXmjV0RFtaRkGxjsyCEi5lFpJ0o5DUXC1ezvb4eTrhbG9LAwcNzva2ONtrcHWww83RFrf6djRxdcDFwRZFUdgTn8YnuxOYtvEUAH5eTkx62I+uLdxo6upAw/p2fLTzLMv3JnLowg0eC2rEv+NSOZqUhaIojH6gBX8b2AadXuHNDSeZvuk0KhV0be7G1+ODsbe1qdF/G7eY9Xj4UVFRMpbOPbDGNoN1tFtXpuf38zc4dOEGkedvcDUznzlPdGRIx8YAXM/RMvrrQ1zKLKSJqwP2Ghuyi0pIzS3mpYdaMfXRNtTT3AyR/GId30Ym8fWBi2QWlKBSwa0U8G5Qj8EdGjO4Q2MSUvP4/NdzpOYWM7hDIz4a2Yn6fwT357+eY8m/Epj/RCDjerWssObfEtKZ8M1RerZqyOsD/Onk44qNWsX6w8l8/us5sgpLDe/1dK5HRn4xinLzzPsvnZrwt4FtaNbQsdw+K/usi0rKWHPwIl/uO09BSRlPdGrCy/1a07bR/XUp3Yler3Do4g08nerhX0nX1ZaTKcyKOE1BSRltvZ0Z2N6bJzo3Kff+Mr3C3/cmcjolh4+f6UQDe9ty+6juWDp1cvA0IaxFfrGOjccuk1lYypgHmuP1R1fIudQ8/vZjNKdTctCoVXRu5opzPTWTvjvOoQstGN+7JRO+OcqN/GLCn3+Anq3cgZsh+P72WL7af4H9Cem0dK/PhYx8kjIKKSnT06+NJ5Mf8adrc1d0eoVinR5HWxvU6ptdON1bNuSprj58feACS/6VQFJGIaufCyYhNZ+l/05geOcmjO3ZotL29GvjycIRHZgZcZrfEw+hUatwcbDlRkEJffzcmfpoW0p1ek5dySH+eh7NGzrSvaUbnZq5Gv6w3C0HOxsm9vdjXK+WFJeW4e5Ur5qfwt1Tq1X0bl3xt6dbnujclL5+HhSWlN32x+sWG7WK1x/xN0aJt5HAF8LELmcWsub3JH48dpn8Yh0qFaz87TyjH2iBp3M9lv47gfp2Nnwa2plHA71xtNNwKiaWn5Nh1YGLrDucTAN7W9a/2JPOzVwN+3Wws+GDJzvwcFsv5m+LJSEtj1YeTjzc1ovHOzQu915bG1WFFzDtbW14bYA/gU1deP27Ewxf8TslOj1tvJxZMKJDlRcWnwluxsAAb45fyuJYchZJGQU8E9yM/m09Dds+8McfqJrgVE9zW/eRqbk71aPmWnh/zOsnI4QFUBSFazlazlzN5ez1XDr4uNKvjedt70vKKGDF3kQiTqSgAoZ2bMz/9fHF1dGWZXsSWRuZRJleISTAm4UjOuDp/N+zVlsbFe8MCaBnK3e+OZjErMEBBDRuUGE9Ie29CWnvfV9teritFxtf7cXz3xyjoKyML8Z0xdHu7uLDrb4djwR480jA/dUg7p9RAz8hIYGJEyfy3HPPMWbMmHKvHTp0iCVLlqBWq/H19eWDDz5ArZbHAkTddjmzkP/75iiJafnl1g9o58Xsoe1p5GLPvrNpbDl5lV1nrmNro2Zcrxa8/FBrGrn8926Wj5/uxOsD/LiWo+UB34aVnknXZpC2a9SAX954kHytjiauDrVyTFGzjBb4hYWFhIWF0atXrwpfnz17Nt9++y2NGjVi8uTJHDhwgH79+hmrHCGM7kJ6PqO/PkxhSRlzh7Wng48Lfp7O/HjsMp/9eo5Hl+7HTqMmv1iHh5MdLzzYihce9MXLueLbFlu416eFe/1absWdNbC3ve2ioqg7jBb4dnZ2rFq1ilWrVlX4ekREBE5OTgA0bNiQrKwsY5UihNElpObx7KrDKIrCDy/1LNe98uJDrXiicxOW7UmktEzP0I5N6NmqoeHebCFqi9Fvy1y2bBlubm63denckpaWxpgxY9iwYQNubm7lXouKisLRseIr21XRarXY29/7Ax91mTW2GWq33SVleuz+FNQ6vcL2s7msO5mFnY2KhY82prmrndHrkM/aelSnzYWFheZ3W+aNGzd45ZVXmD179m1hf0t176+2hnuz/5c1thmM3+4yvcK+s2msO5TMvoR0fNwc6OvnSbtGzqw9mMSFjAL6+Lmz8MmONHev3gnKvZLP2nrU5Jy2Jgv8/Px8XnzxRaZMmULfvn1NVYYQd3Q0KZO//XiSy5lFeDeox/N9fEnOLGRb9FW+P6KjtWd9/vFcMA+39TLK2CdC1CSTBf6iRYsYP368XKgVZmvLyRTe/ukUTd0c+GJ0V0LaexvuVdeV6bmYUUBLj/oV3r8uhDkyWuDHxMSwePFiUlJS0Gg07Nq1iwEDBuDj40Pfvn3ZvHkzycnJbNy4EYChQ4cSGhpqrHKEuGuKorBsTyJL/pXAA74NWTm2G66O5fvlNTbqSh+nF8JcGS3wg4KCCA8Pr/T1mJgYYx1aiGo7czWHeVtjOZKUyYiuTVk0oiN2GjmDF5ZBnrQVVktRFK5kFZGeX8yN/BL2nU3j+yOXcHW0Y9GIDoR2byb98sKiSOALq1RYouONH06yOzbVsM5GrWJ875a88UgbXBzl4SJheSTwhdW5nqPlhW+PEns1lzdC/OnUzBWP+vVo4mpfK6MsCmEqEvjCahTryvjtbDqzt5whT1vK6vHdebidl6nLEqLWSOALi3f2eh5fH7jAzjPXydPqaOrqwMZXe1c6uqQQlkoCX1i0A+fSeSU8CpVKxaOB3gzr1IS+fh5y77ywShL4wmJtPpHCWz9F4+flxNoJPao1mbYQlkQCX1iEMr3Cgh1xHLmYiaOdDbY2av6TmEGvVu6sHNdNhvQVApDvtaLOUxSFdzfHsPo/F3Gws0EBcopKGf1Ac76Z0F3CXog/yBm+qNMUReEfUZlsPJPDxP6tmfZYO1OXJITZksAXdZKiKJxLy+e7w5fYeCaHsT1b8PagtqYuSwizJoEv6hS9XuGzX8+x6fgVrmQVATDI35l5fwmUYRCEqIIEvqgz9HqFdzaf5vsjl+nXxpOJ/f0Y0M6LrKsXUasl7IWoigS+qBMURWHuz2f4/shlXnvYj6mPtjGc0WddNXFxQtQRcpeOqBMW7Ijj28hkXnqoVbmwF0LcPQl8YfZ2nbnOqgMXGderBTMfbydhL0Q1SeALs5ZVUMI7/4whsEkD3hvaXsJeiPsgffjCrM3ZeoacohLCn+8h498IcZ/kN0iYrZ0x19gafZXXB/jLyJZC1AA5wxdm5ez1PA5duMGJS1n8GpdGUNMGvNq/tanLEsIiSOALs7Hh6CWmbzoNgJdzPfr6ezDtsXbSlSNEDZHAF2ZhZ8x1Zkac5kF/DxY91ZEmLvZygVaIGmbUU6eEhARCQkJYt27dba8dPHiQkSNHEhoayooVK4xZhjBzkedvMPmHE3T0cWXl2G40dXWQsBfCCIwW+IWFhYSFhdGrV68KX3///fdZtmwZ33//PQcOHCAxMdFYpQgzdvZ6Hi9+e4wWDR1Z81x3HO3kS6cQxmK0wLezs2PVqlV4ed0+SfTly5dxcXGhcePGqNVq+vXrR2RkpLFKEWYqu7CEF789hqOdDd8+3wO3+namLkkIi2a0wNdoNNjbVzylXHp6Og0bNjQse3h4kJ6ebqxShBnSlel5/fsTXM/R8uXYbjR2cTB1SUJYPJN8f1YU5bZ1lfXZxsXFVesYWq222tvWVXWpzauP3eDAuRym9PLAoeA6cXHXq72vutTummKNbQbrbHdNttkkge/t7U1GRoZhOTU1FU9PzwrfGxAQUK1jxMXFVXvbuqoutFlRFP6+77xh0pI3nwi6733WhXbXNGtsM1hnu6vT5qioqArXm+QGZx8fH/Lz87ly5Qo6nY69e/fSp08fU5QiapG2tIw3Npzko11neaJzE94b2t7UJQlhVYx2hh8TE8PixYtJSUlBo9Gwa9cuBgwYgI+PDwMHDmTu3LlMnToVgMGDB+Pr62usUoQZSMvV8mJ4FNGXs3l7UFsm9m8tt14KUcuMFvhBQUGEh4dX+nr37t3ZsGGDsQ4vzEjctVye/+Yo2UWlrBzbjUGBjUxdkhBWSW56Fka172war313gvr1bPjx5V4ENXUxdUlCWC0JfGE0Px67zMyI07T1dmb1c8Fy66UQJiaBL4zi28gkZm85w4P+HnwxphtO9eSfmhCmJr+Fosat2n+BD3bEERLgzYrRXainsTF1SUIIJPBFDVv523kW/hLPkA6N+XRUZxnaWAgzIoEvasz3Ry6x8Jd4hnZszKehndFI2AthVuQ3UtSIX05f451/nqZfG0+WPCNhL4Q5kt9Kcd8OJmYw5YeTdGnuxpdjumGnkX9WQpgj+c0U9+V6jpbXvj9BC3dH/jG+Ow52coFWCHMlgS+qTVemZ/IPJ9CWlvHFmG64ONqauiQhxB3IRVtRbZ/9eo4jFzNZ8kwn/LycTF2OEKIKcoYvquW3hHSW703k6W4+jOjqY+pyhBB3Qc7wxT1RFIW1B5P4YEcc/l5OzHsi0NQlCSHukgS+uGv5xTpmbDrFtlPXeKSdF0ue6SyTjgtRh8hvq7gredpS/rrqELFXc3l7UFte7dcatVrGsxeiLpHAF1XSlpbx4rfHiL+Wx6pxwTwS4G3qkoQQ1SCBL+5IV6Znyg8nOHQhk89GdZawF6IOk7t0xB3N2XqGXWdSmTusPU90bmrqcoQQ90ECX1Tql9PXWH/4Ei/3a8VzfWTOYSHqOgl8UaHUXC0z/3majj4uvPVoW1OXI4SoARL44jaKovD2xlNoS8tYGipj2gthKeQ3Wdzm28hk9iek887gAFp7ypAJQlgKo96ls2DBAqKjo1GpVMyaNYuOHTsaXlu/fj1bt25FrVYTFBTEO++8Y8xSxF26mFHAwl/i6NfGkzE9W5i6HCFEDTJa4B85coTk5GQ2bNhAYmIiM2fO5KeffgIgPz+f1atXs3v3bjQaDRMmTODkyZN07tzZWOWIu6DXK0zbGI2tjZrFT3VEpZIHq4SwJEbr0omMjCQkJAQAPz8/cnNzyc/PB8DW1hZbW1sKCwvR6XQUFRXh4uJirFLEXVobmcTRpCxmD21PIxd7U5cjhKhhRgv8jIwM3NzcDMvu7u6kp6cDUK9ePSZNmkRISAgDBgygc+fO+PrKbX+mlJRRwOKd8fRv68nIbjL6pRCWyGhdOoqi3LZ8q4sgPz+flStXsnPnTpycnBg/fjzx8fG0a9futv3ExcVV6/harbba29ZV1W2zXlGYvusaNsCEDg7Ex8fXfHFGJJ+19bDGdtdkm40W+N7e3mRkZBiW09LS8PDwAOD8+fM0a9aMhg0bAhAcHExMTEyFgR8QEFCt48fFxVV727qqum1eezCJmFQtHz7VkYeCmxmhMuOSz9p6WGO7q9PmqKioCtcbrUunT58+7Nq1C4DY2Fi8vLxwcrp5i1/Tpk05f/48Wq0WRVGIiYmhZcuWxipF3MGlG4Us+iWefm08eTpYunKEsGRGO8Pv2rUrgYGBjBo1CpVKxZw5c4iIiMDZ2ZmBAwfy/PPPM27cOGxsbOjSpQvBwcHGKkVUQq9XmLYpGo1axcIRHeSuHCEsnFHvw3/rrbfKLf+5y2bUqFGMGjXKmIcXVVh/OJlDFzJZ/FQHmrg6mLocIYSRVRn4p06dYvv27eTl5ZW7ELtw4UKjFiaM60pWIQt/iedBfw+eqYP99kKIe1dl4L/99tu8+OKLhguuou5TFIVZ/4wBkK4cIaxIlYHfqlUrnnrqKQkFC7L5ZAr7E9KZO6w9Pm6Opi5HCFFLqgz8oUOHMnz4cNq2bYuNjY1hvXTp1E038ouZ/3MsXZu7MrZXS1OXI4SoRVUG/qeffspLL72Ep6dnbdQjjGzez7HkF+tY9FRHbGQSciGsSpWB37p1a55++unaqEUY2caoK2yNvsobIf608XY2dTlCiFpWZeC7ubkxevRogoKCynXpTJs2zaiFiZp1/FIWsyJO07u1O6897GfqcoQQJlBl4Pfo0YMePXrURi3CSK7naHk5PIpGLvaseLYrGpnBSgirdFcPXskdOnVXsa6Ml8OPUVisY93zD+BW387UJQkhTKTKwE9ISDD8v06nIzo6Gn9/f4YPH27UwkTN+GLfeaKv5PDlmG60bST99kJYsyoDf/r06eWWy8rKmDx5stEKEjUnKaOAv+87z9COjXksqJGpyxFCmFiVgV9UVFRuOT09nQsXLhitIFEzFEXhvS0x2NmoeW9oe1OXI4QwA1UG/pAhQ1CpVIYJTJydnZkwYUJt1Cbuw7ZT1zhwLoO5w9rj3UCmKxRC3EXg79mzx/D/er0etVru8DB3udpSwrbFEtS0gTxNK4QwqDS94+Pjyw1vPHPmTB544AEGDx7M6dOna6U4UT1hP8eSkV/MB8M7yNO0QgiDSgN//vz5PPPMMwDs37+f6Oho9u/fz1dffcXHH39cawWKe7P7zHV+irrCq/1b06mZq6nLEUKYkUq7dGxsbAwPXP36668MHz4cBwcHfHx85L58M5VdVMbM7acJbNKAKY+0MXU5QggzU+kZfklJCXDzNswDBw7Qp08fw2tardb4lYl7oigKn0Wmk1esY2loZ+w0cq1FCFFepWf4vXv35pVXXqGoqIiWLVsSGBiITqdj+fLltGrVqjZrFHch4ngKhy4X8u6QABkYTQhRoUoDf8qUKRw9epTc3FwefPBBAMMdOnPmzKmd6sRdycgvJmx7LO096zGhj6+pyxFCmKk73pbZvXv3cstqtZo33njDqAWJezf/51gKinVMHtgUtdyVI4SohHT01nF749PYGn2VSQ/70cJVBkYTQlTOqIG/YMECQkNDGTVqFKdOnSr32rVr1/jrX//KyJEjmT17tjHLsFgFxTre3RyDn5cTr/ZvbepyhBBmrtIunfXr199xw9GjR9/x9SNHjpCcnMyGDRtITExk5syZ/PTTT4bXFy1axIQJExg4cCDz5s3j6tWrNGnS5B7Lt24f7z5LSnYRG1/pRT2NTdUbCCGsWqWBn5WVdV87joyMJCQkBAA/Pz9yc3PJz8/HyckJvV5PVFQUS5YsAeQicHWcvJzNNweTGNOzOcEtG5q6HCFEHVBp4L/22muG/y8oKCAnJwe4eX/+/Pnzq9xxRkYGgYGBhmV3d3fS09NxcnIiMzMTJycnPv/8c6KioujSpQt/+9vf5IGuu1RapmfGplN4Oddj2mPtTF2OEKKOqHLwtBUrVhAREUF2djZNmjTh6tWrhIaGVrljRVFuW74V6IqikJqaylNPPcXkyZN56aWX+O233+jfv/9t+4mLi7vLppSn1Wqrva25+/F0NvHX83jvYW9SLiaS8sd6S27znVhju62xzWCd7a7JNlcZ+Pv37+fXX39l7NixhIeHc+bMGXbu3Fnljr29vcnIyDAsp6Wl4eHhAdycGL1x48Y0b94cgF69enHu3LkKAz8gIOBu21JOXFxctbc1Z0kZBXx3KolBgd48Pyi43GuW2uaqWGO7rbHNYJ3trk6bo6KiKlxf5V06t8bCLysrQ6vVEhgYWOnO/qxPnz7s2rULgNjYWLy8vHBycgJAo9HQrFkzkpKSADhz5gy+vvLAUFUURWFGxCnsbNTM+0uQqcsRQtQxVZ7hDxo0iLVr1zJs2DCeeOIJ3N3dcXBwqHLHXbt2JTAwkFGjRqFSqZgzZw4RERE4OzszcOBAZs2axZw5cyguLsbf358BAwbUSIMs2Q9HL3PoQiYLR3SgkYtMaiKEuDdVBv7jjz9Oo0Y350Pt168fWVlZ2NvfXdj8eTx9gHbt/nuBsUWLFnzzzTf3UKp1u5ZTxILtcfRq5c6o7s1MXY4Qog6qtEsnMzOTc+fO8frrr3P+/HkSExMpLCwEYOLEibVWoLjZlfPuP2Mo1etZ9FQHuZtJCFEtlZ7hX7hwgU2bNpGUlMTcuXMN69VqNcOGDauN2sQffj51jV/j03h3SAAt3OubuhwhRB1VaeAHBwcTHBzMsGHD6N27Nzk5OdjY2BguvIrakffH/LQdfVz4PxkJUwhxH6rsw1cUhUGDBqHRaAyTmM+fP59u3brVRn1W7/Nfz5GRX8zX44JlflohxH2pMvCXLVtGeHg4Xl5ewM1Bz6ZOncp3331n9OKs3bnUPNb8nkRocDOZn1YIcd+qvA/f1tbWEPYAjRs3RqOp8u+EuE+KojBn6xnq19PI8AlCiBpRZXL7+Pgwb948evTogaIoHD582PCErDCe7aevcfD8DcKGB9GwvoxzL4S4f5We4U+ePBmAsLAwOnXqRFRUFCdOnKBbt27Mmzev1gq0RkUlZSzYHkf7xg14tof8cRVC1IxKz/Czs7NvvkGjYfjw4QwfPrzWirJ2X+2/wNUcLZ+O6iIXaoUQNabSwL906RIffvhhpRtOmzbNKAVZu6vZRXzxWyJDOjamh6+Mcy+EqDmVBr6DgwP+/v61WYsAFv0Sj6LAzMflQq0QomZVGvgeHh48+eSTtVmL1TuWlMnW6KtMHuCHj5ujqcsRQliYSi/aBgXJ8Lu1SVemZ87WMzRqYM8rMiG5EMIIKg386dOn12YdVm9tZDJnruby3tD2ONrJcw5CiJpX5YNXwviu5RSxZPdZ+rf1ZHCHRqYuRwhhoSTwzcDcrWfQ6RXCngiSoY+FEEYjgW9i/45NZdeZVCY/4k+zhnKhVghhPBL4JpSaq2XmP0/j7+XEiw+2MnU5QggLJ1cHTURbWsbL4VEUFOtY9/wD2Gnkb68Qwrgk8E1AURRmb4nh5OVsvhzTlbaNnE1dkhDCCshppQmEH0rmx2NXmDzAj8eCGpu6HCGElZDAr2U38otZuCOe/m09eSOkjanLEUJYEaMG/oIFCwgNDWXUqFGcOnWqwvd88sknjB071phlmJV//H4Rra6Md4e0Ry0jYQohapHR+vCPHDlCcnIyGzZsIDExkZkzZ/LTTz+Ve09iYiJHjx7F1tbWWGWYlZyiUr49mMzgoMb4eclk8EKI2mW0M/zIyEhCQkIA8PPzIzc3l/z8/HLvWbRoEW+++aaxSjA73x5MIq9Yx6SH/UxdihDCChkt8DMyMnBzczMsu7u7k56ebliOiIigR48eNG3a1FglmJWCYh2rf79ISIAX7Zs0MHU5QggrZLQuHUVRblu+NWxAdnY2ERERrFmzhtTU1DvuJy4urlrH12q11d7WGDbGZJNdWMoQXx3UKMIAABTlSURBVI3R6jK3NtcWa2y3NbYZrLPdNdlmowW+t7c3GRkZhuW0tDQ8PDwAOHToEJmZmYwePZqSkhIuXbrEggULmDVr1m37CQgIqNbx4+Liqr1tTcvVlrJl0z76+nnw5ENdjHYcc2pzbbLGdltjm8E6212dNkdFRVW43mhdOn369GHXrl0AxMbG4uXlhZPTzQuVjz32GDt27ODHH39k+fLlBAYGVhj2luKjnWfJLChh+mMyi5UQwnSMdobftWtXAgMDGTVqFCqVijlz5hAREYGzszMDBw401mHNzvFLWaw7nMz/9falg4+LqcsRQlgxow6t8NZbb5Vbbtfu9jNcHx8fwsPDjVmGyZSW6ZkVcZpGDez526PykJUQwrRkLB0jWvP7ReKv57FybDec6smPWghhWjK0gpGkZBex9F/nGNjem0GBMouVEML0JPCN5IPtsSgozBnW3tSlCCEEIIFvFP85l8GO09d57WE/fNxkFishhHmQwK9hJTo9s7fG0MLdkRdkFishhBmRK4k1bM3vF7mQXsCa57pjb2tj6nKEEMJAzvBrUGquls9/PUdIgDcPt/MydTlCCFGOBH4NWrwzntIyhdlD5UKtEML8SODXkBOXsog4nsILD/rS3F0u1AohzI8Efg3Q6xXm/RyLl3M9JspY90IIMyWBXwM2n0zh5OVspj/WTp6oFUKYLQn8+5RfrGPRL/F0aubKk12sYzIXIUTdJIF/nz7cGU96fjHz/hIok5ILIcyaBP59OJqUybeRN4c+7tzM1dTlCCHEHUngV5O2tIzpm07h4+bAW4Nk6GMhhPmTK4zVtGzPOS6kF7Du+QdwtJMfoxDC/MkZfjXEXcvly98u8EywD339PUxdjhBC3BUJ/HukKArvbY7BxcGWWYOtazJlIUTdJoF/jzYdT+FYchYzHm+Hq6OdqcsRQoi7JoF/D3KKSlm4I46uzV0Z2dXH1OUIIcQ9kauN92DJ7rNkFZawdkIPuedeCFHnyBn+XYpJySH8UDJjerYgqKmLqcsRQoh7JoF/F8r0CjMjTtOwfj2mDmxr6nKEEKJajNqls2DBAqKjo1GpVMyaNYuOHTsaXjt06BBLlixBrVbj6+vLBx98gFptnn9/vjmYxOmUHJb9tQsujramLkcIIarFaAl75MgRkpOT2bBhA++//z5hYWHlXp89ezaff/45P/zwAwUFBRw4cMBYpdyXlOwiPtl9lv5tPRnasbGpyxFCiGozWuBHRkYSEhICgJ+fH7m5ueTn5xtej4iIoFGjRgA0bNiQrKwsY5VSbYqiMHtzDIoCYU8EoVLJhVohRN1ltC6djIwMAgMDDcvu7u6kp6fj5OQEYPhvWloaBw8eZMqUKRXuJy4urlrH12q11d72lt+TC/g1Po0XghuSn5pMXOp97c7oaqLNdZE1ttsa2wzW2e6abLPRAl9RlNuW//cM+caNG7zyyivMnj0bNze3CvcTEFC9p1nj4uKqvS3cHBzthS2/0dbbmRlPPoDGxjyvL/zZ/ba5rrLGdltjm8E6212dNkdFRVW43miB7+3tTUZGhmE5LS0ND4//jjuTn5/Piy++yJQpU+jbt6+xyqi2r/ZfICW7iO9erBthL4QQVTFakvXp04ddu3YBEBsbi5eXl6EbB2DRokWMHz+efv36GauEakvJLuLv+xIZ3KERvVvL4GhCCMtgtDP8rl27EhgYyKhRo1CpVMyZM4eIiAicnZ3p27cvmzdvJjk5mY0bNwIwdOhQQkNDjVXOPVm4Iw5FQQZHE0JYFKPeh//WW2+VW27Xrp3h/2NiYox56Go7fOEG205dY/Ij/vi4OZq6HCGEqDHSOf0npWV6Zm85Q1NXB17t19rU5QghRI2SwP+TtQeTOJuax+xh7XGwszF1OUIIUaMk8P9wPUfL0n8lMKCdF4+29zZ1OUIIUeMk8P8Qtj0WnV5h7rBAeaJWCGGRJPCBA+fS2X7qGpMe9qO5u1yoFUJYJqsP/GJdGbO3nKGluyMvPdTK1OUIIYTRWP2MVyt/u8DFjAK+ndADe1u5UCuEsFxWfYaffKOA5XsTGdKxMQ+18TR1OUIIYVRWG/iKojB36xls1SreG9Le1OUIIYTRWW3g7zpznb1n03lzYBsaudibuhwhhDA6qwz8whId83+OpV0jZ57r3dLU5QghRK2wysBftieRqzla3h8eJEMfCyGshtWlXWJaPl8fuMDIbj4Et2xo6nKEEKLWWFXgK4rCnK0xONjaMOPxdlVvIIQQFsSqAn/bqWv8nniDtwe1xcOpnqnLEUKIWmU1gZ9frOP97bEENW3Asw+0MHU5QghR66zmSdtP/5VAWl4xK8cGY6OWwdGEENbHKs7w467lsuZgEn/t0ZzOzVxNXY4QQpiExQe+Xq/w3uYYXBxsmTaoranLEUIIk7H4wN90/ArHkrOY8Xg7XB3tTF2OEEKYjEUHfnZhCQt/iadbCzdGdvUxdTlCCGFSFh34H+06S3ZhCWFPBKGWC7VCCCtn1MBfsGABoaGhjBo1ilOnTpV77eDBg4wcOZLQ0FBWrFhR48c+m6HluyOXGN+7Je2bNKjx/QshRF1jtMA/cuQIycnJbNiwgffff5+wsLByr7///vssW7aM77//ngMHDpCYmFhjxy7TK6w4dANPp3r8bWCbGtuvEELUZUYL/MjISEJCQgDw8/MjNzeX/Px8AC5fvoyLiwuNGzdGrVbTr18/IiMja+zY3x+5xLkbxbwzJABne9sa268QQtRlRnvwKiMjg8DAQMOyu7s76enpODk5kZ6eTsOG/x24zMPDg8uXL1e4n7i4uHs+9rao63RpVA9/uxzi4nLvvfg6SqvVVuvnVddZY7utsc1gne2uyTYbLfAVRbltWaVSVfgaYHjtfwUEBNzzsde0bsP5c2dp3966ZrKKi4ur1s+rrrPGdltjm8E6212dNkdFRVW43mhdOt7e3mRkZBiW09LS8PDwqPC11NRUPD1rbk5ZBzsbGT5BCCH+h9ECv0+fPuzatQuA2NhYvLy8cHJyAsDHx4f8/HyuXLmCTqdj79699OnTx1ilCCGEwIhdOl27diUwMJBRo0ahUqmYM2cOERERODs7M3DgQObOncvUqVMBGDx4ML6+vsYqRQghBEYeLfOtt94qt9yu3X8nHenevTsbNmww5uGFEEL8iUU/aSuEEOK/JPCFEMJKSOALIYSVkMAXQggroVIqegrKTFT28IAQQog769at223rzDrwhRBC1Bzp0hFCCCshgS+EEFbCqA9emcKCBQuIjo5GpVIxa9YsOnbsaOqSjObDDz8kKioKnU7Hyy+/TIcOHZg2bRplZWV4enry0UcfYWdnmfP4arVahgwZwqRJk+jVq5fFt3vr1q18/fXXaDQapkyZQps2bSy+zQUFBUyfPp2cnBxKS0uZNGkSnp6ezJ07F4C2bdsyb9480xZZgxISEpg4cSLPPfccY8aM4dq1axV+xlu3bmXt2rWo1WpCQ0MZOXLk3R9EsSCHDx9WXnrpJUVRFOXcuXPKyJEjTVyR8URGRiovvPCCoiiKkpmZqfTr10+ZMWOGsmPHDkVRFGXx4sXK+vXrTVmiUS1ZskQZMWKEsmnTJotvd2ZmpvLoo48qeXl5SmpqqvLuu+9afJsVRVHCw8OVjz/+WFEURbl+/boyaNAgZcyYMUp0dLSiKIoyefJkZd++faYsscYUFBQoY8aMUd59910lPDxcURSlws+4oKBAefTRR5Xc3FylqKhIGTRokJKVlXXXx7GoLp07Tbpiabp3785nn30GgIuLC0VFRRw+fJhHHnkEgEceeaRGJ5UxJ+fPnycxMZH+/fsDWHy7IyMj6dWrF05OTnh5eREWFmbxbQZwc3MjOzsbgNzcXFxdXUlJSTF8a7ekdtvZ2bFq1Sq8vLwM6yr6jKOjo+nQoQPOzs7Y29sTHBzM8ePH7/o4FhX4GRkZuLm5GZZvTbpiiWxsbHB0dATgp59+4qGHHqKoqMjwtd7T09Ni27548WJmzJhhWLb0dl+5cgVFUXjjjTd49tlniYyMtPg2AwwZMoSrV68ycOBAxowZw7Rp02jQ4L/zU1tSuzUaDfb29uXWVfQZZ2Rk3DZ51L38DCyqD1+5w6Qrlurf//43Gzdu5B//+AeDBg0yrP/fn4Wl2Lx5M507d6ZZs2aGdX/+jC213ampqSxfvpyrV68ybtw4q2jzli1baNKkCatXryY+Pp7JkycbTnLActt9S0Wf8f1mnEUF/p0mXbFEBw4c4Msvv+Trr7/G2dkZBwcHtFot9vb2pKamlvt6aCn27dvH5cuX2bdvH9evX8fOzs7i2+3u7k6XLl3QaDQ0b96c+vXrY2NjY9FtBjh+/Dh9+/YFbo60W1hYSGFhoeF1S233LRX9u/b29mbfvn2G96SlpdG5c+e73qdFdencadIVS5OXl8eHH37IypUrcXV1BaB3796G9u/evZsHH3zQlCUaxaeffsqmTZv48ccfefrpp5k4caLFt7tv374cOnQIvV5PZmYmhYWFFt9mgBYtWhAdHQ1ASkoK9evXp02bNhw7dgyw3HbfUtFn3KlTJ06fPk1ubi4FBQUcP36c4ODgu96nxT1p+/HHH3Ps2DHDpCt/HoPfkmzYsIFly5aVmzhm0aJFvPvuuxQXF9OkSRMWLlyIra2tCas0rmXLltG0aVP69u3L9OnTLbrdP/zwA9u3b6eoqIhXX32VDh06WHybCwoKmDVrFjdu3ECn0zFlyhQ8PT2ZPXs2er2eTp06MXPmTFOXWSNiYmJYvHgxKSkpaDQavL29+fjjj5kxY8Ztn/HOnTtZvXo1KpWKMWPG8Je//OWuj2NxgS+EEKJiFtWlI4QQonIS+EIIYSUk8IUQwkpI4AshhJWQwBdCCCshgS9qVVJSEi+99BIjR45kxIgRhIWFUVJSwpUrVxgxYgQAb775Jlqt9p72+9VXX3HixIlq1xUREUG/fv0YO3YsY8aMYcKECeUe4qspBQUFDBgw4Lb18+fP58knn6zW2E/Hjx+nbdu2xMfH10SJwoJJ4ItaU1ZWxuuvv84LL7zAxo0b2bRpEwArVqwo976lS5feNq5IVV566SW6dOlyX/UNHjyY8PBw1q1bR9euXQ311YbffvuNtWvXVutBwW3btuHr68u2bduMUJmwJBY1tIIwb7///jutWrWiR48ewM2xQt5++23UajVpaWmG9w0YMICff/6ZsLAwHB0duXDhAllZWSxcuJAGDRowZcoUfH19uXjxIh06dGDu3LnMmDGDQYMGkZWVRVRUFJmZmVy8eJHnn3+ep59+ms2bN7N69WqaNGmCt7c3nTt3NnyjqMiNGzfo1KkTcHPUwqVLlxoeiFm4cCHbtm3j3LlzTJ8+nYKCAoYNG8aePXsYOHAgoaGh7N27l5KSEtasWQPA66+/DlDh/Axff/01aWlpvPLKK6xcuZIvvviC48ePU1ZWxujRoxk+fDhjx47F398fgNmzZxu2LSsrY/fu3SxdupTp06fz1ltvARAfH8+MGTNwdname/fuXL16lUWLFrF+/Xp+/vln1Go1ISEhTJgw4X4+UlHHyBm+qDUXLlwgICCg3Dp7e/s7Ttyh0+n45ptvmDJliuGbwNmzZ5k6dSobN27k9OnTt3VlJCQksHz5clasWMG6devQ6/UsWbKENWvWsHTpUv7zn/9UeKwdO3YwduxYhg4dSmxsrGEwujlz5rB06VLWrVuHi4sLP//8c6X1lpWV0apVK9avX4+Pjw+HDh1iy5Yt+Pv7s2bNGtq2bXvbNi+88AKenp6sWrWK+Ph4zp07xw8//MDatWtZvny5oZvH39+/XNjDzT+irVu3pnv37ri6uhq6tVasWMGkSZMIDw/n/PnzAFy+fJmdO3fy/fffs379enbv3s3Vq1crbYuwPBL4olaVlZXd0/t79+4NQOfOnbl48SIALVu2pHHjxqhUKjp16sSFCxfKbdO5c2dsbGxo1KgReXl5ZGVl4ezsjIeHB46OjvTs2bPCY93q0tm2bRt//etfmT17NtnZ2ahUKho3bgxAcHAwcXFxd6z51tgmt45//vx5Q3fTrW83lYmJiaF79+4AODo60rJlS5KTk4GKvx1s27aNoUOHAjBs2DC2b98O3JwzoFu3bgA8/PDDAJw+fZrk5GTGjRvHuHHjKCgoICUl5Y71CMsiXTqi1rRu3Zr169eXW1dSUkJSUlK5YW//TK/XG/7/1jCwf15X0fCwGk35f9b/+x61uurznEGDBvHZZ5+hUqnKDUmr1+tRqVTl9qfT6cpta2NjU+7YiqIYjvnn2ivyv23587b/O1aOVqtl7969xMbGsm7dOkpLS8nNzWXWrFnlav7z9v3792f+/PlVtl9YJjnDF7WmT58+pKSksGfPHuBm+H300Ufs2LGj0m1uzeZz4sQJWrduDcClS5dIS0tDr9cTHR2Nn5/fHY/r6upKdnY2OTk5aLVajhw5UmWt0dHR+Pr64uLigkqlMnR9HDlyhKCgIJycnAzXHaKiou64L19fX2JiYoCb1wPuJCgoyPCegoICLl26RIsWLSp87549e+jZsyfbtm1jy5Yt7Nixg1atWnH48GGaN29uOOb+/fsBCAwM5PDhwxQVFaEoCu+///493w0l6jY5wxe1Rq1Ws3r1ambPns3y5cuxs7Ojd+/evPbaa5X2JWu1Wl5++WWuX7/Ohx9+CNwM0KVLl5KYmEjXrl0NFzMro9FoePXVVxk9ejQtWrQgKCio3Fn4LTt27DCEJGCYLDssLIypU6ei0Wjw8fFhyJAhaLVavvjiC8aOHUu/fv3uOAnF8OHDmTRpEuPHjzd0s1QmODiYoKAgRo8ejU6nY+rUqZV++9m+ffttE1iPGDGC7du38+qrr/Luu++ydu1a/Pz8yM/Pp0mTJowbN47Ro0djY2NDSEjIPd8NJeo2GS1TmK1bd97c6oOGm9P9TZ48mYiIiHva186dO+nZsyeurq48//zzvPbaa/d9G6c5O3nyJPb29rRr146VK1cC8PLLL5u4KmFqcoYvrEJRURHjx4/HwcGBgIAAiw57uNlf/84772Bvb4+9vT2ffPKJqUsSZkDO8IUQwkrIRVshhLASEvhCCGElJPCFEMJKSOALIYSVkMAXQggrIYEvhBBW4v8B4MmNJpTUKQ8AAAAASUVORK5CYII=\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"def laplace_mech(v, sensitivity, epsilon):\n",
|
|
" return v + np.random.laplace(loc=0, scale=sensitivity/epsilon)\n",
|
|
"\n",
|
|
"epsilon_i = .01\n",
|
|
"plt.plot([laplace_mech(adult['Age'].clip(lower=0, upper=i).sum(), i, epsilon_i) for i in range(100)])\n",
|
|
"plt.xlabel('Clipping Bound for Age')\n",
|
|
"plt.ylabel('Total Sum');"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"The total privacy cost for building this plot is $\\epsilon = 1$ by sequential composition, since we do 100 queries each with $\\epsilon_i = 0.01$. It's clear that the results level off around a value of `upper = 80`, so this is a good choice for the clipping bound.\n",
|
|
"\n",
|
|
"We can use the same approach for data attributes from any numerical domain, but it helps to know something about the scale of the data in advance. For example, trying clipping values between 0 and 100 for yearly incomes would not work very well - we wouldn't even come close to finding a reasonable upper bound.\n",
|
|
"\n",
|
|
"One refinement that can work well when the scale of the data is not known is to test upper bounds according to a logarithmic scale. "
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 22,
|
|
"metadata": {
|
|
"tags": [
|
|
"hide-input"
|
|
]
|
|
},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAESCAYAAADwnNLKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deVzUdf7A8dcMw8AgCMOpNoggJghqireWrWm0nWa1unmkZm1pZa22re6mJmXZWrbbsZW15bXlT9dqNTds10pLTUVFUZDwQMDkHOQcjpn5/UGMotwwB8z7+XjwcOb7nc/3+56PMO/5fr6fQ2E2m80IIYRwWkp7ByCEEMK+JBEIIYSTk0QghBBOThKBEEI4OUkEQgjh5CQRCCGEk+uwiSA1NZXx48ezYcOGRl+XkpLCpEmTmDRpEu+8846NohNCiI6jQyaCsrIy4uLiGDlyZJOvff7554mLi2PLli2kpaVRXl5ugwiFEKLj6JCJQK1Ws2bNGgIDAy3b0tLSmDFjBg899BBz586lqKiIvLw8ysrKiIqKQqlU8vrrr6PRaOwYuRBCOJ4OmQhUKhXu7u51tsXFxbF8+XLWrl3L6NGj2bhxI1lZWfj5+fHCCy/w4IMP8vHHH9snYCGEcGAqewfQXo4dO8bzzz8PQGVlJf3798dsNnPu3Dn++te/4u7uzuTJkxk1ahTXX3+9naMVQgjH0WkSgUajYd26dSgUCsu2jIwM+vTpg1arBSAmJoa0tDRJBEIIcYUO2TRUn4iICHbv3g3Al19+yb59+wgODqa0tJTCwkJMJhPJycmEhYXZOVIhhHAsio44+2hSUhIrV64kKysLlUpFUFAQTz/9NK+99hpKpRI3Nzdee+01fHx8SExMZNWqVVRUVHDjjTfy5JNP2jt8IYRwKB0yEQghhGg/naZpSAghROtIIhBCCCfX4XoNJSQk2DsEIYTokGJiYurd3uESATT8ZpojOTmZyMjIdoymc5H6aZzUT+Okfhpnz/pp7Eu0NA0JIYSTk0QghBBOThKBEEI4OUkEQgjh5CQRCCGEk5NEIIQQTk4SgRBCODlJBEII0Y6+SrpI3PaT9g6jRSQRCCFEO1q//xxr957DaOo483lKIhBCiHZSbTRx9Hwh1SYzP18qt3c4zSaJQAgh2smp7GJKK40AZBRIIhBCCKeTkK63PM7Ul9kxkpaRRCCEEO0kIV2Pv6capQIy9HJFIIQQTichXc+wUF+6dXUns0CuCIQQwqlkFxnI1JczuKcWndaDTLkiEEII51J7f2BIL190vhoy5B6BEEI4l4R0PW4qJf26dyVY68HFIgMV1UZ7h9UskgiEEKIdHErXM1Dng1qlRKfVYDbDz4UGe4fVLJIIhBCijQxVRk5kXSKmlxaAYF8PgA7TPCSJQAgh2uhY5iWqTWZiel6VCDrIoDKrLl6fmprK3LlzmTlzJtOmTbNsz87OZuHChZbnGRkZLFiwALVazapVq+jWrRsAo0aN4vHHH7dmiEII0WaH0gsAGBxSkwi6dXVHpVR0mEFlVksEZWVlxMXFMXLkyGv2BQUFsX79egCqq6uZPn0648aNY+fOnUydOpWZM2daKywhhGh3h9P1hAV0wbeLGgAXpYIePpoOM6jMak1DarWaNWvWEBgY2OjrPvvsM2JjY+nSpQulpaXWCkcIIazCbDaTkK63NAvVCvbVkNFBBpVZLRGoVCrc3d2bfN3mzZu5//77gZqriK+//prZs2cza9YsUlJSrBWeEEK0izN5pejLqogJqZsIdD4dZ1CZVe8RNOXIkSOEhYXh6ekJwIgRIxgwYAAjRozg0KFDPPvss2zbtu2acsnJya0+p8FgaFP5zk7qp3FSP41zxvrZ+VMxAFqjnuTky60absYS8koqOHL8BO6qmu/cjlo/dk0E3377bZ17CAMGDLA8HjJkCAUFBRiNRlxcXOqUi4yMbPU5k5OT21S+s5P6aZzUT+OcsX7WnjyGt8aVCcMHoFQqLNtjKrJYd0SPZ2BP+gR5Afatn4SEhAb32bX76PHjx4mIiLA8f/vtt4mPjwdqehz5+vpekwSEEMKRHErXM7inT50kAKDTagA6RPOQ1a4IkpKSWLlyJVlZWahUKuLj4xk3bhw6nY4JEyYAkJubi5+fn6XMPffcw6JFi1i/fj3V1dW89NJL1gpPCCHarLCskrScEu4ddN01+4K1HWdQmdUSQXR0tKWLaEOubv/X6XRNlhFCCEdx5HwhAIOv6jEEEODlhptK2SF6DsnIYiGEaKVD6QW4KBUMDPa+Zp9CoeA6raZDNA1JIhBCiFZKSNcT1aMrHur6G1eCtR4domlIEoEQQrRCldFEYsalepuFatUMKpMrAiGE6JSSfy6ivMp4zUCyK+m0Hlwqr6LIUGXDyFpOEoEQQrTC5RXJGrki+KXnUKaDXxVIIhBCiFZISNfTw9ud7t6aBl8T7Fuzz9HvE0giEEKIVkhI11umnW6IrvaKwMF7DkkiEEKIFrpQWM7PlwwMaSIRaD1c6aJ2cfixBJIIhBCihWrvD8SE+Db6OoVCQbCvh8MvUCOJQAghWighXY/G1YWI7l5NvlbXAQaVSSIQQogWSkjXc0OwD64uTX+E6rQeZBSUYTabbRBZ60giEEKIFiirrObkz0WNjh+4UrCvB6WVRvRljjuWQBKBEEK0wNGMQowmc7MTweXpqB33PoEkAiGEaIHDv9wobmxqiStZpqN24EFlkgiEEKIFEtL19An0xNvDtVmv13WAQWWSCIQQoplMJjMJ6fpmNwsBdHV3xVvjKk1DQgjRGZzOLaHIUN2iRACOPwupVRevT01NZe7cucycOZNp06bV2Tdx4kS8vC73wV21ahVBQUGsWLGCxMREFAoFixcvrrOgvRBC2NPlgWQtTARaD05lFwNNjzuwB6slgrKyMuLi4hg5cmSDr7l6WcoDBw6Qnp7Opk2bSEtLY9GiRWzevNlaIQohRIscStfj20VNqH+XFpXTaTXsSslx2LEEVmsaUqvVrFmzhsDAwHr3l5aWXrNt3759jB8/HoDw8HCKioooKSmxVohCCNEih9P1DO6pRaFQtKhcsK8HFdUm9OVGK0XWNlZLBCqVCnd39wb3FxYWsmDBAqZMmcLq1asxm83k5eWh1V6+5PLz8yM3N9daIQohRLMVlFZyJq+0xc1CcLkL6cWS6vYOq11Y9R5BY5555hnuvvtu3NzcmDt3Ljt37rzmsslsNtebeZOTk1t9XoPB0KbynZ3UT+OkfhrXmetn//maVgx/LrX4PVYUVgI1g8ocsX7slggefPBBy+Obb76ZU6dOERQURF5enmV7Tk4O/v7+15SNjIxs9XmTk5PbVL6zk/ppnNRP4zpz/XxxLgVXlxzuGj0Qd1eXFpXtVWmELzLJr1DYrX4SEhIa3GeX7qMFBQU88sgjVFXVzL1x8OBB+vTpw+jRo4mPjwfg5MmTBAYG4unpaY8QhRCijsPpeqJ6eLc4CQBo1C74e7qR7WxNQ0lJSaxcuZKsrCxUKhXx8fGMGzcOnU7HhAkTGD58OJMnT0atVtOvXz9iY2NRKpVERUUxZcoUFAoFS5cutVZ4QgjRbJXVJhIzC5k2IqTVxwj21ZBdUtGOUbUfqyWC6Ojoa7qHXmnOnDnMmTPnmu0LFy60VkhCCNEqJy5coqLa1OSKZI3RaT04eNoxe0HKyGIhhGhC7UCyptYobkywVkNuaTVGk+ONJZBEIIQQTUhI16PTagjq2nCX+KYE+3pgNMPPlxxvqglJBEII0Qiz2cyhdH2bmoXgynUJJBEIIUSHkqkvJ7e4olUDya50eV0Cx5uFVBKBEEI0oj3uDwD08NGgADLkikAIITqWhHQ9XdQu9A1q28yhapUSPw8Xh1yXQBKBEEI0IiFdz6CeWlQubf+47ObpSqYDrksgiUAIIRpQUlFNysWiNjcL1QryUjnkkpWSCIQQogFHzxdiMrd8IZqGBHmquFhkoLLa1C7Hay+SCIQQogEJ6XoUChjU06ddjtfN0xWzGS4UOlbzkCQCIYRowKH0AvoGedHV3bVdjhfkWTOrj6M1D0kiEEKIehhNZo6eL2y3+wNwORE42qAySQRCCFGPn3KKKa6obvOI4iv5e6hQKRUON6hMEoEQQtSjdiBZe90oBnBRKujho3G4QWWSCIQQoh4J5/T4e6rp6evRrsfVaTUON6hMEoEQQtQj4byemBBtveumt0Ww1oMMBxtUJolACCGukltcQXp+Wbs2C9UK9tWQV1JBeaWx3Y/dWpIIhBDiKta4P1BL98sspFmFjtM8ZNVEkJqayvjx49mwYcM1+/bv389vfvMbpkyZwqJFizCZTCQlJXHTTTcxffp0pk+fTlxcnDXDE0KIeh0+r0ftoiT6Ou92P3awb826BI7UPGS1NYvLysqIi4tj5MiR9e5fsmQJ69ato1u3bjz11FPs2bMHjUZDbGwsf/rTn6wVlhCtYqgykltcQXKOgfOmi+QWV9T8lFSQ98u/+SWVaFxd8PdS4+/phr+nG36eNY8Drnju56nGTeVi77ckGpGQrqe/ztsq/0+WdQkc6Iax1RKBWq1mzZo1rFmzpt79W7duxdPTEwBfX1/0ej0mk2PNvyE6tyqjifySyl8+0A3kFleQV/v8l5+8kpp/iyuqryh5wfLIr8svH/RebvTs6UFZpZH8kgqOnC8kr6SCsgbagbu6q/D3cvslYVxOHNckDy81Hmqr/Zm2O6PJTElFNV3ULu1+k9VWDFVGjmdeYuboXlY5vr+nG2qV0qEGlVntN0ylUqFSNXz42iSQk5PD3r17mT9/Pnv37iUhIYE5c+ZQXl7Ok08+yYgRI6wVonBCZrOZF7ad5IujWejLqup9jZe7igCvmg/ifj26Wj7oA7zcMOhzGNwvnAAvN3y7qHFtYmrissrqmmTzy5VDXkkleSUV5Jf8knRKKki5WEx+ST6XyuuPx8fDlVmjQnnkplCHTQpGk5nNhzJY+Z/z6MvP4qJU0NVdhbfGla4a17r/utf8W7NNdfmx++XXuSjtl0ROXLhEpdFklfsDAEqlAp1W41CDyuz6W5Wfn89jjz3GkiVL0Gq1REREMG/ePG655RbOnj3LrFmz2LlzJ2q1uk655OTkVp/TYDC0qXxn19nrZ0dqER/vy2NUTw/CfD3RalzQurvU/PvLj7reD/cqoAqDixKXogsUFEFBC86rAYKVENwV6Arg8stWzeUzGM0UGow1P+U1/+rLjaTkGlj931Q+/uE002/Qcmu4l10/KK9kNps5lFXOPxLyOVdYRYSfK/f286a00kRJpYnSSiMllZXk6A2cza7ZVlJppKnJNzWuCrzULniqlXi6Kenq5sLk/j6E+7lZ/T3tSCoEwKsil+Rkfbseu/bvy1dt4qcLBQ7zt2a3RFBSUsIjjzzC/PnzGTNmDAC9e/emd+/eAISGhuLv7092djbBwcF1ykZGRrb6vMnJyW0q39l15vpJzS7m/YPfc2Mff9bOGoayFR+m9qqfhPQCVuxI4W/78vjPaQN//HUE4yIC7dr8cuLCJV7ekcL3aXmE+Hnw96n9CXHR069fv0bLmc1mDFUmLpVXWX6KrnxsuHJ7NUXlVZzIKeaFb3P5bO5ogtt5gNfVMg4eIsTPg1GD+7f7sWt/f/qmVPPl8Z9t+ruUkJDQ4D67JYJXXnmFhx56iLFjx1q2bdmyhbKyMmbMmEFubi75+fkEBQXZK0TRiRiqjDz5zyN4uat47TcDW5UE7CkmxJctj40k/kQ2r36VwsNrDzE81JfFt0cyMLh9pkhurp8vlfPazlT+dTgTb40rS+/qx9ThIahVSpKTC5ssr1Ao0Khd0Khd6Obt3qxzpuWUMOmdH5j98UG2PD4Kb037zAZ6NbPZzOHzem66PsAqx68V7OtBYVkVxYYqvNppZtO2sFoiSEpKYuXKlWRlZaFSqYiPj2fcuHHodDrGjBnD559/Tnp6Olu2bAHgzjvv5LbbbmPhwoXEx8dTWVnJsmXLrmkWEqI1XvzyJKeyi1k7exiBXs378HE0CoWC26K7cUtkIJ8ezOCv/03lnrd/4K6BPXj21r709LPuN+WSimre++40a/acwWSCR28MY+6vwq32oXyl8EBP3ps+hBn/+JF5Gw/z0ayhTd6faY3zBWXklVRa7f5ALUvPoYJy+vXoxIkgOjqa9evXN7g/KSmp3u0N9TISorW+SvqZDfvP8+hNYYy18jc9W3B1UTJ9RAj3DrqO9787zZo9Z/kq6Wemj+jFk+PC0XZp3y9P1UYTnx7M4I3/ppJXUsndA3vwbGxfqzfRXG1kbz9enjSAhZsT+fNnSbxyX/92bxo7dM56A8mupNPW3BvK1JfRr0dXq56rORyzC4IQ7SSrsJw/bDnGQJ03C2/ta+9w2pWnm4rf39qXqSNCWP11Kh/vPcvmhAzm/SqcmaN64e7atj7wZrOZXSk5vPyfFNJyShjWy5cPHorkBhs3RV3p/hgd5/NL+duuNEL8PZh7c3i7Hj/hvB4vNxXXB3q163GvVptEHWUWUkkEotOqNpqY/8kRTGb4228HoVZ1zhlVgrq688p9A5g9JpSV/0nhlf+ksG7vORbG9mXiDde16n5IUtYlXvoymX1n8gn178J702O4tV+QQ4wNeGbC9aQXlPHqV6fo6evBnQN6tNuxD6frGRSitfo9JK2HK13ULg7ThVQSgei0/rYrjUPpet6YfAMhfl3sHY7VXR/kxYczh7L3dB4v70jh9/+XyAd7zrLo9ghu7NO8JrGswnJeiz/F1iNZ+HZR88LdUTw4vKdV2uNbS6FQ8Or9A7hQWM7v/y+R7t6admnKuVRexansYn4d3b0domycQqFAp/VwmEFljvO/K0Q72n8mn7d2/cR9g3VMHHSdvcOxqVG9/fli3mj+9ttBFBmqmP7hAWb84wAnLxQ1WKbIUMXKr1L41apv2X78Zx6/uTffPnszD43q5VBJoJabyoX3pw/hOh8Nj6w7RHp+aZuPeTSjELMZhvSy7v2BWsG+jrMuQZNXBMeOHePLL7+kuLgYs9ls2f7yyy9bNTAhWktfWsnTnx6ll18Xlt8TZe9w7EKpVHD3wB7ERgWxfl86b+5K44439zBpkI4Ft15PD5+am5VVRhOfHDjPG//9iYLSSu4ddB0Lbr3eMkOmI9N2UfOPmUO5950fmPXxQbY+Pgofj9bfKE9I16NUYLPuuDqtB/tO52M2m+3e5NZkInj22Wd55JFH8Pf3t0U8QrSJ2Wzm2S2JFJRW8sFDo+ji5tytn24qF+bcGMYDMcG8820aH+09x/ZjF5g9JpToHt68tvMUZ/JKGRHmy59u70d/XfvPtmlNof5deH/6EKZ98CO/W5/A+oeHt/peUEJ6ARHduuJpo98ZnVZDaaWRwrKqdu/p1VJNvuOwsDDuu+8+u2csIZpj7d5z/Dc5hyV39rPKFMIdlbeHK4tuj2T6yBBe35nK3789DUDvgC58MGMIt0Tad5RyWwwL9eUvDwxg/qdH+ePWY7z2wMAWv5dqo4mj5wu5L0ZnpSivdbnnUJnjJ4I777yTiRMn0rdvX1xcLndHk6Yh4WhOXLjEih0pjIsIZJaVZo7s6HRaD16ffAMP3xhKRkEZ4yODUDngPYCWuueG6zifX8ZrX6fSy68LT93Sp0XlT2UXU1pptPr4gStdOahsgM5+XXKhGYngjTfe4NFHHyUgoOMPxBGdV1llNU9+cgQfD1f+cv+ADvvt1laiengT1aNzXTE9MS6cc/llvP51Kj19PVrUSaB2RbLBPW2XCHS+lweV2VuTiaB379488MADtohFiFZb9u8TnM0rZeOc4fh5Wn+GSuF4FAoFL0/qT1ZhGX/YcowePhqGhfo2q2xCup6grm6WEb+2UDvttiMsUNNkItBqtUydOpXo6Og6TUN/+MMfrBqYEM3178QL/N+hTJ74VTijekunBmemVil5b9oQ7v37Dzy6/hBbHx9FWIBnk+US0vXEhGhtfiUZ7KtxiCUrm2wcHDZsGPfffz8RERH06dPH8iOEIzifX8birccZ3NOH+ePl91LU3Bj/eOYwXBQKZn98kILSykZfn11kIFNfbtNmoVrBWg+HaBpq1l0ihUJxzY8Q9lZlNPHkp0dQKOCvUwY55MAnYR89/Tx4f8YQLlwy8Oi6Qxiq6l8yFC7fHxjSq3nNSO1Jp9WQqS+vM0bLHpr8y0lNTeXUqVOcOnWKEydO8M9//pODBw/aIjYhGvXazlQSMwpZed8Am8+EKRxfTIiW1b+5gUPpev6w5ViDH7YJ6XrcVEr6dbf9LKDBvh5UVJvILa6w+bmv1OQ9gueee67Oc6PRyFNPPWW1gIRojj0/5fLud6f57bCe3N7f+nPDiI7pjgHdOV8QwcqvUgjx82BBPTPQHkrXM1DnY5dJCS1dSPXlBHa13zoZTb7z8vLyOj9ZWVmcOXPGFrEJUa/c4gqe2ZRIn0BPltzZ+LKIQjw2NowpQ4N5c1camw9l1NlnqDJyIusSMTaaX+hqV65LYE9NXhHccccdKBQKy3wYXl5ezJ492xaxCXENk8nMgs2JFBuq2DBnGBp12+bcF52fQqEgbmI0mfpyFm09znU+GkaF1/QuO5Z5iWqTmRg73CgGLHM62Xs66iYTwa5duyyPTSYTSmXzL59SU1OZO3cuM2fOZNq0aXX27d27l9dffx0XFxduuukm5s2bB8CKFStITExEoVCwePFiBgwY0Ozzic7vg+/PsDs1l7iJ0UR0s//KTqJjcHVR8s60wdz3zl4e25DA1rmjCA/04lB6AQCDbTii+EoatQv+nm5270La4Kd6SkoKCxcutDxftGgRw4cP5/bbb+f48eNNHrisrIy4uDhGjhxZ7/4XX3yRN998k08++YQ9e/aQlpbGgQMHSE9PZ9OmTbz44ovExcW14i2Jzioxo5BXvzpFbFQQ04b3tHc4ooPp6u7KR7OGola5MOvjg+SVVHA4XU9YQBd87TjXj06rIbPQvlcEDSaC5cuX85vf/AaA3bt3k5iYyO7du3n//fdZtWpVkwdWq9WsWbOGwMDAa/ZlZGTg7e1N9+7dUSqVjB07ln379rFv3z7Gjx8PQHh4OEVFRZSUlLT2vYlOpNhQxVOfHiHQy42V98kUEqJ1dFoPPnxoCLnFFTyy7lDNQDI7NQvVCvb1cNwrAhcXF4YNGwbA//73PyZOnIhGo0Gn0zXrj1ClUuHuXv9d8NzcXHx9L/fZ9ff3Jzc3l7y8PLTay/8pfn5+5ObmNvvNiM7JbDbz/OdJZBSU8dffDmrTnPNCDAz24Y3JgziaUYi+rMqmE83VJ1ir4UJhOUaT/cYSNHiPoLKyZjSe0Whkz549vPnmm5Z9BoOhTSetrz9v7Q3pq19XX9JJTk5u9bkNBkObynd2jlg/X6cV8/nRXKbdoMWzPJvk5Gy7xeKI9eNIOkr9hLjAIzG+/ONwAYEUkpzc9hXOmqO++nGpKKLaZOb7hCQCPe2zfkaDZx01ahSPPfYY5eXl9OrVi6ioKKqrq3nrrbcICwtr00mDgoLIy8uzPM/OziYgIACVSlVne05OTr0L4kRGRrb63MnJyW0q39k5Wv2cyS3h3U++Z3ioLy/8ZgQuVl5UvCmOVj+OpiPVT2Qk/P4eI+6utut5Vl/95KlyYV8ebn49iAzzs9q5ExISGtzXYNPQ/Pnzefjhh5kxYwbvvvtuzYt/6TG0dOnSNgWk0+koKSkhMzOT6upqvvnmG0aPHs3o0aOJj48H4OTJkwQGBuLp2fSEUaJzqqg28uQnR1CrlLwx5Qa7JwHR+dgyCTQk2AG6kDZ6HTJ06NA6z5VKJU8//XSzDpyUlMTKlSvJyspCpVIRHx/PuHHj0Ol0TJgwgWXLlrFgwQIAbr/9dkJDQwkNDSUqKoopU6agUCjanHBEx/but2c4caGINTOG0N3bdtMDC2FL3X3cUSggU2+/G8ZWa5CKjo5m/fr1De4fOnQomzZtumb7lV1WhfMqq6zm471nGR8ZyIR+QfYORwircVO50K2ru13XJZDpGoVD2nwoE31ZFY+N7W3vUISwumCtB5l27ELa4BXBxo0bGy04derUdg9GCKhZSHzNnjPEhGjtMjWwELam02rYfybfbudvMBHo9XpbxiGExX+SLpKpL+d5mVBOOAmdrwc/H82istpkl1lQG0wETzzxhOVxaWkply5dAmrGFyxfvtz6kQmnZDabeW/3acL8uzAhUu4NCOcQrNVgNsOFwnJ6+Xex+fmbvFn89ttvs3XrVgoLC+nRowcXLlxg8uTJtohNOKG9p/NJyiri5Un9UUp3UeEkamchzdTbJxE0eQ2ye/du/ve//9GvXz+2bdvGunXr6ixiL0R7em/3Gfw93bh30HX2DkUImwn2rekeba+eQ00mgtqpH4xGIwaDgaioqEZHqAnRWicvFLE7NZdZo3s5xEAfIWylW1d3VEqF3QaVNdk0FBsby9q1a7nrrru455578PPzQ6ORwT2i/b2/+zQeahemDQ+xdyhC2JTKRUl3H3e7DSprMhH8+te/plu3bgCMHTsWvV7f4KyiQrRWpr6Mbcd+ZuaoXnh7uNo7HCFsLljr4XhNQwUFBfz00088+eSTnD59mrS0NMrKaoKcO3euzQIUzuEf359DAcweE2rvUISwi2Ct/dYlaPCK4MyZM/zrX//i3LlzLFu2zLJdqVRy11132SI24SQulVXx6cHz3DWwB9f5SLOjcE46rYa8kgoMVbadERUaSQRDhgxhyJAh3HXXXYwaNYpLly7h4uIis4GKdrfhx3TKKo08elPbpjcXoiML9q3tQlpGeKCXTc/d5D0Cs9lMbGwsKpXKsnj98uXLiYmJsUV8opMzVBn56IdzjL0+gMjushi9cF6WLqQF5Y6XCN58803Wr19vWXv4559/ZsGCBfzzn/+0enCi8/vsSBZ5JRX8Tq4GhJO7PKjM9jeMmxxH4OrqWmcB+u7du6NS2Wc5NdG5GE1m1uw+Q//rvBnZ23orMwnREQR4uqFWKcmwQxfSJj/RdTodL7zwAsOGDcNsNvPjjz/Ss2dPW8QmOrmvT2ZzJq+UN387qN61qYVwJkqlAp1WY5dBZQ1eETz11FMAxMXFMXDgQL0hNVoAABtOSURBVBISEjhy5AgxMTG88MILNgtQdE61k8sF+2r4dXQ3e4cjhEPQaT3sMqiswSuCwsLCmheoVEycOJGJEyfaLCjR+R1K13PkfCHL74lC5SLrIwkBNbOQHssstPl5G0wE58+f59VXX22w4B/+8IcmD75ixQoSExNRKBQsXryYAQMGAJCdnV1nScqMjAwWLFiAWq1m1apVlpHMo0aN4vHHH2/2mxEdx3vfnUbr4coDMcH2DkUIhxHs60FhWRXFhiq83G03wr7BRKDRaOjTp0+rD3zgwAHS09PZtGkTaWlpLFq0iM2bNwMQFBRkWc+4urqa6dOnM27cOHbu3MnUqVOZOXNmq88rHF9aTjH/Tc5h/i190Khlcjkhaum0NV1IM/XlRHZ3gETg7+/Pvffe2+oD79u3j/HjxwMQHh5OUVERJSUl1wxI++yzz4iNjaVLly6Ulpa2+nyi43h/9xncXZXMGCmTywlxpeBfupBmFJTZdFxNg42z0dHRbTpwXl4eWq3W8tzPz4/c3NxrXrd582buv/9+AMrKyvj666+ZPXs2s2bNIiUlpU0xCMeTXWTgsyNZPBATjJ+nm73DEcKh1I4utnUX0gavCJ577rk2HdhsNl/z/OougkeOHCEsLMxylTBixAgGDBjAiBEjOHToEM8++yzbtm275tjJycmtjstgMLSpfGdn7fr5R0I+RpOZm7ubOuT/g/z+NE7qp3FN1Y/ZbMZdpeDY6UyS/Qw2i8tqI8OCgoLIy8uzPM/JycHf37/Oa7799ltGjhxpeV57Mxlq5joqKCjAaDResyJaZGRkq+NKTk5uU/nOzpr1U2yo4qtPz/Pr/t25ZfiApgs4IPn9aZzUT+OaUz8hfnmUomn3emxsQTGr9dsbPXo08fHxAJw8eZLAwMBr7g8cP36ciIgIy/O3337bUiY1NRVfX19ZFrMT+eTAeYorqmU6CSEaEeyrsfk0E1a7Ihg8eDBRUVFMmTIFhULB0qVL2bp1K15eXkyYMAGA3Nxc/PwuTy1wzz33sGjRItavX091dTUvvfSStcITNlZZbeIf359jZJgfA3Q+9g5HCIel03qw/0xBvc3p1mLVSYOuHCsA1Pn2D1zT/q/T6SzdSkXn8u/EC1wsMvDyff3tHYoQDk2n1VBSUU1hWRXaLmqbnFOGdAqrM5vNvL/7NH2DvLj5+gB7hyOEQ7vcc8h2zUOSCITVfXsql9TsEn43NkwmlxOiCcGW6aht14VUEoGwune/O013b3fuGtjD3qEI4fB0lgVq5IpAdBJHMwr58WwBD48JxVUmlxOiSV3dXfHWuErTkOg83t99Gi93FVOGyRoWQjRXTRdSaRoSncC5vFL+k3SRaSNC8HSTVe2EaC6dj4c0DYnO4YPvz+CqVDJrVC97hyJEh1J7RXD1VD3WIolAWEVeSQWbD2UyafB1BHZ1t3c4QnQowb4eVFSbyC2usMn5JBEIq1i39xwV1Sbm3CjTSQjRUrXrEthqFlJJBKLdlVVWs25/OhP6BREe6Nl0ASFEHZfHEtjmPoEkAtHu/u9gBoVlVTK5nBCtpLtigRpbkEQg2lW10cQH358lJkTLkF6+9g5HiA5Jo3bB31Ntsy6kkghEu9qRdJFMfblcDQjRRjqth80GlUkiEO3GbDbz3nenCQvowvjIIHuHI0SHFuzrQUaBXBGIDuaHtHxOXCji0RvDUCplcjkh2kKn1XChsByjyfpjCSQRiHbz3u7TBHi5MXHQdfYORYgOL1jrQbXJzMUi669dLIlAtIsTFy6x56c8Zo7qhburLC8qRFsF23AWUqtOALNixQoSExNRKBQsXry4zuL0EydOxMvLy/J81apVBAUFNVpGOK41u8/QRe3CtOEh9g5FiE5BZ8N1CayWCA4cOEB6ejqbNm0iLS2NRYsWsXnz5jqvuXpZyuaUEY7nTG4J2479zKxRvfD2cLV3OEJ0Cj183FEobHNFYLWmoX379jF+/HgAwsPDKSoqoqSkxLK/tLS0xWWE4ymrrGbuxsN4uatkOgkh2pGbyoVuXd1t0oXUaokgLy8PrVZree7n50dubq7leWFhIQsWLGDKlCmsXr0as9ncZBnhWMxmM8/96zinsov525RBdPOWyeWEaE86rW3WJbBa09DV06eazeY669U+88wz3H333bi5uTF37lx27tzZZJlaycnJrY7LYDC0qXxn15L6+exEIdsSC5g5WEuAMY/k5DwrR2d/8vvTOKmfxrW0froqqziWXW71OrVaIggKCiIv7/IHQ05ODv7+/pbnDz74oOXxzTffzKlTp5osUysyMrLVcSUnJ7epfGfX3PrZdzqfDw+fJTYqiKUPxDjNovTy+9M4qZ/GtbR+ojJd2HX2J3r36Yta1bYGnISEhAb3Wa1paPTo0cTHxwNw8uRJAgMD8fSsmYmyoKCARx55hKqqKgAOHjxInz59Gi0jHMfPl8p54p+HCfHzYNUDA50mCQhhazqtBrO55m/Omqx2RTB48GCioqKYMmUKCoWCpUuXsnXrVry8vJgwYQLDhw9n8uTJqNVq+vXrR2xsLEql8poywrFUVBt5bMNhDFVGNk0fgZe79BISwlqCLbOQlhPi18Vq57HqOIKFCxfWeR4REWF5PGfOHObMmdNkGeFYlv37JIkZhbw7bTDhgV5NFxBCtJplUJmVew7JyGLRbJ8eOM8nB87z+M29uS26u73DEaLT69bVHRelwuoL1EgiEM1yNKOQJV+cYEy4Pwtv7WvvcIRwCioXJT183K0+C6kkAtGkvJIKHt+QQICXG2/+dhAuMrOoEDYTbIN1CSQRiEZVG008+c8jFJRW8t70GLRd1PYOSQinYotBZZIIRKNejT/FvjP5vHRvf6Kv87Z3OEI4nWCtB7nFFRiqjFY7hyQC0aDtxy7w/u4zzBgZwv0xOnuHI4RTCvatnYXUes1DkghEvU5dLOYPW44RE6Llz3f0s3c4Qjgtnba2C6n1mockEYhrXCqv4rENCXRxU/HO1MFtHtouhGg9yxWBFaejlr9wUYfJbOb3m46SUVDGO1MHE9RVZhQVwp4CPN1Qq5RyRSBs55NjhfwvJYfn7+zH0F6+9g5HCKenVCrQ+WjkHoGwjV0p2Ww8qmfSoOuYMVKWnBTCUeh8Paw6qEwSgQDgXF4pT396lFBfNS/d219mFBXCgQRrNVYdVCaJQFBWWc1jGxJQKBQ8f3MQGrWLvUMSQlxBp/WgsKyKYkOVVY4vicDJmc1m/li73ORvB9HNS6aVFsLR1M5Caq0RxpIInNw/fjjHvxMvsPDWvoy9PsDe4Qgh6nF5XQLrNA9JInBi+8/ks2JHMrFRQcy9ube9wxFCNODy6GK5IhDtSJabFKLj0Hq44qF2sdoNY0kETqii2sjjGw5TXmnk/ekxstykEA5OoVAwqrc/vh7Wmf3XqktVrlixgsTERBQKBYsXL2bAgAGWffv37+f1119HqVQSGhrKSy+9xMmTJ5k7dy4hITV92K+//nqef/55a4bolF7YdpKjstykEB3KBw8NsdqxrZYIDhw4QHp6Ops2bSItLY1FixaxefNmy/4lS5awbt06unXrxlNPPcWePXvQaDTExsbypz/9yVphOb1NB8/zzx9luUkhxGVWSwT79u1j/PjxAISHh1NUVERJSQmenp4AbN261fLY19cXvV6PyWSyVjhOr6SimlXxp1i775wsNymEqMNq9wjy8vLQarWW535+fuTm5lqe1yaBnJwc9u7dy9ixYykrKyMhIYE5c+YwdepU9u/fb63wnMrXJ7OZ8Pp3rN13jmnDQ3h3eowsNymEsLDaFYHZbL7m+dU9U/Lz83nsscdYsmQJWq2WiIgI5s2bxy233MLZs2eZNWsWO3fuRK2ue4MkOTm51XEZDIY2le9I8suq+fuBfH5IL6WXjyuv3daDyEAXMs781GAZZ6qf1pD6aZzUT+MctX6slgiCgoLIy8uzPM/JycHf39/yvKSkhEceeYT58+czZswYAHr37k3v3jX92UNDQ/H39yc7O5vg4OA6x46MjGx1XMnJyW0q3xGYTGY2HjjPq/9JocJo4tnYvjxyY1iz1hVwhvppC6mfxkn9NM6e9ZOQkNDgPqs1DY0ePZr4+HgATp48SWBgoKU5COCVV17hoYceYuzYsZZtW7ZsYd26dQDk5uaSn59PUFCQtULslE5dLOaB9/bx/OdJ9Nd5E//0Tcz7VbgsLiOEaJDVrggGDx5MVFQUU6ZMQaFQsHTpUrZu3YqXlxdjxozh888/Jz09nS1btgBw5513ctttt7Fw4ULi4+OprKxk2bJl1zQLifoZqoy8tSuNd787jZe7itceGMikwdfJQDEhRJOsOo5g4cKFdZ5HRERYHiclJdVbZs2aNdYMqVPam5bH4s+Ocy6/jEmDr+PPd/TDt4skUCFE81g1EQjr0pdW8tKOZLYkZBLi58GGh4czpo9/0wWFEOIKkgg6ILPZzOdHs4jbnkxReRVzb+7NU7f0wd1V1hEQQrScJIIOJj2/lD9/nsSen/K4IdiHlyf1J7J7V3uHJYTowCQRdBBVRhMf7DnLG/9NxdVFyfJ7opg6PEQGhgkh2kwSQQdw5LyeRVuPk3KxmNioIJbdHUV3b429wxJCdBKSCBxYsaGK13amsnbfOQK93Hh3Wgy3RXezd1hCiE5GEoEDKiyrJP7ERVZ//RPZxQZmjAhhYWxfWTdACGEVkggcRHaRgZ0nLvLViYvsP1OA0WQmopsX70wbzOCe2qYPIIQQrSSJwI7O5ZUS/8uH/5HzhQCE+XfhdzeFERvVjQE6bxkZLISwOkkENmQ2m0m5WMxXSReJP3GRlIvFAERf15WFt15PbFQ3wgM95cNfCGFTkgiszGQycySjsOabf9JFzheUoVDA0BBfnr+zH7f2CyLY18PeYQohnJgkAiuoMpr48UwBX534mZ0nsskprsDVpWbx6cdv7s34yCACvNzsHaYQQgCSCNqNocrI7tRcvjpxkf8l53CpvAqNqws39w3gtuhu/CoikK7S60cI4YAkEbSC0WQmS19OWm4xp3NKOXxez7encimvMtLVXcX4fkHcFtWNm64PkPl/hBAOTxJBI8orjZzJKyEtp4TTuaWczi3hdE4JZ/JKqaw2WV7Xras798Vcx21R3Rke5ouriywCI4ToOJw+EZjNZvJLKzmdU0Jabgmnc2o+8NNySsgqLLe8TqmAYF8PwgM8uen6AHoHdCE80JMwf0+0Mve/EKIDc6pE8POlcn7MKGVP7uk63/ILy6osr9G4uhAW0IUhvbRMDgimd4An4YGehPh5SDOPEKJTsmoiWLFiBYmJiSgUChYvXsyAAQMs+/bu3cvrr7+Oi4sLN910E/PmzWuyTFvkl1Qw9tVvqTSagGz8Pd3oHdCFO/p3p3eAJ70Daz7wu3d1RykzegohnIjVEsGBAwdIT09n06ZNpKWlsWjRIjZv3mzZ/+KLL/Lhhx8SFBTEgw8+SGxsLAUFBY2WaQuth5p/zBxKzoUMbhkajbeH9OARQgiwYiLYt28f48ePByA8PJyioiJKSkrw9PQkIyMDb29vunfvDsDYsWPZt28fBQUFDZZpK6VSwZg+/iRX50oSEEKIK1ite0teXh5a7eXJ0vz8/MjNzQUgNzcXX19fyz5/f39yc3MbLSOEEMI6rHZFYDabr3leO4fO1fsAFApFo2WulJyc3Oq4DAZDm8p3dlI/jZP6aZzUT+MctX6slgiCgoLIy8uzPM/JycHf37/efdnZ2QQEBKBSqRosc6XIyMhWx5WcnNym8p2d1E/jpH4aJ/XTOHvWT0JCQoP7rNY0NHr0aOLj4wE4efIkgYGBlrZ+nU5HSUkJmZmZVFdX88033zB69OhGywghhLAOq10RDB48mKioKKZMmYJCoWDp0qVs3boVLy8vJkyYwLJly1iwYAEAt99+O6GhoYSGhl5TRgghhHVZdRzBwoUL6zyPiIiwPB46dCibNm1qsowQQgjrkklxhBDCyUkiEEIIJ6cw19eX04E1dudbCCFEw2JiYurd3uESgRBCiPYlTUNCCOHkJBEIIYSTk0QghBBOzqkTwbFjx1i8eDGLFi0iKyvL3uE4nJycHObPn99uU4F3RkeOHGHx4sU899xzJCUl2Tsch5OQkMCzzz7L008/zfHjx+0djkPKzc1lzJgxVFdX2y2GTpkIUlNTGT9+PBs2bLBsW7FiBZMnT2bKlCkcO3YMgM2bN7Ns2TLmzp3LZ599Zq9wba659aNUKpk8ebK9wrSr5taRRqNh6dKlzJw5k0OHDtkrXJtrbv14enry4osvMnv2bA4cOGCvcG2uufUD8NFHHzF06FB7hGnR6RJBWVkZcXFxjBw50rLtykVyXnzxReLi4gCorKxErVYTEBBATk6OvUK2qZbUj7+/Py4uzrc8Z0vqKCIigqqqKj755BMmTpxor5BtqiX107dvX/bv38+qVauYMGGCvUK2qZbUzxdffMGtt96Km5ubvcIFOmEiUKvVrFmzhsDAQMu2hhbJcXNzo6KigosXL9KtWzd7hWxTLakfZ9WSOiouLubVV1/lmWeewcfHx14h21RL6icxMZGbbrqJN954g48//thOEdtWS+tnz549JCcn8+WXX9or5M63eL1KpUKlqvu28vLyiIqKsjyvXfBm8uTJLFu2DKPRyO9//3tbh2oXLamf48eP88knn1BcXIyPj4/TfKNrSR199tlnlJaW8s477zBkyBBiY2NtHa7NtaR+CgsLWbJkCWVlZdx99922DtUuWlI/S5YsASArK4s77rjDpnFeqdMlgvo0tOBNVFQUL7/8sp2ichwN1c/IkSPrXN46s4bqyFm+QDSlofoZO3YsY8eOtVNUjqOpRbdeeeUVW4dUR6drGqpPY4vkCKmf5pA6apzUT+McvX6cIhHIgjeNk/ppmtRR46R+Gufo9dPpmoaSkpJYuXIlWVlZqFQq4uPjefPNN2XBm19I/TRN6qhxUj+N64j1I5POCSGEk3OKpiEhhBANk0QghBBOThKBEEI4OUkEQgjh5CQRCCGEk5NEIIQQTk4SgXAI586d49FHH+X+++9n0qRJxMXFUVlZSWZmJpMmTQLgmWeewWAwtOi477//PkeOHGl1XFu3bmXs2LFMnz6dadOmMXv27DojRNtLaWkp48aNu2b78uXLuffee1s1CeDhw4fp27cvKSkp7RGi6MQkEQi7MxqNPPnkk8yZM4ctW7bwr3/9C4C33367zutWr16Nu7t7i4796KOPMmjQoDbFd/vtt7N+/Xo2bNjA4MGDLfHZwnfffcfatWtbNQp1+/bthIaGsn37ditEJjqTTjeyWHQ8P/zwA2FhYQwbNgwAhULBs88+i1KprLNOxLhx49i2bRtxcXF4eHhw5swZ9Ho9L7/8Ml27dmX+/PmEhoZy9uxZ+vfvz7Jly/jjH/9IbGwser2ehIQECgoKOHv2LA8//DAPPPAAn3/+OR9++CE9evQgKCiIG264wXIFUp/8/HwGDhwIwI8//sjq1atRqVQEBQXx8ssvs337dn766Seee+45SktLueuuu9i1axcTJkxg8uTJfPPNN1RWVvLRRx8B8OSTTwIwYMCAa871wQcfkJOTw2OPPcZ7773H3//+dw4fPozRaGTq1KlMnDiR6dOn06dPHwDLTJZQk1x37tzJ6tWree6551i4cCEAKSkp/PGPf8TLy4uhQ4dy4cIFXnnlFTZu3Mi2bdtQKpWMHz+e2bNnt+W/VHQwckUg7O7MmTNERkbW2ebu7o5arW6wTHV1NR9//DHz58+3XDmcOnWKBQsWsGXLFo4fP35Nk0hqaipvvfUWb7/9Nhs2bMBkMvH666/z0UcfsXr1ar7//vt6z7Vjxw6mT5/OnXfeycmTJy1TTS9dupTVq1ezYcMGvL292bZtW4PxGo1GwsLC2LhxIzqdjv379/PFF1/Qp08fPvroI/r27XtNmTlz5hAQEMCaNWtISUnhp59+4tNPP2Xt2rW89dZbluaiPn361EkCUJNce/fuzdChQ/Hx8bE0j7399tvMmzeP9evXc/r0aQAyMjL46quv+OSTT9i4cSM7d+7kwoULDb4X0flIIhAOwWg0tuj1o0aNAuCGG27g7NmzAPTq1Yvu3bujUCgYOHAgZ86cqVPmhhtuwMXFhW7dulFcXIxer8fLywt/f388PDwYMWJEveeqbRravn07v/3tb1myZAmFhYUoFAq6d+8OwJAhQ0hOTm405iFDhgBYzn/69GlLs1Xt1VBDkpKSLMsZenh40KtXL9LT04H6rya2b9/OnXfeCcBdd91lWfTk9OnTxMTEAPCrX/0KgOPHj5Oens6MGTOYMWMGpaWlsoa3k5GmIWF3vXv3ZuPGjXW2VVZWcu7cOTw8POotYzKZLI9r53W/ctvV870D1ywWcvVrlMqmvxfFxsby17/+FYVCUWeOeZPJhEKhqHO8qxcjv3LZT7PZjNlstpzzytjrc/V7ubKsq6trnX0Gg4FvvvmGkydPsmHDBqqqqigqKmLx4sV1Yr6y/M0338zy5cubfP+ic5IrAmF3o0ePJisri127dgE1H4p/+ctf2LFjR4NlDh8+DMCRI0fo3bs3AOfPnycnJweTyURiYiLh4eGNntfHx4fCwkIuXbqEwWBo1uLqiYmJhIaG4u3tjUKhsDShHDhwgOjoaDw9PS33NRISEho9VmhoKElJSUDN/YbGREdHW15TWlrK+fPnCQkJqfe1u3btYsSIEWzfvp0vvviCHTt2EBYWxo8//kjPnj0t59y9ezcAUVFR/Pjjj5SXl2M2m3nxxRdb3DtLdGxyRSDsTqlU8uGHH7JkyRLeeust1Go1o0aN4oknnmiwrdpgMPC73/2Oixcv8uqrrwI1H6yrV68mLS2NwYMHW26iNkSlUvH4448zdepUQkJCiI6OrvOtvdaOHTssH54Ay5YtAyAuLo4FCxagUqnQ6XTccccdGAwG/v73vzN9+nTGjh17zTf5K02cOJF58+bx0EMPWZprGjJkyBCio6OZOnUq1dXVLFiwoMGrpS+//JL777+/zrZJkybx5Zdf8vjjj/PnP/+ZtWvXEh4eTklJCT169GDGjBlMnToVFxcXxo8f3+LeWaJjk2moRYdT2xOoto0bIDMzk6eeeoqtW7e26FhfffUVI0aMwMfHh4cffpgnnniizd1NHdnRo0dxd3cnIiKC9957D4Df/e53do5K2JtcEQinVl5ezkMPPYRGoyEyMrJTJwGouR/wpz/9CXd3d9zd3XnttdfsHZJwAHJFIIQQTk5uFgshhJOTRCCEEE5OEoEQQjg5SQRCCOHkJBEIIYSTk0QghBBO7v8BZQvHvqnNTnwAAAAASUVORK5CYII=\n",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 1 Axes>"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"xs = [2**i for i in range(15)]\n",
|
|
"plt.plot(xs, [laplace_mech(adult['Age'].clip(lower=0, upper=i).sum(), i, epsilon_i) for i in xs])\n",
|
|
"plt.xscale('log')\n",
|
|
"plt.xlabel('Clipping Bound for Age')\n",
|
|
"plt.ylabel('Total Sum');"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "markdown",
|
|
"metadata": {},
|
|
"source": [
|
|
"This approach allows us to test a huge range of possible bounds with a small number of queries, but at the expense of less precision in determining the perfect bound. As the upper bound gets really large, the noise will start to overwhelm the signal - notice how the sum fluctuates wildly for the largest clipping parameters! The key is to look for a region of the graph which is relatively smooth (meaning low noise) and also not increasing (meaning the clipping bound is sufficient). Here, this occurs at roughly $2^8 = 256$, which is a reasonable approximation of the upper bound we derived earlier."
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
}
|
|
],
|
|
"metadata": {
|
|
"celltoolbar": "Tags",
|
|
"kernelspec": {
|
|
"display_name": "Python 3",
|
|
"language": "python",
|
|
"name": "python3"
|
|
},
|
|
"language_info": {
|
|
"codemirror_mode": {
|
|
"name": "ipython",
|
|
"version": 3
|
|
},
|
|
"file_extension": ".py",
|
|
"mimetype": "text/x-python",
|
|
"name": "python",
|
|
"nbconvert_exporter": "python",
|
|
"pygments_lexer": "ipython3",
|
|
"version": "3.9.9"
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 2
|
|
}
|