截面门槛:thresholdtest Y X X1 X2 X3,q(门槛处理变量Xn) trim_per(0.15) rep (500) //测算有无门槛
其中核心自变量X可以和门槛处理变量Xn相同,0.15为默认值 rep(500)是指在样本中抽样500次,可以视样本量修改,我1800样本抽样200次要将近30分钟。 输出结果关注bootstrap P_value 的值 显著时可以认为存在门槛
第一个门限:thresholdreg Y X1 X2 X3 ,q(门限变量) h(1) 。。。。。。。。。。假设门限是1998
第二个门限:drop if 门限变量Xn<=1998
重复第一个thresholdtest命令