做随机前沿回归,回归一直跑,不出结果,总是显示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
有没有大佬知道怎么解决啊