我帮你用MAXDEA测考虑非期望的DEA效率值

ata = csvread('C:\Users\Administrator\Desktop\论文\deaml\data.csv');

% （id,xm,year)
X = zeros(6,3,19);
for i = 1:114
     X(fix((i-0.1)/19)+1,:,mod(i-1,19)+1) = data(i,1:3);
end

Y = zeros(6,1,19);
for i = 1:114
     Y(fix((i-0.1)/19)+1,:,mod(i-1,19)+1) = data(i,4);
end

Yu = zeros(6,1,19);
for i = 1:114
     Yu(fix((i-0.1)/19)+1,:,mod(i-1,19)+1) = data(i,5);
end
%%

ml = deamalmluen(X, Y, Yu,'orient','ddf')
deadisp(ml);  

T = dea2table(ml);
writetable(T, 'C:\Users\Administrator\Desktop\论文\deaml\dea_results.csv')