运用SAS程序实现层次分析法的综合评价

1、运用sas程序实现层次分析法的综合评价【摘要】目的：用sas程序帮助实现层次分析法的综合评价，从 而简化运算过程，提高计算的准确度。方法：对层次分析法中权重系 数计算和检验过程编sas程序。结果：将原始判断优势矩阵输入sas 程序中的相应数组后，运行程序便可很快得到归一化权重系数和随机 一致性比率，并得出各评价对象综合评分指数和排序结果。结论：笔 者编写的层次分析法综合评价sas程序具有通用、简便、实用和可操 作性强的特点，在实际综合评价工作实践中具有很强的实用价值。【关键词】层次分析法；权重系数；矩阵；综合评价；sas层次分析法(analytic hierarchy process,简称ah

2、p)由美国科 学家t. l. saaty于20世纪70年代提出，是用系统分析的方法，对评 价对象依评价目的所确定的总评价目标进行连续性分解，得到各级(各 层)评价目标，并以最下层作为衡量目标达到程度的评价指标。然后依 据这些指标计算出一综合评分指数对评价对象的总评价目标进行评 价，依其大小确定评价对象的优劣等级1。层次分析法是综合评价中 较为重要的方法,是分析一类相互关联、相互制约众多因素构成复杂系 统的系统分析方法。目前该法多用于卫生事业管理方面，如医院工作 质量的评价。虽然层次分析法优点较多，目前的运用变得越来越普遍，但 是其计算过程较为烦锁，表现在计算公式多，且需要重复多次计算。 多步骤

3、的四舍五入而累积的计算误差可能导致结果产生较大的偏甚至 错误。因此为层次分析法的计算及检验过程偏写sas程序以简化运算 并提高计算的准确性就显得很有必要。以下就以医院工作质量综合评 价为例来演示运用sas程序实现层次分析法的综合评价。1层次分析法综合评价法步骤及计算公式示例1：对医院工作质量进行综合评价的基本步骤为：1.1建立目标图对总评价冃标进行连续性分解以得到不同层次的评价冃标， 将各层评价目标用图有机地标示出来，即建立目标树图。例如对某个 综合医院（评价对象）的工作质量（总评价目标）进行评估，按系统分析 的方法，医院工作质量这个总评价目标可通过医疗工作，护理工作、 膳食供应等三个次级目标

4、来反映，而这三个次级目标又可通过各自的 次级冃标来反映，像医疗工作可通过医疗制度、医疗质量、病床使用 等三个次级目标反映。如此分解下去便可建立一个医院工作质量的目 标图，如图1所示。图1医院工作质量评价口标树图1.2计算权重系数在同一层评价0标依据各个评价b标对总评价0标作用价值 的大小分别赋予一定的系数。例如医院工作质量而言，医疗工作最为 重要，应赋予较大的权数，护理工作次z,膳食供应更次z,因而权 重也小一点。权重计算步骤为：1.2.1对目标树自上而下分层一一对比打分，建立成对比较 判断优选矩阵，评分标准见表1。表1冃标树图各层次评分标准以口标树中第1个子ii标分层为例，3个评价口标成对比

5、较 判断优选矩阵见表2。表2第层了目标成对比较判断优选矩阵1. 2. 2按公式wz i=mail, ai2, , aim计算初始权重系数 wz i得wz 1=31 x 3 x 5=2.4662 ,同理得 2=1.0000 , wz 3=0. 40551.2.3按公式wi二射i /mi二i计算归一化权重系数wi得wl=2.46622. 4662+1. 0000+0. 4055=0. 6370同理得 w2 =0. 2583 , w3=0. 1047用同样方法可获得其它分层屮项目指数的权重系数。1. 3对权重系数是否符合逻辑进行检验在确定权系数过程中，靠人的主观判断而给出的判断矩阵， 还必须通过一致

