《多元统计分析》（第6版）课件 第6-11章：因子分析理论与应用 -多维标度法理论与应用

因子分析理论与应用CONTENTS目录01因子分析概述02因子分析的基本思想与理论03因子载荷的求解方法04因子旋转与因子得分CONTENTS目录05因子分析与主成分分析的区别06因子分析的步骤与逻辑框图07因子分析的上机实现08案例分析因子分析概述01学习目标

01理解因子分析方法的思想领会因子分析通过降维将复杂变量归结为少数公共因子的核心思路，把握其对变量相关性分组的本质。

02了解因子分析的基本理论掌握因子分析模型的构成，包括公共因子、特殊因子、因子载荷等关键概念及相关数学表达。

03掌握求解因子的方法步骤熟悉从确定因子载荷、因子旋转到计算因子得分的完整流程，明确各环节的实现方式。

04分辨因子分析与主成分分析的异同对比两种方法在模型设定、变量关系、假设条件等方面的区别与联系，理解各自适用场景。

05能够用SPSS软件进行因子分析并理解结果学会运用SPSS的FactorAnalysis模块执行因子分析，正确解读输出的载荷矩阵、方差解释表等结果。2026/5/14因子分析的定位与核心思想作为主成分分析的推广，因子分析基于降维思想，从原始变量相关矩阵内部依赖关系出发，将错综复杂的变量归结为少数综合因子。因子分析的特点与主成分分析相比，更倾向于描述原始变量之间的相关关系，出发点为原始变量的相关矩阵，通过公共因子和特殊因子构建变量分解模型。因子分析的发展与应用领域思想始于1904年查尔斯·斯皮尔曼对学生考试成绩的研究，如今已成功应用于心理学、医学、气象、地质、经济学等多个领域，理论方法不断丰富。因子分析的概念因子分析的基本思想与理论02基本思想

变量分组原则根据相关性大小对原始变量分组，同组内变量相关性较高，不同组变量相关性较低，每组代表一个基本结构即公共因子。

原始变量构成原始变量可分解为两部分：少数不可观测公共因子的线性函数与该变量特有的特殊因子，特殊因子与公共因子相互独立。

物价变动案例反映物价变动无需调查所有商品价格，通过提取“综合商品”价格这一公共因子，即可代表某类商品物价变动，体现降维与信息浓缩思想。2026/5/14斯皮尔曼的例子

研究背景与数据1904年斯皮尔曼研究33名学生古典语（C）、法语（F）等6门考试成绩，得到相关矩阵，发现非对角元素大致成比例规律。

单公共因子模型提出模型xi=aiF+ei，xi为标准化成绩（均值0、方差1），F为公共因子（一般智力，均值0、方差1），ei为特殊因子且与F独立。

因子载荷与共同度ai为因子载荷，其平方a²i称为共同度，代表公共因子解释xi方差的比例，满足1=a²i+var(ei)，var(ei)为特殊因子方差（特殊度）。2026/5/14一般因子分析模型模型数学形式xi=ai1F1+ai2F2+…+aimFm+ei，xi为标准化变量（均值0、方差1），F1…Fm为独立公共因子（均值0、方差1），ei为特殊因子（均值0、与F独立）。基本假设条件原始变量X均值向量E(X)=0、协方差矩阵Σ=相关阵R；公共因子F协方差矩阵为单位阵I；特殊因子ε协方差矩阵为对角阵且与F独立。矩阵表达形式模型矩阵形式为X=AF+ε，其中X=(X1…Xp)为可观测变量向量，A为p×m因子载荷矩阵，F=(F1…Fm)为公共因子向量，ε=(ε1…εp)为特殊因子向量。2026/5/14因子载荷aij的含义aij是xi与Fj的协方差，因xi和Fj均为标准化变量（均值0、方差1），故aij同时也是两者的相关系数，其绝对值反映xi与Fj的相依程度。变量共同度h²ih²i=a²i1+a²i2+…+a²im，代表所有公共因子解释xi方差的比例，满足var(xi)=1=h²i+var(εi)，h²i越大，因子分析效果越好。公共因子方差贡献g²jg²j=a²1j+a²2j+…+a²pj，是公共因子Fj对所有原始变量方差贡献的总和，用于衡量Fj的相对重要性，g²j越大，Fj对X的影响越显著。载荷矩阵的统计意义因子载荷的求解方法03主成分法基本思路先对数据进行主成分分析，将前m个主成分作为未旋转的公共因子，通过对主成分标准化处理得到因子模型。载荷矩阵与共同度载荷矩阵A=(√λ1γ1,...,√λmγm)，其中λ为相关阵特征根，γ为标准正交特征向量；共同度h²i=Σ(aij)²，反映公共因子对变量方差的解释比例。优缺点及适用情况优点：简单直观，易于实现；缺点：特殊因子不独立，不完全符合因子模型假设。适用于共同度较大、特殊因子影响可忽略的场景。2026/5/14主轴因子法

核心原理以调整相关矩阵R*=R-Σε（主对角线为共同度h²i）为出发点，求解其特征根与特征向量，得到因子载荷矩阵A=√λ*γ*。

与主成分法的区别主成分法基于原始相关阵，解释全部方差；主轴因子法基于调整相关阵，仅解释公共因子方差，更符合因子模型“部分方差解释”的假设。

共同度初始估计通常先通过主成分分析得到初始共同度估计，再迭代优化调整相关矩阵，直至结果稳定。2026/5/14极大似然法

假设前提假定公共因子F和特殊因子ε均服从正态分布，即F~N(0,I)，ε~N(0,Σε)，且两者相互独立。

估计方法通过极大化似然函数估计因子载荷A和特殊因子方差Σε，需添加唯一性条件A'Σε⁻¹A=Λ（Λ为对角阵）以确定唯一解。

特点理论上更严谨，依赖正态分布假设，在大样本下估计效果较好，适用于对模型假设有严格要求的分析场景。2026/5/14因子旋转与因子得分04因子旋转的目的与方法因子旋转的核心目的解决初始因子解意义模糊问题，通过线性组合使各主因子表达式中变量系数差异增大，让公共因子实际意义更明确，便于对实际问题分析。正交旋转的特点与方法保持公共因子彼此独立，由初始载荷矩阵右乘正交阵实现。常用方差最大正交旋转，目标是使各列元素平方的相对方差之和最大，使载荷系数接近0或±1。斜交旋转的特点与方法放弃因子间独立限制，可能得到更简洁形式，实际意义更易解释。常用最优斜交旋转（promax方法），输出结果含因子载荷矩阵（patternMatrix）和相关阵（structureMatrix）。2026/5/14因子得分的定义指公共因子F1,F2,…,Fm在每一个样品点上的得分，用于反映样品在公共因子上的取值，可据此对样品性质及相互关系进行分析。因子得分的计算方法基于回归思想，以公共因子为因变量、原始变量为自变量建立回归方程，在最小二乘意义下得到估计值公式：F=A'R⁻¹X，其中A为因子载荷矩阵，R为原始变量相关阵，X为原始变量向量。因子得分的应用可用于样本点比较分析、聚类分析，当因子数较少时，能将样本点在因子构成的空间中标示，直观描述样本分布，还可代替原始数据进行后续回归等分析。因子得分的概念与计算因子分析与主成分分析的区别05模型与目的差异

因子分析的模型与目的因子分析将变量表示为公共因子和特殊因子的线性组合，目的是探寻对变量起解释作用的公共因子和特殊因子，以及它们的组合系数，从数据中提取潜在结构。主成分分析的模型与目的主成分分析把主成分表示为各变量的线性组合，目的是从空间生成角度寻找能解释诸多变量绝大部分变异的几组彼此不相关的新变量（主成分），实现数据降维。2026/5/14假设与提取方法不同

