在数模中,确定权重方法有很多种,比如主成分分析法、层次分析法、熵权法、相关系数作为权重等。 网上很多教程都是用spss计算权重,这里主要讲利用python通过主成分分析法确定权重。
主成分分析法概述
主成分分析法是一种线性的降维算法,通过将N维特征进行正交变换,得到相互独立的k维(k<N)数据,通过分析得出“主成分”,利用“主成分”确定影响权重。
数据归一化
为了消除不同变量的量纲的影响,在对数据进行主成分分析前,首先要对数据进行标准化,之所以要对数据进行标准化。 涉及m个样本,n个指标,第j个样本的第i个指标值为,将各个归一化值按如下公式进行归一化为:
from sklearn.preprocessing import MinMaxScaler
mms = MinMaxScaler()
data = mms.fit_transform(data)
pca降维
from sklearn.decomposition import PCA
from sklearn import preprocessing
pca = PCA()
pca.fit(data)
先大致看一下各主成分的累计贡献值,确定k个主成分
pca.components_
pca.explained_variance_
pca.explained_variance_ratio_
pca = PCA(5)
pca.fit(data)
low_d = pca.transform(data)
确定权重
k1_spss = pca.components_ / np.sqrt(pca.explained_variance_.reshape(-1, 1))
j = 0
Weights = []
for j in range(len(k1_spss)):
for i in range(len(pca.explained_variance_)):
Weights_coefficient = np.sum(100 * (pca.explained_variance_ratio_[i]) * (k1_spss[i][j])) / np.sum(
pca.explained_variance_ratio_)
j = j + 1
Weights.append(np.float(Weights_coefficient))
print('Weights',Weights)
权重归一化
Weights=pd.DataFrame(Weights)
Weights1 = preprocessing.MinMaxScaler().fit(Weights)
Weights2 = Weights1.transform(Weights)
print('Weights2',Weights2)
'''
Weights2 [[0.56203305]
[1. ]
[0.63623309]
[0. ]
[0.07484027]]
'''
关注茶谈大数据回复【主成分代码】获取完整的代码
|