下一个包，FEAR，在这：http://www.economics.clemson.edu/faculty/wilson/Software/FEAR
看你具体要做什么，有时候还需要nonparaeff
给你一个例子，参考一下——
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(ccrx1,ccrx2,ccrx3,ccrx4,ccr$x5),
nrow=70,ncol=5))
y=t(matrix(c(ccry1,ccry2,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
扫码加好友，拉您进群
全部版块
我的主页
收藏
