如何调用回归估计结果中某一个变量的系数估计值和相应的·伴随概率？

我现在做了一个回归，比如说 reg y x1 x2，我想调用变量x2的估计系数，和它的相应的伴随概率。最好是可以写成scalar的形式，我要通过程序控制，请问通过什么语句可以实现呀？大侠快来帮忙呀！

sysuse auto, clear
rename price y
rename wei   x1
rename mpg   x2
reg y x1 x2
scalar b1 = _b[x1]
mat v = e(V)
scalar t_b1 = b1/sqrt(v[1,1])
local df = e(N) - e(df_m)
scalar p_b1 = ttail(`df' , t_b1)*2

dis b1
dis t_b1
dis p_b1