因子分析的假设与提取方法因子分析需假设公共因子之间不相关、特殊因子之间不相关、公共因子和特殊因子之间不相关；提取方法多样，包括主成分法、主轴因子法、极大似然法等。主成分分析的假设与提取方法主成分分析不需要专门假设；仅用主成分法提取，主成分数量固定（一般有几个变量就有几个主成分），且主成分固定，因子分析中因子可旋转得到不同结果。2026/5/14解释与应用场景区别

因子分析的解释与应用场景因子分析可使用旋转技术帮助解释因子，在解释方面更有优势，适用于需要明确因子实际意义的场景，如探究影响学生成绩的潜在能力因子等。

主成分分析的解释与应用场景主成分分析适合将现有变量变成少数几个几乎带有原来所有变量信息的新变量，用于后续分析；实际中也可通过计算因子得分处理类似场景，但区分并非绝对。2026/5/14因子分析的步骤与逻辑框图06因子分析的步骤选取原始变量根据研究问题确定分析对象，选择具有相关性的指标变量，如衡量企业经济效益的多项财务指标。数据标准化与相关阵分析对原始变量进行标准化处理（均值为0，方差为1），计算相关矩阵并分析变量间相关性，通过KMO检验（>0.7适合分析）和Bartlett球形检验（拒绝单位阵假设）判断适用性。求解初始公共因子及载荷矩阵采用主成分法、主轴因子法等提取公共因子，确定因子载荷矩阵，如主成分法通过相关阵特征根与特征向量求解，保留特征值>1的因子。因子旋转对初始因子进行正交（如方差最大旋转）或斜交旋转，使载荷系数向0或±1集中，便于解释因子意义，旋转后共同度不变但载荷矩阵更新。计算因子得分通过回归方法建立公共因子与原始变量的线性关系（如F=A'R⁻¹X），得到因子得分矩阵，用于样品评分及后续分析。因子得分的进一步分析利用因子得分进行样本比较、聚类分析或绘制因子得分散点图，直观展示样品分布特征，辅助决策。2026/5/14因子分析的逻辑框图数据输入与预处理阶段

输入原始变量数据，完成标准化转换，计算相关矩阵并进行KMO和Bartlett检验，确保数据适合因子分析。因子提取与载荷矩阵构建阶段

基于相关阵或调整相关阵，通过主成分法等提取公共因子，求解初始因子载荷矩阵，确定因子数量（如特征值>1准则）。因子旋转与解释优化阶段

对初始因子载荷矩阵进行旋转（正交/斜交），得到结构更清晰的旋转后载荷矩阵，明确各公共因子的实际意义。因子得分计算与应用阶段

通过回归模型计算因子得分，将得分用于样本排序、分类或可视化（如散点图），实现对原始数据的降维分析与解释。步骤间逻辑关系

各阶段依次衔接：预处理为因子提取提供数据基础，旋转优化因子解释性，得分计算将抽象因子转化为可应用的量化指标，共同构成完整分析流程。2026/5/14因子分析的上机实现07SPSS操作步骤

01模块选择与变量导入依次点选Analyze→DimensionReduction→Factor，进入FactorAnalysis对话框，将待分析的指标变量选入Variables框中。

02提取方法与参数设置点击Extraction按钮，在Method选项框选择提取方法（默认主成分法），Analyze选项框默认从相关阵出发，Extract选项框可通过特征值大于1或固定因子数目（如输入2）确定因子个数，完成后点击Continue。

03因子得分设置点击Scores按钮，选中Displayfactorscorecoefficientmatrix选项以输出因子得分系数矩阵，点击Continue后返回主对话框，点击OK运行分析。2026/5/14输出结果解读（一）01共同度表（Communalities）展示各变量的初始共同度（均为1.000）和提取共同度，提取共同度反映变量被公共因子解释的程度，如X2提取共同度为0.993，表明其方差的99.3%可由公共因子解释。02总方差解释表（TotalVarianceExplained）包含初始特征值、提取平方和载荷，显示各因子的方差贡献及累计解释率，如前2个因子累计解释96.809%的方差，说明其能较好概括原始变量信息。03因子载荷阵（ComponentMatrix）呈现各变量在公共因子上的载荷系数，如X1在因子1上载荷0.913、因子2上0.320，表明X1与因子1相关性较强，载荷绝对值越大，变量与因子关系越密切。2026/5/14因子得分系数矩阵给出公共因子用标准化原始变量表示的线性系数，如F1=0.155X1+0.168X2+…-0.097X7，可直接代入变量值计算因子得分。因子载荷与得分系数的关系因子得分系数等于因子载荷除以对应因子的特征根，如X1在因子1的得分系数0.155=0.913/5.910（因子1特征根为5.910），因公共因子需标准化为方差1。主成分法下的因子模型保留前m个主成分作为公共因子，剩余部分为特殊因子，如X1=0.913F1+0.320F2+特殊因子，特殊因子方差为1-累计解释率（如3.191%）。输出结果解读（二）变量相关性检验KMO检验用于判断变量间相关性和偏相关性，KMO值越接近1效果越好，0.7以上适合因子分析，案例中KMO值为0.718，表明适合进行因子分析。Bartlett球形检验原假设为相关阵是单位阵，若拒绝原假设则变量相关，案例中Bartlett统计量为645.127，sig.=0.000，拒绝原假设，说明变量间存在较强相关性。检验操作路径在FactorAnalysis对话框中点击Descriptives，勾选CorrelationMatrix下的KMOandBartlett'stestofsphericity，可输出检验结果。2026/5/14因子旋转结果分析正交旋转（Varimax）旋转后因子载荷矩阵元素更倾向0或±1，如旋转后X1在因子1载荷0.964、因子2载荷0.082，变量与因子关系更清晰，且公共因子仍保持独立，旋转3次收敛。斜交旋转（Promax）放弃因子独立性约束，输出PatternMatrix（因子载荷阵）和StructureMatrix（因子与变量相关阵），如X7在因子2载荷-1.052，因子间相关系数为0.546，解释更简洁。旋转效果对比旋转不改变共同度和累计方差解释率，但载荷矩阵和得分系数矩阵变化，正交旋转适合因子独立场景，斜交旋转适合追求因子实际意义的场景。2026/5/14因子得分的保存与应用

因子得分保存在FactorScores对话框中勾选Saveasvariables，系统默认用回归法计算得分，数据窗口新增FAC1-1、FAC2-1等变量，分别对应各公共因子得分。

得分标准化验证通过Descriptives分析得分变量，均值为0、标准差为1，如案例中FAC1-1均值0.000、标准差1.000，符合标准化要求。

得分应用场景可用于样品比较（如企业经济效益排名）、聚类分析（划分企业类型）、绘制散点图（直观展示样品在因子空间的分布，如广东、江苏在因子1得分较高）。2026/5/14案例分析08餐饮企业经济效益因子分析

数据来源与指标选择沿用衡量限额以上餐饮企业经济效益的7项指标数据，通过SPSS的FactorAnalysis模块进行因子分析，以建立因子分析模型。

因子分析模型建立在FactorAnalysis对话框中，将7个指标变量选入variables框，提取特征值大于1的因子，默认用主成分法提取，设置提取因子数目为2，并要求输出因子得分矩阵。

