因子分析论文.doc

1、关于2000年美国各州犯罪数量的因子分析【摘要】本文采用因子分析方法对美国各州的犯罪情况州际差异及结构差异进行了探讨,研究表明人口与犯罪数量有密切相关，其主要为因贪图财色而进行犯罪，人口越多，犯罪数量越多，尤其是因财物而引起的抢劫、入室行窃、盗窃、自盗数量越多，而大部分的攻击性犯罪发生在美国南部，可能是由于历史原因。 【关键字】美国 各州 犯罪数量 因子提取 正交旋转 因子分析 因子得分【正文】一、 因子分析原理因子分析是一种将多变量化简的多元统计方法，它可以看作是主成份分析的推广。因子分析的目的是分解原始变量，从中归纳出潜在的“类别”，相关性较强的变量归为一类，不同类间的变量的相关性则较低。

2、每类变量代表了一个“共同因子”，即一种内在结构（联系）。因子分析就是寻找这种内在结构（联系）的方法。从全部计算过程来看作R型因子分析与作Q型因子分析都是一样的，只不过出发点不同，R型从相关系数矩阵出发，Q型从相似系数阵出发都是对同一批观测数据，可以根据其所要求的目的决定用哪一类型的因子分析。（一）模型主要模型形式：(2)矩阵型式（二）相关概念解释1、因子载荷aij称为因子载荷（实际上是权数）。因子载荷的统计意义：就是第i个变量与第j个公共因子的相关系数，即表示变量xi依赖于Fj的份量（比重），心理学家将它称为载荷。2、变量共同度二、 主要计算方法及步骤（一）方法说明1、因子载荷矩阵估计方法因子

3、载荷的求解方法主要有主成分法，主轴因子旋转法和极大似然法。主成分法指在进行因子分析之前先对数据进行主成分分析，把前几个主成分作为未旋转的公因子，但是此种方法得到的特殊因子间并不相互独立，当变量的共同度较大时，特殊因子所起的作用较小，它们之间的相关性可以忽略。主轴因子法与主成分分析方法类似，都是都分析矩阵的结构入手，主轴因子法的不同之处在于，其假定m个公因子只能解释原始变量的部分方差，利用变量共同度来代替相关矩阵中对角元素1，并以新矩阵为出发点求解特征值和特征向量。极大似然估计法假定公因子与特殊因子服从正态分布，通过构造似然函数求因子载荷和特殊因子方差的极大似然估计。2、因子旋转因子分析的目的不

4、仅是找出主因子，更重要的是知道每个主因子的意义。主因子的意义是根据主因子与可观测变量Xi的关系来确定的。因此希望主因子Fj对Xi（i=1,2,p）的载荷平方,有的值很大,有的值很小,(向0,1两极分化)，因子载荷矩阵的这种特征称“因子简单结构”。 但是用上述方法所求出的主因子解，初始因子载荷矩阵并不满足“简单结构准则”，各因子的典型代表变量不很突出，因而容易使因子的意义含糊不清，不便于对因子进行解释。为此须对因子载荷矩阵施行旋转，因子轴方差最大正交旋转的目的即使因子载荷矩阵成为“简单结构”的因子载荷矩阵。使得因子载荷的平方按列向0和1两极转化，较大的载荷值只集中在少数X变量上,达到其结构简化的

5、目的。易于因子命名。经过旋转后，主因子对Xi的方差贡献（变量共同度）并不改变，但各主因子的方差贡献可能有较大的改变，即不再与原来相同，因此，可以通过适当的旋转求得令人满意的主因子。为了对公因子F能够更好的解释，可通过因子旋转的方法得到一个好解释的公因子。所谓对公因子更好解释，就是使每个变量仅在一个公因子上有较大的载荷，而在其余的公因子上的载荷比较小。这种变换因子载荷的方法称为因子轴的旋转。因子旋转的方法很多，常用的为方差最大正交旋转。3、因子得分在分析中，人们往往更愿意用公共因子反映原始变量，这样更有利于描述研究对象的特征。因而往往将公共因子表示为变量（或样品）的线性组合，即：称上式为因子得分

