主成分分析與因子分析及SPSS實現 - 程式前沿

主成分分析與因子分析及SPSS實現一、主成分分析（1）問題提出在問題研究中，為了不遺漏和準確起見，往往會面面俱到，取得大量的指標來進行分析。

程式語言前端開發IOS開發Android開發雲端運算人工智慧伺服器搜尋資料庫軟體開發工具主成分分析與因子分析及SPSS實現2018.08.03程式語言SPSS主成分分析,svm主成分分析,主成分分析,主成分分析中心化,主成分分析怎麼看,主成分分析視覺化HOME程式語言主成分分析與因子分析及SPSS實現Advertisement主成分分析與因子分析及SPSS實現一、主成分分析（1）問題提出在問題研究中，為了不遺漏和準確起見，往往會面面俱到，取得大量的指標來進行分析。

比如為了研究某種疾病的影響因素，我們可能會收集患者的人口學資料、病史、體徵、化驗檢查等等數十項指標。

如果將這些指標直接納入多元統計分析，不僅會使模型變得複雜不穩定，而且還有可能因為變數之間的多重共線性引起較大的誤差。

有沒有一種辦法能對資訊進行濃縮，減少變數的個數，同時消除多重共線性？這時，主成分分析隆重登場。

（2）主成分分析的原理主成分分析的本質是座標的旋轉變換，將原始的n個變數進行重新的線性組合，生成n個新的變數，他們之間互不相關，稱為n個“成分”。

同時按照方差最大化的原則，保證第一個成分的方差最大，然後依次遞減。

這n個成分是按照方差從大到小排列的，其中前m個成分可能就包含了原始變數的大部分方差（及變異資訊）。

那麼這m個成分就成為原始變數的“主成分”，他們包含了原始變數的大部分資訊。

注意得到的主成分不是原始變數篩選後的剩餘變數，而是原始變數經過重新組合後的“綜合變數”。

我們以最簡單的二維資料來直觀的解釋主成分分析的原理。

假設現在有兩個變數X1、X2，在座標上畫出散點圖如下： 可見，他們之間存在相關關係，如果我們將座標軸整體逆時針旋轉45°，變成新的座標系Y1、Y2，如下圖：根據座標變化的原理，我們可以算出：Y1=sqrt(2)/2*X1sqrt(2)/2*X2Y2=sqrt(2)/2*X1–sqrt(2)/2*X2其中sqrt(x)為x的平方根。

通過對X1、X2的重新進行線性組合，得到了兩個新的變數Y1、Y2。

此時，Y1、Y2變得不再相關，而且Y1方向變異（方差）較大，Y2方向的變異（方差）較小，這時我們可以提取Y1作為X1、X2的主成分，參與後續的統計分析，因為它攜帶了原始變數的大部分資訊。

至此我們解決了兩個問題：降維和消除共線性。

對於二維以上的資料，就不能用上面的幾何圖形直觀的表示了，只能通過矩陣變換求解，但是本質思想是一樣的。

二、因子分析（一）原理和方法：因子分析是主成分分析的擴充套件。

在主成分分析過程中，新變數是原始變數的線性組合，即將多個原始變數經過線性（座標）變換得到新的變數。

因子分析中，是對原始變數間的內在相關結構進行分組，相關性強的分在一組，組間相關性較弱，這樣各組變數代表一個基本要素（公共因子）。

通過原始變數之間的複雜關係對原始變數進行分解，得到公共因子和特殊因子。

將原始變數表示成公共因子的線性組合。

其中公共因子是所有原始變數中所共同具有的特徵，而特殊因子則是原始變數所特有的部分。

因子分析強調對新變數（因子）的實際意義的解釋。

