E14
动力学财富交换模型中税收再分配规则的比较：以对 WID 实测基尼系数与顶端份额的联合拟合为判据
1 · 研究问题
在动力学财富交换模型（kinetic wealth exchange model）中加入不同的税收再分配规则（比例税 vs 按实际税制标定的累进税），哪一种能在同一组参数下同时拟合世界不平等数据库（WID）中 2–3 个国家的基尼系数与前 10%/前 1% 财富份额？拟合所需的再分配强度参数与该国实测宏观税负的差距有多大？
2 · 研究背景与空白
动力学交换模型把经济体视为大量随机配对交换财富的主体，是统计物理与经济学交叉的经典 ABM：纯随机交换收敛到指数分布（Boltzmann–Gibbs，Drăgulescu & Yakovenko 2000，有解析解）；引入储蓄倾向得到类 Gamma 分布（Chakraborti & Chakrabarti 2000）；Bisi、Toscani 等证明了含再分配算子的模型在连续交易极限下收敛到 Gamma 密度且形状因子由再分配权重决定（又一个解析锚点）。
已有工作：税收与再分配对 Pareto 指数的影响已被多篇论文研究（如 Entropy 2023 关于极端不平等模型稳定性的数值研究；Toscani 等的再分配理论）。但这些工作绝大多数停留在"改规则→看分布形状变化"的定性层面。
空白在评价维度：很少有工作把模型输出与具体国家的实测不平等指标做联合定量校准，并检验"拟合出的再分配参数是否与该国真实税负同量级"——这是模型与现实的双向对表。解析极限完备、数据免费、算力极小，且"模型无法联合拟合"同样是有效结论（说明该类模型缺失关键机制），适合一年期课题。
3 · 可检验假设
- H1：含累进再分配的模型可在单一参数组下同时把基尼系数拟合到实测值 ±0.03、前 10% 份额 ±5 个百分点，而比例税模型无法同时满足两者（至少对高不平等国家如美国失败）。
- H2：拟合出的有效再分配率与该国实测"税收+转移支付占 GDP 比重"（OECD/WID 口径）符号一致且同一数量级（比值在 0.3–3 之间）；若数量级不符，该失配量化本身是主结论。
4 · 量化验收标准
- 方法学校验（硬门槛）：（i）零储蓄、无税模型的稳态分布与解析指数分布做 KS 检验不拒绝（N=10^4 主体、10^4 次/主体交换后，p>0.1，多种子稳定）；（ii）复现 Drăgulescu–Yakovenko 2000 图 1 的指数温标 T=⟨m⟩，偏差 ≤2%；（iii）复现含再分配模型的解析 Gamma 形状因子，偏差 ≤5%。任一不过关，后续全部结论无效。
- 模拟统计口径预先写死：每个参数点 ≥20 个随机种子，报基尼与份额的均值 ± 标准差与 95% 置信区间；稳态判定用滑动窗口基尼漂移 <0.001/10^3 步的显式收敛判据。
- 参数敏感性用全局方法：对（储蓄倾向、税率、再分配靶向参数）做 Sobol 一阶与总指数（Saltelli 采样 ≥1,024 组）；不得只做 OAT。
- 校准目标预先写死：以 WID 2022 年（或最新可得年）美国、法国、（中国，可得性需核实）三国的基尼 + 前 10% + 前 1% 份额为目标向量，拟合用同一损失函数（标准化 L2），禁止逐国换损失函数。
- 有限规模效应检查：N=10^3/10^4/10^5 三档规模扫描，主指标随 N 的漂移 <1% 方可下结论。
- 全部代码开源（NumPy/Numba），单参数点运行时长与总机时写入附录。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 实测不平等指标 | World Inequality Database（wid.world，免费 CSV/API：财富与收入基尼、顶端份额，分国分年） |
| 宏观税负对照 | OECD Revenue Statistics / World Bank WDI（免费；中国税收转移口径可比性需核实，不可比则换英国/德国） |
| 累进税制标定 | 各国法定税率表（OECD Taxing Wages 免费出版物） |
| 模拟框架 | Python + NumPy + Numba（自写核心循环，约 200 行；不依赖重型 ABM 框架） |
| 解析对照 | Drăgulescu & Yakovenko (2000)、Toscani 等再分配论文——仅用于校验与对比，不计入本项目贡献 |
算力：N=10^4 主体 × 10^8 次交换的单次运行在 Numba 下约 1–3 分钟（笔记本实测为准）；Sobol 1,024 组 × 20 种子 ≈ 2×10^4 次运行，总机时约 3–7 天可分批跑，纯 CPU 可行。N=10^6 全国尺度超出范围。降规模版本：Sobol 降为 256 组 + 10 种子（<1 天）。
6 · 方法路径
- 实现基础交换模型，通过三项解析校验（第 1 条验收）。
- 加入储蓄倾向与两种税收再分配算子（比例 / 按税率表分段累进），单元测试财富守恒（每步总量漂移 <10^-10）。
- 下载 WID 目标数据，写统一损失函数与网格 + 局部细化的拟合流程。
- 核实中国口径可比性与 WID 财富（非收入）序列覆盖（第 6 周核实，不边做边发现）。
- 对三国分别拟合两种规则，输出可行参数域与最优拟合及置信区间。
- 做 Sobol 全局敏感性与有限规模扫描。
- 独立交叉校验：对最优参数点用纯 Python（非 Numba）慢速实现复算基尼，偏差 <0.5%。
7 · 新颖性边界
本课题不声称提出新模型（交换核与再分配算子均已发表），不声称模型机制是现实财富分布的因果解释（它是最简代理模型，须在结论中明示）。已有工作：指数/Gamma 稳态、税收改变 Pareto 指数、高税率下双峰分布等均已发表。本项目贡献是评价维度的转换：前人评估"规则改变分布形状"，本项目评估"规则 + 单一参数组能否联合命中实测多指标，及拟合参数与真实税负的数量级对表"——把模型从定性演示推到可证伪的定量校准。丘奖历届 ABM/统计物理类获奖极少（2022 铜奖为三方演化博弈、2024 铜奖为元宇宙交互数学框架，均非财富分布），差异清晰。"两种规则都拟合不上"是有效结论，须报告最小可达损失及其置信区间以证明不是优化失败。
8 · 决策门槛（go / no-go）
- 第 4 周末：三项解析校验全部通过。KS 检验反复失败多为稳态未达或抽样相关性问题→限期 2 周排查；第 6 周仍未过则降级 A：缩减到纯复现 + 收敛性与有限规模研究（可复现性贡献定位），保留"模拟-解析对照"框架。
- 第 10 周末：单次运行时长实测 ≤5 分钟。若 >5 分钟→N 降到 3×10^3 并加大种子数，Sobol 降为 256 组（降规模版本，主结论保留）。
- 第 22 周末：至少一国完成双规则拟合。若两规则均严重失配（损失下界远超阈值）→失配即主结论，补一节"缺失机制讨论"（收入流、资产回报异质性），框架不变。
- 已知困难：WID 的财富份额本身是估算值——把"目标数据不确定性传播到拟合参数"写成研究内容（对 WID 值 ±1 标准误扰动重拟合）。
- 预算裁剪顺序：先砍第三国，再砍累进税的分段精细度（降为两段），最后砍 Sobol 规模；主结论（双规则联合拟合对比）在裁剪下成立。