劉 剛
(山東工商學院,山東 煙臺 264005)
近些年來,我國煤礦災害事故頻繁發生,給國家和人民帶來巨大的人身和財產損失。在這種情況下,開展礦山安全評價具有現實意義。在煤礦企業的安全管理中,需要對礦井的安全生產狀況進行分析和評價。針對這一情況,本文采用因子分析,對礦井安全評價結果進行統計分析,增強安全管理的有效性。
因子分析的基本思想,是根據相關性大小把原始變量分組,使得同組內的變量之間相關性較高,而不同組的變量間的相關性則較低。每組變量代表一個基本結構,并用一個不可觀測的綜合變量表示,這個基本結構就稱為公共因子。對于所研究的某一具體問題,原始變量就可以分解成兩部分之和的形式,一部分是少數幾個不可測的所謂的公共因子的線性函數,另一部分是與公共因子無關的特殊因子。少數幾個公共因子中,每一個主要因子就代表經濟變量間相互依賴的一種經濟作用。抓住這些主要因子,就可以幫助我們對復雜的經濟問題進行分析和解釋。
因子分析的一般模型為
式中,X1、X2、…、Xm為實測變量;aij(i=1,2,…,m;j=1,2,…,n)為因子載荷;Fi(i=1,2,…m)為公共因子;εi(i=1,2,…m)為特殊因子。
因子載荷aij是第i個變量在第j個主因子上的載荷,或者說,第i個變量與第j個因子的相關系數。載荷較大,則說明第i個變量與第j個因子的關系越密切;載荷較小,則說明第i個變量與第j個因子的關系越疏遠。因子載荷矩陣中各行數值的平方和,稱為各變量對應的共同度。
公共因子是在各個變量中共同出現的因子,在高維空間中,它們是相互垂直的坐標軸。特殊因子,實際上就是實測變量與估計值之間的殘差值。如果特殊因子為零,則稱為主成分分析。
為了使找到的主因子更易于解釋,往往需要對因子荷載矩陣進行旋轉。旋轉的方法有很多,最常用的是最大方差旋轉法(Varimax)。進行因子旋轉的目的,就是要使因子荷載矩陣中因子荷載的平方值向0和1兩個方向分化,使大的荷載更大,使小的荷載更小。因子旋轉過程中,如果因子對應軸相互正交,則稱為正交旋轉;如果因子對應軸相互間不是正交的,則稱為斜交旋轉。常用的斜交旋轉方法有Promax法等。
將因子表示為變量的線性組合時,所得到的計算結果稱為因子得分,它是對公共因子的估計值。利用它,可以作進一步的分析。
①根據研究問題選取原始變量,并使指標正向化;②對原始數據進行標準化并求其相關矩陣,分析變量之間的相關性;③求解初始公共因子及因子載荷矩陣;④因子旋轉(采用主成分法時可以省略);⑤因子得分;⑥根據因子得分值做進一步分析。
某煤礦有8個礦井,安全管理人員對各礦井的評價結果見表1。

表1 礦井安全狀況原始數據
使用SPSS 17.0軟件,比較各個礦井的安全狀況,評價優劣。
(1)將原始數據標準化(SPSS自動完成),求得相關系數矩陣,見表2。

表2 相關矩陣
從矩陣可以看出,各指標之間具有一定的相關關系。
(2)計算特征值、貢獻率和累計貢獻率。
采用主成分法提取公因子,得到總方差解釋表,如表3所示。

表3 總方差解釋表
從表3可以看出,前兩個特征根都大于1,并且累計貢獻率達到了77%,再結合表4中變量不出現丟失,所以提取因子個數m=2。

表4 成分矩陣
由旋轉后的因子載荷矩陣可以看出:因子1主要反映了與礦井爆炸事故有關的因素;因子2主要反映了礦井非爆炸事故影響因素。
(3)因子得分
SPSS直接給出礦井的單因子得分,再利用綜合因子得分公式:Z=(53.596Z1+23.759Z2)/77.355計算出八個礦井的綜合得分,見表5。

表5 各礦井的綜合得分及排序
在安全管理人員的評價結果基礎上,利用SPSS軟件進行因子分析,并得到最終的排序結果,相對比較客觀,增加了評價結果的可信度。
對于不同的煤礦企業,安全文化、技術水平、組織管理等諸多方面不盡相同,提取的因子個數和意義也會有所不同,并不影響方法的實用性。
[1] 何曉群.多元統計分析[M].北京:中國人民大學出版社,2004: 167-181,66-98.
[2] 張文彤. SPSS11統計分析教程(基礎篇)[M]. 北京:希望電子出版社,2002:56-68.
[3] 任英華,王亞雄. 湖南城市發展的因子分析及政策建議[J].統計教育,2003(2).
[4] 王芳.主成分分析與因子分析的異同比較及應用[J]. 統計教育, 2003(5).
[5] 林海明,林敏子,丁潔花.主成分分析法與因子分析法應用辨析[J]. 數量經濟技術經濟研究,2004,21(9).
[6] 王旭,霍德利.主成分聚類分析法在煤礦安全評價中的應用[J].中國礦業, 2009,18(2).
[7] 曾康生.我國煤礦安全評價方法的應用現狀及其改進[J].中國煤炭, 2008,34(11).