【C0126】《管理世界》复现！CL-SW-GEKS法测算购买力平价（PPP）

4. 复现流程

第一步：准备原始数据

将以下原始文件放入 data/ 目录：

• • p2021fi.xlsx：价格数据。
• • q2021fi.xlsx：数量数据。
• • tdata.xlsx：供 Stata 执行 QUAIDS 估计使用的整理后数据。

当前项目中已包含上述三个 Excel 文件。

第二步：运行 Stata 脚本生成 QUAIDS 结果

需要先修改do文件中的数据路径。

运行文件：

• • QUAIDS估计+双曲正弦变换.do

该脚本完成的工作包括：

• • 从 data/tdata.xlsx 导入数据。
• • 计算总支出 expgdp。
• • 计算各类支出份额。
• • 对份额做双曲正弦变换。
• • 估计 QUAIDS 模型。
• • 生成预测结果。
• • 导出 data/quaidsresu.csv。

注意：quaidsresu.csv 是后续 R 脚本运行的必要输入，由于精度问题（原作者选用float精度，这里提供double精度，能精确到更多的小数位数，以供特殊需求），会与原作者的quaidsresu.csv文件有细微差距。

第三步：运行 R 脚本计算 CL-SW-GEKS 及对比方法

运行文件：

• • 表1 各汇总方法PPP测算结果.R

该脚本会：

• • 读取 data/p2021fi.xlsx、data/q2021fi.xlsx、data/quaidsresu.csv。
• • 调用 CL-SW-GEKS方法.R 计算 CL-SW-GEKSu、CL-SW-GEKSs 和 GEKS。
• • 同步计算 CPD、IDB、GK、GAIA、SS、MPCP、Blended、Own Share、Van Yzeren、Gerardi 等方法。
• • 整理为对比表并输出 Excel 文件。
• •

第四步：查看输出结果

输出目录为：

• • output/

当前脚本会输出两个 Excel 结果文件：

• • 表1 各种测算PPP方法的对比（原始数据）.xlsx
• • 表1 各种测算PPP方法的对比.xlsx

其中一份是原始结果对比表，另一份包含标准化后的展示结果。

左边为原文展示的结果，右边为复现代码产出的结果。

注意：
由于第二步引入的精度偏差会被第三步的计算脚本放大，导致某些国家的结果会与原文的结果有一定差距，属于正常现象。若希望完全复现可以自行修改。