[拼音]：duobianliang tongji fenxi

[外文]：multivariate statistical analysis

统计资料中有多个变量（或称因素、指标）同时存在时的统计分析，是统计学的重要分支，是单变量统计的发展。例如对630名炊事员高血压病进行调查，检查项目中除血压外，尚有年龄、性别、体重、体胖等15个项目（变量）。如果用单变量统计分析法考察超重与血压的关系，一般是把数据做成表1的形式。从表1可见，超重组与不超重组相比，高血压患病率高出一倍以上。但如果把资料按体胖者与不体胖者划分成两组，再考察每组内的超重与高血压患病率的关系，就未能发现超重与高血压患病率有任何明显的联系。也就是说，单变量统计分析忽视了另外因素（如此例中的体胖及年龄等）的影响。对于有多个变量客观存在而又相互影响的资料，采用简单的单变量统计分析是不合理的。多变量统计分析就能把变量间的内在联系和相互影响考虑在内。

统计学中的多变量统计分析起源于医学和心理学。1930年代它在理论上发展很快，但由于计算复杂，实际应用很少。1970年代以来由于计算机的蓬勃发展和普及，多变量统计分析已渗入到几乎所有的学科。到80年代后期，计算机软件包已很普遍，使用也方便，因此多变量分析方法也更为普及。

多变量统计的理论基础和工具是数学中的概率论和矩阵。但对于实际应用者而言，只要有合适的计算机和软件包以及掌握一些初步的多变量统计知识就可以使用它来解决实际问题。多变量统计的内容很多，但从实际应用角度看，主要包括回归分析、判别分析、因子分析、主成分分析、聚类分析、生存分析等六个大的分支。

当多个变量x₁，x₂，…，x_m（称为回归变量或自变量、独立变量）同时影响某个指标 y（称为因变量或依赖变量）时，可进行回归分析，回归分析的第一个任务就是求回归变量对指标 y的影响的统计规律性（也称回归关系）；第二个任务是寻找众多的回归变量中哪一些能对指标 y产生影响（常称为因素分析或变量的筛选）；第三个任务(也称相关分析)是在固定（或称消除）其他变量的影响后，考察每一个回归变量对指标 y的相关程度（称为偏相关系数）。上述三个任务常是相互联系，可以同时完成。

回归变量x₁，x₂，…，x_m与因变量y之间最常见的统计关系有两大类型：线性模型和非线性模型。线性模型中假定y 的主要部分（记为），可由x₁，x₂，x_m线性表示为

其中b₀，b₁，b₂，…，b_m是未知常数，需用样本去估计，ε 是用取代y后的误差。这是最常用的模型，称为多重线性回归或多元线性回归。用样本估计线性回归模型中未知常数的方法也很多，经典的方法为最小二乘法，它的理论较为完善，此法较适用于回归变量之间的相关性不很大时。其他求未知常数b₀，b₁，b₂，…，b_m的方法还有岭回归、特征根回归、主成分回归等，它们常用于回归变量之间相关性很大时。

非线性回归模型中y的主要部分与x₁，x₂，…，x_m的关系为非线性函数：

其中┃ 的形式已知，未知常数 α₁，α₂；…用样本去估计。医学中最常见的非线性回归是logistic回归，它常用于疾病对照研究以及生长发育问题中。

在前述的炊事员高血压病调查中，使用线性模型和最小二乘法求出未知常数，再用逐步回归选取变量，可求得15个变量中有 7个变量对炊事员舒张压有显著的影响，它们按偏相关系数大小排列为：年龄 (0.297)，体胖程度 (0.253)，肾炎史(0.162)，性别(0.117)，工作类别(0.081)，高血压家族史(0.061)，嗜咸程度(0.052)。从相关性大小看，体胖对舒张压的影响与年龄的影响大体相当。另外还可看出：工种，家族史和嗜咸对舒张压虽有影响，但影响甚小。

