有谁知道CCR对偶模型的lingo程序？

DEA模型中CCR对偶模型的lingo程序。。。

lingo处理DEA的最基本的CCR模型，不是一个理想的选择，我建议还是用DEAP比较好。至于CCR的对偶模型，跟CCR模型有什么本质的区别吗？

CCR对偶模型可以得出系数值lamda在最优下的和，可以判断规模效率。DEAP我有尝试，但不知道dta为何物，如何数据输入。。。所以作罢。能指点一二么

关注。。。。。。。。。。。。。。

MODEL: ! Data Envelope Analysis of Decision Maker Efficiency ;
SETS:
  DMU/1..6/: ! Six schools;
     SCORE; ! Each decision making unit has a;! score to be computed;
  FACTOR/COST RICH WRIT SCIN/;! There is a set of factors, input & output;
  DXF( DMU, FACTOR):  F; ! F( I, J) = Jth factor of DMU I;
ENDSETS
DATA: ! Inputs are spending/pupil, % not low income;
! Outputs are Writing score and Science score;
  NINPUTS = 2;  ! The first NINPUTS factors are inputs;
!      The inputs,    the outputs;
F= 89.39        64.31        25.2        223
86.25        99        28.2        287
108.13        99.6        29.4        317
106.38        96        26.4        291
62.4        96.2        27.2        295
47.19        79.9        25.2        222;

ENDDATA
!----------------------------------------------------------;
! The Model;
SETS:
   ! Weights used to compute DMU I's score;
   DXFXD(DMU,FACTOR) : W;
ENDSETS
! Try to make everyone's score as high as possible;
  MAX = @SUM( DMU: SCORE);
! The LP for each DMU to get its score;
  @FOR( DMU( I):
   SCORE( I) = @SUM( FACTOR(J)|J #GT# NINPUTS:
    F(I, J)* W(I, J));
! Sum of inputs(denominator) = 1;
   @SUM( FACTOR( J)| J #LE# NINPUTS:
    F( I, J)* W( I, J)) = 1;
! Using DMU I's weights, no DMU can score better than 1;
   @FOR( DMU( K):
    @SUM( FACTOR( J)| J #GT# NINPUTS: F( K, J) * W( I, J))
     <= @SUM( FACTOR( J)| J #LE# NINPUTS:
      F( K, J) * W( I, J))
   )
  );
! The weights must be greater than zero;
  @FOR( DXFXD( I, J): @BND( .00001, X, 100000));
END

5# deman007
谢啦，不过这好像是元模型不是对偶模型的代码诶。。。您还知道对偶的代码么？

6# laurence2009
原模型代码应该是下面的吧：
MODEL:
sets:
DMU/1..6/:S,T,P;  ! Decision Making Unit;
II/1..2/:w;       ! Input Index;
OI/1..2/:u;       ! Output Index;
IV(II,DMU):X;     ! Input Variable;
OV(OI,DMU):Y;     ! Output Variable;
endsets
data:
p=?;
X=89.39 86.25 108.13 106.38 62.40 47.19
    64.31 99 99.6 96 96.2 79.9;
Y=25.2 28.2 29.4 26.4 27.2 25.2
   223 287 317 291 295 222;
enddata
[email=max=@sum(DMU:P*T]max=@sum(DMU:P*T[/email]);
@for(DMU(j):
S(j)=@sum(II(i):w(i)*X(i,j));
T(j)=@sum(OI(i):u(i)*Y(i,j));
S(j)>=T(j));
@sum(DMU: P*S)=1;
END

这些程序我也是从网上看到的，楼主能帮我具体解释一下这个DEA的规划模型吗？还有就是LINGO运行出来的结果是什么意思？lambda值是哪些？还有为什么会有48个松弛变量？万分感谢！如果解释内容多的话，可以发到我邮箱：qiudebin007@163.com

用matlab编程求解DEA方法中的CCR对偶模型http://www.pudn.com/downloads190/sourcecode/others/detail891919.html