因子分析——建立载荷矩阵

到这里已经学了好多的多元分析方法了,有聚类分析法,有主成分分析法,尤其是主成分分析法,为什么还要讨论因子分析法呢?很多地方都有对主成分分析法和因子分析法的区别比较,这里就不多说了,只记录一下最重要的地方。

主成分分析法:是对原始变量的线性组合,且相互垂直。

因子分析法:研究众多变量之间的内部依赖关系,潜在的假想变量+随机变量的线性组合。

因子载荷,反映了第 i 个变量与第 j 个公共因子的相关系数,反映了第 i 个变量与第 j 个公共因子的相关重要性。绝对值越大,相关的密切程度越高。

因子载荷矩阵的估计方法:

  • 主成分分析法

  • 主因子法

  • 极大似然估计

这里就只用其中一种方法解决吧,好像大部分都用的是主成分分析法。

举个栗子吧,P243,公司股票有5个因素影响,探讨一下这5个因素的内在联系。

clc,clear;
​
r=[1.000 0.577 0.509 0.387 0.4620.577 1.000 0.599 0.389 0.3220.509 0.599 1.000 0.436 0.4260.387 0.389 0.436 1.000 0.5230.462 0.322 0.426 0.523 1.000];
​
% 利用关系矩阵进行主成分分解
[vec1,val,rate] = pcacov(r);
​
f1 = repmat(sign(sum(vec1)),size(vec1,1),1);
vec2 = vec1.*f1;        % 使得每个特征向量分量和为正
​
% 因子分析公式
f2 = repmat(sqrt(val)',size(vec2,1),1);
​
% 全部的因子的载荷矩阵
a = vec2.*f2;
​
% 提取第一个因子的载荷矩阵
a1 = a(:,1);
% 计算第一个因子的特殊方差
tcha1 = diag(r-a1*a1');
​
% 提取两个因子的载荷矩阵
a2 = a(:,[1,2]);
% 计算两个因子的特殊方差
tcha2 = diag(r-a2*a2');
% 求两个因子时的残差矩阵
ccha2 = r - a2*a2' - diag(tcha2);
​
% 求累积贡献率
con = cumsum(rate);

之前说主成分分析已经是很复杂的多元分析方式了,看来这句话说早咯!!!

首先用主成分分析​对相关系数矩阵操作,求出特征向量,特征值,和主成分贡献率。

然后变形使得每个特征向量的分量和为正。

利用特征值和载荷矩阵公式求出载荷矩阵:

载荷矩阵求出来以后:

然后依然进行需要,选择因子的个数。

比如选择一个 ​ F1 公共因子。

那么就可以表示5个变量(5种股票)的内在关系:

并对此进行分析:

特殊方差是什么呢? 其实可以理解为从原来的变量到公共因子空间的转换效果的一个评价。

公式:

即:

转载于:https://www.cnblogs.com/TreeDream/p/8337748.html

因子分析——建立载荷矩阵相关推荐

  1. 一篇文章带你搞定数学建模中的载荷矩阵、相关系数矩阵、主成分分析(11年土壤重金属污染示例讲解含代码)

    文章目录 一.题目分析 二.基于主成分分析法的重金属污染评价模型 1. 模型建立 2. 模型求解 三.问题求解代码 四.相关系数矩阵的了解 五.载荷矩阵的了解 一.题目分析 通过数据分析,说明重金属污 ...

  2. 如何在R语言中建立六边形矩阵热图heatmap可视化

    原文链接:http://tecdat.cn/?p=18879 这是一个六边形热图可视化程序,主要用到的知识RColorBrewer,fields,也就是R中的可视化绘图库(点击文末"阅读原文 ...

  3. matlab因子载荷矩阵正交旋转,如何对载荷矩阵进行旋转

    用stata做主成分分析,只把累计贡献率算出来了,如下,看起来需要有5个主成分.下一步就不知道怎么做了,特别是"对载荷矩阵进行旋转"得到各变量的特征向量. . screeplot ...

  4. 2021.03.03循环建立储存矩阵,或者是操作矩阵像素

    C++ 记录一下用循环建立储存矩阵,或者是操作矩阵 为了之后的矩阵操作方便,数据类型都用为double,CV_64FC1 Mat x, y; meshgrid(Range(-s, s), Range( ...

  5. 从生态繁荣到市场猛增:财报背后,李彦宏如何建立DuerOS矩阵

    有人说成功的商业项目像风筝,一飞冲天:但事实上不回头的风筝往往会断线,真正好的商业项目,应该像座码头--有进有出,四通八达. 一个商业体的成功,并不是仅仅以增长速度来判断,更重要的是要看这个商业体的产 ...

  6. matlab建立荧光矩阵,基于Matlab的溶解有机物的三维荧光矩阵构建

    基于Matlab 的溶解有机物的三维荧光矩阵构建 钱伟1,2,3,杨柳明1,2,3,王宁3,张莉3 [摘要]摘要:三维荧光光谱(EEM)结合平行因子分析(PARAFAC)能够根据水体中有机物的荧光性质 ...

  7. python建立空矩阵_SciPy创建稀疏矩阵

    3. SciPy创建稀疏矩阵 严格意义上讲ndarray数据类型应属数组而非矩阵,而matrix才是矩阵,这个在NumPy创建matrix一章里有讲述,是最基本的矩阵matrix创建方法,忘记了可以回 ...

  8. python建立空矩阵_创建空矩阵Python

    首先,您应该在最里面的列表中插入一些内容(比如None).其次,当您在最外层列表中使用乘法时,它会将引用复制到内部列表,因此当您更改一个元素时,您也会在所有其他列表中更改此元素:>> pa ...

  9. SPSS因子分析中成分矩阵、旋转成分矩阵和成分得分系数矩阵的区别

最新文章

  1. 国产AI框架再进化!百度Paddle Lite发布:率先支持华为NPU在线编译,全新架构更多硬件支持...
  2. LeetCode 503 Next Greater Element II(stack)
  3. 018-继承-OC笔记
  4. chrome 63 android分类,Chrome 63 Beta新特性介绍
  5. 怎么学ui设计入门?怎样才能少走弯路?
  6. Python数据分析与展示[第二周]
  7. mysql索引操作(转)
  8. 在Mac下连接阿里云服务器
  9. python eval函数简单使用
  10. protel DXP2004批量修改
  11. zblog如何调用HTML,Zblog教程:模板文件与调用标签详细说明
  12. 可靠传输协议 rdt 1.0、rdt 2.0、rdt 2.1、rdt 2.2、rdt3.0
  13. GDI+绘图基础,绘制表格
  14. 【巴比特:区块链是什么】笔记
  15. SQL 给时间字段增加一固定时间
  16. 汇顶2018年 IC校招笔试题目
  17. 面向对象(文字格斗进阶版)
  18. RecyclerView 配合 DiffUtil,RecyclerView局部刷新
  19. error: void value not ignored as it ought to be
  20. 我所认为的产品经理——从读PONY的产品笔记想到的

热门文章

  1. STM32通信接口RS485
  2. 检查域名是否被注册过或被K过
  3. 白嫖彩虹云主机自制图床
  4. 转 3GP文件格式分析
  5. hibernate的sessionfactory
  6. 现在3d建模师的工资高么
  7. msq_table's methods
  8. 我再copy回来。中海真是有心人。
  9. xcode断点看不到调试信息 变量值 字典值
  10. 基于Qt平台利用UDP协议传输大文件(百度网盘免费代码)