使用matlab求解Mlamquist-Luenberger指数出现问题求助

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近本人在学习DEA，找到了一个deatoolbox，但是在使用deamalmluen这个命令进行求解时，出现了如下情况：

1. Warning: DMU 11. First Step. Optimization exit flag: -2
   
2. > In deaund (line 197)
   
3.   In deamalmluen (line 111)
   
4.   In forum (line 141)
   
5. Warning: DMU 11. First Step. Optimization doesn't return a result in
   
6. First Step. Efficiency set to NaN.
   
7. > In deaund (line 202)
   
8.   In deamalmluen (line 111)
   
9.   In forum (line 141)
   
10. Warning: DMU 15. First Step. Optimization exit flag: -2
    
11. > In deaund (line 197)
    
12.   In deamalmluen (line 111)
    
13.   In forum (line 141)
    
14. Warning: DMU 15. First Step. Optimization doesn't return a result in
    
15. First Step. Efficiency set to NaN.
    
16. > In deaund (line 202)
    
17.   In deamalmluen (line 111)
    
18.   In forum (line 141)

复制代码

请问这种情况应该如何解决？
数据和代码如下：

1. x(:,:,1)=[20933        276        18611
   
2. 38338        34801        38298
   
3. 155978        173429        143098
   
4. 17300        641        8204
   
5. 108236        1941        14829
   
6. 531024        575562        486604
   
7. 191323        10146        42640
   
8. 374585        17268        156399
   
9. 45448        153832        43725
   
10. 779585        1060621        1289362
    
11. 332572        2315180        826248
    
12. 574660        214837        566399
    
13. 239210        689488        862193
    
14. 358536        61993        698899
    
15. 620604        827028        1001372
    
16. 200289        6611        373456
    
17. 665667        117225        948715
    
18. 431867        54937        1131910
    
19. 595328        876453        1371681
    
20. 254382        279452        637803
    
21. 62179        244777        240323
    
22. 68560        11176        505511
    
23. 183936        8382        1229134];
    
24. y(:,:,1)=[37763.39063
    
25. 67743.78906
    
26. 174187.5938
    
27. 8030.772949
    
28. 14907.68164
    
29. 951548.875
    
30. 27373.375
    
31. 110279.0078
    
32. 118164.5781
    
33. 1323135.125
    
34. 1755794.125
    
35. 393147.0313
    
36. 1784757
    
37. 369612.0313
    
38. 2286594.5
    
39. 116007.2969
    
40. 518146.125
    
41. 257510.375
    
42. 1802576.5
    
43. 334858.2188
    
44. 286889.2188
    
45. 22058.81055
    
46. 8982.219727];
    
47. yu(:,:,1)=[23744.24414
    
48. 61913.15234
    
49. 136914.2969
    
50. 5186.003906
    
51. 13981.53613
    
52. 550892.4375
    
53. 15922.27051
    
54. 85320.26563
    
55. 112180.0781
    
56. 860442.1875
    
57. 807219.9375
    
58. 322136.3125
    
59. 837553.6875
    
60. 305820.1563
    
61. 857084.5625
    
62. 57816.04688
    
63. 423289.0313
    
64. 248024.3906
    
65. 1046380.688
    
66. 277902.9063
    
67. 47799.63672
    
68. 10660.13184
    
69. 7462.730957];
    
70. x(:,:,2)=[20049        232        17872
    
71. 40389        33846        38338
    
72. 158027        166090        163585
    
73. 18130        726        7629
    
74. 110797        1908        17600
    
75. 560317        576557        493222
    
76. 189612        7642        41589
    
77. 379763        16202        150279
    
78. 42543        152175        43211
    
79. 800246        1102210        1261832
    
80. 323363        2213101        794522
    
81. 589133        224480        585088
    
82. 247793        710336        871298
    
83. 363831        64289        712713
    
84. 667115        788455        1028134
    
85. 203726        7341        383755
    
86. 671462        112748        988551
    
87. 453585        46773        900686
    
88. 599003        854607        1401403
    
89. 259411        283627        636203
    
90. 63517        261647        237353
    
91. 29523        2168        39198
    
92. 64711        1605        13935];
    
