各位大神 我用了循环语句回归 得出系数 想提取p值，b1 和t_bo是可以提取出来的，经过循环之后会存储在stata的表格里，但是每次运行到replace p_bo 就会提示命令无效，同样的命令为什么到了p_bo就不行？
我是用t= b/se 就算出t，再用t调用函数求p的


reg rate3 bondr10y tspread cspread sml hml rmrf if code == `n'
  *estat vif
  replace b1=_b[bondr10y] in `n'
  replace t_bo = _b[bondr10y]/_se[bondr10y] in `n' //求出t值
  replace p_bo = 2*ttail(e(df_r), abs(`t_bo')) in `n'  //计算得出pvalue

******************************详细回归程序如下********************************************

*开头
cap log close
set more off
cd "E:\job\model"
use panalmix2.dta, clear
xtset fundid week
egen code = group(fundid)

//作图发现有好多异常值
graph hbox rate in 1/3080, over(fundid)
graph save "E:\job\model\ratebox40.gph", asis replace //保存图
graph use "E:\job\model\ratebox40.gph" //打开图


//删除离群点
gen rate3 = rate
*br
_pctile rate3, percentile (10 90)
drop if rate3 < r(r1)
drop if rate3 > r(r2)

*进行回归

gen b1=. //用b来记录5Y国债收益率系数
gen b2=. //记录期限利差系数
gen b3=. //记录信用利差系数
gen t_bo=.
gen t_ts=.
gen t_cs=.
gen p_bo=.
gen p_ts=.
gen p_cs=.
replace b1=_b[bondr10y]
replace t_bo = _b[bondr10y]/_se[bondr10y]
local 2*ttail(e(df_r), abs(`t_bo'))
di

forvalues n = 1/114{
  dis _s(8) `n'
  reg rate3 bondr10y tspread cspread sml hml rmrf if code == `n'
  *estat vif
  replace b1=_b[bondr10y] in `n'
  replace t_bo = _b[bondr10y]/_se[bondr10y] in `n' //求出t值
  replace p_bo = 2*ttail(e(df_r), abs(`t_bo')) in `n'  //计算得出pvalue

  replace b2=_b[tspread] in `n'
  replace t_ts= _b[tspread]/_se[tspread] in `n' //求出t值
  replace p_ts= 2*ttail(e(df_r),abs('t_ts')) in `n'  //计算得出pvalue

  replace b3=_b[cspread] in `n'
   local t_ts = _b[tspread]/_se[tspread]
   replace t_cs= _b[cspread]/_se[cspread] in `n' //求出t值
  replace p_cs=2*ttail(e(df_r),abs(`t_cs')) in `n'  //计算得出pvalue

  }
sum b1 b2 b3, d
*tabstat b1 p_bo b2 p_ts b3 p_cs,stats(mean sd p25 med p75 min max) c(s) f(%6.2f)
tabstat b1 p_bo b2 p_ts b3 p_cs,c(s) f(%6.2f)