在coefplot中仅绘制交互作用项

时间:2018-07-23 18:15:16

标签: stata stata-macros coefplot

在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上交叉发布了这个问题。

2 个答案:

答案 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)