6、性检验，以便尽可能消除人的主观判断的不一致性 。通用常用一致性指标ct检验该项目的相对优先顺序有无逻辑混 乱，一般认为,当ci<0.01吋，可能逻辑混乱，即计算得的各项权 重可以接受。其中ci= x max-mm-1, 入 max二 mi=l x i / m , 入 i二 mj二laijwj / wi式中m为受检验层次的子目标数，入max为最大特征根，x i 为该层子冃标成对比较判断优选矩阵的特征根。为了度量不同阶段矩 阵是否具有满意的一致性，我们还需引入判断矩阵的平均随机一致性 指标rt值。通常采用美国运筹学家saaty教授创立的19级标度法 3。对于19阶判断矩阵，ri值见表

7、3。 表3 19阶平均随机一 致性指标ri的取值对于1, 2阶判断矩阵，ri只是形式上的，因为1, 2阶判断 矩阵具有完全一致性。当阶数大于2时，判断矩阵一致性指标ci与同 阶平均随机一致性指标ri z比称为随机一致性比率，记为crocr=ciri ,当cr<0. 1吋，即可认为判断矩阵具有满意的 一致性，否则就需要调整判断矩阵，并使之具有满意的一致性。本例,x 1=(1/ 0. 637+3x0. 2583+5x0. 1047) / 0. 6370=3.0383x 2= (1/3 x 0. 637+1x 0. 2583+3 x 0. 1047) / 0. 2583=3.0298x

8、 3=(1/5 x 0.637+1/3 x 0.2583+1 x 0. 1047)/0. 1047=3. 0298入 max二(3. 0383+3. 0298+3. 0309)/3=3. 033ci=(3. 033-3)/(3-1) =0. 0165<0. 1cr=o. 0165/0. 58=0. 0284<0. 1认为第一层子目标各项判断无逻辑错误。依此类推，运用此 法可检测各层子目标权重系数是否合乎逻辑。1.4计算各评价指标权重系数及进行检验的sas程序/*计算各子目标归一化权重系数并检验是否符合逻辑*/data ahp;%let a-3: /*a为判断矩阵层数

9、*/array s (& a) ; /*定义初始权重系数数组*/array w (&amp: a) : /*定义归一化权重系数数组*/array z (& a, & a)(1 3 5 0.33 1 3 0.2 0. 33 1) ; /*判断优选矩阵原始数据*/array t (& a) ; /*定义各指标特征根数组*/array ri (9) (0 0 0. 58 0.9 1. 12 1.24 1. 32 1.41 1.45);/*19价平均随机一致性指标rt的取值*/do i=l to &amp:a:p二 1；d

