(* Data Envelopment Analysis *);



P = {{2, 4, 2, 3, 2}, {6, 12, 8, 12, 2}, {6, 8, 2, 8, 8}};

ie = IEfficiencyCRS[1, 2, P]; MatrixForm[ie]

ie4 = IEfficiencyCRS[4, 1, 2, P]



Y = TakeRows[P, 1]; X = TakeRows[P, -2];

MatrixForm[OEfficiencyCRS[X, Y]]



(1 - Part[ie4, 1]) Part[Transpose[X], 4]



Y.Part[ie4, Range[2, 6]] - Part[Transpose[Y], 4]