怎么给把standard error刮号加上去？

跑完回归后，我们可以用mat se= vecdiag(cholesky(diag(vecdiag(e(V)))))， 把standard error提出来。这个时候是一个数字。

这样我们就有了两个数字，一个是系数，一个是standard error。但是我们的standard error一般需要加挂号，怎么用写程序的方法把standard error刮号加上去？

econfj 发表于 2012-8-17 20:49
如果要导出了，比如到excel，还有我想在导出到excel之前就只显示两位小数后的数字，不用mata可以吗？

*思路简单，但步骤啰嗦一些（可考虑安装outreg2）：
statsby _b _se,clear: reg y x*
xpose,clear v
split _varname,p("_") g(name)
ren (name3 name1) (name b)
sort name name2
replace b=string(v1,"%9.2f")
replace b="("+b+")" if name2=="se"
replace name="" if name2=="se"
export excel name b using reg, first(var) replace

lu guo

**用mata比较方便（设因变量是y，自变量是x1、x2、x3）

statsby _b _se,clear: reg y x*
mata
v=(strofreal(st_data(.,"_b*"))',J(4,1,"(")+ strofreal(st_data(.,"_se*"))'+J(4,1,")"))
v
end

**分解写法（更一般些）
statsby _b _se,clear: reg y x*
mata
b=strofreal(st_data(.,"_b*"))'
se=strofreal(st_data(.,"_se*"))'
n=rows(b)
v=(b,J(n,1,"(")+se+J(n,1,")"))
v
end

为什么不用outreg2呢？输出结果自然包括括号的。

reg...
outreg2 using out, excel replace dec(2)

你要的都有了。

看不懂

sungmoo 发表于 2012-8-17 16:21
**用mata比较方便（设因变量是y，自变量是x1、x2、x3）

statsby _b _se,clear: reg y x*

谢谢！sungmoo

这样的确在屏幕上能显示刮号，不过如果要导出了，比如到excel，还有我想在导出到excel之前就只显示两位小数后的数字，不用mata可以吗？不熟悉mata，^_^

请问vecdiag的命令在没有装，在哪里能获取？