(手寫)PCA原理及其Python實現(xiàn)圖文詳解
為什么需要降維呢?
因為數(shù)據(jù)個數(shù) N 和每個數(shù)據(jù)的維度 p 不滿足 N >> p,造成了模型結(jié)果的“過擬合”。有兩種方法解決上述問題:
增加N;減小p。
這里我們講解的 PCA 屬于方法2。
2、樣本均值和樣本方差矩陣





''' -*- coding: utf-8 -*- @ Time : 2021/8/15 22:19 @ Author : Raymond @ Email : wanght2316@163.com @ Editor : Pycharm'''from sklearn.datasets import load_digitsfrom sklearn.decomposition import PCAimport pandas as pdimport matplotlib.pyplot as pltdigits = load_digits()print(digits.keys())print('數(shù)據(jù)的形狀為: {}'.format(digits[’data’].shape))# 構(gòu)建模型 - 降到10 dpca = PCA(n_components=10)pca.fit(digits.data)projected=pca.fit_transform(digits.data)print(’降維后主成分的方差值為:’,pca.explained_variance_)print(’降維后主成分的方差值占總方差的比例為:’,pca.explained_variance_ratio_)print(’降維后最大方差的成分為:’,pca.components_)print(’降維后主成分的個數(shù)為:’,pca.n_components_)print(’original shape:’,digits.data.shape)print(’transformed shape:’,projected.shape)s = pca.explained_variance_c_s = pd.DataFrame({’b’: s,’b_sum’: s.cumsum() / s.sum()})c_s[’b_sum’].plot(style= ’--ko’,figsize= (10, 4))plt.rcParams[’font.sans-serif’] = [’SimHei’] # 指定默認字體plt.rcParams[’axes.unicode_minus’] = False # 解決保存圖像是負號’-’顯示為方塊的問題plt.axhline(0.85, color= ’r’,linestyle= ’--’)plt.text(6, c_s[’b_sum’].iloc[6]-0.08, ’第7個成分累計貢獻率超過85%’, color=’b’)plt.title(’PCA 各成分累計占比’)plt.grid()plt.savefig(’./PCA.jpg’)plt.show()
結(jié)果展示:

本篇文章就到這里了,希望能給你帶來幫助,也希望您能夠多多關(guān)注好吧啦網(wǎng)的更多內(nèi)容!
相關(guān)文章:
1. JS中6個對象數(shù)組去重的方法2. Java commons-httpclient如果實現(xiàn)get及post請求3. 資深程序員:給Python軟件開發(fā)測試的25個忠告!4. 一文帶你徹底理解Java序列化和反序列化5. PHP程序員簡單的開展服務(wù)治理架構(gòu)操作詳解(二)6. PHP利用curl發(fā)送HTTP請求的實例代碼7. Python基于requests庫爬取網(wǎng)站信息8. vscode運行php報錯php?not?found解決辦法9. PHP laravel實現(xiàn)導(dǎo)出PDF功能10. python中文本字符處理的簡單方法記錄

網(wǎng)公網(wǎng)安備