我尝试做一个两期DID模型回归，构建了2011-2020的地级市面板。第一期是2014年的第一批试点城市，第二期则推广到2016年全国，后来发现由于TWFE估计量由于负权重问题存在严重的偏误，导致基准回归系数由正变负，因此换用了did2s命令（即二阶段DID估计量）。
但研究后续我想针对政策展开异质性分析，以分析政策在不同地区的异质性效果，我想请问在did2s命令下，应当如何展开异质性分析。
我原先设想是做交乘项，如是，应当这样编写代码吗？
did2s y, first_stage(i.id i.year z z*treat) second_stage(i.treat) treatment(treat) cluster(id) 。z是用于区分异质性的解释变量，例如地区的人均GDP一类。

或者直接进行分组回归，但是分组回归系数不能跨组直接比较的问题会有所限制，我还是想如果能用交乘项最好