R软件做Bootstrap DEA修正

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

R软件怎么做Bootstrap DEA修正？求指教！万分感谢！

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：Bootstrap Bootstra boots boot Trap 软件

我给你报价了。。。

我就不报价了......下一个包，FEAR，在这：http://www.economics.clemson.edu/faculty/wilson/Software/FEAR
给你一个例子，参考一下——
library(FEAR)

data(ccr)

x=matrix(nrow=5,ncol=70)

x[1,]=ccr$x1

x[2,]=ccr$x2

x[3,]=ccr$x3

x[4,]=ccr$x4

x[5,]=ccr$x5

y=matrix(nrow=3,ncol=70)

y[1,]=ccr$y1

y[2,]=ccr$y2

y[3,]=ccr$y3

Alternatively,onemighttype

data(ccr)

x=t(matrix(c(ccr$x1,ccr$x2,ccr$x3,ccr$x4,ccr$x5),

nrow=70,ncol=5))

y=t(matrix(c(ccr$y1,ccr$y2,ccr$y3),nrow=70,ncol=3))

dhat=dea(XOBS=x,YOBS=y)

tmp=boot.sw98(XOBS=x,YOBS=y,DHAT=dhat,NREP=2000)

n=ncol(x)

#number of DMUs

table.in=matrix(nrow=n,ncol=7)

table.in[,1]=c(1:n)

table.in[,2]=dhat

table.in[,3]=dhat-tmp$bias

#bias-corrected

estimate

table.in[,4]=tmp$bias

table.in[,5]=tmp$var

table.in[,6:7]=tmp$conf.int

table.in[1:9,1]=paste("

",table.in[1:9,1],sep="")

table.in[,2]=ifelse(nchar(table.in[,2])==1,paste(table.in[,2],".",sep=""),

table.in[,2])

table.in[,2:7]=paste(table.in[,2:7],"000000",sep="")

table.in[,c(2:3,5:7)]=substr(table.in[,c(2:3,5:7)],1,6)

table.in[,4]=substr(table.in[,4],1,7)

table.in=paste(table.in[,1],"&",table.in[,2],"&",table.in[,3],"&",

table.in[,4],"&",

table.in[,5],"&",

table.in[,6],"&",

table.in[,7],"\\",sep="")

特别详细的一篇文章：http://wenku.baidu.com/link?url= ... DKZQinRo7kJB3_kEl3K