zoukankan      html  css  js  c++  java
  • 特征提取:PCA主成分分析法和NMF非负矩阵分解法

    1 from sklearn.datasets import load_wine
    2 from sklearn.preprocessing import StandardScaler
    3 wine=load_wine()
    4 X,y=wine.data,wine.target
    5 scaler=StandardScaler()
    6 X_scaled=scaler.fit_transform(X)
    1 from sklearn.decomposition import PCA
    2 pca=PCA(n_components=2)
    3 pca.fit(X_scaled)
    4 X_pca=pca.transform(X_scaled)
    5 X0=X_pca[wine.target==0]
    6 X1=X_pca[wine.target==1]
    7 X2=X_pca[wine.target==2]
    1 import matplotlib.pyplot as plt
    2 plt.scatter(X0[:,0],X0[:,1],c='b',s=60,edgecolors='k')
    3 plt.scatter(X1[:,0],X1[:,1],c='g',s=60,edgecolors='k')
    4 plt.scatter(X2[:,0],X2[:,1],c='r',s=60,edgecolors='k')
    5 plt.legend(wine.target_names,loc='best')
    6 plt.xlabel("component 1")
    7 plt.ylabel("component 2")
    8 plt.show()
    1 plt.matshow(pca.components_,cmap='plasma')
    2 plt.yticks([0,1],['component 1','component 2'])
    3 plt.colorbar()
    4 plt.xticks(range(len(wine.feature_names)),wine.feature_names,rotation=20,ha='left')
    5 plt.show()
     1 from sklearn.decomposition import NMF
     2 from sklearn.preprocessing import MinMaxScaler
     3 mms=MinMaxScaler()
     4 nmf=NMF(n_components=2)
     5 X_scaled=mms.fit_transform(X)
     6 nmf.fit(X_scaled)
     7 X_nmf=nmf.transform(X_scaled)
     8 print(X_nmf.shape)
     9 X0=X_nmf[wine.target==0]
    10 X1=X_nmf[wine.target==1]
    11 X2=X_nmf[wine.target==2]
    1 import matplotlib.pyplot as plt
    2 plt.scatter(X0[:,0],X0[:,1],c='b',s=60,edgecolors='k')
    3 plt.scatter(X1[:,0],X1[:,1],c='g',s=60,edgecolors='k')
    4 plt.scatter(X2[:,0],X2[:,1],c='r',s=60,edgecolors='k')
    5 plt.legend(wine.target_names,loc='best')
    6 plt.xlabel("component 1")
    7 plt.ylabel("component 2")
    8 plt.show()
  • 相关阅读:
    BZOJ3578:GTY的人类基因组计划2(集合hash,STL)
    【BZOJ 1022】 [SHOI2008]小约翰的游戏John
    【BZOJ 1295】 [SCOI2009]最长距离
    【BZOJ 1103】 [POI2007]大都市meg
    【BZOJ 3172】 [Tjoi2013]单词
    【BZOJ 1067】 [SCOI2007]降雨量
    【BZOJ 1491】 [NOI2007]社交网络
    【BZOJ 1087】[SCOI2005]互不侵犯King
    【BZOJ 1009】 [HNOI2008]GT考试
    【BZOJ 1053】[HAOI2007]反素数ant
  • 原文地址:https://www.cnblogs.com/St-Lovaer/p/12269600.html
Copyright © 2011-2022 走看看