93. y(:,:,2)=[35187.14453
    
94. 77367.70313
    
95. 208592.3594
    
96. 8699.649414
    
97. 14714.12891
    
98. 1097135.125
    
99. 27209.96094
    
100. 119441.4219
     
101. 134439.4531
     
102. 1599358.5
     
103. 1914448.25
     
104. 440216.9688
     
105. 1899217.25
     
106. 679891.8125
     
107. 2493645
     
108. 125093.8984
     
109. 606324.5
     
110. 279395.7813
     
111. 1799590.875
     
112. 361206.2188
     
113. 62733.65625
     
114. 23071.30273
     
115. 12637.12109];
     
116. yu(:,:,2) = [33742.45703
     
117. 85424.78125
     
118. 244769.6094
     
119. 9342.887695
     
120. 18307.11133
     
121. 1221806
     
122. 29398.92188
     
123. 117318.5313
     
124. 136693.9688
     
125. 1671097.875
     
126. 2089436.875
     
127. 161083.8438
     
128. 682680.1875
     
129. 316933.1875
     
130. 2194546
     
131. 388743.9688
     
132. 357223.9375
     
133. 79161.20313
     
134. 298215.0938
     
135. 24558.83594
     
136. 83439.00781
     
137. 23664.6543
     
138. 15048.27539];
     
139. ml = deamalmluen(x, y, yu);
     
140. deadisp(ml);

复制代码

deamalmluen的语法如下：

1. deamalmluen Data envelopment analysis Malmquist-Luenberger indices
   
2.     Computes data envelopment analysis Malmquist-Luenberger indices
   
3. 
   
4.     out = deamalmluen(X, Y, Yu, Name, Value) computes data envelopment
   
5.     analysis Malmquist-Luenberger indices with inputs X and outputs Y.
   
6.     Model properties are specified using one or more Name ,Value pair
   
7.     arguments.
   
8. 
   
9.     Additional properties:
   
10.     - 'names': DMU names.
    
11.     - 'orient': orientation. Directional distane function with undesirable
    
12.     outputs 'ddf' (Aparicio, Pastor and Zofio, 2013), default. Directional
    
13.     distance function with undesirable outputs 'ddf_ccf' (Chung, Fare and
    
14.     Grosskopf).
    
15.     - 'fixbaset': previous year 0 (default), first year 1.
    
16.     - 'geomean': compute geometric mean for technological change. Default
    
17.     is 1.
    
18. 
    
19.     Example
    
20. 
    
21.        malmluen = deamalmluen(X, Y, Yu);
    
22.        malmluenfixed = deamalmluen(X, Y, Yu, 'fixbaset', 1);

复制代码

附件中为工具箱
求赐教！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

这个模型可以做出来么？
我用MAXDEA专业版做
不知道结果是不是一样的

PX0706 发表于 2018-4-4 15:11
这个模型可以做出来么？
我用MAXDEA专业版做
不知道结果是不是一样的

是可以的，但是只有一种解法，并且容易出现问题。

PX0706 发表于 2018-4-4 15:11
这个模型可以做出来么？
我用MAXDEA专业版做
不知道结果是不是一样的

这个工具箱中的算法是Chung, Fare and Grosskopf (1997) 和Aparicio, Pastor and Zofio, 2013中的算法，即共两种算法。没用过maxdea，不知道这种软件有几种算法。

PX0706 发表于 2018-4-4 15:11
这个模型可以做出来么？
我用MAXDEA专业版做
不知道结果是不是一样的

加一下QQ  1107085694

刘勇超063 发表于 2018-10-31 14:55
加一下QQ  1107085694

好的，可以的，如果你要测算ML指数请加我，我已经加你啦

楼主，请问你的问题解决了吗？我也遇到了相同的问题，然后结果中有部分为NaN，你知道该如何解决吗？

拂去尘缘 发表于 2018-4-3 22:24
最近本人在学习DEA，找到了一个deatoolbox，但是在使用deamalmluen这个命令进行求解时，出现了如下情况：
...

楼主，能加QQ979534884咨询问题吗

不知道LZ后来怎么解决的？
我也是出现这个问题。

拂去尘缘 发表于 2018-4-3 22:24
最近本人在学习DEA，找到了一个deatoolbox，但是在使用deamalmluen这个命令进行求解时，出现了如下情况：
...

请问你们解决了这个问题吗