结果解读与效益评价通过输出的公共因子方差解释比例表、因子载荷阵和因子得分系数矩阵等结果，可得到公共因子对原始变量的解释情况及因子得分表达式，进而对餐饮企业经济效益进行评价，如分析各企业在不同公共因子上的表现。2026/5/14研究背景与数据说明近年来交通运输业上市公司股票投资受关注，选取30家该行业上市公司的8项财务指标数据（来源于巨潮资讯），包括基本每股收益、每股净资产等，用于评价公司整体投资效益。因子分析方法应用使用因子分析方法，对8项财务指标进行分析，提取公共因子，通过因子载荷矩阵等结果明确各公共因子的意义，如反映企业规模、盈利能力等方面。为投资者决策提供依据根据因子得分对上市公司投资效益进行评价，帮助投资者掌握公司经营状况、盈利状况等，以便更好地进行投资决策，例如判断哪些公司投资效益较好。交通运输业上市公司投资效益评价对应分析理论与实践应用对应分析的方法步骤04对应分析的基本理论03列联表及列联表分析02对应分析概述01目录CONTENTS对应分析的应用价值与总结07对应分析案例分析06SPSS软件进行对应分析操作05目录CONTENTS01对应分析概述2026/5/14对应分析的学习目标

理解列联表分析及对应分析的基本思想明确列联表用于描述属性变量状态及关系的作用，掌握对应分析通过降维同时处理行与列数据、揭示变量间关系的核心思路。了解对应分析的基本理论知晓对应分析中涉及的行剖面、列剖面、距离、总惯量等关键概念，理解R型与Q型因子分析的对等关系及原理。掌握对应分析的方法熟悉从列联表数据规格化、计算Z矩阵，到进行因子分析、提取公共因子并在低维图形呈现结果的完整分析流程。能用SPSS软件进行对应分析并正确理解输出结果能够运用SPSS的对应分析模块处理数据，对输出的列联表、总惯量、因子载荷、二维投影图等结果进行合理解释与分析。对应分析的基本思想R型与Q型因子分析的结合对应分析整合了R型因子分析（研究变量间关系）和Q型因子分析（研究样品间关系），弥补单独分析的局限，实现对变量和样品的同时研究。降维思想的应用借助降维手段，将高维数据结构简化，把属性变量不同取值以二维图上的点表示，直观呈现行与列及各状态间的相互关系。行与列的同时处理区别于因子分析仅对行或列单独处理，对应分析同时对数据表的行与列进行分析，以低维图形形式清晰展示两者间的关联。2026/5/14对应分析的发展历程

思想的提出1933年，理查森（Richardson）和库德（Kuder）首次提出对应分析的基本思想，为该方法的发展奠定了基础。

方法的详细论述与发展法国统计学家让·保罗·贝内泽（Jean-paulBenzécri）和日本统计学家林知己夫（Hayashichikio）对对应分析方法进行了深入且详细的论述，推动其进一步发展和完善。02列联表及列联表分析2026/5/14列联表的概念与形式列联表的定义列联表是描述属性变量（定类或定序尺度）各种状态或相关关系的表格形式，通过横栏与纵列交叉位置的频数呈现数据。两变量列联表示例以公司消费者满意度调查为例，横栏为不同职业（如一般工人、管理者等），纵列为评价等级（非常满意、比较满意等），交叉单元格为相应频数，底部和右侧为汇总行/列。经济研究中的列联表应用在企业研究中，可将横栏设为企业规模，纵列设为获利能力水平，通过列联表分析两者关系；也可按上市/非上市、行业等分类，研究盈利能力、偿债能力等指标。2026/5/14列联表的作用与信息提取基础信息获取从列联表可直接获取不同类别下的频数分布，如各职业消费者对产品的评价频数、所有被调查者的整体评价及职业构成情况。类别间相关关系判断通过比较交叉单元格频数（需抵消不同类别在总样本中的比例影响），判断类别相关性。例如，若管理者与“比较满意”交叉单元格频数相对较大，说明两者有较强相关性。综合指标分析支持列联表可容纳简单或综合指标，如将因子分析提取的公共因子按取值范围分类后，研究其与其他属性变量的关系，为多维度分析提供基础。频率意义上的列联表频率列联表的构建将原始列联表中每个元素除以总频数T，得到频率意义上的列联表，其中每个单元格元素pij表示特性A第i状态与特性B第j状态同时出现的概率。边缘概率与符号定义记PⅠ=（p1·,p2·,…,pn·）为行边缘概率向量，PJ=（p·1,p·2,…,p·m）为列边缘概率向量，且满足PⅠ1=1、PJ1=1（1为全1向量），表格右下角合计为1。概率意义的应用频率列联表中，pij具有概率含义，可用于研究属性变量独立性：若A与B独立，则对任意i、j有pij=pi·×p·j，为后续对应分析等方法提供数据基础。03对应分析的基本理论2026/5/14行剖面与列剖面

行剖面的定义当属性变量A的取值固定为i（i=1,2,…,n）时，变量B各个状态相对出现的概率情况，即矩阵P第i行元素除以行边缘概率pi·，坐标为（p_i1/p_i·,p_i2/p_i·,…,p_iq/p_i·），表示条件概率P(B=j|A=i)，且各分量之和为1。

列剖面的定义当属性变量B的取值固定为j（j=1,2,…,q）时，变量A各个状态相对出现的概率情况，即矩阵P第j列元素除以列边缘概率p·j，坐标为（p_1j/p·j,p_2j/p·j,…,p_nj/p·j），表示条件概率P(A=i|B=j)，且各分量之和为1。

欧氏空间表示行剖面集合是n个点在q维欧氏空间超平面x₁+x₂+…+x_q=1上的分布；列剖面集合是q个点在n维欧氏空间超平面x₁+x₂+…+x_n=1上的分布，对应分析通过降维将两类点同时呈现在二维图中。2026/5/14距离与总惯量加权距离公式

为消除变量B各状态边缘概率影响，行剖面k与l的加权距离公式为：d(k,l)=√[Σ(p_kj/p_k·-p_lj/p_l·)²/p·j]，同理列剖面s与t的加权距离为：d(s,t)=√[Σ(p_is/p·s-p_it/p·t)²/p_i·]。行剖面总惯量

行剖面集合各点与其重心（坐标为√p·j）的加权欧氏距离之和，公式为I_r=Σp_i·Σ[(p_ij/p_i·-p·j)²/p·j]，反映行剖面间的总差异，与χ²统计量关系为χ²=T·I_r（T为总频数）。列剖面总惯量

列剖面集合各点与其重心（坐标为√p_i·）的加权欧氏距离之和，公式为I_c=Σp·jΣ[(p_ij/p·j-p_i·)²/p_i·]，且行剖面总惯量等于列剖面总惯量（I_r=I_c），均反映属性变量间的相关关系。2026/5/14R型与Q型因子分析的对等关系

协方差矩阵关系对行剖面数据做变换Z_ij=(p_ij-p_i·p·j)/√(p_i·p·j)，则R型因子分析协方差矩阵Σ_r=Z^TZ，Q型因子分析协方差矩阵Σ_c=ZZ^T，二者存在对等关系。

非零特征根相同矩阵Σ_r与Σ_c具有完全相同的非零特征根λ₁≥λ₂≥…≥λ_r，这些特征根表示各公共因子解释的总惯量份额，且最大维数为min(n,q)-1（n、q分别为行、列状态数）。

因子结果推导若u_j是Σ_r对应特征根λ_j的特征向量，则Z^Tu_j是Σ_c对应λ_j的特征向量，因此可由R型因子分析结果直接推导Q型因子分析结果，实现变量与样品在同一因子平面上的可视化。2026/5/14对应分析应用于定量变量的情况

适用数据类型适用于定距/定比尺度数据，可将每个观测视为一类（行）、每个变量视为一类（列），通过规格化处理转化为类似列联表的结构，分析观测与变量间的相关关系。

数据预处理要求对应分析要求数据阵中所有元素≥0，若存在负值，需先进行加工（如对变量各取值加常数），确保数据非负后再进行分析，例如经济截面数据可通过分类汇总转化为符合要求的格式。

