enjoying~~~
direc.dea <-
function(base = NULL, frontier = NULL, ngood = 1, nbad = 1)
{
b.mat <- as.matrix(base)
f.mat <- as.matrix(frontier)
if(ncol(b.mat) != ncol(f.mat))
stop("Number of columns of the base matrix and the frontier matrixis not identical!")
n<- nrow(b.mat)
n1<- nrow(f.mat)
s<- ngood
m<- nbad
l<- ncol(b.mat) - s - m