Реинжиниринг имитационной модели работы железной дороги на основании данных из Информационного Фонда производственно-экономических показателей

quit;

data EK.pyt ;

set ek.pyt;

drop fact fact_fot ;

run;

proc sql

create table ek.p330_it as

select distinct b.*, c.fact, c.fact_fot

from ek.pyt as b, ek.q322_it as c

where c.stati=b.stati and c.dv_id=b.dv_id and c.var_id_izm=b.var_id_izm;

quit;

data EK.pyt;

set kat.stat_1 ;

if stati='333, 334' then output;

run;

data ek.p2;

set kat.factor_first;

rename fact=fact_p;

d=year(date);

if d=&d1 and (var_id=50040 or var_id=37300 or var_id=50050 or

var_id=37310) then output;

drop d;

run;

proc sort data=ek.p2;

by date dor_id operiod var_id;

run;

proc transpose data=ek.p2

out=ek.p_3;

by date dor_id operiod ;

* var fact_p;

run;

data ek.p21;

set ek.p_3;

rename col1=l_k_g;

rename col2=l_c_p;

rename col3=n;

rename col4=l_g;

drop _name_;

run;

data ek.p21;

set ek.p21;

alpr=l_c_p*0.4+0.05*n;

ak2=(l_g+l_k_g*0.75)/alpr;

ak2=round(ak2,.001);

run;

proc sql;

create table ek.p_330 as

select distinct b.*, c.*

from ek.pyt as b, ek.p21 as c

where c.dor_id=b.dor_id ;

quit;

data ek.p_330;

set ek.p_330;

e_g=fact*ak2;

fe_g=fact_fot*ak2;

keep stati e_g fe_g dv_id;

run;

proc sort data= ek.p_330;

by stati dv_id;

run;

proc transpose data=ek.p_330

out=ek.p330;

by stati dv_id;

var e_g fe_g;

run;

data ek.p330;

set ek.p330;

rename col1=fact;

rename _name_=var;

drop col2 col3 col4;

run;

data ek.p330_f;/*к условно-пост части*/

set ek.p330;

if var='fe_g' then var_id_izm='50000' ;

rename fact=fact_fot;

if var='fe_g' then output;

run;

data ek.p330;

set ek.p330;

if var='e_g' then var_id_izm='50000' ;

if var='e_g' then output;

drop var;

run;

proc sort data=ek.p330_f;

by stati dv_id var_id_izm;

run;

proc sort data=ek.p330 ;

by stati dv_id var_id_izm;

run;

data ek.p330_11 ;

merge ek.p330_f ek.p330 ;

by stati dv_id var_id_izm;

run;

data EK.pyt ;

set ek.pyt;

drop fact fact_fot ;

run;

proc sql

create table ek.p333_it as

select distinct b.*, c.fact, c.fact_fot

from ek.pyt as b, ek.p330_11 as c

where c.stati=b.stati and c.dv_id=b.dv_id and c.var_id_izm=b.var_id_izm;

quit;

proc sort data=ek.p330_it;

by stati dv_id var_id_izm;

run;

proc sort data=ek.p333_it ;

by stati dv_id var_id_izm;

run;

data ek.p330_333 ;

merge ek.p330_it ek.p333_it ;

by stati dv_id var_id_izm;

run;

proc sort data=ek.p330_333 ; /*итоговая */

by stati dv_id var_id_izm;

run;

3.8 Листинг SCLкодов интерфейса 1

init:

t3.text=1.331;

t4.text=1.474;

t5.text=1.992;

t6.text=1.402;

t7.text=1.662;

submit continue

proc sql;

create table work.nazvdor as

select distinct name

from ref.dor;

quit;

endsubmit;

declare list nazvaniya;

nazvaniya=makelist();

tab=open('work.nazvdor','i');

do while(fetch (tab) ne -1);

n=getvarc(tab,1);

rc=insertc(nazvaniya, n);

end;

rc=close (tab);

list1.items=copylist(nazvaniya);

call putlist(nazvaniya);

return;

pushbutton1:

d1=t1.text;

call symputn('d1',t1.text);

d2=t1.text;

call symputn('d2',t2.text);

F=t3.text;

call symputn('F',t3.text);

M=t4.text;

call symputn('M',t4.text);

A=t5.text;

call symputn('A',t5.text);

call symputn('d1',t1.text);

E=t6.text;

call symputn('E',t6.text);

T=t7.text;

call symputn('T',t7.text);

dor1=list1.selecteditem;

dor1="'"||dor1||"'";

submit continue;

%let dor_id=&dor1;

proc sql;

select distinct dor_id into:ww

from ref.dor

where name=&dor1;

quit;

%include 'c:\dip\factor1.sas';

endsubmit;

return;

pushbutton3:

call display ('kat.prog.fact.frame');

return;

pushbutton2:

call display ('kat.prog.pokaz1.frame');

return;

pushbutton4:

call display ('kat.prog.izm.frame');

return;

pushbutton5:

call display ('kat.prog.factors1999.frame');

return;

pushbutton6:

call display ('kat.prog.znach.frame');

return;

pushbutton7:

call display ('анин фрейм');

return;

exit:

_STATUS_='H';

return;

3.9 Листинг SCLкодов интерфейса 2

init:

dcl list name_f value_f;

name_f=makelist();

value_f=makelist();

tab=open('kat.factor_first');

do while (fetch(tab) ne -1);

namenum=varnum(tab,'name');

a=getvarc(tab,namenum);

rc=insertc(name_f,a);

valnum=varnum(tab,'fact');

b=getvarn(tab,valnum);

rc=insertn(value_f,b);

end;

factors.items=name_f;

a1=-1;

rc=close(tab);

return;

factors:

if a1 >0 then do;

rc=SETITEMN(value_f,val_f.text,a1);

end;

a1=factors.SelectedIndex;

c1=getitemn(value_f,a1);

val_f.text=c1;

return;

upd:

if a1 >0 then do;

rc=SETITEMN(value_f,val_f.text,a1);

end;

submit continue;

data kat.data;

set kat.factor_first;

delete;

run;

endsubmit;

tabw=open('kat.data','UN');

aaa=listlen(name_f);

put aaa;

do i= 1 to listlen(name_f);

rc=initrow(tabw);

namenum=varnum(tabw,'name');

a=getitemc(name_f,i);

call putvarc(tabw,NAMENUM,a);

valnum=varnum(tabw,'fact');

b=getitemn(value_f,i);

call putvarn(tabw,valnum,B);

rc=append(tabw,'NOINIT' );

end;

rc=close(tabw);

submit continue;

proc sql;

create table kat.factor_first as

select distinct b.date, b.dor_id , b.operiod, b.var_id, c.*

from kat.factor_first as b, kat.data as c

where c.name= b.name;

quit;

endsubmit;

call display('fact.FRAME');

return;

pushbutton1:

submit continue;

%include 'c:\dip\izmerit.sas';

endsubmit;

return;

exit:

_STATUS_='H';

return;

3.10 Листинг SCLкодов интерфейса 3

pushbutton1:

submit continue;

%include 'c:\dip\statya.sas';

endsubmit;

call display('factors1999.FRAME');

return;

exit:

_STATUS_='H';

return;

Размещено на Allbest.ru