分析思路与属性变量分析步骤一致：计算规格化矩阵→构建Z矩阵→进行R型/Q型因子分析→在低维图中呈现观测与变量的分布，直观揭示变量间依赖关系及观测分类特征。对应分析需要注意的问题二维投影的局限性多维空间中相隔较远的点在二维投影图上可能距离较近，易产生误导，需结合状态对公共因子的贡献度（CTR(i)=p_i·a_ik²/λ_k）判断各点对维度的影响，避免误判。缺乏量化相关统计量对应分析仅通过图形提示变量间关系，无法提供具体统计量（如相关系数）度量相关程度，结论主观性较强，需结合列联表分析、χ²检验等方法交叉验证。对异常值敏感数据中极端值或频数较小的状态可能对总惯量和因子载荷产生较大影响，分析前需检查数据质量，必要时对小样本状态合并或剔除，确保结果稳健性。04对应分析的方法步骤2026/5/14对应分析的步骤01计算规格化概率列联表将原始列联表中每个元素除以总和T，得到频率意义上的列联表，其中元素pij表示特性A第i状态与特性B第j状态同时出现的概率，边缘概率pi●和p●j分别为行、列之和。02计算Z矩阵基于规格化列联表，通过公式zij=(pij-pi●p●j)/√(pi●p●j)计算Z矩阵，该矩阵消除了边缘概率影响，用于后续因子分析。03进行R型或Q型因子分析并推导结果对Z矩阵进行R型因子分析（分析变量）或Q型因子分析（分析样品），利用二者协方差矩阵Σr=ZᵀZ与Σc=ZZᵀ的对等关系，由一种因子分析结果推导出另一种，提取主因子并确定特征根与特征向量。04在二维图上画变量状态并分析相关性取前两个公共因子，将行、列变量的各状态以二维坐标点形式绘制在同一张图中，通过点的距离直观判断变量间及各状态间的相关关系，距离越近表示关联性越强。2026/5/14对应分析的逻辑框图数据输入与预处理输入原始列联表数据，检查数据有效性（如非负性），若为定量数据需先分类汇总，确保符合对应分析数据格式要求。规格化与矩阵转换将原始频数列联表转换为频率列联表，计算边缘概率pi●、p●j，进而构建Z矩阵，完成数据标准化以消除量纲和边缘概率影响。因子分析与结果推导对Z矩阵进行R型或Q型因子分析，求解协方差矩阵特征根与特征向量，利用R型与Q型因子分析的对等关系，共享特征根并推导对应特征向量，确定低维（通常二维）因子空间。图形绘制与结果解释根据因子载荷计算行、列状态的二维坐标，绘制对应分析图，通过点的分布、距离及聚类情况，解释变量间相关性及各状态特征，结合专业知识得出结论。05SPSS软件进行对应分析操作2026/5/14SPSS对应分析模块介绍

模块功能定位SPSS的correspondenceAnalysis模块是专门用于对应分析的功能模块，能够实现对属性变量列联表数据的降维分析，直观展示行与列变量各状态间的关系。

核心分析目标该模块通过对数据的处理和分析，可在二维图上同时呈现两类属性变量的各种状态，帮助研究者揭示变量间及变量各状态间的相关关系，简化数据结构。数据录入与变量设定数据窗口录入以高校星级排名与学校类型为例，将原始数据（如高校名称、星级排名、学校类型及类型取值）依次录入SPSS数据窗口，确保数据准确对应。变量名称设定在variableview窗口，为变量命名，如“university”“rank”“type”，明确变量含义，便于后续分析操作。变量标签设定对分类变量（如type）的取值设置标签，如“1=综合类”“2=理工类等”，使输出结果更具可读性，操作时通过valueLabels对话框完成标签添加。2026/5/14对应分析对话框操作

模块调用路径在SPSS中依次点选Analyze→DimensionReduction→correspondenceAnalysis，打开对应分析主对话框，准备进行变量设置。

行变量选择与范围定义将行变量（如rank）选入Row框，点击DefineRange按钮，设定其取值范围（如4~8），点击Update确认，完成行变量配置。

列变量选择与范围定义按同样方法将列变量（如type）选入Column框，设定取值范围（如1~5），完成后点击OK运行分析，生成初步结果。2026/5/14模型与统计选项设置

模型参数设定在Model对话框中，可规定对应分析的最大维数（默认2）、距离量度方法（如卡方距离或欧氏距离）、标准化方法及正态化方法（如symmetrical），根据研究目的调整参数。

统计量选择输出在Statistics对话框中，可选中Rowprofiles和Columnprofiles选项，输出行剖面与列剖面数据，还可选择其他统计量以检验对应分析效果，丰富分析结果。图形编辑与优化参考线添加在SPSS结果输出窗口双击图形进入编辑窗口，通过Options→XAxisReferenceLine和YAxisReferenceLine，添加X=0和Y=0参考线，清晰展示各点相对位置。标记类型修改选中图形图例中变量对应的标记，双击打开properties对话框，在Marker中选择不同标记类型（如菱形），区分不同变量的状态点，提升图形辨识度。图形可读性优化通过调整标记大小、颜色等外观属性，使二维图上的各状态点分布更清晰，便于直观观察变量间及状态间的关系，辅助分析结论的得出。06对应分析案例分析2026/5/14高校星级排名与学校类型对应分析案例数据与列联表构建选取2025年中国“双一流”大学排行榜中73所高校，以星级排名（4-8星）和学校类型（综合类、理工类等5类）为属性变量，构建5×5列联表，观测总数73，如4星级仅1所综合类高校，8星级含7所综合类和4所理工类。总惯量与显著性检验总惯量为0.409，反映行列关联强度；χ²值29.848（自由度16），Sig=0.019<0.05，拒绝独立性假设，表明星级排名与学校类型显著相关。维度解释与二维投影分析前两维度累计解释90.5%总惯量，二维图显示：7星级集中于理工类，8星级集中于综合类，5-6星级与医药类、文法类等关联，4星级因样本量小无显著类型特征。2026/5/14农村居民人均可支配收入对应分析

01数据背景与标准化方法基于2022年全国31省区农村居民收入数据（工资性、经营净等4类收入），因数据为分类汇总指标，选择“Columntotalsareequalizedandmeansareremoved”标准化方法消除量纲与均数影响。

02维度提取与解释能力总惯量0.384，前两维度累计解释94.8%（第一维81.8%、第二维13.1%），特征根分别为0.314和0.050，可有效降维展示省区与收入类型关系。

03省区与收入类型关联特征二维投影显示：北京、上海等经济发达地区以工资性收入为主，内蒙古、吉林等农业区以经营净收入为主，江西、河南等地转移净收入占比高，财产净收入整体贡献较低。07对应分析的应用价值与总结2026/5/14对应分析的应用价值

直观揭示变量关系通过二维图同时呈现两类属性变量状态，直观展示变量间及状态间相关关系，如高校星级排名与学校类型的关联分析。

高效简化数据结构利用降维思想，在总惯量信息损失最小前提下简化数据，从原始数据中提取较多信息，适用于列联表等复杂数据。

广泛适用于多领域可应用于消费者满意度调查、企业规模与获利能力研究、农村居民收入来源分析等多个领域，解决实际问题。2026/5/14对应分析的局限性

结果解释主观性较强仅通过图形提示变量关系，无法提供具体统计量度量相关程度，易导致研究者得出主观结论。

图形投影存在偏差风险多维空间中相隔较远的点在二维平面投影可能接近，需结合贡献度等指标深入分析，避免误判。

