简单来说,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
#####
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什么意思,以后有需要的时候再说。具体情况看上面那个链接。