面向科研政策评估的完整 DID 工具集：双向固定效应估计、个体层聚类稳健标准误、事件研究法平行趋势检验、出版级可视化。
方案说明
基准模型为双向固定效应 DID：
$$y_{it} = \alpha_i + \lambda_t + \beta \cdot \left(\text{Treat}_i \times \text{Post}t\right) + \varepsilon{it}$$
其中 $\alpha_i$ 为个体固定效应，$\lambda_t$ 为时间固定效应，交互项系数 $\beta$ 即政策处理效应。估计上用稀疏虚拟变量承载固定效应并配合反斜杠求解，避免逐元素循环与 inv；标准误采用聚类稳健 sandwich 形式 $V = (X'X)^{-1} M (X'X)^{-1}$，其中 $M = \sum_g X_g' u_g u_g' X_g$。
平行趋势检验采用事件研究法，以政策前一期（$k = -1$）为基期。若政策前的系数 $\beta_k$（$k < -1$）联合不显著异于 0，则平行趋势假设成立。
文件清单
| 文件 | 作用 |
|---|
estimateDid.m | DID 主估计，返回政策效应、聚类稳健标准误、p 值、置信区间 |
estimateEventStudy.m | 事件研究法，输出动态系数与平行趋势 Wald 检验 |
plotDidResults.m | 绘制均值趋势图与事件研究系数图 |
buildDummyMatrix.m | 工具函数：分组变量转稀疏虚拟变量矩阵 |
computeClusterRobustVariance.m | 工具函数：聚类稳健方差-协方差矩阵 |
五个文件放在同一目录下即可直接调用。
1. estimateDid.m
<EMPTY_PARAGRAPH>
function results = estimateDid(panelTable, options)
