各位大神，帮帮忙，最近在做面板门槛模型，得出单一门槛，我自己按照单一门槛值分样本进行回归，但是结果好像与连老师给出的程序中最后的交乘项的系数估计不一样。为什么呢？
       连老师给出的是关于双重门槛的代码，我自己改成单一门槛的，请大家帮我看看有什么问题，如下：
xttr_graph
dis e(rhat)
global q1 = e(rhat)
dis "$q1"
dropvars d1 lev_x_*
gen d1 = (cash<=$q1)
gen d2 = (cash> $q1)
gen lev_x_cash1 = lev*d1
gen lev_x_cash2 = lev*d2
local x "size prof top1 tang grow beta year2 year3 year4 year5 year6 year7 year8 year9 year10"
xtreg q `x' lev_x_cash* lev, fe
est store fe
xtreg q `x' lev_x_cash* lev, fe robust
est store fe_robust
local m "fe fe_robust"
esttab `m', mtitle(`m') nogap s(r2 r2_w N F) star(* 0.1 ** 0.05 *** 0.01) order(`x' lev_x_cash1 lev_x_cash2 )

       最后得出的结果还是除了lev_x_cash1 和 lev_x_cash2的估计系数，还是会出现lev（虽然没有显示数值），但是这样会有影响吗？为什么这样的结果和我自己按照门槛值分样本进行回归的结果不一样呢？
       求大神帮我看看