标题: [算法] Topsis综合分析方法(SAS) [打印本页] 作者: 挚爱/Sunsgne° 时间: 2013-8-16 16:56 标题: [算法] Topsis综合分析方法(SAS) Topsis法是technique for order preference by similarity to ideal solution的缩写,意为与理想方案相似性的顺序选优技术,是系统工程中有限方案多目标决策分析的一种常用方法,可用于效益评价、卫生决策和卫生事业管理等多个领域[1]。此方法对样本量、指标多少及数据的分布无特殊要求和限制,灵活、方便、实用。但目前几种常用的统计软件包并没有提供现成的计算机程序,而现有已发表的关于Topsis法的SAS程序在应用于不同数据时,需要做较多修改,并且不能进行加权Topsis法的运算。本研究就Topsis法综合分析的SAS程序做了进一步的改进,使之更简单实用,以满足实际工作的需要。% m6 `% J. o2 D& `$ ^
! }' D" |& U+ h0 t" j# i5 }1 o 1 原理与方法" R( K0 g, U9 M) H/ K1 k
; z. N- _2 A$ z 基于归一化后的原始数据矩阵,找出有限方案中的最优方案和最劣方案构成一个空间:待评价的某方案可视为该空间上的一个点,据此可获得该点与最优方案和最劣方案问的距离(常用欧氏距离(Euclidean distance)),从而得出该方案与最优方案的相对接近程度,借此可进行方案优劣的评价[1]。其原理和方法如下:- ]: j. `1 Q% l' z5 {
, b0 A( V4 Z, V- X' N
① 同趋势化:通常采用低优指标转换为高优指标的方式,可采用倒数法和差值法。 : }9 L; p/ J+ g& O' S3 G $ h% {# m7 b+ K$ Z' i
② 无量纲化:为了消除计量单位的影响,需要做归一化处理。设(xij)n×m为同趋化后的指标矩阵,aij=xij/ ?ni=1 x2ij,j=1,2,…,m,建立则归一化矩阵(aij)n×m。6 Z4 W. `$ }, ~8 F) e
0 ^( j0 d' k8 \3 G
③ 求有限方案中的最优方案A+和最劣方案A- 。最优方案为:A+=(ai1+ ,ai2+ ,… ,aim+),最劣方案为:A-=(ai1- ,ai2- ,… ,aim-)。 , w- `4 Y8 ^& G7 \ ; k5 z5 P; J5 R/ ^7 |! U% F% ^4 { ④ 分别计算诸评价对象各指标值与最优方案及最劣方案的距离,通常采用欧氏距离。8 x' [& O3 J6 m( K
0 |( ?* F0 L/ A0 K6 F
D+i=?mj=1(aij-a+ij)2, D-i=?mj=1(aij-a-ij)2 " X( _5 o4 Q `) _4 \& f3 i$ L" | * R! E! m9 | h; U: N9 ?; e ⑤ 计算诸评价对象与最优方案的相对接近程度Ci。4 Q3 ?) E, W% z6 y2 Z6 l0 G
. Z: C1 p. w* v' F4 B Ci=D-i / (D-i+D+i)% ^1 X' b8 H9 e2 F" U) K5 O
6 u2 z" D6 b8 i7 g$ D* c
2 实例分析与SAS宏程序 / D9 q5 E& w. P8 `3 D/ E) V # ` s$ K6 \( k+ W K: Y. ? P
以中文版SF-36为健康调查问卷,于2004年7月~2004年12月间对≥60岁老年人(602人)进行老年人生存质量(QOL)的研究,问卷包含8个领域的指标,分别为x1(生理机能)、x2(生理职能)、x3(躯体疼痛)、x4(一般健康状况)、x5(精力)、x6(社会功能)、x7(情感职能)、x8(精神健康)[2]。得到数据如下:3 x2 F# Z$ g8 T) S4 W3 e
0 N! M8 j1 E5 {7 `! B d5 }' I 表1 7个地区老年人QOL各领域平均分(略), V ~9 O8 C' ` Q9 K5 B4 P+ _
x) |! L x; ]6 p5 Q# C6 ?# P 以下是宏程序:! }& r% h& ~5 n. u5 ^2 T
( c) [9 M+ H( v/ o3 i0 h, A
%macro topsis(data_name=,var_name=,var_keep=&var_name,weight=0); $ V1 i: g! D; X* ^0 F
7 x% ~6 t7 E4 q; n
%let i=1; 8 ?1 U6 o" p5 L" W : Z$ S7 d: ?. m, Q# f! X0 W
%let var&i=%scan(&var_name,&i,' ');: s2 z8 G6 ]! M& F% v- ]
2 P. ?/ c- n) M; m1 \4 M# Z %do %until (&&var&i=); % t8 z7 b5 V9 ~9 |* L' \* k5 N
- W* B& r! r2 k& s, U
%let i=%eval(&i+1); : }: [4 z; s2 K* J% a ( l* Q6 |+ @0 i %let var&i=%scan(&var_name,&i,' '); " ~' l( |+ w" M) d 1 e- w& h I% H" y+ x- @& [
%end; / j6 N$ T+ O9 p0 f8 l% z! G 1 \ r/ C/ ]/ M2 @ %let var_num=%eval(&i-1); *自变量个数;; v1 g0 r4 ~9 C: ~4 R7 Y9 S( }
( P: p [8 L& ^$ p( G! C
*top_1为原始数据,top_2为权重数据。必须把权重数据放在最后一行; 8 B) Y4 i9 I9 K! |( V# } , t2 E3 |+ @, M+ K %if &weight=1 %then %do; 6 X+ F% D7 i4 S+ p + r' b8 _( j& c! Q" `6 u' a) V
data top_1 top_2;/ S; `& g& @% l4 P
* I9 H/ p) w/ R% j
set &data_name end=last;- L8 @" h" n& K1 e# w
a8 {( a- [: N) W+ |) M if last then output top_2;9 ]$ Y3 O. A$ T+ M/ h" ]
& P$ n5 ~' |& K R5 U7 G5 J* \
else output top_1;9 [% ~- a* Q+ A' }+ I- w
' r1 x8 \- ~# H/ J& H run;%end;/ B. }, I9 E b. v0 A! }1 C
" w$ G- N5 }: ?2 [+ a
%else %do;8 ^! B: P+ G! K- Y* o2 V: L) _
: N+ O6 G+ N% g; s* C: p# `1 o5 n data top_1;( w1 r8 m7 b+ F5 U6 R
# P0 ?1 D; {% d! _ set &data_name;' c# d- K( B" S; Z7 L3 P
. J( O/ C/ ?4 {: V data top_2; # q, i5 ^) Z5 \% J : E7 G% s& D. c0 W" H
%do i=1 %to &var_num;7 p9 P' z4 V( y* V; p) i4 c- I
$ @" ]; j/ _7 q1 |: g &&var&i=1; & a/ H# }2 q6 D8 o- x 1 o" P. r& b* n7 f; ^- l, c %end;% g4 D8 I. g" E1 A) V
0 {) \8 k9 g$ z) e& |
run;%end; 5 ]2 u; O9 b8 v4 P% R $ m7 }. @/ J4 N# O1 Q7 W) V data _null_;*样本个数;: A. \* M) l0 h
( Q6 Y) V7 U3 ?! u& n( b
set top_1 end=last; 8 o" R9 L) E* c! T; O4 o 4 d8 J$ k# _2 v0 N0 g if last then call symput('obj_num',compress(_n_)); 5 R. w$ u. F- I7 [/ r b0 p: e5 K N5 x, s8 K. Q& e* T" J %let sqr_str=; , G$ s- E/ `3 U9 y$ \ 9 s! k" Y1 Y t
data top_a; *求各指标值的平方; ' U. s5 J* I' C 1 o0 A: x W4 C1 \/ M6 y6 O/ S2 R set top_1 end=last;1 ?3 G" `) w4 r0 O+ n
3 M2 L5 R. H- X! ~ %do i=1 %to &var_num; / D( S% k7 i2 K5 q2 ~/ q N4 H3 u 1 d- h& L- ~; O) D7 D &&var&i.._2=&&var&i ** 2;* T! K. x4 v3 ] J0 j) w# B$ A
* X- \& a2 g& H2 ?% \
%let sqr_str=&sqr_str &&var&i.._2;%end; 1 K m y' g# \9 ^ ( A" B- J( c- e g& U proc transpose data=top_a(keep=&sqr_str)prefix=top_m out=top_b; *对平方值矩阵转置; ! O0 O) P! _0 y7 p( ] 6 D3 {* W% l6 J( |5 r% y2 W# R9 a; a proc transpose data=top_a(keep=&var_name)prefix=top_n out=top_C; *对原始数据矩阵转置;; ?7 p% N, |6 n) F0 j. q
% r" k w, y. D( _
proc transpose data=top_2(keep=&var_name)prefix=top_w out=top_w; *对权重数据矩阵转置; 3 {" \. x3 P7 S1 h2 r ! N; P* F: ^/ W/ ~9 R/ A. n, X data top_d; *求各指标值的平方和; " |6 U; L7 i. @% S6 e( U 2 A; \, j( F# m5 y+ L
set top_b;) z, q+ }' M2 [0 h
/ B" u& W# Y! z$ q9 _- `
z=sum(of top_m1-top_m&obj_num); ! F( c- a0 N% C# H# I . Z0 Y! j* p+ N T, J% D9 i. o %let min_str=;%let max_str=;0 ?* ?/ X, L# D8 [/ g
+ c( B" `# R' F8 C
data top_e; ! n' }: n3 J* Y. R4 C, { 3 e# U6 ?3 H" p/ z6 d/ g& O merge top_c top_d top_w; 6 q" e& ?% ?% X. g ' t3 q5 o+ `0 j& `0 ^3 i
%do i=1 %to &obj_num; *求归一化矩阵的转置矩阵; 5 U8 F: G$ x4 H2 l $ S% N0 |/ P; e! V1 ^ z&i=top_n&i/sqrt(z);%end;+ g2 a( h: v. Z4 ]3 s
6 z6 I" ?; p- h( W& m
max=max(of z1-z&obj_num); *求最优方案和最劣方案;8 d" n, v8 z! {4 }0 G- {! s
# Z6 X7 q! V7 `! e5 s min=min(of z1-z&obj_num);$ k& d R+ J" ?6 h0 p9 t2 t. [
2 y: V2 z- Q( y& K %do i=1 %to &obj_num;# _ O* g& o" J) i+ h
) S. ]+ V/ M/ L d_max&i=top_w1*(z&i-max)**2;2 _7 R5 W& N+ C
, K m, u3 I- j, _$ H d_min&i=top_w1*(z&i-min)**2;, M- ?1 }% V2 |/ ~8 {
+ p. ^( V' c# O' \$ k- X %let max_str=&max_str d_max&i;5 U* ]" Q4 K( S; }2 Y. h: j2 [! X, x4 Q
3 x6 X2 r7 ?2 ^. k k; L3 E
%let min_str=&min_str d_min&i;%end; P$ @6 ~' E# Z1 w# ?$ L / L: y# J& E" w2 {2 y
proc transpose data=top_e(keep=&max_str) prefix=max out=top_f; ( S9 I3 y& I* k4 z# f 0 n5 S: T4 G, {9 k4 V* L proc transpose data=top_e(keep=&min_str) prefix=min out=top_g;*求各评价对象与最优方案的相对接近程度C,并排序;( l& H& X1 t( U; N. g+ \ r
" J8 d: Y$ ]; i data top_h(keep=&var_keep d_max d_min c);: ^2 `7 M, w$ n/ X, V
" x# q* K5 O8 S9 `& P' _
merge top_1 top_f top_g ;7 z$ w) R: }& D' P0 H3 C
/ Q* i. X9 J. l& l max_sum=sum(of max1-max&var_num);4 o) X' ]; Y6 _- D% m
# [ k2 S+ g/ S2 R- ` min_sum=sum(of min1-min&var_num); * w7 H7 i, o5 X# o0 W7 B7 s P( A/ G8 [# B6 A d_max=sqrt(max_sum); I% D% \" Q4 T1 d& y
9 v! e" ]& a. T6 h5 i: D4 J d_min=sqrt(min_sum);6 X4 T% g ?/ h1 [! A
" _& P5 [9 t/ r* a8 I C=d_min/(d_max+d_min);( M- F+ H; @0 G/ S6 q$ p
+ ?+ e9 D1 a% ?' x- N6 m proc rank data=top_h descending out=topsis_result; : c0 m$ F5 E( y& _ " i4 m/ r" g; j/ A7 }- l$ X
ranks order; / I7 Y( P4 N' Z7 ~. v% Y. g! h , B% l& g k+ m
var c; ' z) } t( g6 B6 U: P. \ 4 U( l& m" m6 k& w. s9 z proc print;run ; ' R& T% }( Z1 C$ F) Q# c1 \* A2 |9 W: K) n: t. {
%mend topsis; ( t$ {$ i0 P* ]" T1 G- M* D1 D + @' e4 e6 R1 t& Y 设数据集名为topsis,地区变量名为city,var_keep为需要打印出的变量,weight为是否进行加权Topsis法运算,默认为0(否)。运行宏程序如下: 6 E- s' z+ l' h v+ s/ |1 Y' ]) p$ i9 Y. c
%topsis(data_name=topsis,var_name=x1 x2 x3 x4 x5 x6 x7 x8,var_keep=city,weight=0); * ?0 l4 Z+ A) _ ) @) f4 H; J8 y! M
运行结果如下: # T3 X8 T1 Q9 a9 E) v3 \6 D" C* J* x3 ]6 t$ m9 S
Obs CITY d_maxd_minCorder1长沙0.186680.161690.4641442苏州0.178900.115590.3925173沈阳0.195440.142200.4211564广州 " x- q Z3 f1 V9 h& {( ?" ]( n 0.128570.175710.5774625成都0.093750.238970.7182316天津0.205630.164550.4445157上海0.144700.194060.572863 " a8 b! A/ e4 c7 r- l$ ~! D' |1 M7 t' ~0 K8 f5 [/ q. H. h( c6 Q( i
3 讨论- W+ I% B) k# m9 p# O+ [8 a3 P
0 ~' }- W- \" L( N2 z
Topsis法对原始数据进行同趋势和归一化处理,从同一指标与最优指标的距离之比计算中体现出同一指标间的接近度,这在一定程度上反映了散程度,排序结果充分利用原始数据信息,能定量反映不同评价单元的优劣程度,直观、可靠。 " A* N6 K2 i. N4 Y4 X6 u9 h " I0 y8 i: C- X 相对于之前已发表的Topsis法的SAS程序[3,4],本程序的特点是对不同的数据,只需做最少的修改(数据集名称和变量名)即可完成Topsis法的运算。另外,本程序还能进行加权Topsis法的运算。确定权重的常用方法有专家评价法、因子分析法、模糊综合评价的逆估计法、统计计算法等。只需将确定的权重加入到数据集的最后一行,并将参数weight=1即可完成加权Topsis法的评价。 9 K5 w' [: @. h作者: 亦久亦旧。 时间: 2013-8-17 17:13
赞一个