求问，如何用R语言的Benchmarking做DEA

小弟想请教一下如何用R语言中Benchimarking做DEA数据包络分析，我用这个数据去分析，给出的e是最后的评分吗，为什么跟我用MATLAB分析得到的结果不一样？想知道做DEA的完整过程以及得到的结果的解释，谢谢！
> a
   年份         y1   y2    x1   x2    x3
1  1990    3621.00 0.00 14.40 0.65 31.30
2  1991    3943.00 0.09 16.90 0.72 32.20
3  1992 4086367.00 0.07 15.53 0.72 31.87
4  1993 4904367.00 0.13 15.40 0.76 32.23
5  1994    6311.67 0.37 14.17 0.76 32.40
6  1995    8173.33 0.59 13.33 0.69 30.77
7  1996   10236.00 0.51 12.83 0.61 29.23
8  1997   12094.33 0.44 13.00 0.63 28.20
9  1998   13603.33 0.58 13.40 0.75 28.80
10 1999   14841.00 1.00 14.00 0.84 29.10
> x <- with(a, cbind(x1,x2,x3))
> y <- with(a, cbind(y1,y2))
> e<-dea(x,y)
> e
[1] 0.9385 0.8758 1.0000 1.0000 0.9054 0.9768 1.0000 1.0000 0.9881 1.0000

dea_result <- Benchmarking::dea(x, y, RTS="crs", XREF=x, YREF=y, ORIENTATION="out", SLACK=TRUE);
efficiency  <- dea_result$eff

efficiency 即为效率值。如果和matlab不一样，请检查RTS的设定（可以为“crs”或“vrs”）及ORIENTATION的设定（“out”或“in”）

halleyg 发表于 2019-8-20 08:24
dea_result

谢谢您，我回头试试

halleyg 发表于 2019-8-20 08:24
dea_result

您好，我尝试了，也变换了你说的两个参数，依然跟matlab计算结果不一致，请问是什么情况

halleyg 发表于 2019-8-20 08:24
dea_result

我感觉matlab的结果更加真实

你可以用DEAP2.1先算一遍