数据要求有一定限制要求数据阵中数值非负，对含负值数据需先加工处理；分类汇总数据需选择合适标准化方法，增加操作复杂度。对应分析的未来发展展望拓展多领域应用场景有望在大数据分析、人工智能等领域发挥作用，如结合机器学习算法处理高维复杂数据，挖掘潜在关联。与其他统计方法融合加强与因子分析、聚类分析等方法的结合，形成更综合的数据分析体系，提升对数据的解释力和预测能力。软件功能与可视化优化进一步优化SPSS等软件的对应分析模块，提供更丰富的标准化方法和可视化工具，增强结果的准确性与可读性。谢谢观看THEEND典型相关分析理论与应用典型相关分析的步骤04典型相关分析模型与求解03典型相关分析的统计思想02典型相关分析概述01目录CONTENTS总结与展望08典型相关分析的应用与拓展07实例分析：城市经济发展与空气质量关系06实例分析：生理指标与运动指标关系05目录CONTENTS01典型相关分析概述典型相关分析的定义核心研究对象典型相关分析是研究两组变量之间整体线性相关关系的多元分析方法，将每组变量作为一个整体进行研究，而非分析组内单个变量。基本思想借用主成分分析降维思想，分别对两组变量提取主成分，使从两组提取的主成分间相关程度达到最大，同一组内部提取的各主成分互不相关。衡量指标用从两组分别提取的主成分（典型相关变量）的简单相关系数（典型相关系数）来描述两组变量整体的线性相关关系。2026/5/14典型相关分析的发展历程思想提出典型相关分析的思想首先由霍特林于1936年提出，为研究两组变量相关性提供了理论基础。发展契机计算机的发展解决了典型相关分析在应用中计算方面的困难，使其从理论走向实际应用。应用现状目前已成为普遍应用的两组变量之间相关性分析的技术，在多个领域得到广泛使用。2026/5/14学习目标理解思想理解典型相关分析研究两组变量整体线性相关关系的核心思想，包括其降维和提取典型相关变量的思路。了解理论方法了解典型相关分析的基本理论，如总体与样本典型相关、典型变量的性质等，以及具体的分析方法步骤。掌握实现与解释掌握利用SPSS或R语言实现典型相关分析的方法，并能正确理解和解释分析过程中产生的各种输出结果，如典型相关系数、典型载荷等。02典型相关分析的统计思想2026/5/14研究对象与要求

研究对象：两组变量整体线性相关关系典型相关分析将每组变量作为整体研究，而非关注组内单个变量间的关系，旨在揭示两组变量间的综合线性关联。

变量尺度要求：至少为间隔尺度两组变量需满足间隔尺度以上测量标准，以确保变量间线性关系的可分析性与结果的有效性。

变量地位：自变量与因变量或同等地位两组变量可存在自变量与因变量的对应关系，也可处于同等分析地位，适用范围广泛。2026/5/14典型相关变量与系数

典型相关变量：线性组合生成的综合变量通过对每组变量构建线性组合（如Ui=aᵢ¹X₁+...+aᵢₚXₚ，Vi=bᵢ¹Y₁+...+bᵢqYq），生成能代表原始变量主要信息的综合变量。

典型相关系数：衡量相关程度的指标典型相关变量间的简单相关系数称为典型相关系数，用于量化两组变量整体的相关强度，按绝对值大小排序。

典型相关变量特性：各对间互不相关不同对典型相关变量之间互不相关，保证每对变量反映两组变量间独立的相关关系，避免信息重叠。典型相关分析的作用降维研究：简化两组变量相关关系通过少数几对典型相关变量代替原两组变量的研究，减少分析维度，更易抓住问题本质与核心关联。显著性检验：筛选有效典型相关变量对典型相关系数进行显著性检验，剔除不显著的典型相关变量，确保分析结果的可靠性与代表性。结合定性分析：深化实际问题解释依据典型相关变量的相关程度及原始变量系数大小，结合研究背景的定性分析，可对实际问题给出深刻的解释与结论。03典型相关分析模型与求解2026/5/14总体典型相关和典型变量随机向量与协方差矩阵

设随机向量(x=(X_1,X_2,…,X_p)^T)，(y=(Y_1,Y_2,…,Y_q)^T)，协方差矩阵(Sigma=begin{pmatrix}Sigma_{11}&Sigma_{12}Sigma_{21}&Sigma_{22}end{pmatrix})，其中(Sigma_{11})、(Sigma_{22})为组内协方差阵，(Sigma_{12}=Sigma_{21}^T)为组间协方差阵，且(Sigma)正定时(Sigma_{12})、(Sigma_{21})正定。线性组合与目标函数

考虑线性组合(U=a^Tx)，(V=b^Ty)，在(var(U)=1)、(var(V)=1)约束下，目标是最大化(corr(U,V)=a^TSigma_{12}b)，通过拉格朗日乘数法转化为求解极值问题。典型相关变量与系数推导

利用拉格朗日乘数法得到方程组(Sigma_{12}b=lambdaSigma_{11}a)、(Sigma_{21}a=lambdaSigma_{22}b)，推导得出(lambda^2)是(Sigma_{11}^{-1}Sigma_{12}Sigma_{22}^{-1}Sigma_{21})的特征根，(a)、(b)为对应特征向量，(lambda)为典型相关系数，(U)、(V)为典型相关变量。总体典型相关和典型变量典型变量的性质同一组典型变量互不相关且方差为1，即(cov(u_i,u_j)=0)（(i≠j)），(var(u_i)=1)；不同对典型变量(u_i)与(v_j)（(i≠j)）互不相关，同一对相关系数为(lambda_i)，即(cov(u_i,v_i)=lambda_i)，(cov(u_i,v_j)=0)（(i≠j)）。2026/5/14样本典型相关和典型变量

总体协方差阵的估计实际中总体协方差阵(Sigma)未知，设样本((x_i,y_i)(i=1,2,…,n))来自正态总体，极大似然估计(hat{Sigma}=frac{1}{n}A)，其中(A)为样本离差阵，当(n>p+q)时(hat{Sigma})以概率1正定。

样本典型相关系数与变量用(hat{Sigma})替代(Sigma)，计算(hat{Sigma}_{11}^{-1}hat{Sigma}_{12}hat{Sigma}_{22}^{-1}hat{Sigma}_{21})的非零特征根(hat{lambda}_1^2geqhat{lambda}_2^2geq…geqhat{lambda}_k^2)（(k=min(rank(hat{Sigma}_{11}),rank(hat{Sigma}_{22})))），(hat{lambda}_i)为样本典型相关系数，对应特征向量(hat{a}_i)、(hat{b}_i)构成样本典型变量(hat{u}_i=hat{a}_i^Tx)、(hat{v}_i=hat{b}_i^Ty)。

典型变量得分及应用将样本数据代入典型变量表达式得到得分，可绘制散点图用于样品分类研究，如通过典型变量得分分布特征对观测对象进行分组分析。2026/5/14典型相关系数的显著性检验检验原理与原假设若两组变量独立，则(Sigma_{12}=0)，典型相关系数(lambda_i=0)。检验第(j)个及后续典型相关系数是否为0，原假设(H_0:lambda_j=lambda_{j+1}=…=lambda_k=0)，备择假设(H_1)：至少(lambda_j≠0)。巴特莱特X²检验统计量计算统计量(Q_j=-(n-1-frac{p+q+1}{2})ln(prod_{i=j}^k(1-hat{lambda}_i^2)))，其中(n)为样本量，(p)、(q)为两组变量个数。(Q_j)近似服从自由度(df=(p-j+1)(q-j+1))的(chi^2)分布。检验步骤与判断标准1.对第一对典型相关系数检验：(j=1)，(df=pq)，若(Q_1>chi_{alpha}^2(pq))，拒绝(H_0)；2.若前(j-1)个显著，检验第(j)个：(df=(p-j+1)(q-j+1))，若(Q_j>chi_{alpha}^2(df))，拒绝(H_0)，否则停止检验。04典型相关分析的步骤确定分析目标判断变量组独立性或关系强弱明确两组变量是相互独立，还是存在关联及关联程度，为后续分析提供基础方向。推导最优权重生成典型变量为每组变量推导权重，使生成的线性组合（典型变量）间相关程度最大，且后续组合与前面组合独立。解释变量对典型函数的贡献通过分析每个变量在典型函数中的相对贡献，揭示自变量组与因变量组间的内在联系。2026/5/14设计分析方案

