需要利用stata返回的 e(b) 和 e(V) 两个矩阵自己计算,存储于矩阵中后再行引用。
*---------------命 令-------------
sysuse auto, clear
reg price wei mpg turn
mat b = e(b)
mat Var_b = e(V)
mat se_b = cholesky(diag(vecdiag(Var_b)))
mat b0 = diag(b)
mat inv_se_b = inv(se_b)
mat t = vecdiag(hadamard(b0, inv_se_b))
mat list t
*---------------结 果-------------
. sysuse auto, clear
(1978 Automobile Data)
. reg price wei mpg turn
Source | SS df MS Number of obs = 74
-------------+------------------------------ F( 3, 70) = 14.59
Model | 244284435 3 81428145.1 Prob > F = 0.0000
Residual | 390780961 70 5582585.15 R-squared = 0.3847
-------------+------------------------------ Adj R-squared = 0.3583
Total | 635065396 73 8699525.97 Root MSE = 2362.7
------------------------------------------------------------------------------
price | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
weight | 3.524339 .8171395 4.31 0.000 1.894605 5.154072
mpg | -72.86501 81.29538 -0.90 0.373 -235.0035 89.27351
turn | -395.1902 122.6444 -3.22 0.002 -639.7968 -150.5836
_cons | 12744.24 4760.102 2.68 0.009 3250.519 22237.97
------------------------------------------------------------------------------
.
. mat b = e(b)
. mat Var_b = e(V)
.
. mat se_b = cholesky(diag(vecdiag(Var_b)))
.
. mat b0 = diag(b)
. mat inv_se_b = inv(se_b)
. mat t = vecdiag(hadamard(b0, inv_se_b))
.
. mat list t
t[1,4]
weight mpg turn _cons
r1 4.3130197 -.89629955 -3.2222433 2.6773047
.
[此贴子已经被作者于2007-11-15 17:20:09编辑过]