执行xtthres y , thres(q) dthres(x) bs1(500) bs2(500) bs3(500)之后就有报告对应的固定效应回归结果，为什么要有如下程序？

呈现估计结果

    local q1 = e(rhat21)       // 取出第二轮搜索的第一个门槛值。

         local q2 = e(rhat22)      // 取出第二轮搜索的第二个门槛值

         gen d1 = (q<=`q1')         // 生成虚拟变量

         gen d2 = (q> `q2')

         gen xd1 = x*d1             // 交乘项

         gen xd2 = x*d2

         xtreg y x xd1 xd2, fe      // 常规标准误

         est store fe

         xtreg y x xd1 xd2, 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)