因子分析是用少数几个公共因子去提取研究对象的绝大部分信息,既减少了因子的数目,又把握住了研究对象的相互关系.在因子分析中根据研究对象的不同,分为R 型和Q 型,如果研究变量的相互关系时则采用R 型因子分析;如果研究样品间相互关系时则采用Q 型因子分析.但无论是R 型或Q 型都未能很好地揭示变量和样品间的双重关系,另一方面当样品容量n 很大(如n>1 000),进行Q型因子分析时,计算n 阶方阵的特......
2023-11-18
将原始数据矩阵A 变换为B 矩阵后,记SR=BT B,SQ=BBT,SR 和SQ 这两个矩阵存在明显的简单的对应关系,而且将原始数据aij 变换为bij 后,bij 关于i,j 是对等的,即bij 对变量和样品是对等的.
为了进一步研究R 型与Q 型因子分析,我们利用矩阵代数的一些结论.
引理11.2.2 设SR=BT B,SQ=BBT,则SR 和SQ 的非零特征值相同.
引理11.2.3 若v是BT B 相应于特征值λ 的特征向量,则u=Bv是BBT 相应于特征值λ 的特征向量.
定义11.2.1 (矩阵的奇异值分解)设B 为n×p 矩阵,且
rank(B)=m ≤min(n-1,p-1),
BT B 的非零特征值为λ1 ≥λ2 ≥…≥λm>0,令,则称di 为B 的奇异值.
如果存在分解式(www.chuimin.cn)
其中,U 为n×n 正交矩阵,V 为p×p 正交矩阵,Λ=,这里Λm=diag(d1,d2,…,dm),则称分解式B=UΛVT 为矩阵B 的奇异值分解.
记U=(U1⋮U2),V=(V1⋮V2),Λm=diag(d1,d2,…,dm),其中U1 为m×n 的列正交矩阵,V1 为p×m 的列正交矩阵,则奇异值分解式(11.2.4)等价于
引理11.2.4 任意非零矩阵B 的奇异值分解必存在.
引理11.2.4的证明就是具体求出矩阵B 的奇异值分解式(高惠璇,统计计算(1995)).从证明过程中可以看出:列正交矩阵V1 的m 个列向量分别是BT B 的非零征值为λ1,λ2,…,λm 对应的特征向量;而列正交矩阵U1 的m 个列向量分别是BBT 的非零征值为λ1,λ2,…,λm 对应的特征向量,且.
矩阵代数的这几个结论为我们建立了因子分析中R 型与Q 型的关系.借助以上引理11.2.2和引理11.2.3,我们从R 型因子分析出发可以直接得到Q 型因子分析的结果.
由于SR 和SQ 有相同的非零特征值,而这些非零特征值又表示各个公共因子所提供的方差,因此变量空间Rp 中的第一公共因子、第二公共因子…,直到第m个公共因子,它们与样本空间Rp 中对应的各个公共因子在总方差中所占的百分比全部相同.
从几何的意义上看,即Rp 中诸样品点与Rp 中各因子轴的距离平方和,以及Rp 中诸变量点与Rp 中相对应的各因子轴的距离平方和是完全相同的.因此可以把变量点和样品点同时反映在同一因子轴所确定的平面上(即取同一个坐标系),根据接近程度,可以对变量点和样品点同时考虑进行分类.
有关应用多元统计分析:基于R的实验的文章
因子分析是用少数几个公共因子去提取研究对象的绝大部分信息,既减少了因子的数目,又把握住了研究对象的相互关系.在因子分析中根据研究对象的不同,分为R 型和Q 型,如果研究变量的相互关系时则采用R 型因子分析;如果研究样品间相互关系时则采用Q 型因子分析.但无论是R 型或Q 型都未能很好地揭示变量和样品间的双重关系,另一方面当样品容量n 很大(如n>1 000),进行Q型因子分析时,计算n 阶方阵的特......
2023-11-18
对应分析的具体计算步骤如下:由原始数据矩阵A 出发计算对应矩阵P 和对应变换后的新数据矩阵B,计算公式见式和式.计算行轮廓分布,记R 矩阵由A 矩阵的每一行除以行和得到,其目的在于消除行点出现“概率”不同的影响.记N={Ri,i=1,2,…,m)是B 的奇异值.式给出Q 的分解式,第i个因子(i=1,2,…......
2023-11-18
对于广义线性模型,除了以上介绍的Logistic回归模型外,还有其他的模型,如Poisson模型,这里就不作详细介绍.以下简要介绍R 软件中“glm( )”关于这些模型的使用方法.Poisson分布族模型和拟Poisson分布族模型的使用方法如下:fm<-glmfm<-glm其直观意义是ln[E]=β0+β1x1+β2x2+…+βp xp,即E=exp(β0+β1 x1+β2 x2+…......
2023-11-18
,xn 是从一元总体中抽取的.一元数据的数字特征主要有:均值,方差s2=,标准差s=,等等.对于多元数据,除分析各分量的取值特征外,还要分析各分量之间的相关关系.由于多元统计分析中的符号多而杂,因此需要说明:在一元统计学中一般用大写和小写字母分别来区分随机变量及其观测值,在本书后面的章节里,由于其他复杂的符号,我们可能不再遵守此约定,请读者注意一个符号在每一章中的意义.......
2023-11-18
以下我们来讨论两个总体的距离判别,分别讨论两个总体协方差矩阵相同和不同的情况.设总体X1 和X2 的均值向量分别为μ1 和μ2,协方差矩阵分别为Σ1 和Σ2.给定一个样本x,要判断x 来自哪个总体.首先考虑两个总体X1 和X2 的协方差矩阵相同的情况,即μ1 ≠μ2,Σ1=Σ2=Σ.要判断x 来自哪个总体,需要计算x 到总体X1 和X2 的马氏距离的平方d2(x,X1)和d2(x,X2),然后进行......
2023-11-18
,F5 的线性组合表示出来Xi=μi+ai1 F1+ai2 F2+…,F5 的值却是未知的,有关参数的意义也有很大的差异.因子分析的首要任务就是估计因子载荷aij 和方差σ2i,然后给因子Fi 一个合理的解释,若难以进行合理的解释,则需要进一步作因子旋转,希望旋转后能发现比较合理的解释.特别需要说明的是这里的因子和试验设计里的因子(或因素)是不同的,它比较抽象和概括,往往是不可以单独测量的.......
2023-11-18
多元统计分析是统计学中应用性很强的一个分支,它的应用范围十分广泛.多元统计分析可以应用于几乎所有的领域,主要包括经济学、农业、地质学、医学、工业、气象学、金融、精算、物理学、地理学、军事科学、文学、法律、环境科学、考古学、体育科学、遗传学、教育学、生物学、管理科学、水文学等,还有一些交叉学科或方向等.多元统计分析的应用实在是难以一一罗列,以下简要地介绍一下多元统计分析在文学、数据挖掘(作为交叉学科......
2023-11-18
那就是这些被选的主成分所代表的主轴的长度之和占了主轴长度总和的大部分.有些文献建议,所选的主轴总长度占所有主轴长度之和的大约80%即可.其实,这只是一个大体的说法;具体选几个,要看实际情况而定.但如果所有涉及的变量都不那么相关,就很难降维.不相关的变量就只有自己代表自己了.......
2023-11-18
相关推荐