1、2019/6/2,西南科技大学生命科学与工程学院周海廷制作,1,第十一讲,因子分析,Factor Analysis（2）,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,2,4.4 用FACTOR过程实现因子分析,例4.1某医学院测得20例肝病患者的项肝功能指标SGPT(转氨酶)、肝大指数、ZnT(硫酸锌浊度)和AFP(胎甲球)，依次用X1至X4表示， 试进行因子分析。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,3,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,4,程序的第部分写了个过程步，每个过程步中都含有个重要的选择，即 提取初始公因子的主要方法
2、（即在选择项METHOD=或M=后填个选择项） PRINCIPAL 主成分分析法 PRINIT 迭代主因子分析法 ML 最大似然法 给各变量的公因子方差h2i赋初值(即在PRIORS后填个选择项，中括号内为 缩写形式),2019/6/2,西南科技大学生命科学与工程学院周海廷制作,5,通常在用主成分分析法时，才让PRIORS=ONE。 PRIORS=后其他选择项及其含义如下： ASMCA 置h2i正比于xi与其他指标之间的全相关系数的平均; MAXM 置h2i为xi与其他指标中某一个的最大绝对相关系数; ONEO 置所有的h2i ; RANDOMR 置h2i为在与之间服从均匀的伪随机数; SMC
3、S 置h2i为xi与其他指标之间全相关系数的平均。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,6,注意：上述个选择项通常是配套发挥作用的，如：当M=PRIN PRIORS=ONE时则指定用 主成分分析法； 若MPRIN PRIORS其他内容时则指定用无迭代的主因子分析法； 当MPRINIT，无论PRIORS=什么，则指定用迭代的主因子分析法； 当M=ML，无论PRIORS=什么，则指定用最大似然法。M=后还可写其他一些方法(从略)。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,7,规定因子轴的旋转法(即在选择项ROTATE=或R=后填个选择项) VARIMAX
4、V 最大方差旋转法; ORTHOMAX 正交最大方差旋转法; PROMAX 在正交最大方差旋转基础上进行斜交旋转。 SCORE表示要求算出因子得分模型中的系数。 另外，用户还可事先在过程步语句中增加选择项NFACTOR=n，要求保留n个公因子,否则,只保留特征值大于的那些公因子。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,8,由以上的解释可知：本程序的个过程分别用主成分分析法、无迭代主因子分析法、有迭代主因子分析法、最大似然法提取初始因子，都用最大方差旋转法旋转、并求出得分因子模型。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,9,输出结果及其解释,这是用主成
5、分分析法提取初始公因子的第部分结果，相关矩阵的特征值总和为4(指标数)，前个特征值1.718252和1.093536都大于,下面将根据这个较大的特征值提取个相应的初始公因子。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,10,含有个公因子的初始公因子模型为: x1=0.91753F1+0.09935F2 x2=0.90420F1-0.29662F2 x3=0.11527F1+0.94550F2 x4=0.21337F1+0.31893F2,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,11,第和第公因子能解释的方差分别为1.718252和1.093536。个标准化
6、变量最后的共性估计值之和为2.811787， 与x1， x2， x3， x4相对应的共性估计值依次为：0.851727，0.905562，0.907256，0.147242。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,12,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,13,经最大方差旋转法旋转后的因子模型为: x1= 0.87226G1+0.30149G2 x2= 0.94758G1-0.08748G2 x3=-0.09851G1+0.94739G2 x4= 0.13687G1+0.35848G2 旋转后的第和第公因子能解释的方差分别为1.687177和1.1
7、24611；个标准化指标共性之和以及它们各自的共性估计值与旋转前相同。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,14,因子得分模型为: G1=0.50027x1+0.57347x2-0.12744x3+0.05600x4G2=0.20766x1-0.14705x2+0.85781x3+0.31200x4,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,15,下面分别是由第个过程步产生的最后部分结果，即因子得分模型。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,16,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,17,下面仅根据第个过
8、程步给出的完整结果，对公因子进行命名，并加以应用。,从旋转后的因子模型看，第公因子G1基本上支配着x1和x2；第公因子G2基本上支配着x3。这揭示(以专业知识为依据):,第公因子反映肝脏慢性炎症，可简称为慢性炎症因子。,第公因子反映肝脏急性炎症，可简称为急性炎症因子；,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,18,设现有一肝病病人的项肝功能指标测量值为: X1=50、 X2=2.0、 X3=31、 X4=45 把它们分别代入用原指标表达的公因子得分模型，求得: G1=-0.922464， G2=1.824968 该病人在反映慢性炎症的第公因子上的得分明显大于反映急性炎症的第
9、公因子上的得分，故可诊断该病人很可能为肝脏慢性炎症患者。,将标准化变量xj还原成原变量Xj，得到用原指标表达的公因子得分模型:,G1=0.005628X1+0.543702X2-0.017176X3+0.002560X4-1.874012 G2=0.002336X1-0.139417X2+0.115612X3+0.014260X4-2.238670,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,19,例2 以10名游泳选手为受试者，观测了气力、耐力、速度与协调性四项指标，观察值如下表。试分析这四项体能的共同因素。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,20,
10、【SAS 程序】 OPTIONS NODATE PS=40 LS=80 NONUMBER； DATA FACTOR； INPUT SUBJECT $ A B C D； 资料集定义 LABEL A=气力 B=耐力 C=速度 D=协调； CARDS；,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,21,PROC FACTOR DATA=FACTOR OUTSTAT=FACTOUT OUT=FOU METHOD=P PRIORS=SMC ROTATE=V NFACTOR=2 因子分析(主轴法) SIMPLE CORR SCORE； VAR A B C D； PROC PRINT DAT
11、A=FACTOUT； 打印因子分析的统计输出 PROC PRINT DATA=FOUT； 打印因子分数的资料输出 RUN；,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,22,结果输出,(1)平均数、标准差、及相关矩阵 Means and Standard Deviations from 10 observations,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,23,Correlations(相关性),由上表得知气力和耐力之间有高度正相关，而速度与协调之间也呈现高度正相关。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,24,(2)抽取共同因子：主
12、轴法(PRINCIPAL FACTOR METHOD) (未转轴) 估计共同性的方法：复相关平方估计法(SMC)Initial Factor Method：Principal FactorsPrior Communality Estimates：SMC(共同性估计值),Eigenvalues of the Reduced Correlation Matrix：Total = 2.80800128 Average = 0.70200032,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,25,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,26,2 factors wil
13、l be retained by the NFACTOR criterion.,因子組型(未转轴),Factor Pattern,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,27,解释量 Variance explained by each factor,Initial Factor Method: Principal Factors,共同性 Final Communality Estimates: Total = 2.997940,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,28,由以上报表所知，第一因子解释性较第二因子強，且未转轴時，以 C (速度)因子共通
14、性較強，解释83.69 。而由相关矩阵得知，A、B具有高度正相关，C、D 也具有高度正相关。我們大致可区分成A、B一类，C、D一类，即气力和耐力、速度和协调分类。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,29,(3)转轴法：正交转轴法(最大变异法，VARIMAX ROTATION) Rotation Method：Varimax,转换矩阵 Orthogonal Transformation Matrix,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,30,因素組型(已转轴) Rotated Factor Pattern,2019/6/2,西南科技大学生命科学与
15、工程学院周海廷制作,31,解释量 Variance explained by each factor,共同性 Final Communality Estimates：Total = 2.997940,经过转轴后，A、B、C、D之共同性依然不变。,2019/6/2,西南科技大学生命科学与工程学院周海廷制作,32,(4)因子分数系数Scoring Coefficients Estimated by Regression Squared Multiple Correlations of the Variables with each FactorFACTOR1 FACTOR20.855688 0.799598 Rotation Method：VarimaxStandardized Scoring Coefficients,(标准化后之因子得分),2019/6/2,西南科技大学生命科学与工程学院周海廷制作,33,经过转轴后，将因子与因子之差异拉开，可明显区分成气力与耐力、速度与协调二类。【结论】 由以上报表我们可得知，对于影响游泳选手的因子气力、耐力、速度、协调，经过转轴后，可将气力、耐力归为一类，将其命名为体能因子；而速度、协调可将其归为技巧因子，所以影响选手测验成绩我們可大致知道有二类共同因素。,