舉個例子：比如在市場調查中我們收集了食品的五項指標（x1-x5）:味道、價格、風味、是否快餐、能量，經過因子分析，我們發現了：x1=0.02*z10.99*z2e1x2=0.94*z1–0.01*z2e2x3=0.13*z10.98*z2e3x4=0.84*z10.42*z2e4x5=0.97*z1–0.02*z2e1（以上的數字代表實際為變數間的相關係數，值越大，相關性越大）第一個公因子z1主要與價格、是否快餐、能量有關，代表“價格與營養”第二個公因子z2主要與味道、風味有關，代表“口味”e1-5是特殊因子，是公因子中無法解釋的，在分析中一般略去。

同時，我們也可以將公因子z1、z2表示成原始變數的線性組合，用於後續分析。

（二）使用條件：（1）樣本量足夠大。

通常要求樣本量是變數數目的5倍以上，且大於100例。

（2）原始變數之間具有相關性。

如果變數之間彼此獨立，無法使用因子分析。

在SPSS中可用KMO檢驗和Bartlett球形檢驗來判斷。

（3）生成的公因子要有實際的意義，必要時可通過因子旋轉（座標變化）來達到。

三、主成分分析和因子分析的聯絡與區別聯絡：兩者都是降維和資訊濃縮的方法。

生成的新變數均代表了原始變數的大部分資訊且互相獨立，都可以用於後續的迴歸分析、判別分析、聚類分析等等。

區別：（1）主成分分析是按照方差最大化的方法生成的新變數，強調新變數貢獻了多大比例的方差，不關心新變數是否有明確的實際意義。

（2）因子分析著重要求新變數具有實際的意義，能解釋原始變數間的內在結構。

SPSS沒有提供單獨的主成分分析方法，而是混在因子分析當中，下面通過一個例子來討論主成分分析與因子分析的實現方法及相關問題。

 一、問題提出 男子十項全能比賽包含100米跑、跳遠、跳高、撐杆跳、鉛球、鐵餅、標槍、400米跑、1500米跑、110米跨欄十個專案，總分為各個專案得分之和。

為了分析十項全能主要考察哪些方面的能力，以便有針對性的進行訓練，研究者收集了134個頂級運動員的十項全能成績單，將通過因子分析來達到分析目的。

 二、分析過程 變數檢視： 資料檢視（部分）：選單選擇（分析->降維->因子分析）：開啟因子分析的主介面，將十項成績選入”變數“框中（不要包含總分），如下：點選”描述“按鈕，開啟對話方塊，選中”係數“和”KMO和Bartlett球形度檢驗“：上圖相關解釋：”係數“：為變數之間的相關係數陣列，可以直觀的分析相關性。

”KMO和Bartlett球形度檢驗“：用於定量的檢驗變數之間是否具有相關性。

點選”繼續“，回到主介面，點選”抽取“，開啟對話方塊。

”方法“=>”主成分“，”輸出“=>”未旋轉的因子解“和”碎石圖“，”抽取“=>”基於特徵值“，其餘選擇預設。

解釋：①因子抽取的方法：選取預設的主成分法即可，其餘方法的計算結果可能有所差異。

②輸出：”未旋轉的因子解”極為主成分分析結果。

碎石圖有助於我們判斷因子的重要性（詳細介紹見後面）。

③抽取：為抽取主成分（因子）的方法，一般是基於特徵值大於1，預設即可。

點選”繼續“，回到主介面，點選”確定“，進入分析。

輸出的主要表格如下：（1）相關性檢驗因子分析要求變數之間有相關性，所以首先要進行相關性檢驗。

首先輸出的是變數之間的相關係數矩陣：可以直觀的看到，變數之間有相關性。

但需要檢驗，接著輸出的是相關性檢驗：上圖有兩個指標：第一個是KMO值，一般大於0.7就說明不了之間有相關性了。

第二個是Bartlett球形度檢驗，P值<0.001。

綜合兩個指標，說明變數之間存在相關性，可以進行因子分析。

否則，不能進行因子分析。

（2）提取主成分和公因子接下來輸出主成分結果：這就是主成分分析的結果，表中第一列為10個成分；第二列為對應的”特徵值“，表示所解釋的方差的大小；第三列為對應的成分所包含的方差佔總方差的百分比；第四列為累計的百分比。

