用极大似然法估计因子载荷矩阵_因子分析的方法和公式 - Minitab
极大似然法估计在假定数据遵循多变量正态分布时的因子载荷。由其名称便可以看出,此方法通过最大化与多变量正态模型相关的似然函数找到因子载荷和唯一方差的估计值。这也可以通过最小化涉及残差方差的表达式达到同样的效果。算法会进行迭代,直到找到最小值或达到最大的指定迭代数(默认为 25)。
Minitab 使用基于 Joreskog1、2的算法通过一些调整加强收敛。我们在这里对该算法做了简要小结。
假定我们有 p 个变量并且想要使用 m 个因子拟合模型。设 R 作为变量的 p × p 相关矩阵,设 L 作为因子载荷的 p × m 矩阵,设 Ψ 作为对角元素是唯一方差的 p × p 对角矩阵 Ψi。然后,我们需要找到最大化似然函数 f(L,Ψ) 的 L 和 Ψ 值。这涉及两个步骤,首先为 Ψ 找一个值,然后为 L 找一个值。
您可以间接地指定 Ψ 的初始值。在“因子分析 - 选项”子对话框中,输入包含使用初始公因子方差估计于中公因子方差初始值的列。Minitab 之后将 Ψ 对角元素计算为(1 − 公因子方差)。
对于 Ψ 的固定值,我们要相对于 L 最大化 f(L,Ψ)。这是一种简单的矩阵计算。之后将 L 的值代入 f(L,Ψ)。现在 f 可视为 Ψ 的函数。此函数的简单变换提供
其中,λ1 < λ2 < ... λp 是 Ψ R- 1Ψ 的特征值。我们之后通过 Newton-Raphson 过程最小化 g(Ψ)。这得出 Ψ 的估计值,该值之后代入似然 f(L,Ψ)。然后,似然又一次相对于 L 进行最大化。计算出 g(Ψ) 的新值,以此类推。默认情况下,如果未达到收敛,迭代会继续多达 25 个步长。如果算法在 25 个步长中未达到收敛,您可能要在“选项”子对话框中更改默认的最大迭代次数。
如果下面任一项为真,则在第 n 个步长达到收敛:
函数 g(Ψ) 在连续性步长之间不会发生太大变化。具体地说,如果:
| [第 n 个步长的 g(Ψ)] − [第 (n − 1) 个步长的 g(Ψ)] | < 10-6
所有唯一方差在连续性步长之间都不会发生太大变化。具体地说,如果:
| ln(第 n 个步长的 Ψi)− ln(第 n − 1 个步长的 Ψi) | < K2,
对于所有 i = 1, ... , p,其中 Ψ 的第 i 个对角元素 Ψi 是对应于变量 i 的唯一方差。
K2 的值可在“选项”子对话框中的收敛中指定。默认情况下,该值为 0.005。
在“结果”子对话框中选择全部和最大似然提取迭代可显示每次迭代的信息。会显示目标函数 g(Ψ) 的值,之后在 ln(Ψi) 中发生最大变化。如果在一次迭代中,g(Ψ) 的值未减少,那么将采用较小(一半尺寸)的步长。继续半步直到 g(Ψ) 减少或采用 25 个半步。此时显示半步的数量。如果 g(Ψ) 在 25 个半步中未减少,则算法停止,并显示一条消息。
第二个导数的矩阵用于 g(Ψ) 的最小化。该矩阵并非总为正定。如果不是,则使用近似。当 Minitab 使用准确的矩阵时,结果中会出现星号。
最小化函数 g(Ψ) 时,可能发现 Ψ 的对角元素的值为 0 或为负。为防止出现这种情况,Minitab 的算法将 Ψ 的对角元素限制在 0 之外。具体地说,如果唯一方差 Ψi 小于 K2,则其设定为等于 K2。K2 是在“选项”子对话框中的“收敛”中设定的值。
当算法收敛时,会对唯一方差执行最终检查。如果任何一个唯一方差小于 K2,则它们设定为等于 0。对应的公因子方差则等于 1。此结果称为 Heywood 案例,并且 Minitab 会显示一条消息告知用户此结果。优化算法,比如用于极大似然因子分析的算法,可提供不同的答案,在输入中略有变化。例如,如果您更改了一些数据点、在使用初始公因子方差估计于中更改了初始值,或者在收敛中更改了收敛标准,您可能会发现因子分析结果出现差异。这对于解出现在极大似然曲面上相对平坦区域的情况尤为如此。
用极大似然法估计因子载荷矩阵_因子分析的方法和公式 - Minitab相关推荐
- 用极大似然法估计因子载荷矩阵_第7章-因子分析.ppt
7.3 因子正交旋转 在第7.1节我们已经看到,满足方差结构Σ =AAT+Ф的因子模型并不惟一,模型的公因子与载荷矩阵不惟一.如果F是模型的公因子,A是相应的载荷矩阵,而T是m×m正交矩阵,则F*=T ...
- 用极大似然法估计因子载荷矩阵_应用统计学与R语言实现学习笔记(十三)——因子分析...
Chapter 13 Factor Analysis 本篇是第十三章,内容是因子分析. 这篇博客的完整内容包含各类数学表达.可以见我CSDN和hexo搭的个人博客. CSDN博客 1 因子分析概念 因 ...
- matlab求因子载荷矩阵,Matlab因子分析
% 从相关系数矩阵出发进行因子分析 %***************************定义相关系数矩阵PHO***************************** PHO = [1 ...
- 因子分析 factor analysis (三) : 因子载荷矩阵的估计方法
因子分析系列博文: 因子分析 factor analysis (一 ):模型的理论推导 因子分析 factor analysis (二 ) : 因子分析模型 因子分析 factor analysis ...
- python主成分得分系数矩阵_用spss得出的成分矩阵就是初始因子载荷矩阵吗?成分得分系数矩阵又是什么...
用spss得出的成分矩阵就是初始因子载荷矩阵吗?成分得分系数矩阵又是什么呢? 第一问是的,译法不同而已:第二问,成分得分系数矩阵是用来求成分得分的,用标准化后的原始数据矩阵乘以成分得分系数矩阵就可以得 ...
- R语言使用psych包的fa函数对指定数据集进行因子分析(输入数据为相关性矩阵)、使用rotate参数指定进行斜交旋转提取因子、编写自定义函数通过因子模式矩阵与因子相关性矩阵相乘计算因子载荷矩阵
R语言使用psych包的fa函数对指定数据集进行因子分析(输入数据为相关性矩阵).使用rotate参数指定进行斜交旋转提取因子.编写自定义函数通过因子模式矩阵与因子相关性矩阵相乘计算因子载荷矩阵 目录
- 主成分与因子分析异同_因子分析
1 导言 因子分析是将多个实测变量简化为较少变量的方法. 利用因子分析能简化分析,更能反映事物本质. 学习因子分析能确定哪些变量应该保留,哪些应该剔除,并得到主要成分的表达式,这些表达式是原有变量的线 ...
- 主成分分析法案例_因子分析案例及操作解析
分析一个省的科技创新能力受哪些潜在因素的影响?(本数据来源于网络收集,x1~x15代表各省市相关经济指标) 第一步:数据输入(如下图) 第二步:操作步骤:分析--降维--因子分析(F). 得出下图:( ...
- 因子分析累计贡献率_累计方差贡献率_spss累计方差贡献率_因子分析方差贡献率...
spss中因子方差贡献率–怎样用SPSS求方差贡献率和方差累计贡献率 analyze 下面选择 diemnsion reducation ,再选择factor 因子分析,把你需要计算方差贡献率的变量放 ...
- matlab因子载荷矩阵正交旋转,因素分析中的矩阵旋转
因素分析中的矩阵旋转 因素分析法 因素分析是一种统计技术,目的是从众多的可观测的"变量"中,概括和推论少数"因素".用最少数的"因素"来概括 ...
最新文章
- 符合RESTful的接口规范
- r语言提取列名_R语言基础(2)向量|矩阵|数组|数据框|数据IO|实用函数
- 【Paper】2021_Robust Near-Optimal Coordination in Uncertain Multiagent Networks With Motion Const
- SQL SERVER 子查询的用法
- 盘点Chrome 灵魂插件!爱了爱了!
- sqlserver快速查找所有存储过程中是否包含某字符
- Dubbo简单介绍及实例
- ffmpeg编译(生成Windows或Win32平台dll, lib)
- Java使用递归实现全排列的代码
- 【重磅】关于本人CSDN的博客答疑
- com app.html,downloadApp.html
- 淘宝美工设计师细说何为天猫透明背景
- android禁止录屏后键盘,怎样取消华为按键录屏功能 | 手游网游页游攻略大全
- mysql 空格键和回车键查询
- pycharm 学生资格过期,提示No suitable licenses associated with account解决方案
- AI笔记: 数学基础之方向导数的计算和梯度
- NeuroAI作者,冷泉港主任独家专访:忽视神经科学给我们的启示是愚蠢的
- (寒假开黑gym)2018 ACM-ICPC, Syrian Collegiate Programming Contest(爽题)
- 网页回拨-Web CallBack
- 万年历(C语言代码实现)
热门文章
- xps数据怎么导出为txt_如何处理XPS原始数据
- c语言中cmp的用法,cMP+cplay值得一试!(附下载地址和简要使用说明)
- FPGA实现FIR滤波器
- 嵌入式从入门到精通的学习攻略
- 计算机电子设计论文,计算机毕业设计电子信息毕业设计电子信息毕业论文
- 【数据结构 严蔚敏版】 链式栈基本操作
- Finetune时更改tokenizer词表
- 实践与交流:采取“硬刷”方式完美激活Windows 7旗舰版的具体方法步骤(图文详解教程)...
- Android+FFmpeg音视频学习笔记
- 软件测试安装包,安装包的测试