在Stata中进行回归之后,我试图仅绘制交互项的系数。
我无法使用社区贡献的命令coefplot
来执行此操作。
这是一个可复制的示例,也是我尝试的解决方法:
sysuse auto, clear
reg price foreign i.turn foreign#i.turn
*this plots all coefficients:
coefplot,
*this drops _cons and foreign but not i.turn
coefplot, drop(i.turn _cons foreign )
*variations with keep also do not work
coefplot, keep(foreign#i.turn )
还有其他方法吗?
我已经在Statalist上交叉发布了这个问题。
答案 0 :(得分:3)
您只需要指定交互:
sysuse auto, clear
reg price foreign i.turn foreign#i.turn, coeflegend noheader
local coefinter 1.foreign#33.turn 1.foreign#34.turn 1.foreign#35.turn ///
1.foreign#36.turn 1.foreign#37.turn
coefplot, keep(`coefinter')
编辑:
您还可以按如下方式获取所有非零系数:
sysuse auto, clear
reg price foreign i.turn i.foreign#i.turn, coeflegend noheader
matrix A = e(b)
local namecol "`: colnames A'"
tokenize `namecol'
forvalues i = 1 / `=colsof(matrix(A))' {
local mv = A[1,`i']
if `mv' != 0 & strmatch("``i''" , "*#*") {
local coefinter `coefinter' ``i''
}
}
coefplot, keep(`coefinter')
答案 1 :(得分:2)
一种比已经提出的更简单的解决方案是使用通配符(“*”或“?”)。
sysuse auto, clear
reg price foreign i.turn foreign#i.turn
coefplot, keep(*.foreign#*.turn)