我正在尝试使用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编码是不可能的?
答案 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时,这些都会给你带来外国来源对价格的边际影响。