# Convert mu = 1 GDP to (epsilon, delta)-DP with epsilon = 1
dp_guarantee <- gdp_to_epsdelta(mu = 1.0, epsilon = 1.0)
dp_guarantee
# Round delta to 6 decimal places for reporting
dp_rounded <- gdp_to_epsdelta(mu = 1.0, epsilon = 1.0, dp = 6)
dp_rounded
# Compare the original GDP with its (epsilon, delta)-DP representation
fdp(gdp(1.0),
gdp_to_epsdelta(mu = 1.0, epsilon = 1.0),
.legend = "Privacy Mechanism")
# Explore how delta varies with epsilon for a fixed mu
mu_fixed <- 1.0
epsilons <- c(0.1, 0.5, 1.0, 2.0)
res <- fdp(gdp(mu_fixed))
for (eps in epsilons) {
res <- res+fdp(gdp_to_epsdelta(mu = mu_fixed, epsilon = eps))
}
res
Run the code above in your browser using DataLab