EC3133
# True tax burden (progressive with income)
true_burden <- 0.15 * income + rnorm(n, 0, 100)
# Flat poll tax (same for everyone)
poll_tax <- mean(true_burden)
# Calculate distributional impact
impact <- data.frame(
Income_Quintile = cut(income, breaks = 5, labels = c("Lowest", "Lower", "Middle", "Higher", "Highest")),
True_Burden = true_burden,
Poll_Tax = poll_tax,
Change = poll_tax - true_burden
)##
##
## |Income_Quintile | Avg_Change| Pct_Change|
## |:---------------|----------:|----------:|
## |Lowest | 1660.26| 91.28|
## |Lower | 126.16| 5.78|
## |Middle | -1821.38| -32.07|
## |Higher | -3644.28| -49.16|
## |Highest | -5693.74| -60.23|