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="")