你可以使用lincom与Stata中用##编码的交互模型吗?

时间:2016-08-24 21:29:04

标签: stata

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

guard

回归效果很好,但lincom给出了以下错误:

  

曝光##性无效名称

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

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

然后我得到

  

无效的矩阵条纹;   c.l10_mep_i_sg2_pg ##性别

我做错了什么,或者这种#interactive编码是不可能的?

1 个答案:

答案 0 :(得分:2)

尝试在回归命令的末尾添加, 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时,这些都会给你带来外国来源对价格的边际影响。