满足多元分析共同要求需考虑测量误差影响、变量类型适配性及必要的数据变换，与多元回归、判别分析等方法要求类似。

确保充足样本量为避免数据“过度拟合”，建议每个变量至少对应10个观测，小样本可能掩盖有意义的相关关系。

合理选择变量组构成研究者需谨慎确定自变量组和因变量组的变量数量，避免因变量过多而忽视样本量的实际含义。2026/5/14检验基本假定线性假定的两方面影响一是变量间相关系数基于线性关系，非线性关系需对变量进行变换；二是典型相关衡量变量间线性关系，非线性关系可能无法被捕捉。正态性假定的要求与处理虽不严格要求所有变量正态分布，但变量分布（如高度偏态）不应削弱与其他变量的相关性，建议检验单变量正态性，必要时进行变换。多元正态性的统计检验典型函数的多元正态性检验有必要，若多元检验不可行，流行准则是保证每个单变量的正态性以满足分析条件。2026/5/14推导典型函数与评价拟合

典型函数的推导过程类似无旋转因子分析，先提取第一对最大相关典型变量，后续变量基于剩余残差提取，且各对典型变量正交独立，最大提取数为最小变量组的变量数。

典型根与共同方差典型相关系数的平方称为典型根（特征根），表示一个典型变量通过另一个典型变量解释的方差，即两者间的共同方差。

冗余指数的概念与计算冗余指数是一组变量方差被另一组变量方差解释的比例，通过三步计算：共同方差比例、解释方差比例（典型相关系数平方）、两者乘积，克服典型根的潜在偏误。2026/5/14解释典型变量

典型权重（标准化系数）反映原始变量对典型变量的贡献大小和方向，但受变量间相关性影响，解释时需谨慎，不能单独作为变量重要性的依据。

典型载荷（结构系数）是原始变量与典型变量的简单线性相关系数，类似因子载荷，能反映原始变量与典型变量的共同方差，是解释典型函数的基础。

典型交叉载荷使原始因变量与自变量典型变量直接相关，提供更直接测量因变量组与自变量组关系的指标，作为典型载荷的有效替代。2026/5/14验证模型

子样本分析比较在样本量允许时构造两个子样本，分别进行分析，比较典型函数相似性、典型载荷等，确保结果代表总体而非单个样本。

变量剔除灵敏度测量通过测量剔除一个因变量或自变量后结果的变化，评估典型权重和典型载荷的稳定性，验证模型的可靠性。

典型相关分析的局限性包括反映变量组线性组合共享方差而非提取方差、典型权重不稳定性、解释困难、难以识别变量子集间有意义关系等，需在应用中注意。05实例分析：生理指标与运动指标关系数据介绍样本信息研究对象为20名中年男性，测量3项生理指标与3项运动指标，数据来源于《多元统计分析（第6版）》例8-1。生理指标变量包括体重（kg）、腰围（cm）、脉搏（次/分钟），反映身体基本生理状态。运动指标变量包括引体向上（次）、仰卧起坐（次）、跳跃次数（次），衡量运动能力与身体素质。数据样例如序号1：体重191kg、腰围36cm、脉搏50次/分钟；引体向上5次、仰卧起坐162次、跳跃次数60次。2026/5/14R语言实现步骤

01数据读取与变量分组使用read.csv函数读取“例8_1.csv”数据，将前3列设为生理指标（PHY），后3列设为运动指标（EXER）。

02安装与加载CCA包通过install.packages("CCA")安装典型相关分析包，用library(CCA)加载包以调用分析函数。

03相关矩阵计算调用matCor(PHY,EXER)函数，输出两组变量的自相关矩阵及交叉相关矩阵，初步观察变量关系。

04典型相关分析与结果输出使用cc(PHY,EXER)执行分析，通过CC1[1]输出典型相关系数，CC1[3:4]输出原始典型系数，CC1[5]输出得分及载荷矩阵。

05标准化系数与显著性检验计算并输出标准化典型系数（典型权重），编写代码进行巴特莱特X²检验，判断典型相关系数的显著性。2026/5/14结果分析01原始变量相关阵解读体重与腰围正相关（0.870），腰围与引体向上、仰卧起坐负相关（-0.552、-0.646），表明肥胖可能降低运动能力。02典型相关系数与显著性第一典型相关系数为0.796（P=0.062），在0.1水平显著；第二、三系数分别为0.201、0.073，均不显著，故主要关注第一对典型变量。03典型变量表达式生理第一典型变量U1=0.0314体重-0.4932腰围+0.0082脉搏，运动第一典型变量V1=0.0661引体向上+0.0168仰卧起坐-0.0140跳跃次数。04典型载荷与交叉载荷分析生理U1与腰围负相关（-0.925），反映肥胖程度；运动V1与仰卧起坐、引体向上正相关（0.818、0.728），反映运动能力。交叉载荷显示腰围与V1负相关，表明腰围越大运动能力越差。05结论生理指标中的腰围是影响运动能力的关键因素，肥胖（腰围大）与引体向上、仰卧起坐能力呈显著负相关，运动对体形有重要影响。06实例分析：城市经济发展与空气质量关系2026/5/14数据介绍经济发展水平指标包含6项指标：地区生产总值（X1）、第二产业增加值（X2）、第三产业增加值（X3）、住户存款余额（X4）、社会消费品零售总额（X5）、医院数（X6），数据来源于2023年《中国统计年鉴》。空气质量指标包含7项指标：细颗粒物（PM2.5）年平均浓度（Y1）、可吸入颗粒物（PM10）年平均浓度（Y2）、二氧化硫年平均浓度（Y3）、一氧化碳日均值第95百分位浓度（Y4）、二氧化氮年平均浓度（Y5）、臭氧（O3）日最大8小时第90百分位浓度（Y6）、空气质量达到或好于二级的天数（Y7）。数据覆盖范围选取我国31个城市2022年的相关数据，涵盖北京、上海、广州等直辖市及省会城市，具有广泛的区域代表性。2026/5/14SPSS操作步骤

数据准备将经济发展水平指标（X1-X6）和空气质量指标（Y1-Y7）数据合并，确保数据格式正确，无缺失值。

菜单选择依次点选“Analyze”→“Correlate”→“CanonicalCorrelation”，打开典型相关分析对话框。

变量设置将X1-X6选入“Set1”下方的框中，将Y1-Y7选入“Set2”下方的框中。

选项配置点击“Options”，在“Display”下勾选“Pairwisecorrelations”“Loading”“Varianceproportions”“Coefficients”，点击“Continue”返回主界面，再点击“OK”运行分析。2026/5/14结果分析

相关系数矩阵经济发展水平各指标间相关性极显著；空气质量指标中Y3（二氧化硫）与部分指标相关性不显著，其余指标间多呈极显著相关，适合进行典型相关分析。

典型相关系数第一对典型变量相关系数为0.8314，在0.1334的显著性水平下显著；第二对相关系数为0.7127，后续典型相关系数逐渐减小，表明第一对典型变量解释了主要相关关系。

典型系数与载荷阵标准化典型系数显示经济发展水平第一典型变量（U1）与X1、X2、X3等呈负相关，反映总体经济水平；空气质量第一典型变量（V1）与各污染物指标相关性符号不同，需结合载荷阵分析。典型载荷阵和交叉载荷阵表明经济发展水平高的城市NO2、O3、PM2.5含量较高，PM10、SO2、CO含量及优良天数较低。结果分析经济与空气质量关系结论城市经济发展水平与空气质量存在一定相关关系，经济水平高的城市因汽车尾气、工业规模等因素，NO2、O3、PM2.5等污染物含量较高；经济水平较低的城市因燃煤、自然环境等因素，PM10、SO2、CO等污染物含量较高。07典型相关分析的应用与拓展2026/5/14应用领域

