Question

我正在尝试获取分类变量的所有效果差异组合到表中。

到目前为止，我已经接近我想要的：

sysuse auto
reg gear_ratio i.rep78

foreach i in 1 2 3 4 5 { 
reg gear_ratio b`i'.rep78
contrast r.rep78
mat A`i'=r(table)
}

mat B = (A1\A2\A3\A4\A5)
esttab matrix(B)

是否可以将列作为行？我只需要系数。

我也尝试过：

mat C = (A1,A2,A3,A4,A5)

但是，结果不是我想要的。

Answer 1

您只需要转置每个矩阵AX（我还对您的代码进行了一些改进）：

sysuse auto, clear
regress gear_ratio i.rep78

forvalues i = 1 / 5 {
    regress gear_ratio b`i'.rep78
    contrast r.rep78
    matrix A`i' = r(table)' // transpose
    matrix A`i' = A`i'[1...,1] // select column
    local all `all' A`i'
}

local all = subinstr("`all'", " ", "\", .)
matrix B = `all'

esttab matrix(B), collabels(none)

但是，您也可以通过简单地使用pwcompare的{{1}}选项来完成此操作：

margins

创建效果差异组合表

1 个答案: