一、引言在进行综合评价的过程中权重的确定确实是很重要的,对最终的结果会起着决定性的影响。现在大多数采用的都是主客观结合的组合赋权法。主观赋权用的最多的是AHP,客观赋权用的多的是熵值法、粗糙集的方法,正如楼主所讲的,客观赋权完全依赖于样本数据,当样本数据变化时,权重也会发生变化,从统计规律来讲,随着样本容量的增加,权重的变化应该越来越小,最终趋于一个稳定的值,但在我们实际的评价过程中不可能让样本数达到足够大,因此我们实际还是要把整个评价系统看作是一个不确定性的系统,运用已知的信息来最大限度的挖掘系统的规律,所以我们在有限样本下求出的只能是近似值。主观权重方法简单,但人为因素太强;客观权重又过于依赖样本,这两种方法都存在着信息的损失,采用组合赋权就是最大限度的减少信息的损失,使赋权的结果尽可能的与实际结果接近。现在组合赋权的核心问题在于如何确定两种方法的权重分配,目前这方面的研究比较多,但个人认为大部分方法的数学推导过于繁琐,应用性比较差,基本上不具备可操作性,至于如何准确的
进行组合赋权,现在似乎还没有一套很好的方法。 二、主观赋权:层次分析法AHP2.1 优缺点优点:它完全依靠主观评价做出方案的优劣排序,所需数据量少,决策花费的时间很短。从整体上看,AHP在复杂决策过程中引入定量分析,并充分利用决策者在两两比较中给出的偏好信息进行分析与决策支持,既有效地吸收了定性分析的结果,又发挥了定量分析的优势,从而使决策过程具有很强的条理性和科学性,特别适合在社会经济系统的决策分析中使用。 缺点:用AHP进行决策主观成分很大。当决策者的判断过多地受其主观偏好影响,而产生某种对客观规律的歪曲时,AHP的结果显然就靠不住了。 适用范围:尤其适合于人的定性判断起重要作用的、对决策结果难于直接准确计量的场合。要使AHP的决策结论尽可能符合客观规律,决策者必须对所面临的问题有比较深入和全面的认识。另外,当遇到因素众多,规模较大的评价问题时,该模型容易出现问题,它要求评价者对问题的本质、包含的要素及其相互之间的逻辑关系能掌握得十分透彻,否则评价结果就不可靠和准确。 改进方法: (1) 成对比较矩阵可以采用德尔菲法获得。 (2) 如果评价指标个数过多(一般超过9个),利用层次分析法所得到的权重就有一定的偏差,继而组合评价模型的结果就不再可靠。可以根据评价对象的实际情况和特点,利用一定的方法,将各原始指标分层和归类,使得每层各类中的指标数少于9个。 2.2 基本步骤步骤一、建立递阶层次结构模型:目的层,准则层,方案层
/ t4 C6 y6 O' D
步骤二、构造出各层次中的所有判断矩阵;判断矩阵为两两指标重要性的判断,一般判别矩阵构造如下。
- y: `! k3 H: ~3 T# j k2 N/ R/ n6 Z( `
, [7 z9 |6 k5 p5 l6 t: p% f) d
步骤三:层次单排序及一致性检验;
$ @/ j) _: t. ^% f ^
& O6 v9 U% U% ^) {( A步骤四:层次总排序及一致性检验整套算法实际上是用了两次重要性权重。 准则层,从准则的重要性矩阵(nxn矩阵)中,抽取重要性权重。它的现实意义是 每个准则的重要程度,也就是步骤三。 方案层,对每个准则,m个方案都有个mxm矩阵(总共是n个mxm矩阵)。也就是说,对每个准则,都可以算出m个方案的重要性权重。 然后n个重要性权重组合起来,与准则层的重要性权重相乘。就得到了每个方案的重要性权重,也就是层次总排序。 2.3 程序实现
8 `# G6 x2 p4 M2 eclass AHP:4 m* Z8 l! h7 o" T3 x' q' z
def __init__(self, criteria, b):5 Y b( a0 L+ J: {( N' W& _5 Z
self.RI = (0, 0, 0.58, 0.9, 1.12, 1.24, 1.32, 1.41, 1.45, 1.49)
- O, _ A1 f: c self.criteria = criteria
% r: Y& W6 M& s( Z I self.b = b
# T6 b" b4 ^& h. T* E" {! X% y" [ self.num_criteria = criteria.shape[0]
2 E, R/ f- b, r self.num_project = b[0].shape[0]) m$ ~/ `$ B1 v9 j; T& ?8 x
/ v& c$ N9 |3 K- d5 F5 ^- k" F* a
def cal_weights(self, input_matrix):
( b5 z Z X; x9 u, ^) [% s0 G% q input_matrix = np.array(input_matrix)$ G! ]( e* p, H9 ^$ o5 I! [& O
n, n1 = input_matrix.shape
% \0 L! i6 U9 `2 r6 u8 V0 Q5 a assert n == n1, '不是一个方阵'7 M. d% `# ~" q3 b
for i in range(n):
6 a) c, v( D1 c# F for j in range(n):7 X$ M) O; s# v. v' q# q3 _ I
if np.abs(input_matrix[i, j] * input_matrix[j, i] - 1) > 1e-7:5 R1 u( l8 ]3 u8 P9 \7 O6 G
raise ValueError('不是反互对称矩阵')
4 d7 g2 `- {* H* s% K3 c: \! B+ P
eigenvalues, eigenvectors = np.linalg.eig(input_matrix)
4 n$ ^$ U0 R( l* Y2 l; P5 l2 b+ M/ ?9 X/ W4 }- \0 s2 @
max_idx = np.argmax(eigenvalues)
1 z" k) o& w! P( E% B max_eigen = eigenvalues[max_idx].real
3 ]- l' N& b# X/ u0 M3 a) g5 l; g0 H eigen = eigenvectors[:, max_idx].real$ R' I' `. d, A6 F+ z8 ^) N
eigen = eigen / eigen.sum()
9 U* k6 C6 `( \% G) X/ _$ l" R( Q" h& I2 C, K7 b V
if n > 9:
; x: X7 ]5 U6 k v5 ?$ | CR = None0 L; R& f" W* r" |7 o; U
warnings.warn('无法判断一致性')
4 a) y# }+ @5 e! o else:
; c0 G0 Z6 m1 W CI = (max_eigen - n) / (n - 1)
2 i0 \* j1 u) O) V L/ x% f CR = CI / self.RI[n]3 \" k; H: d- I' C4 T
return max_eigen, CR, eigen
, I( @- ?8 @5 _! |
0 v/ V* K" i+ t# Y def run(self):( n& k0 m! l- m) M/ b* J
max_eigen, CR, criteria_eigen = self.cal_weights(self.criteria)( r8 B1 b7 P0 C8 p, {: g
print('准则层:最大特征值{:<5f},CR={:<5f},检验{}通过'.format(max_eigen, CR, '' if CR < 0.1 else '不')), K' P9 A: S! T d) }9 q
print('准则层权重={}\n'.format(criteria_eigen))2 Y+ N6 t$ L: Z* t5 V& s6 j
/ q1 e, z, s" I. ^+ Q, ? max_eigen_list, CR_list, eigen_list = [], [], []* p, J: r7 o- |2 e
for i in self.b:
! i# c/ y; ]/ h' ? max_eigen, CR, eigen = self.cal_weights(i)
) B; q- y( n0 i v: K. } max_eigen_list.append(max_eigen)2 I; o, _# j# @: B# D4 N" q
CR_list.append(CR)
" i: P9 z: P. w8 j' G+ X9 _ eigen_list.append(eigen) l3 z% M' u# u; X) F
8 q0 w5 F8 E; S0 q0 O
pd_print = pd.DataFrame(eigen_list,: e% L! J6 e8 D7 m$ J
index=['准则' + str(i) for i in range(self.num_criteria)],
% D9 w0 |! A t) r columns=['方案' + str(i) for i in range(self.num_project)],# K3 Q1 ^, N* B, |2 \# I; f) k
)- N) b9 p' b4 F" |5 d
pd_print.loc[:, '最大特征值'] = max_eigen_list+ s) B- f# S* o% Q0 q8 J
pd_print.loc[:, 'CR'] = CR_list8 d: K* u; J; Q6 B
pd_print.loc[:, '一致性检验'] = pd_print.loc[:, 'CR'] < 0.1! h) p2 |) }; f6 K; t' H% T& {/ n
print('方案层')% A! }: _# V Y4 Z8 J8 d
print(pd_print)+ t% ]8 a, i( a6 Q& ^8 W6 U" E
& s' ]2 U( {& v" B& h6 ` # 目标层
% `% j) {9 c$ A3 T6 _ ^9 _3 Q obj = np.dot(criteria_eigen.reshape(1, -1), np.array(eigen_list))( ]) J# r2 x0 p. g- [
print('\n目标层', obj)
' N# h, c& H/ y. C! F print('最优选择是方案{}'.format(np.argmax(obj)))
& U; z N) ]9 ~# c5 d+ A2 `' _ return obj
u0 N* |0 B, Y5 O) b# Y; q0 `* d4 I$ `; U0 ? @
二、客观赋权因子分析、熵权法等,不再赘述,见yys , ? T) L* ?! ~5 I" ~6 E$ Q" b+ \ u/ d
三、组合赋权常见的主要有:方差最大化赋权法、最佳协调赋权法、组合目标规划法,组合最小二乘法。
- D9 V6 l C) P9 c* Y" @
# v" a! P' r0 K
. u6 z/ x1 N; V" g1 _) K3 I |