6、函数，用它可计算每个样品的公因子得分。估计因子得分的方法很多。（一）、背景介绍犯罪是危害他人、社会或者国家的行爲，订立法律在于保护人民所享有的法律权益，凡是违反法律，破坏他人受到法律所保障的权益行为，都是犯罪行为。本文将以美国2000年各类犯罪数据为样本，研究美国各州各类犯罪的差异性和相似性，收集到2000年美国50个州7类犯罪（谋杀、强奸、抢劫、攻击、入室行窃、盗窃、自盗）的犯罪数量数据(数据来源：美国司法统计局Bureau of Justice Statistics，BJS)。由于涉及变量较多，直接进行地区间的比较分析较为繁琐，因此首先考虑因子分析方法减少变量个数，之后在进行比较和综合分析

7、。（二）、分析过程详解1、 数据标准化2、 考察原有变量是否适合进行因子分析利用SPSS软件，correlate功能计算相关系数矩阵，计算皮尔森相关系数并进行卡方双尾检验，可以看出变量间存在这很大的相关性。表2-1相关性谋杀强奸抢劫攻击入室行窃盗窃自盗谋杀Pearson 相关性1.520*.294*.813*.167.072.110显著性（双侧）.000.038.000.245.618.448N50505050505050强奸Pearson 相关性.520*1.478*.696*.679*.585*.441*显著性（双侧）.000.000.000.000.000.001N50505050505

8、050 抢劫Pearson 相关性.294*.478*1.545*.549*.398*.572*显著性（双侧）.038.000.000.000.004.000N50505050505050攻击Pearson 相关性.813*.696*.545*1.468*.318*.330*显著性（双侧）.000.000.000.001.024.019N50505050505050入室行窃Pearson 相关性.167.679*.549*.468*1.748*.678*显著性（双侧）.245.000.000.001.000.000N50505050505050盗窃Pearson 相关性.072.585*.39

9、8*.318*.748*1.529*显著性（双侧）.618.000.004.024.000.000N50505050505050自盗Pearson 相关性.110.441*.572*.330*.678*.529*1显著性（双侧）.448.001.000.019.000.000N50505050505050*. 在 .01 水平（双侧）上显著相关。*. 在 0.05 水平（双侧）上显著相关。表2-1是原有变量的相关系数矩阵。可以看到：大部分的相关系数都较高，除了少数变量间的相关性较低外，其他变量的相关性都大于0.3，各变量呈交较强的线性关系，能够从中提取公因子，适合进行分析。表2-2 KMO 和

10、 Bartlett 的检验取样足够度的 Kaiser-Meyer-Olkin 度量。.780Bartlett 的球形度检验近似卡方218.365df21Sig.000KMO值：0.9以上非常好；0.8以上好；0.7一般；0.6差；0.5很差；0.5以下不能接受。通过观察上面的计算结果，可以知道，巴特利球体检验的统计量的观测值为218.365，相应的概率P接近0，则应拒绝零假设，认为相关系数阵与单位阵有显著的差异。同时KOM值为0.780，接近比较好的范围内，根据Kaiser给出的KOM度量标准可知原有变量适合进行因子分析。3、提取因子表3-1公因子方差初始提取Zscore(谋杀)1.000.9

11、10Zscore(强奸)1.000.748Zscore(抢劫)1.000.548Zscore(攻击)1.000.919Zscore(入室行窃)1.000.850Zscore(盗窃)1.000.719Zscore(自盗)1.000.680提取方法：主成份分析。表3-1是因子分析的初始解，显示了所有变量的共同度数据。第一列是因子分析初始解下的变量共同度，它表明：对原有10个变量如果采用主成分分析法提取所有特征根，那么原有变量的所有方差都可被解释，变量的共同度都为1（原有变量标准化后的方差为1）。事实上，因子个数小于原有变量个数才是因子分析的目标，所以不可以提取全部特征根，第二列是在指定提取条件（特

12、征根大于1）提取特征根时的共同度。可以看到，谋杀、强奸、攻击、入室行窃、盗窃、自盗等变量的绝大部分信息可被因子解释，因子提取的总体效果较好。表3-2解释的总方差成份初始特征值提取平方和载入旋转平方和载入合计方差的 %累积 %合计方差的 %累积 %合计方差的 %累积 %13.91455.90955.9093.91455.90955.9093.04843.54943.54921.45920.83976.7481.45920.83976.7482.32433.19876.7483.6799.70486.4514.3785.40691.8575.2623.74295.5996.1932.76098.3