根据样本的某些指标来决定样本归属的类别。例如在医疗诊断中，要确定一个病人是否患有急性阑尾炎，这就是一个判别问题。为了回答这一问题往往需要对病人进行多项指标（变量）的检测，然后根据各项指标的观测值将病人归入患有急性阑尾炎或不患有急性阑尾炎的类中。判别分析通常是先确立一个判别函数，将各指标的观测值代入相应的变量，再根据某判别规则（如函数值大于某值）作出判断、鉴别或决策。例如，为了研究亚硝酸基盐化合物与胃癌的关系，有人曾对胃癌组（记为H₁）、萎缩性胃炎组(H₂)、浅表性胃炎组(H₃)3组病人测量以下6个指标(变量)：性别（x₁，男为1，女为0）、年龄(x₂)、胃液中pH值(x₃)、唾液中亚硝酸盐浓度(x₄)、胃液中亚硝酸盐浓度(x₅)、胃液中二甲基亚硝酸胺浓度 (x₆)。用判别分析法，可求出6个指标（变量）在 3个疾病组中分布有显著不同的是x₁，x₂，x₄，x₆；其余两个指标在不同组中的分布大体相同。对应于每个疾病组可建立以下的判别函数：u₁=-11.48+2.68x₁+0.37x₂+0.04x₄+0.90x₆ (H₁)u₂=-14.06+3.79x₁+0.35x₂+0.50x₄+1.82x₆ (H₂)u₃=-6.36+1.84x₁+0.27x₂+0.34x₄+0.84x₆ (H₃)

在判别分析时，可将测得的病例值(x₁，x₂，x₄，x₆)代入判别函数，求得一组函数值u₁，u₂，u₃。这里的判别规则是：如果u₁最大，则病例判属疾病组H₁；如果u₂最大，则判属H₂；如果u₃最大，则属H₃。这样，诊断就变成了数据的处理及分析，现代化医院自动诊断的原理就基于此。通常说的把医生的经验和知识存入计算机，也就是在计算机中建立诊断的经验方式──判别函数。判别函数中变量前的系数含有重要的信息。上列中变量x₃，x₅前的系数都为0；x₁前的3个系数(2.68，3.79，1.84)说明相对于女性(x₁=0)而言，男性 (x₁=1)更容易得萎缩性胃炎(3.79)或胃癌(2.68)； x₂前的3个系数说明相同年龄者得胃癌、萎缩性胃炎、浅表性胃炎的比为0.37:0.35:0.27；等等。

也称因素分析。医学、生物学及一切社会和自然现象中各变量（或事物）之间常存在有相关性或相似性。这是因为变量（或事物）之间往往存在有共性因素（称为公因子或共性因子），这些共性因子同时影响不同的变量（或事物）。因子分析的根本任务就是从众多的变量（或事物）中由表及里找出隐含于它们内部的公因子，指出公因子的主要特点，并用由实际测量到的变量（或事物）构造公因子。因子分析有R型及Q型之分，用于变量之间时称为 R型因子分析，用于事物之间时称为Q型。

以R型因子分析为例，设样本中的变量为x₁，x₂，…，x_m，隐藏的公因子为┃₁，┃₂，…，┃_k。这时每个变量在理论上常可写成下面的形式：

上式右边的第一部分是变量中公因子(┃₁，┃₂，…，┃_k)起作用的部分，后一部分是与公因子无关的部分（称为独立性部分）。因子分析的根本任务就是用样本求出┃₁，┃₂，…，┃_k及其系数｛α_嗎｝，系数α_嗎 称为权或负荷系数，当样本是标准化数据且假定公因子之间彼此不相关时，则权α_嗎就是公因子┃j与变量x_i间的相关系数。利用因子分析方法可以从所观测到的变量中推断出少数因子，用最少的因子来解释所观测到的变量，从而揭示事物之间内在的联系。对因子的实际解释必须结合专业知识并由实践检验。例如中国学者梁月华、孙尚拱曾用因子分析法找出隐含在6个易测量的生理指标（收缩压、 舒张压、呼吸、心率、体温及唾液量）内部的公因子┃₁，并用实验判定┃₁可很好地代表交感神经的平衡状态，最后用 ┃₁ 论证了中医的“寒热”其本质就是交感神经的抑制或兴奋。

