stata在计算过程中好像没有直接存储置信区间的数值我们知道回归系数的置信区间是根据点估计值与标准差计算出来的,以stata自带数据集为例,运行下列语句
sysuse auto
regress mpg weight c.weight#c.weight foreign
mat B=e(b) /*将回归得到的系数结果保存在矩阵B中*/
mat SE=e(V) /*将回归得到的协方差矩阵保存在矩阵SE中*/
此时如果我们输入
di B[1,1]
di sqrt(SE[1,1])
则会发现,所得到的值正好是第一个自变量weight的回归系数和标准差,从而我们可以计算其95%可信区间
local lci=B[1,1]-invttail(e(df_r),0.025)*sqrt(SE[1,1]) /*计算95%可信下限*/
local uci=B[1,1]+invttail(e(df_r),0.025)*sqrt(SE[1,1]) /*计算95%可信上限*/
di `lci'
di `uci'
此时我们通过调用局部宏就可以存储置信区间的信息了
回归过程中存储的更多中间变量请直接help regress查询
|