* 1. 读取数据
use your_dataset.dta, clear

* 2. 计算除当年度外的 FDI 年度均值作为工具变量
bysort region year: egen mINCEN = mean(FDI) if year != year[_n]

* 3. 执行两阶段最小二乘法回归
ivregress 2sls digitization_transform (FDI = mINCEN) control_variables, robust
* 说明： - `digitization_transform` 是被解释变量（数字化转型）
* - `FDI` 是内生解释变量
* - `mINCEN` 是工具变量
* - `control_variables` 是其他控制变量（如有需要）

* 4. 检验工具变量的有效性
estat firststage  // 检查第一阶段回归的结果
estat overid      // 检查工具变量的过识别限制检验（如果有多个工具变量）

* 5. 解释 Kleibergen-Paap rk Wald F 统计量。从回归结果中获取 Kleibergen-Paap rk Wald F 统计量。
di "Kleibergen-Paap rk Wald F statistic: " e(rkstat)