是研究如何把彼此相关的变量综合成一个（或少数几个）综合指标（或称主成分），而该综合指标应能最大程度地反映观测变量所提供的信息。如记(x₁，x₂，…，x_m)为观测变量，欲求的综合指标Z一般可写成

实际上Z往往只能吸收m个变量中相关最大的一部分信息（类似于因子分析中┃₁），此当观测变量间彼此很少有相关性时，使用主成分分析是不合适的。如果观测变量间相关性可以分成几组而各组间又很少相关，这时不能用一个主成分综合全体变量，而应多取几个主成分。

实际使用时，由于主成分分析与因子分析极为相似，所以不少统计学家常把两种分析不加区别，名称也相互套用。

主成分分析在医学研究中有很多应用，例如有人把5个易测量的老化征 （白斑、老年斑、闭目单腿直立时间、老年环、脱齿数）综合成一个指标Z，计算表明综合值Z可以吸收5个老化征全体信息中的43％，能综合地反映出形体老化的程度。

也称为分类学。经典的分类学诞生于几百年前，比如化石分类、植物标本分类等。过去的分类多依靠一些特异性指标。如果对于所需分类的事物，不存在或难以使用特异性指标时就只能采用多变量统计分析法。把数学方法引进分类学并称之为“聚类分析”是60年代的事。此后聚类分析发展很快，并取得广泛应用，但还不大成熟。

聚类分析也可分为 R型及Q型，对变量作分类称R型，对样品（观察单元、事物）作分类称 Q型。分类的基础是相似性或距离。如果两个变量（或样品）彼此相似或距离很近，自然就分在同一类。因此在进行聚类分析时必须先定义相似性或距离。相似性或距离的定义法种类繁多。例如，常用变量间的相关系数代表变量间的相似性，以几何中两点间的欧氏距离（先要去量纲）代表两个样品间的距离。然后选用分类的数学公式，对它们的分类作出判别。这些公式也是种类繁多。至今没有一种公式是最优的。实际工作者常选用多种方法试算，再结合专业知识确定分类的结果。

生存分析起源于寿命表。生物的生存时间除了受健康的影响外，同时还受社会因素，生活条件等影响。生存分析研究哪些因素对“寿命”有显著影响，它的风险程度如何。20世纪末生存分析已不仅用于研究人的寿命问题，还用于一切广义的“寿命”或有关“死亡”的问题，比如发动机的寿命，病人手术后的生存时间，两种疗效的对比分析等。生存分析有多种模型，最常用的有Cox回归模型，它的特点是：m个变量联合作用的相对风险可以表示成每个变量单独作用时相对风险的乘积（故也称为乘法模型）。另外常用的模型为可加性模型，它的特点是：m 个变量联合作用的相对风险可表示为每个变量单独作用之和。究竟应使用什么样的模型应在具体问题中结合专业知识确定。

多变量统计分析除了上述六个大的分支外，通径分析和典则相关分析也很常用。一般回归分析只能计算每一个变量（在固定其他变量时）对指标 y的直接作用大小，而通径分析可同时计算每一个变量对指标 y的间接作用（即通过与它相关的变量作用于 y）。通径分析在流行病的遗传研究中已有不少应用。典则相关分析也是回归分析的进一步发展。对每个事物同时测量多个指标(y₁，y₂，…)和多个自变量(x₁，x₂，…)，分析指标的综合与自变量的综合是如何相关时多使用典则相关分析。

参考书目

1. 孙尚拱：《实用多变量统计方法与计算程序》，北京医科大学与中国协和医科大学联合出版社，北京，1990。