10、o j二1 to &a;p二p*z(i, j);s(i)=p*(l/&a) ; /*计算初始权重系数*/end; end;do i二1 to &a;w(i)=s(i)/sum (of s(*); /*各子目标归一化权重系数*/ do i=l to &a;num二0;do j二1 to &a;num二num+z(i, j)*w(j)/w(i);end;t(i)=num; /*该层子目标成对比较判断优选矩阵的特征根*/end;tmax=sum (of t (*)/& a; /*最大特征根 */ci= (

11、tmax-& a) / (& a-1) ; /*一致性指标 ci*/ cr=ci/ri (&a) ; /* 随机一致性比率 cr*/run;2计算综合评分指数，对评价对象的总评价冃标进行综合评估2. 1计算综合评分指数gi并对其进行排序首先手工计算各评价指标组合权重系数c* ,它等于各评价 指标（即最后子冃标）与其所属各上级子冃标归一化权重系数的连乘 积。综合评分指数gi按下列公式求得。gi二micipi ,式屮pi为第 i个评价指标的测量值，m为评价指标的个数。某市6所医院综合评价数据见表4。表4某市6所医院工作 质量的综合评价表评价指标pi组

12、合权重ci医0. 189295. 092. 094. 895. 689. 177. 40. 229288. 191. 290. 094. 093. 692. 20. 114515.48. 37. 93. 19. 53. 7 医0. 104174. 753. 461. 950. 061. 967. 10. 258354. 720. 726. 120. 027. 435. 30. 104741. 341. 422. 820. 034. 030. 3治疗有效率症收容率疗制度执行优良率护理制度执行优良率膳食供应优良率综合评i分指数院abcdef病床使用率gi66. 254. 555. 052. 556

13、. 555. 5 工作质量顺位 1546232. 2综合评分指数gi计算及排序sas程序/*计算工作质量评分指数gi*/data score;%let b=6; /*b为评价单位数*/%let c=6; /*c为评价指标数*/array score(&b);array cw (& b) (0. 1892 0. 2292 0. 1145 0. 1041 0. 25830. 1047 );/*最后子目标与其所属各上级子目标归一化权重系数的连乘积*/array percent(&b，&c)(95 92 94. 8 95. 6 89. 1

14、 77. 488. 1 91.2 90 94 93.6 92.215.4 8. 3 7. 9 3. 1 9. 5 3. 774. 7 53.4 61.9 50 61.9 67. 154. 7 20. 7 26. 1 20 27. 4 35. 541. 3 41.4 22.8 20 34 30. 3)；apercent数组为各评价对象各评价指标的原始数据*/ array rank(&b);do i二1 to &b;num二0;do j=l to &amp:c:num二num+cw(j)*percent(j, i);score(i)=num; ascore数

15、组为各评价对象综合评分指数gi*/end ;end;/*对gi指数进彳亍排序*/do i二1 to &b; rank(i)=l;do j二1 to &b;if score (i)=score (j) thenif i ne j then rank(i)=rank(i)+ 0.5;if score (i) end;end;/*rank数组内各元素为各评价对象对应的综合评价名次，遇 相同秩次时取平均秩*/run;3讨论3. 1层次分析法的优点缺点评价层次分析法较其他方法更为合理地解决了权重的计算，各项 指标对比等级划分较细,在综合评价吋充分显示权重的重要；对原始数

16、据直接加权计算出综合评价指数，充分利用原始信息量，使评价结果 更具切实性、合理性;建立层次结构是定性指标定量化，解决指标性质 不一致难以比较的问题;运用一致性指标(ci)检验以判断确定各项指 标的权重有无逻辑上的混乱，具较强的科学性4。但是层次分析法在 应用过程中确实也存在一些不足和弊端。主观判断难免会存在各种偏 差，特别在指标较多的情况下，偏差可能越大，判断所得的判断矩阵 就很难保证具有满意的一致性。另外，层次分析法的信息基础是判断 矩阵，判断矩阵的一致性要求是层次分析法应用的关键前提，当评价 者所得判断矩阵一致性较差时，需要结合研究目的来确定该判断矩阵 是舍弃还是进行必要的检验与修正5。3

17、.2笔者编写的层次分析法综合评价sas程序具有通用、简 便、实用和可操作性强的特点，具有一定的实用价值鉴于层次分析方法的诸多优点及综合评价的运用越来越普 遍，如何简化层次分析法综合评价计算及检验过程，提高其通用性、 实用性和可操作性成为一个值得研究的问题。由于目前较为流行和权 威的spss统计软件包和sas统计软件包都没有提供层次分析法计算和 检验的专用模块或是过程，手工计算量很大，费时费力而且容易出现 错误，因此笔者对层次分析法的计算和检验过程编写了 sas程序，并 在程序中提供了详细的注释。通过运用sas程序对层次分析法综合评 价实例进行计算和检验,可知整个评价sas运算过程非常简单和快速, 无需建立sas数据集，只需更改sas程序中判断矩阵层次数，再将原 始数据输入程序中的相应数组，运行程序便可很快得到归一化权重系 数及一致性检验结果。之后把最后评价指标与其所属各上级子目标归 一化权重系数连乘得到各评价指标的组合权重系数，再将组合权重系 数输入后续sas程序中的和应数组，运行程序便可得到各评价对象综 合评分指数和排序结果。整个sas程序简明、流畅，可读性强，得出 的计算结果详细、准确，并很好地避免了手工计算过程中因四舍五入 导致的结果偏差甚至错误。综上所述，笔