一般來說，選擇”特徵值“大於1的成分作為主成分，這也是SPSS預設的選擇。

在本例中，成分1和2的特徵值大於1，他們合計能解釋71.034%的方差，還算不錯。

所以我們可以提取1和2作為主成分，抓住了主要矛盾，其餘成分包含的資訊較少，故棄去。

下面，輸出碎石圖，如下：碎石圖來源於地質學的概念。

在岩層斜坡下方往往有很多小的碎石，其地質學意義不大。

碎石圖以特徵值為縱軸，成分為橫軸。

前面陡峭的部分特徵值大，包含的資訊多，後面平坦的部分特徵值小，包含的資訊也小。

由圖直觀的看出，成分1和2包含了大部分資訊，從3開始就進入平臺了。

接下來，輸出提取的成分矩陣：上表中的數值為公因子與原始變數之間的相關係數，絕對值越大，說明關係越密切。

公因子1和9個運動專案都正相關（注意跑步運動運動的計分方式，時間越短，分數越高），看來只能稱為“綜合運動”因子了。

公因子2與鐵餅、鉛球正相關，與1500米跑、400米跑負相關，這究竟代表什麼意思呢？看來只能成為“不知所云”因子了。

（三）因子旋轉前面提取的兩個公因子一個是大而全的“綜合因子”，一個不知所云，得到這樣的結果，無疑是分析的失敗。

不過，不要灰心，我們可以通過因子的旋轉來獲得更好的解釋。

在主介面中點選“旋轉”按鈕，開啟對話方塊，“方法”=>“最大方差法”，“輸出”=>“旋轉解”。

點選“繼續”，回到主介面點選“確認”進行分析。

輸出結果如下：這是選擇後的成分矩陣。

經過旋轉，可以看出：公因子1得分越高，所有的跑步和跨欄成績越差，而跳遠、撐杆跳等需要助跑類專案的成績也越差，所以公因子1代表的是奔跑能力的反向指標，可稱為“奔跑能力”。

公因子2與鐵餅和鉛球的正相關性很高，與標槍、撐杆跳等需要上肢力量的專案也正相關，所以該因子可以成為“上肢力量”。

經過旋轉，可以看出公因子有了更合理的解釋。

（四）結果的儲存在最後，我們還要將公因子儲存下來供後續使用。

點選“得分”按鈕，開啟對話方塊，選中“儲存為變數”，方法採用預設的“迴歸”方法，同時選中“顯示因子得分系數矩陣”。

SPSS會自動生成2個新變數，分別為公因子的取值，放在資料的最後。

同時會輸出一個因子係數表格：由上圖，我們可以寫出公因子的表示式（用F1、F2代表兩個公因子，Z1~Z10分別代表原始變數）：F1=-0.16*Z10.161*Z20.145*Z30.199*Z4-0.131*Z5-0.167*Z60.137*Z70.174*Z80.131*Z9-0.037*Z10F2同理，略去。

注意，這裡的變數Z1~Z10，F1、F2不再是原始變數，而是標準正態變換後的變數。

Advertisement写评论取消回覆很抱歉，必須登入網站才能發佈留言。

近期文章Spark入門（一）用SparkShell初嘗Spark滋味2019.12.08Spark入門（二）如何用Idea運行我們的Spark項目2019.12.08Spark入門（三）Spark經典的單詞統計2019.12.08Spark入門（四）Spark的map、flatMap、mapToPair2019.12.08Spark入門（五）Spark的reduce和reduceByKey2019.12.08Spark入門（六）Spark的combineByKey、sortBykey2019.12.08Spark入門（七）Spark的intersection、subtract、union和distinct2019.12.08Spark實戰尋找5億次訪問中，訪問次數最多的人2019.12.08Spark實戰搭建我們的Spark分佈式架構2019.12.08【機器學習】深度學習開發環境搭建2019.12.08AdvertisementAdvertisement