就是正常的用一个direc.dea函数,值如下,运行时却出现了“ In re[i] <- sol : 被替换的项目不是替换值长度的倍数”,想请教一下各位大神,我是错在哪了?
还请指教!谢谢!!!
> base
yg yb x1 x2 x3
10 339503.91 25.526392 13.875834 25598 134831.70
5 89598.61 39.560303 20.753479 4665 23009.78
3 305970.10 29.734445 11.909268 11468 112985.32
8 989531.81 182.252293 98.961262 26147 528322.90
1 759930.11 27.812514 13.003372 83486 418878.67
6 3618128.51 408.045368 223.906764 61835 1155138.94
7 584683.85 45.826740 24.786893 14015 321347.36
2 2217939.69 161.981028 119.153212 68633 2295364.97
9 447141.87 6.023103 3.151996 4037 54394.32
4 110277.38 7.999352 3.835442 7940 50227.01
> frontier
yg yb x1 x2 x3
20 417280.3 24.762803 13.428558 18909 109485.08
15 283143.6 56.691582 28.572816 5805 55126.08
13 375563.4 52.606288 21.388228 12758 165021.17
18 1489330.0 249.374400 135.556138 27155 647519.73
11 1163804.5 45.549334 21.620174 86706 508700.67
16 5098905.0 479.056704 262.082436 79541 1688268.53
17 974376.6 56.454885 30.507996 17987 399284.89
12 2919789.5 261.707596 139.784376 57085 3285062.56
19 557806.5 7.154884 3.686734 4166 74187.68
14 204380.2 12.663109 6.029145 10101 81336.86
> direc.dea(base, frontier, ngood = 1, nbad = 1)
[1] 0.7085282 0.0000000 0.7668014 0.8697832 0.4809645 0.7957591 0.7187270 0.7012121 0.0000000 0.6994839
Warning messages:
1: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
2: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
3: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
4: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
5: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
6: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
7: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
8: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
9: In re[i] <- sol : 被替换的项目不是替换值长度的倍数
10: In re[i] <- sol : 被替换的项目不是替换值长度的倍数