Topsis法在SAS宏中的统计实现_曾邦伟.pdf
(130.5 KB)
%macro topsis(database=,vara=,vark=,weight=0);/*database:被分析变量所在的数据集
vara:要分析的变量
vark:结果中保留的变量
weight:权重值可取0和1,默认去0*/
%let i=1;%let var&i.=%scan(&vara.,&i.,'');
%do %until(& &var&i.=);
%let i=%eval(&i.+1);
%let var&i.=%scan(&vara.,&i.,'');
%end;
%let var_num=%eval(&i.-1);
%if &weight. %then %do;
data top_1 top_2;
set &database. end=last;
if last then output top_2;
else output top_1;
run;
%end;
%else %do;
data top_1;
set &database.;
data top_2;
%do i=1 %to &var_num.; & &var. &i.; %end; run; %end;
data _null_;set top_1 end=last;
if last then call symput('obj_num',compress(_n_));
%let sqr_str=;
data top_a;set top_1 end=last;%do i=1 %to &var_num.;
& &var. &i.._2=& &var.&i.**2;
%let sqr_str=&sqr_str & &var&i..2;%end;
proc transpose data=top_a(keep=&sqr_str.) prefix=top_m out=top_b;run;
proc transpose data=top_a(keep=&vara.) prefix=top_n out=top_c;run;
proc transpose data=top_2(keep=&vara.) prefix=top_w out=top_w;run;
data top_d;set top_b;z=sum(of top_m1-top_m&obj_num.);
%let min_str=;%let max_str=;
data top_e;merge top_c top_d top_w;%do i=1 %to &obi_num.;z&i.=top_n&i./sqrt(z);
%end;
max=max(of z1-z&obj_num.);min=min(of z1 z&obj_num.);
%do i=1 %to %obj_num.;d_max&i.=top_w1*(z&i.-max)**2;d_min&i.=top_w1*(z&i.-min)**2;
%let max_str=&max_str d_max&i;%let min_str=&max_str d_min&i;%end;
proc transpose data=top_e(keep=&max_str.) prefix=max out=top_f;run;
proc transpose data=top_e(keep=&min_str.) prefix=min out=top_g;run;
data top_h(keep=&vark. d_max d_min c);merge top_1 top_f top_g;
max_sum=sum(of max1-max&var_num.);
min_sum=sum(of min1-min&var_num.);
d_max=sqrt(max_sum);d_min=sqrt(min_sum);C=d_min/(d_max+d_min);
proc rank data=top_h descending out=topsis_result;ranks order;var c;run;
proc print;title "The Topsis Report of &vark.";run;
%mend topsis;
data SF236;input city $ x1-x8;
format city $ 5. x1-x8 6.2;
cards;
CS 72.07 64.08 75.40 46.85 69.10 78.83 81.78 75.63
SZ 74.58 63.17 68.04 53.07 70.87 76.57 71.86 68.91
SY 77.79 54.77 72.32 61.83 69.35 73.65 63.41 71.75
GZ 82.01 71.03 73.10 58.93 68.83 73.00 72.27 72.64
CD 79.78 79.06 78.27 56.77 70.85 82.57 82.96 77.54
TJ 66.37 59.91 73.16 69.27 79.09 67.92 68.14 57.12
SH 67.11 76.15 72.37 54.63 68.13 77.22 77.71 78.63
;
run;
proc format;
value $ cityfmt CS='ChangSha' SZ='SuZhou' SY='ShenYang' GZ='GuangZhou' CD='ChengDu' TJ='TianJin' SH='ShangHai';run;
proc print data=SF236;format city $cityfmt.;title "Data of SF236";run;
%topsis(database=SF236,vara=x1-x8,vark=city,weight=0);



雷达卡





京公网安备 11010802022788号







