第11讲因子分析(2).pptVIP

西南科技大学生命科学与工程学院周海廷制作西南科技大学生命科学与工程学院周海廷制作第11讲因子分析(2)*西南科技大学生命科学与工程学院周海廷制作*第1页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*4.4用FACTOR过程实现因子分析例4.1某医学院测得20例肝病患者的４项肝功能指标∶SGPT(转氨酶)、肝大指数、ZnT(硫酸锌浊度)和AFP(胎甲球)，依次用X1至X4表示，试进行因子分析。第2页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*第3页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*程序的第２部分写了４个过程步，每个过程步中都含有４个重要的选择，即∶

１．提取初始公因子的主要方法（即在选择项METHOD=或M=后填１个选择项）

①PRINCIPAL主成分分析法

②PRINIT迭代主因子分析法

③ML最大似然法２．给各变量的公因子方差h2i赋初值(即在PRIORS＝后填１个选择项，中括号内为

缩写形式)第4页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*通常在用主成分分析法时，才让PRIORS=ONE。PRIORS=后其他选择项及其含义如下：

①ASMC[A]置h2i正比于xi与其他指标之间的全相关系数的平均;

②MAX[M]置h2i为xi与其他指标中某一个的最大绝对相关系数;

③ONE[O]置所有的h2i＝１;

④RANDOM[R]置h2i为在０与１之间服从均匀的伪随机数;

⑤SMC[S]置h2i为xi与其他指标之间全相关系数的平均。第5页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*注意：上述２个选择项通常是配套发挥作用的，如：当M=PRINPRIORS=ONE时则指定用

主成分分析法；若M＝PRINPRIORS＝其他内容时则指定用无迭代的主因子分析法；当M＝PRINIT，无论PRIORS=什么，则指定用迭代的主因子分析法；当M=ML，无论PRIORS=什么，则指定用最大似然法。M=后还可写其他一些方法(从略)。第6页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*３．规定因子轴的旋转法(即在选择项ROTATE=或R=后填１个选择项)

①VARIMAX[V]最大方差旋转法;

②ORTHOMAX正交最大方差旋转法;

③PROMAX在正交最大方差旋转基础上进行斜交旋转。４．SCORE表示要求算出因子得分模型中的系数。

另外，用户还可事先在过程步语句中增加选择项NFACTOR=n，要求保留n个公因子,否则,只保留特征值大于１的那些公因子。第7页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*由以上的解释可知：本程序的４个过程分别用主成分分析法、无迭代主因子分析法、有迭代主因子分析法、最大似然法提取初始因子，都用最大方差旋转法旋转、并求出得分因子模型。第8页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*输出结果及其解释这是用主成分分析法提取初始公因子的第１部分结果，相关矩阵的特征值总和为4(指标数)，前２个特征值1.718252和1.093536都大于１,下面将根据这２个较大的特征值提取２个相应的初始公因子。第9页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*含有２个公因子的初始公因子模型为:

x1=0.91753F1+0.09935F2x2=0.90420F1-0.29662F2

x3=0.11527F1+0.94550F2x4=0.21337F1+0.31893F2第10页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*第１和第２公因子能解释的方差分别为1.718252和1.093536。４个标准化变量最后的共性估计值之和为2.811787，与x1，x2，x3，x4相对应的共性估计值依次为：0.851727，0.905562，0.907256，0.147242。第11页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*第12页，共33页，星期日，2025年，2月5日*西南科技大学生命科学与工程学院周海廷制作*经最大方差旋转法旋转后