用 R做 DEA 工作记录

简单来说，dea算TFP是这个意思。

因为没有找到python算dea的很多资料，不知道怎么下手，用了网上一个人写的代码结果也是错的，还有别的方法太复杂。查到一个 pydea 包，也没看到任何人用过。

最后查了一下，发现R有成熟的包，函数看起来也很简单，就用了一下。

弄了半天，男朋友又不要算这个指标了，他根本不知道这个指标是什么意思，就让我算。好想杀他。

总觉得上课的时候老师讲过一点DEA，总觉得说过，但是什么都想不起来了。

###

library(foreign)                    #可以直接读取stata12以下格式的
df  <- read.dta("auto.dta")
library(readstata13)              #可以读取stata13.stata14格式的
df  <- read.dta13("auto.dta")



x=firmlevel[,c(6,3)] #提取第6列和第3列

library(readstata13) 
library(Benchmarking)
firmlevel=read.dta13("E:/文件/zcy数据/200107 firm level.dta", convert.factors=F)
x=firmlevel[,c(6,3)]
y=firmlevel[,7]
str(firmlevel)
x<- data.frame(a,b)

a<- c(1,2,3)
b<- c(1,2,3)
x<- data.frame(a,b)
y<-c(1,1,1)
e<-dea(x,y)             #全部结果，是一个列表，不需要设置什么orientation参数，设置了就不对了，就默认的就行
f<-eff(e)                #效率，是一串向量，结果是   1,0.5,0.3333333

扫描二维码关注公众号，回复： 8476770 查看本文章

#####

R 先分组，再用几个指标，进行一个多参数的函数的计算：

https://bbs.pinggu.org/thread-6870668-1-1.html

by(data[c("dbh", "tht", "bolt", "dob")], list(data$no), 
   function(df) {
     f(dfdbh,dfdbh, dftht, dfbolt,dfbolt, dfdob)
   })  


#没有试验过，先记下来，还不知道那个df什么意思，以后有需要的时候再说。具体情况看上面那个链接。