经济领域：经济发展与空气质量关系研究选取我国31个城市的经济发展水平指标（如地区生产总值、第二产业增加值等）与空气质量指标（如PM2.5浓度、空气质量达标天数等），通过典型相关分析揭示两者间的相关关系，发现经济发展水平较高的城市在部分污染物含量上呈现特定规律。

医学领域：生理指标与运动能力关系研究对20名中年男性的生理指标（体重、腰围、脉搏）和运动指标（引体向上、仰卧起坐、跳跃次数）进行典型相关分析，结果显示第一典型变量主要反映体形肥胖程度与运动能力，腰围大的人体形较胖，运动能力较差。

社会科学领域：多变量组间关系探讨可用于研究教育投入指标（如教育经费、教师数量）与学生成绩指标（如平均分、及格率）等不同变量组之间的整体线性相关关系，帮助分析教育资源配置对学生成绩的综合影响。2026/5/14软件实现对比R语言：功能齐全且灵活R语言可自行编写程序，能实现典型相关分析的各项功能，如计算典型相关系数、典型载荷、进行显著性检验等，适合有编程基础的用户，能满足个性化分析需求。SPSS：操作简便易上手通过菜单操作即可完成典型相关分析，无需编程，能输出相关系数矩阵、典型载荷、交叉载荷等结果，适合非专业编程人员快速进行分析，但灵活性相对较低。SAS：输出结果全面SAS软件在典型相关分析中输出的结果较为全面，包含多种统计量和检验结果，但操作相对复杂，需要一定的学习成本，在专业统计分析领域应用广泛。冗余分析补充冗余指数的计算步骤第一步计算共同方差的比例，通过平方每个因变量的典型载荷并简单平均得到；第二步计算解释的方差比例，即自变量典型变量与因变量典型变量间相关系数的平方；第三步将共同方差比例乘以典型相关系数平方得到冗余指数。冗余指数的意义冗余指数用于衡量一组变量的方差能被另一组变量的方差解释的比例，它像多元回归的R²统计量，为评估两组变量间的解释能力提供了综合测量指标，有助于更全面地理解典型相关分析结果。鼓励读者尝试实现冗余指数的计算在R语言中容易实现，读者可基于典型相关分析的结果，按照上述步骤编写代码进行计算，以深入探究变量组间的方差解释情况，完善分析结论。08总结与展望2026/5/14主要内容总结

典型相关分析的核心思想借鉴主成分分析降维思想，通过提取两组变量的典型相关变量（线性组合），使各对典型变量间相关程度最大且组内各典型变量互不相关，以典型相关系数衡量两组变量整体线性相关关系。

基本理论与方法要点包括总体与样本典型相关变量及系数的定义、模型求解（基于协方差矩阵特征根与特征向量）、典型相关系数显著性检验（巴特莱特X²检验），且要求变量为间隔尺度。

分析步骤与关键环节需依次完成确定目标、设计分析（样本量建议每个变量至少10个观测）、检验假设（线性关系、变量正态性）、推导函数与评价拟合（典型根、冗余指数）、解释变量（权重、载荷、交叉载荷）、验证模型（子样本分析等）。

实例应用与价值体现通过中年男性生理与运动指标、城市经济发展与空气质量的典型相关分析案例，展示了其在揭示两组变量复杂关联（如肥胖程度与运动能力负相关、经济水平与污染物含量关系）中的实际应用价值。2026/5/14未来发展方向与其他多元分析方法的融合探索与因子分析、回归分析、结构方程模型等结合，如将典型相关变量作为因子分析输入或纳入回归模型，以更全面挖掘变量间多层次关系，提升分析深度与解释力。大数据时代的应用拓展针对高维、海量数据，研究改进典型相关分析算法以提高计算效率，拓展其在基因组学（基因与疾病指标关联）、社交媒体（用户行为与内容特征关系）等大数据领域的应用。方法局限性的改进研究针对典型权重不稳定性、变量解释困难等局限，开发更稳健的估计方法（如正则化典型相关分析），引入变量旋转技术辅助解释，或提出新的冗余指数检验方法，增强方法可靠性与实用性。跨学科应用与理论创新推动典型相关分析在心理学（认知与行为变量关系）、环境科学（多污染物与影响因素关联）等更多学科的应用，同时结合非参数方法、机器学习理论，创新典型相关分析理论框架，适应复杂数据结构分析需求。谢谢观看THEEND定性数据的建模分析：对数线性模型与Logistic回归Logistic回归模型基础04对数线性模型的SPSS实现与案例分析03对数线性模型基础02引言与学习目标01目录CONTENTS总结与思考08Logistic回归的方法步骤与扩展07未分组数据的Logistic回归模型06分组数据的Logistic回归模型05目录CONTENTS01引言与学习目标定性数据与建模需求定性数据的特点定性数据是离散型数据，常以列联表形式呈现，描述分类变量间的关系，如收入水平与产品满意度的交叉分组数据。列联表分析的局限性列联表分析无法系统评价变量间联系，不能估计交互作用大小，难以深入揭示变量间复杂关系。建模方法的必要性对数线性模型可将方差分析方法应用于列联表，分析变量间交互效应；Logistic回归适用于因变量为二分类的情况，克服判别分析的严格假设，二者为定性数据建模提供有效工具。本章核心价值通过学习对数线性模型和Logistic回归，掌握定性数据的建模分析方法，提升对离散型数据关系的描述、解释和预测能力，为实际问题分析提供支持。2026/5/14本章学习目标01掌握对数线性模型的基本原理理解对数线性模型将概率取对数后分解为因素主效应和交互效应的原理，以及模型饱和与非饱和的判定条件。02掌握对数线性模型的建模方法学会基于列联表数据，通过参数估计和检验，分析变量间主效应和交互效应，如利用SPSS软件实现模型分析。03掌握Logistic回归分析结果的解释能够解释Logistic回归中的系数含义、优势比率，以及模型拟合优度检验结果，如通过wald统计量判断自变量显著性。04理解判别分析与Logistic回归的优缺点判别分析依赖多元正态性和等协方差阵假设，Logistic回归无此限制且结果稳定，对自变量类型要求宽松，更具实用性。05掌握通过SPSS软件实现Logistic回归熟练操作SPSS软件进行数据准备、模型参数设置、结果输出与解读，如完成分组或未分组数据的Logistic回归分析。02对数线性模型基础2026/5/14对数线性模型的基本原理

对数线性模型的定义对数线性模型是用于离散型数据或交叉列联表数据的统计分析工具，可将方差分析和线性模型方法应用于交叉列联表，以描述和分析定性变量间的关系。

模型的核心作用弥补列联表分析的不足，能系统评价变量间联系，估计变量间交互作用大小，是处理定性变量关系问题的有效方法。

概率对数变换与分解思想先将概率取对数，再进行分解处理，通过建立对数线性关系，将复杂的定性变量关系转化为可分析的线性模型结构。2×2维交叉列联表与模型结构2×2维频数表与概率表频数表展示交叉分类下的观测频数，如A因素与B因素不同水平组合的n11、n12等；概率表则对应各组合的概率πij，边际概率πi·、π·j及总和为1。对数线性模型公式推导通过对概率取对数并分解，得到模型ηij=Ai+Bj+(AB)ij，其中Ai对应lnπi·，Bj对应lnπ·j，(AB)ij对应ln(πij/(πi·π·j))。与双因素方差分析模型类比模型结构类似有交互效应的双因素方差分析，αi、βj为A、B因素的主效应，γij为交互效应，满足各效应和为零的约束条件。2026/5/14模型参数估计与效应判断

参数估计方法以交叉列联表频率为概率的无偏估计，将频率取对数后代入模型，通过公式计算交互效应γij的估计值，调整i、j取值上限可适用于非双水平因素。

