求助-求助DEA分析代码（面板数据的数据包络分析）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

（悬赏1000论坛币）求助-求助DEA分析代码（面板数据的数据包络分析），最好是matlab，有demo的那种，感激不尽

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：数据包络分析 DEA分析 数据包络 面板数据 数据包

太史公 发表于 2017-9-14 08:34
（悬赏1000论坛币）求助-求助DEA分析代码（面板数据的数据包络分析），最好是matlab，有demo的那种，感激不 ...

面板分析不了了，只能够一年分析。

您說的是麥氏DEA,可惜,我只會使用R實現。其實還有一套不錯的免費軟體DEAP,它是專門用來實現麥氏DEA的。

cheng334 发表于 2017-9-14 16:14
您說的是麥氏DEA,可惜,我只會使用R實現。其實還有一套不錯的免費軟體DEAP,它是專門用來實現麥氏DEA的。

那请问您有r的关于面板dea的乘数模型代码吗~

简单的用DEAP做，复杂的，带非期望的，可以用MATLAB或者MAXDEA做，有问题私信

library(Benchmarking)

> A=array(c(50,40,20,30,40,35,20,30,20,30,12,12,15,28,10,18,
+           50,48,25,40,35,50,20,50,30,32,13,25,20,30,11,30,
+           45,50,30,40,40,45,28,45,32,35,15,30,25,30,12,38),dim=c(4,4,3))
>
> tse_1=dea(A[,c(1,2),1],A[,c(3,4),1],RTS="crs",ORIENTATION="in")
> A             #頭兩列是投入,後兩列是產出
, , 1

     [,1] [,2] [,3] [,4]
[1,]   50   40   20   15
[2,]   40   35   30   28
[3,]   20   20   12   10
[4,]   30   30   12   18

, , 2

     [,1] [,2] [,3] [,4]
[1,]   50   35   30   20
[2,]   48   50   32   30
[3,]   25   20   13   11
[4,]   40   50   25   30

, , 3

     [,1] [,2] [,3] [,4]
[1,]   45   40   32   25
[2,]   50   45   35   30
[3,]   30   28   15   12
[4,]   40   45   30   38

> tse_2=dea(A[,c(1,2),2],A[,c(3,4),2],RTS="crs",ORIENTATION="in")
> tse_3=dea(A[,c(1,2),3],A[,c(3,4),3],RTS="crs",ORIENTATION="in")
> TSE=cbind(tse_1$eff,tse_2$eff,tse_3$eff);TSE
          [,1]      [,2]      [,3]
[1,] 0.5833333 1.0000000 1.0000000
[2,] 1.0000000 1.0000000 1.0000000
[3,] 0.8000000 0.9288462 0.6959607
[4,] 0.8571429 1.0000000 1.0000000
>
>
> iei_1_2=dea(A[,c(1,2),1],A[,c(3,4),1],XREF=A[,c(1,2),2],YREF=A[,c(3,4),2],RTS="crs",ORIENTATION="in")
> iei_1_3=dea(A[,c(1,2),1],A[,c(3,4),1],XREF=A[,c(1,2),3],YREF=A[,c(3,4),3],RTS="crs",ORIENTATION="in")
> iei_2_3=dea(A[,c(1,2),2],A[,c(3,4),2],XREF=A[,c(1,2),3],YREF=A[,c(3,4),3],RTS="crs",ORIENTATION="in")
> IEI=cbind(iei_1_2$eff,iei_1_3$eff,iei_2_3$eff);IEI
          [,1]      [,2]      [,3]
[1,] 0.6456693 0.6250000 1.0714286
[2,] 1.3340659 1.1362126 0.9054594
[3,] 0.9111111 0.8225806 0.8323529
[4,] 1.0000000 0.7105263 0.8333333

上述是我之前用來練習的麥氏DEA的代碼,投入產出資料是A矩陣,共有三個時期,四個DMU,使用r套件為Benchmarking,未處理到MI的計算,但MI就是TSE及IEI的比值罷了,查查公式就能計算。給您做參考。