求助一下关于门槛回归的问题

     关于门限模型这一块，我用的是王群勇老师的xthreg命令，我了解了一天，但还是有几个点不明白：在stata上运行code后，得出运行结果来，在门槛估计量上的输出中得到门限值，但是门限值的个数该怎么确定，我用了王老师原论文上的画图命令得出命令结果图，只能画出王老师论文里的图，但之后再进行三门槛的命令运行后就画不出图来了，只是 不断报错-------------由此我的疑问是是不是可以根据似然比的门槛结果图来确定门限值个数，如果是，画图的命令该怎么改进（报错结果为“matrix should have at least 2 columns”），画完图之后我该怎么确定门限值个数
这是根据命令画出来的图，这种图该怎么解释


问的问题有点细碎，但是查了好久的资料，就是没找到解决这几个问题的方法，希望好心人给解决一下，解决一部分也行