做随机前沿回归，回归一直跑，不出结果，总是显示not concave或者跑完出现convergence not achieved,或者是discontinuous region encountered           cannot compute an improvement.
运算代码如下
xtset id t
sort id t
egen iid= tag(id)  
global xxx x1 x2
regress y $xxx
matrix b1 = e(b)  
version 10.1  
sfmodel y, prod dist(h) frontier($xxx) usigmas(t) vsigmas()  
sf_init, frontier(b1) usigmas(0 0) vsigmas(0)  
ml max, difficult gtol(1e-5) nrtol(1e-5)  

数据形式如下
id        y        x1        x2        t
1        0        -0.886666667        136.47        1
1        0        -0.686666667        -22.93        2
1        0        -0.486666667        -121.83        3
1        63.40629449        0.313333333        201.17        4
1        72.90032096        -0.086666667        61.57        5
1        62.99606867        -0.686666667        151.67        6
1        53.85679615        0.113333333        -120.93        7
1        83.46912918        0.413333333        175.87        8
1        0        0.313333333        -179.13        9
1        0        -0.186666667        -104.43        10
1        0        0.013333333        -94.93        11
1        0        0.313333333        -102.33        12
1        0        -0.086666667        -60.03        13
1        0        0.113333333        12.67        14
1        0        -0.186666667        -78.53        15
1        0        0.113333333        -79.23        16
1        0        0.613333333        -33.13        17
1        0        -0.086666667        99.17        18
1        0        0.013333333        -108.23        19
1        0        -0.886666667        -34.13        20

有没有大佬知道怎么解决啊