13、597.1151.641100.000提取方法：主成份分析。提取2个公因子，从方差贡献率可以看出，其中第一个公因子解释了总体方差的43.549%，第二个公因子解释了总体方差的33.198，两个公因子的累计方差贡献率为76.748%,可以较好的解释总体方差。图3-1 通过此图可以明显看出前三个因子可以解释大部分的方差，到第四个因子以后，线逐渐平缓，解释能力不强。但提取三个因子后，因子含义不如提取两因子明晰，因此我们提取两个公因子。表3-3成份矩阵a成份12Zscore(谋杀).539.787Zscore(强奸).853.140Zscore(抢劫).738-.055Zscore(攻击).787.5

14、48Zscore(入室行窃).844-.370Zscore(盗窃).717-.452Zscore(自盗).710-.419提取方法 :主成份。a. 已提取了 2 个成份。通过因子载荷矩阵可以看出因子的意义并不是十分明确，为了对因子进行解释与说明，进行因子旋转，选取方差最大因子旋转方法，并保留因子得分。4、 因子旋转表4-1旋转成份矩阵a成份12Zscore(谋杀)-.033.953Zscore(强奸).604.619Zscore(抢劫).627.394Zscore(攻击).308.908Zscore(入室行窃).899.204Zscore(盗窃).845.062Zscore(自盗).820.0

15、84提取方法 :主成份。 旋转法 :具有 Kaiser 标准化的正交旋转法。a. 旋转在 3 次迭代后收敛。通过上表中旋转的结果，我们可以看出第一个公因子在强奸、抢劫、入室行窃、盗窃、自盗上有较大的载荷，主要由于求财及求色，可以将其命名为财色欲望因子；第二个因子在谋杀、强奸、攻击上有较大的载荷，可以理解为攻击欲望因子； 最后，通过上面的因子我们可以通过两个主要的因子来表示七个犯罪变量，分别为财色欲望因子和攻击欲望因子。接下来计算各因子得分。图4-15、 因子得分表5-1成份得分系数矩阵a成份12Zscore(谋杀)-.209.516Zscore(强奸).119.206Zscore(抢劫).17

16、4.082Zscore(攻击)-.061.422Zscore(入室行窃).324-.076Zscore(盗窃).331-.141Zscore(自盗).317-.124提取方法 :主成份。 旋转法 :具有 Kaiser 标准化的正交旋转法。 构成得分。a. 系数已被标准化。表5-1中，得出因子得分系数矩阵，根据因子得分系数和原始变量的标准化值，可以观测各犯罪类型的各因子的得分数，旋转后的因子表达式可以写成：F1=-0.209谋杀+0.119强奸+0.174抢劫-0.061攻击+0.324入室行窃+0.331盗窃+0.317自盗F1=0.516谋杀+0.206强奸+0.082抢劫+0.442攻击-

17、0.076入室行窃-0.141盗窃-0.124自盗因子表达式中各变量均为标准化后的变量。表5-2成份得分协方差矩阵成份1211.000.0002.0001.000提取方法 :主成份。 旋转法 :具有 Kaiser 标准化的正交旋转法。 构成得分。 表5-2为因子间的相关矩阵，可以看出旋转后F1和F2是完全不相关的，这也是由于正交旋转后因子仍然正交。6、 各州的综合评价表6-1排名财色欲望因子得分攻击欲望因子综合得分人口1加利福尼亚州 北卡罗来纳州 加利福尼亚州 加利福尼亚州 2内华达州 佛罗里达州 密歇根州 纽约州 3密歇根州 路易斯安那州 佛罗里达州 得克萨斯州 4科罗拉多州 南卡罗来纳州

18、马里兰州 宾夕法尼亚州 5纽约州 阿拉巴马州 阿利桑那州 伊利诺斯州 6夏威夷州 乔治亚州 纽约州 佛罗里达州 7阿利桑那州 马里兰州 科罗拉多州 俄亥俄州 8华盛顿州 得克萨斯州 内华达州 密歇根州 由上表可知，加利福尼亚州人口最多，而其在财色欲望因子上得分最高，在综合得分中也是最高，其犯罪数量最多。在综合得分中得分较高的州其人口较多。而在攻击欲望因子中得分较高的八个州均是位于美国南部的州，这可能与黑人和白人的矛盾，以及南部各州的经济有关系，也可能与其历史有一定的原因。可知，人口与犯罪数量有密切相关，其主要为因贪图财色而进行犯罪，人口越多，犯罪数量越多，尤其是因财物而引起的抢劫、入室行窃、盗

