如何用Stata的reghdfe或R实现配对堆叠数据的双重差分回归？

我尝试在堆叠数据集上运行**双重差分（difference in difference）**回归。此处的堆叠是指每个处理组公司都会与一个控制组公司配对，这些配对组被垂直堆叠。数据的大致结构如下：

共有50组配对，除结果变量外还有其他控制变量。

我想要运行的回归模型如下：

其中T表示配对组d中的处理单元s，P表示配对组d中的处理后时期t。该模型可将时间和个体固定效应与配对标识符结合。

我不确定如何在Stata中使用reghdfe或在R中正确设置该模型，也不清楚如何指定此处的多水平固定效应与交互项。

我曾尝试在Stata的reghdfe中加入配对ID与个体标识符的交互项、配对ID与日期的交互项作为固定效应，但输出结果存在问题。

非常感谢您的帮助。

内容的提问来源于stack exchange，提问作者ECF