在Stata中用GLM实现差分估计

我正在尝试用Stata 13.0实现一个GLM模型。我感兴趣的参数是预期值相对于二元治疗组指标T和二元治疗后时期指标S相互作用的导数（T#S，而不是相对于T的全导数）。这一方法在《论国家主义》的结尾部分进行了解释。这是我的代码：

glm y i.T##i.S, exposure(e) cluster(user_id) link(log) family(poisson) robust
preserve
    replace e = 30
    margins rb0.T#rb0.S
restore

保留/替换/恢复步骤是必要的，因为页边距不允许将

at（）

选项用于曝光变量

两个问题

如何获得该效果的p值

是否有可能以半弹性形式获得这种效果，或者以某种方式使用eydx（）的边距