19、窃、自盗数量越多，而大部分的攻击性犯罪发生在美国南部，可能是由于历史原因。参考文献：1于秀林，任雪松 多元统计分析中国统计出版社 2SPSS统计分析从基础到实践（第2版）:北京:电子工业出版社，20093何晓群.多元统计分析M.中国人民大学出版社，20044茆诗松,周纪芗.概率论与数理统计M.高等统计出版社，19995高惠璇应用多元统计分析M.北京大学出版社20056薛微 基于SPSS的数据分析中国人民大学出版社原始数据：该数据集包括了7个变量的50个观测值度量，它报告了美国2000年50个州的7种罪行的犯罪数量。其中，地区一栏中，1代表东北部，2代表中西部，3代表南部，4代表西部。州名人口谋

20、杀强奸抢劫攻击入室行窃盗窃自盗地区缅因州1,1641.57.012.6625621,0551461新罕布什尔州9982.06.012.1365669291721佛蒙特州5351.310.37.6557319691241马萨诸塞州5,8223.512.099.5881,1341,5318781罗得岛州9683.23.678.31201,0192,1868591康涅狄格州3,1743.59.170.4871,0841,7514841纽约州17,7837.915.5443.32091,4142,0256821新泽西州7,5625.712.9169.4901,0411,6895571宾夕法尼亚州11,

21、8535.311.3106.090594113401俄亥俄州10,7446.616.0145.91168541,9444932印第安纳州5,4994.817.9107.5958601,7914292伊利诺斯州11,5359.620.4251.11877652,0285182密歇根州9,0889.427.1346.61931,5712,8974642威斯康辛州4,7752.06.733.1445391,8602182明尼苏达州4,1932.09.789.1518021,9023462爱荷华州2,8841.96.228.6485071,7431752密苏里州5,02910.727.42.81671

22、,1872,0745382北达科他州6850.56.26.5212861,295912南达科他州7083.811.117.1604711,396942内布拉斯加州1,6063.09.357.31155051,5722922堪萨斯州2,4504.814.575.11088822,3022572特拉华州6227.718.6105.51961,0562,3205593马里兰州4,3929.223.9388.62531,0512,4175483弗吉尼亚州5,7068.415.492.01438061,9802973西弗吉尼亚州1,9366.26.727.384389774923北卡罗来纳州6,2551

23、1.812.953.02937661,3381693南卡罗来纳州3,34714.618.160.11931,0251,5092563乔治亚州5,97615.310.195.8177901,8693093佛罗里达州11,36612.722.2186.12771,5622,8613973肯塔基州3,72611.113.772.81237041,2123463田纳西州4,7628.815.582.01698071,0252893阿拉巴马州4,02111.718.550.32157631,1252233密西西比州2,61311.58.919.0140351694783阿肯色州2,35910.117.1

24、45.61508851,2111093路易斯安那州4,48111.723.1140.82388901,6283853俄克拉荷马州3,3015.915.654.91278411,6612803得克萨斯州16,37011.621.0134.11951,1512,1833943蒙大拿州8263.210.522.3755941,9562224爱达荷州154.612.320.5866742,2141444怀俄明州5095.712.322.0736462,0491654科罗拉多州3,2316.236.0129.11851,3812,9925884新墨西哥州1,4509.421.766.11961,1422

25、,4083924阿利桑那州3,1879.527.0120.22141,4933,5505014犹他州1,6453.410.953.1709152,8333164内华达州9368.819.6188.41821,6613,0446614华盛顿州4,4093.518.093.51061,4412,8533624俄勒冈州2,6874.618.0102.51321,2732,8253334加利福尼亚州26,3656.935.1206.92261,7533,4226894阿拉斯加州52112.226.171.81687902,1835514夏威夷州1,0543.611.863.3431,4563,1065

26、814标准化数据：州名Z谋杀Z强奸Z抢劫Z攻击Z入室行窃Z盗窃Z自盗缅因州-1.39242-1.17253-0.91416-1.07702-0.92758-1.15943-1.11147新罕布什尔州-1.26248-1.30862-0.91948-1.45842-0.91706-1.32757-0.98122佛蒙特州-1.44439-0.72344-0.96734-1.17970-0.48300-1.27420-1.22169马萨诸塞州-0.87266-0.492090.01010-0.695620.57717-0.524232.55569罗得岛州-0.95063-1.63523-0.2153