交互效应γij的判断γij>0表明两因素对应水平存在正效应，γij<0为负效应；γij=0时因素相互独立。如例9.1中，高收入与满意的γij=-0.260为负效应。

饱和模型与非饱和模型若所有γij=0，模型为非饱和模型（因素独立）；否则为饱和模型（因素有交互效应）。饱和模型中实际频数与期望频数相同，残差为0。03对数线性模型的SPSS实现与案例分析2026/5/14数据准备与SPSS操作步骤

数据输入格式以3×2维交叉列联表（收入情况：高/中/低；满意情况：满意/不满意）为例，数据需按“频数-收入情况-满意情况”三列输入SPSS，如收入1（高）、满意1对应频数53，收入1、满意2对应频数38等，共6行数据。

激活频数操作通过Data→WeightCases菜单，选中“Weightcasesby”，将“频数”变量选入FrequencyVariable框，点击OK完成激活，确保分析时以频数为权重，否则结果会出错。

Loglinear模块参数设置依次选择Analyze→Loglinear→ModelSelection，将“收入情况”“满意情况”选入Factor(s)框，分别定义范围（收入1-3，满意1-2）；“频数”选入CellWeights框；Options中勾选SaturatedModel的ParameterEstimates，点击OK执行分析。2026/5/14输出结果解读（一）：模型信息与拟合优度

模型基本信息输出显示有效样本量792，分为6类（3×2交叉），收入情况3水平、满意情况2水平，分析效应包含主效应及交互效应（收入情况※满意情况）。

收敛情况系统经1次迭代后达到收敛，相邻两次估计之差≤0.001，表明模型迭代收敛速度快，结果稳定。

拟合优度检验饱和模型下，观测频数与期望频数完全一致（如高收入满意组均为53.5），残差及标准化残差均为0；似然比卡方与Pearson卡方值均为0，自由度0，提示模型完全拟合数据。2026/5/14输出结果解读（二）：效应检验与参数估计K-way效应检验

K-way及Higher-orderEffects检验显示，1维及以上效应（主效应+交互效应）的似然比卡方=662.843（Sig=0.000），2维交互效应卡方=22.087（Sig=0.000），均拒绝原假设，表明主效应与交互效应显著。参数估计结果

收入情况主效应：高收入=-0.683，中收入=0.883，低收入=-0.200；满意情况主效应：满意=0.425，不满意=-0.425；交互效应：高收入满意=-0.260，中收入满意=0.269，低收入满意=-0.009。效应方向与大小判断

参数值为正表示正效应，负为负效应。如中收入满意交互效应为正（0.269），表明中等收入对满意度有促进作用；高收入满意为负（-0.260），表明高收入对满意度有抑制作用。案例结论与应用启示案例分析结论顾客满意度整体偏向满意（β满意=0.425）；收入阶层满意度差异显著，中收入最高（α中收入=0.883），高收入最低（α高收入=-0.683）；交互效应显示中等收入对满意度有正效应，高、低收入有负效应，企业核心消费群体为中等收入者。模型应用价值对数线性模型可系统分析定性变量间的主效应与交互效应，克服列联表分析无法量化交互作用的局限，适用于市场调研、社会调查等领域，帮助识别关键影响因素及群体差异。实践注意事项应用时需确保数据格式正确并激活频数，根据研究目的选择饱和或非饱和模型；通过参数估计值的正负与大小可直观判断效应方向及强度，为决策提供量化依据。04Logistic回归模型基础2026/5/14Logit变换与模型构建直接处理概率的局限性因变量为二元变量时，概率π的取值范围为[0,1]，难以用线性模型描述其与自变量的关系；且当π接近0或1时，微小变化不易被普通方法捕捉。Logit变换的定义引入Logit变换Q=ln(π/(1-π))，将π从[0,1]区间映射到(-∞,+∞)，解决了线性模型的取值限制问题，对π在0或1附近的微小变化敏感。Logistic线性回归模型表达式假设Q与自变量线性相关，即ln(π/(1-π))=β₀+β₁x₁+...+βₖxₖ，反变换后得到Logistic回归模型π=1/(1+e^-(β₀+β₁x₁+...+βₖxₖ))，克服了线性模型的局限。Logistic回归与判别分析的比较假设条件对比判别分析依赖严格的多元正态性和等协方差阵假设，实际中常难以满足；Logistic回归无此假设，结果在假设不满足时仍稳定。适用场景差异判别分析适用于被解释变量为多组的情况；Logistic回归在被解释变量为两组时更受欢迎，且对自变量类型（度量/非度量）无限制。优缺点总结Logistic回归优势：无需严格假设、结果稳定、类似回归分析易解释、可包含非线性效果和诊断；判别分析在满足假设时分类精度可能较高，但假设限制严格。2026/5/14模型估计与拟合度评价

极大似然估计法通过迭代寻找使似然函数最大化的参数估计值，不同于多元回归的最小二乘法，适用于Logistic变换的非线性特征。

拟合度指标-2LL用-2倍似然比值的自然对数（-2LL）衡量模型拟合度，值越小拟合越好；完全拟合时-2LL=0，实际中通过比较不同模型的-2LL值判断改进效果。

Wald统计量的作用用于检验回归系数的显著性，Wald值越大（sig.值越小），系数越显著；不同于多元回归的t检验，专门适用于Logistic回归的参数检验。05分组数据的Logistic回归模型2026/5/14分组数据模型的改进与Logit变换

01回归模型的两方面改进针对0-1型因变量问题，回归函数改用限制在[0,1]区间内的连续曲线，如Logistic函数；因变量用yi等于1的比例代替本身作为因变量，因E(yi)=πi是自变量为xi时yi等于1的比例。

02Logistic函数形式与图形Logistic函数形式为特定表达式，其图形呈S形。图9-1和图9-2分别展示了不同Logistic函数的图形，体现其在[0,1]区间内单调递增的特征。

03购房意向数据的Logit变换计算以例9-2购房意向数据为例，计算实际购房比例πi=mi/ni，再进行Logit变换π9i=ln(πi/(1-πi))。如序号1，年家庭收入1.5（10万元），签订意向书人数25，实际购房人数8，πi=0.32，π9i=-0.75377。普通最小二乘估计与结果解读建立线性回归模型对例9-2数据，令π9i=β0+β1xi+εi，构建普通一元线性回归模型。以Logit变换后的π9i为因变量，年家庭收入x为自变量。经验回归方程与检验结果得到经验回归方程：9=-0.886+0.156x，判定系数r²=0.9243，显著性检验P值≈0，模型高度显著，表明年家庭收入与Logit变换后的购房比例线性关系强。购房比例预测应用利用模型预测，当x0=8（年家庭收入80万元）时，代入方程得π9=-0.886+0.156×8=0.362，再还原为购房比例，预计实际购房比例为59%，即签订意向书的该收入家庭购房概率为59%。2026/5/14加权最小二乘估计与异方差处理分组数据模型的异方差性分组数据的Logistic回归模型存在异方差性，式（9.7）的回归模型不是等方差的，当ni较大时，πi的近似方差为特定形式，需采用加权最小二乘估计处理。加权最小二乘估计的权数选择选取权数wi=niπi(1-πi)，其中ni为签订意向书人数，πi为实际购房比例。如例9-2中序号1，wi=25×0.32×(1-0.32)=5.440，通过权数调整消除异方差影响。SPSS加权回归操作及结果对比在SPSS中依次选择Analyze→Regression→Linear，将逻辑变换、年家庭收入、权重分别选入对应框。输出结果显示，加权后模型R²=0.881，调整后R²=0.864，对x0=8预测购房比例为58.5%，与未加权结果接近，提升了模型精度。06未分组数据的Logistic回归模型2026/5/14未分组数据模型的极大似然估计