您可以将 lincom 与在 Stata 中使用 ## 编码的交互模型一起使用吗？

Question

我正在尝试使用 lincom 对使用 ## 进行交互编码的 Stata 模型中的回归系数求和。曝光和体重指数是连续变量。性别是二元的。

regr bmi c.exposure##sex covar1 covar2 covar3 i.covar4 i.covar5
lincom chemical + chemical#sex

回归工作正常，但 lincom 给出以下错误：

exposure##sex 无效的名字

或者，如果我将第二行编码为

lincom chemical + c.chemical##sex

然后我得到

无效的矩阵条纹； c.l10_mep_i_sg2_pg##sex

是我做错了什么还是#interaction coding不可能做到这一点？

Answer 1

尝试在回归命令的末尾添加, coeflegend 选项。这将允许您查看 Stata 对每个系数的称呼。

这是一个可重现的示例：

sysuse auto
reg price i.foreign##c.mpg, coeflegend
lincom 1.foreign+ 1.foreign#c.mpg*25

或者，使用margins 通常可以更轻松地完成这类事情：

margins, dydx(foreign) at(mpg=25)

当每加仑英里数为 25 时，这两者都会给您提供外国来源对价格的边际效应。