27、8-0.226200.274640.349842.46050康涅狄格州-0.87266-0.88675-0.29940-0.710290.44564-0.230650.58184纽约州0.27079-0.015793.666761.079371.313770.134991.57377新泽西州-0.30094-0.369610.75356-0.666280.33252-0.313380.94755宾夕法尼亚州-0.40489-0.587350.07924-0.66628-0.84340-2.55261-0.13957俄亥俄州-0.067050.052260.50361-0.28488-0.159

28、420.026900.62692印第安纳州-0.534830.310820.09519-0.59293-0.14364-0.177270.30630伊利诺斯州0.712580.651041.622520.75664-0.393550.139000.75217密歇根州0.660611.562832.638260.844661.726791.298640.48164威斯康辛州-1.26248-1.21336-0.69613-1.34107-0.98809-0.08519-0.75077明尼苏达州-1.26248-0.80509-0.10051-1.23838-0.29622-0.02914-0.1

29、0951爱荷华州-1.28847-1.28140-0.74399-1.28239-1.07227-0.24132-0.96619密苏里州0.998451.60366-1.018400.463260.716600.200380.85236北达科他州-1.65229-1.28140-0.97904-1.67846-1.65365-0.83916-1.38701南达科他州-0.79470-0.61457-0.86630-1.10636-1.16698-0.70438-1.37198内布拉斯加州-1.00260-0.85953-0.43873-0.29955-1.07753-0.46952-0.380

30、04堪萨斯州-0.53483-0.15187-0.24941-0.40223-0.085760.50464-0.55538特拉华州0.218820.406090.073920.888660.371980.528660.95757马里兰州0.608631.127353.084971.724810.358830.658100.90246弗吉尼亚州0.40073-0.02939-0.069670.11119-0.285690.07494-0.35499西弗吉尼亚州-0.17100-1.21336-0.75781-0.75429-1.38269-1.53442-1.38200北卡罗来纳州1.28431

31、-0.36961-0.484472.31158-0.39092-0.78178-0.99625南卡罗来纳州2.011960.33804-0.408950.844660.29043-0.55359-0.56039乔治亚州2.19388-0.75066-0.029250.60995-2.16927-0.07318-0.29488佛罗里达州1.518200.896000.931182.076881.703111.250600.14599肯塔基州1.10240-0.26074-0.27388-0.18219-0.55402-0.94992-0.10951田纳西州0.50468-0.01579-0.17

32、6030.49259-0.28306-1.19947-0.39507阿拉巴马州1.258320.39248-0.513191.16738-0.39881-1.06602-0.72572密西西比州1.20635-0.91397-0.846090.06719-1.48266-1.64117-1.45214阿肯色州0.842520.20195-0.563180.21388-0.07787-0.95126-1.29683路易斯安那州1.258321.018480.449371.50477-0.06472-0.394790.08587俄克拉荷马州-0.24896-0.00218-0.46426-0.12

33、352-0.19362-0.35075-0.44016得克萨斯州1.232330.732700.378110.874000.621900.345840.13096蒙大拿州-0.95063-0.69622-0.81099-0.88632-0.843400.04292-0.73073爱达荷州-0.58680-0.45127-0.83014-0.72496-0.632940.38721-1.12149怀俄明州-0.30094-0.45127-0.81419-0.91566-0.706600.16702-1.01628科罗拉多州-0.171002.774010.324930.727301.226961

34、.425421.10285新墨西哥州0.660610.82796-0.345140.888660.598220.646090.12094阿利桑那州0.686591.549220.230271.152711.521592.170050.66700犹他州-0.89865-0.64179-0.48341-0.959660.001051.21324-0.25981内华达州0.504680.542170.955640.683291.963551.494811.46857华盛顿州-0.872660.32443-0.05371-0.431571.384801.23993-0.02936俄勒冈州-0.586800.324430.04201-0.050170.942841.20256-0.17464加利福尼亚州0.010912.651531.152411.328742.205571.999241.60884阿拉斯加州1.388261.42674-0.284510.47793-0.327780.345840.91749夏威夷州-0.84668-0.51931-0.37492-1.355731.424261.577551.0677814