储存门槛值ereturn list mat list e(Thrss)
mat s=e(Thrss)
global q1 = s[2,1] （连老师示例为双重门槛，我想做单一门槛，这的命令该怎么没？这个s[2.1]是什么意思？）
lobal q2 = s[3,1]  

  估计结果
local q1 = $q1
local q2 = $q2
dis "q1=" in y `q1' in g "  q2=" in y `q2' （这个命令什么意思，如果做成单一门槛该怎么改？）
dropvars d1 d2 xd1 xd2
gen d1 = (q>`q1'&q<`q2') /
gen d2 = (q>=`q2') 817   gen xd1 = x*d1
local q1s = substr("`q1'",1,4)
ocal q2s = substr("`q2'",1,4)
label var xd1 "xd1(q>`q1s'&q<`q2s')"
gen xd2 = x*d2
label var xd2 "xd2(q>=`q2s')"
xtreg y x xd1 xd2 z1 z2, fe