stata - Stata で GLM を使用して差差推定器を実装する

Question

Stata 13.0 の GLM モデルを使用して差分推定器を実装しようとしています。私が興味を持っているパラメータは、バイナリ治療群指標 T とバイナリ治療後期間指標 Sのみの相互作用に関する期待値の導関数です(T に関する完全な導関数ではなく、T#S)。このアプローチについては、Statalistのこのスレッドの後半で説明されています。これは私のコードです：

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()マージンでは露出変数でオプションを使用できないため、保存/置換/復元の手順が必要です。

2 つの質問です。

1. この効果の p 値を取得するにはどうすればよいですか?
2. おそらく何らかの方法で eydx() でマージンを使用することにより、半弾性形式で効果を得ることができますか?