Model state 8 ?!

130 views
Skip to first unread message

Saeid Abbasi Parizi

unread,
Apr 25, 2016, 2:47:11 AM4/25/16
to gamsworld
Hi dear friends

I will glad if anyone can help me about model stat 8. what is means?
it means that a feasible integer solution is founded for the model or optimal integer solution can be achieved?

kind regard
Saeid

Vafa Mousavi

unread,
Apr 26, 2016, 8:07:30 AM4/26/16
to gamsworld
سلام مهندس
به این معنی است که مدل شما با اینکه
MLP
هست ولی به طور اتفاقی جواب شما عدد صحیح شده است که این جواب عدد صحیح ممکن است با توجه به مقادیری که برای
optcr و  optca
داده اید این طوری عدد صحیح شده است.
حالا شما به بنده جواب بدهید
بنده مدل تحلیل پوششی داده ها را نوشتم ولی کارایی همه واحدها را یکسان می دهد و ارور هم نمی دهد و وقتی هم ورودی ها و یا خروجی ها را دستکاری می کنم کارایی همه واحدها با هم یکسان تغییر می کنند و باز هم ارور نمی دهد
مدل اصلی  که کارایی همه واحدها را یک می دهد و ارور نمی دهد به شرح زیر است:
sets
i   "inputs"   /i1  "doctors",i2  "nurses"/
r   "outputs"  /o1  "outpatients",o2 "inpatients"/
j   "units"    /dmu01*dmu12/
is(j) "selected units"
scalar epsilon  "non-archimedean value"  /0.005/;
parameters
*let dmu01  be under evaluation
   xo(i)   "inputs of under evaluation dmu"  /i1 20,i2  151/
   yo(r)   "outputs of under evaluation dmu" /o1 100,o2 90/;

table x(i,j)
    dmu01 dmu02 dmu03 dmu04 dmu05 dmu06 dmu07 dmu08 dmu09 dmu10 dmu11 dmu12
i1   20    19    25    27    22     55    33    31    30    50    53    38
i2  151   131  160  168  158    255   235  206  244  268  306   284;
table y(r,j)
    dmu01 dmu02 dmu03 dmu04 dmu05 dmu06 dmu07 dmu08 dmu09 dmu10 dmu11 dmu12
o1   100   150   160   180   94    230   220   152   190   250   260   250
o2   90     50     55     72    66    90     88     80    100   100   147   120;
variables
v(i) "output weights"
u(r) "input weights"
z    "efficiency";
positive variables
v
u;
v.lo(i)=epsilon;
u.lo(r)=epsilon;
equations
objective
const(j);
objective..z=e=(sum(r,u(r)*yo(r)))/sum(i,v(i)*xo(i));
const(j)..sum(r,u(r)*y(r,j))/sum(i,v(i)*x(i,j))=l=1;
model   fractionalccr_model /all/;
solve    fractionalccr_model using nlp maximizing z;
set jj(j) set of units to analyze /dmu01*dmu12/;
parameters efficiency summary report;
loop(jj,is(jj)=yes;
solve fractionalccr_model using nlp max z;
 efficiency(jj)=z.l;
is(jj)=no;);
option decimals=4;
display efficiency ;
مدلی که تغییر دادم فقط در جدول ورودی ها جدول ایکس واحد 1 را شامل 21 دکتر و 100 پرستار کردم و همه کارایی ها را برابر 0.7875 کرده است و ارور هم نمی دهد.
مشکل کار کجاست که همه کارایی ها یکسان می دهد و ارور هم نمی دهد؟

Ali Reza Bahari

unread,
Apr 29, 2016, 2:48:25 AM4/29/16
to gams...@googlegroups.com
Hi My friends Saeid and Vafa
mr. saeid please set the following statements befor solve statemnet to get optimal soloton 
option optcr=0;
option optca=0;


خانم وفا با توجه به تعریف xo و yo مدل شما کارایی واحد اول را برای همه واحدها محاسبه میکند. در واقع با عوض شدن واحد تحت ارزیابی داده های تابع هدف متناسب با آن عوض نشده و کاریی فقط برای اول محاسبه میشود برنامه اصلاح شده به صورت زیر است
Sets
i   "inputs"   /i1  "doctors",i2  "nurses"/
r   "outputs"  /o1  "outpatients",o2 "inpatients"/
j   "units"    /dmu01*dmu12/
is(j) "selected units";
is(j)=yes;
Alias(j,jj);

scalar epsilon  "non-archimedean value"  /0.005/;

table x(i,j)
    dmu01 dmu02 dmu03 dmu04 dmu05 dmu06 dmu07 dmu08 dmu09 dmu10 dmu11 dmu12
i1   20    19    25    27    22     55    33    31    30    50    53    38
i2  151    131   160   168   158    255   235   206   244   268   306   284;
table y(r,j)
    dmu01 dmu02 dmu03 dmu04 dmu05 dmu06 dmu07 dmu08 dmu09 dmu10 dmu11 dmu12
o1   100   150   160   180   94    230   220   152   190   250   260   250
o2   90    50    55    72    66    90     88   80    100   100   147   120;
variables
v(i) "output weights"
u(r) "input weights"
z    "efficiency";
positive variables
v
u;
v.lo(i)=epsilon;
u.lo(r)=epsilon;
equations
objective
const(j);
objective(is) ..z=e=sum(r,u(r)*y(r,is))/sum(i,v(i)*x(i,is));
const(j)      ..sum(r,u(r)*y(r,j))/sum(i,v(i)*x(i,j))=l=1;
model   fractionalccr_model /all/;
solve    fractionalccr_model using nlp maximizing z;
parameters efficiency summary report;

is(j)=no;

loop(jj,
  is(jj)=yes;

solve fractionalccr_model using nlp max z;
 efficiency(jj)=z.l;

is(jj)=no);

option decimals=4;
display efficiency ;




--
You received this message because you are subscribed to the Google Groups "gamsworld" group.
To unsubscribe from this group and stop receiving emails from it, send an email to gamsworld+...@googlegroups.com.
To post to this group, send email to gams...@googlegroups.com.
Visit this group at https://groups.google.com/group/gamsworld.
For more options, visit https://groups.google.com/d/optout.

Reply all
Reply to author
Forward
0 new messages