stata - 从 Stata 中的自定义程序发布 e(b) 向量

Question

我编写了一个计算加权回归的程序，现在我希望将我的估计结果存储为 e(b) 向量，以便引导命令可以轻松访问结果，但我一直收到错误消息。我的程序看起来像:

capture program drop mytest
program mytest, eclass
version 13
syntax varlist [if]
marksample touse
// mata subroutine creates matrix `b', such as mata: bla("`varlist'", "`touse'")
tempname b
matrix `b' = (1\2\3)
ereturn post `b'
end 

mytest town_id
ereturn list

但我不断收到 conformability error r(503);在运行脚本时。当我改为发布一个普通矩阵时，例如 ereturn matrix x = b ，一切正常，但我想将我的系数“正确”存储在 e(b) 中向量。

我检查了 Stata 的文档，但无法找出为什么这不起作用。他们的建议是编码

tempname b V
// produce coefficient vector `b' and variance–covariance matrix `V'
ereturn post `b' `V', obs(`nobs') depname(`depn') esample(`touse')

ereturn post的选项都是可选的。谁能告诉我我在这里缺少什么？谢谢!

Answer 1

使用“行”向量而不是“列”向量。例如，如果您检查 regress 的存储结果，您会看到这是预期的结果。

capture program drop mytest
program mytest, eclass
version 13
syntax varlist [if]
marksample touse
// mata subroutine creates matrix `b', such as mata: bla("`varlist'", "`touse'")
tempname b
matrix `b' = (1,2,3)
ereturn post `b'
end 

*----- tests -----

clear
sysuse auto

// mytest test
mytest mpg weight
ereturn list
matrix list e(b)

// regress example
regress price weight mpg
ereturn list
matrix list e(b)