快捷導(dǎo)航

Python Numpy實(shí)現(xiàn)計(jì)算矩陣的均值和標(biāo)準(zhǔn)差詳解

更新時(shí)間：2021年11月22日 09:55:16 作者：葉庭云

NumPy（Numerical Python）是Python的一種開源的數(shù)值計(jì)算擴(kuò)展。這種工具可用來存儲(chǔ)和處理大型矩陣，比Python自身的嵌套列表結(jié)構(gòu)要高效的多。本文主要介紹用NumPy實(shí)現(xiàn)計(jì)算矩陣的均值和標(biāo)準(zhǔn)差，感興趣的小伙伴可以了解一下

一、前言

CRITIC權(quán)重法是一種比熵權(quán)法和標(biāo)準(zhǔn)離差法更好的客觀賦權(quán)法：

它是基于評(píng)價(jià)指標(biāo)的對(duì)比強(qiáng)度和指標(biāo)之間的沖突性來綜合衡量指標(biāo)的客觀權(quán)重?？紤]指標(biāo)變異性大小的同時(shí)兼顧指標(biāo)之間的相關(guān)性，并非數(shù)字越大就說明越重要，完全利用數(shù)據(jù)自身的客觀屬性進(jìn)行科學(xué)評(píng)價(jià)。
對(duì)比強(qiáng)度是指同一個(gè)指標(biāo)各個(gè)評(píng)價(jià)方案之間取值差距的大小，以標(biāo)準(zhǔn)差的形式來表現(xiàn)。標(biāo)準(zhǔn)差越大，說明波動(dòng)越大，即各方案之間的取值差距越大，權(quán)重會(huì)越高；

指標(biāo)之間的沖突性，用相關(guān)系數(shù)進(jìn)行表示，若兩個(gè)指標(biāo)之間具有較強(qiáng)的正相關(guān)，說明其沖突性越小，權(quán)重會(huì)越低。

對(duì)于 CRITIC 權(quán)重法而言，在標(biāo)準(zhǔn)差一定時(shí)，指標(biāo)間沖突性越小，權(quán)重也越小；沖突性越大，權(quán)重也越大；另外，當(dāng)兩個(gè)指標(biāo)間的正相關(guān)程度越大時(shí)，（相關(guān)系數(shù)越接近1），沖突性越小，這表明這兩個(gè)指標(biāo)在評(píng)價(jià)方案的優(yōu)劣上反映的信息有較大的相似性。

在用 Python 復(fù)現(xiàn) CRITIC 權(quán)重法時(shí)，需要計(jì)算變異系數(shù)，以標(biāo)準(zhǔn)差的形式來表現(xiàn)，如下所示：

Sj表示第 j 個(gè)指標(biāo)的標(biāo)準(zhǔn)差，在 CRITIC 權(quán)重法中使用標(biāo)準(zhǔn)差來表示各指標(biāo)的內(nèi)取值的差異波動(dòng)情況，標(biāo)準(zhǔn)差越大表示該指標(biāo)的數(shù)值差異越大，越能放映出更多的信息，該指標(biāo)本身的評(píng)價(jià)強(qiáng)度也就越強(qiáng)，應(yīng)該給該指標(biāo)分配更多的權(quán)重。

研究收集到湖南省某醫(yī)院 2011 年 5 個(gè)科室的數(shù)據(jù)，共有 6 個(gè)指標(biāo)，當(dāng)前希望通過已有數(shù)據(jù)分析各個(gè)指標(biāo)的權(quán)重情況如何，便于醫(yī)院對(duì)各個(gè)指標(biāo)設(shè)立權(quán)重進(jìn)行后續(xù)的綜合評(píng)價(jià)，用于各個(gè)科室的綜合比較等。數(shù)據(jù)如下：

二、詳解計(jì)算均值和標(biāo)準(zhǔn)差

初始化一個(gè)簡單的矩陣：

a = np.array([
    [1, 2, 3],
    [4, 5, 6],
    [7, 8, 9]
    ])
a

分別計(jì)算整體的均值、每一列的均值和每一行的均值：

print("整體的均值：", np.mean(a))              # 整體的均值
print("每一列的均值：", np.mean(a, axis=0))    # 每一列的均值
print("每一行的均值：", np.mean(a, axis=1))    # 每一行的均值

分別計(jì)算整體的標(biāo)準(zhǔn)差、每一列的標(biāo)準(zhǔn)差和每一行的標(biāo)準(zhǔn)差：

print("整體的方差：", np.std(a))              # 整體的標(biāo)準(zhǔn)差
print("每一列的方差：", np.std(a, axis=0))    # 每一列的標(biāo)準(zhǔn)差
print("每一列的方差：", np.std(a, axis=1))    # 每一行的標(biāo)準(zhǔn)差

結(jié)果如下：

三、實(shí)踐：CRITIC權(quán)重法計(jì)算變異系數(shù)

導(dǎo)入需要的依賴庫：

import numpy as np
import pandas as pd

提取數(shù)據(jù)：

df = pd.read_excel("./datas/result03.xlsx")
df

datas = df.iloc[:, 1:]
datas

如下所示：

數(shù)據(jù)正向和逆向化處理：

X = datas.values
xmin = X.min(axis=0)
xmax = X.max(axis=0)
xmaxmin = xmax - xmin
n, m = X.shape
print(m, n)
for i in range(n):
    for j in range(m):
        if j == 5:
            X[i, j] = (xmax[j] - X[i, j]) / xmaxmin[j]   # 越小越好
        else:
            X[i, j] = (X[i, j] - xmin[j]) / xmaxmin[j]   # 越大越好

X = np.round(X, 5)
print(X)

如下所示：