偏相关分析MATLAB代码实现
偏相关分析
概念:偏相关分析,用于分析当两个变量都与第三个变量相关时,将第三个变量的影响剔除,只分析另外两个变量之间的相关程度。
准备:首先获得每个变量之间的简单相关系数,即相关性分析(协方差/标准差之积)得到的系数,数学公式如下:
变量x的标准差:
变量y的标准差:
变量x、y的协方差:
最后得出变量x、y的相关系数:
注意,如果是要进行偏相关分析,需要所有变量之间的相关系数(即一阶相关性系数),可通过二重循环得到相关系数矩阵。
接下来就进行偏相关分析的流程:
下面先给出偏相关系数的公式:
三个变量的偏相关分析(剔除变量h的影响):
在计算式里面rij 是变量xi 与xj 的简单相关系数,rih 是变量xj 与xh 的简单相关系数,rjh 是变量 xj 与 xh 的简单相关系数.
四个变量的偏相关分析:
一次类推到n的变量的偏相关分析:
有上述公式可知,偏相关分析的计算公式是随着变量的增多而发生改变,且高阶偏相关分析(n阶)的计算公式的右方是由n-1阶的偏相关系数通过计算得到的,所以我们可以采用递归的方式去实现n阶偏相关分析的计算,代码如下:
function f = Pc_analysis(co_array, subscript)if length(subscript) == 2f = co_array(subscript(1),subscript(2));elseif length(subscript) > 2f_ = Pc_analysis(co_array, subscript(1:end-1));r = Pc_analysis(co_array, [subscript(1), subscript(end), subscript(3:end-1)]);r_ = Pc_analysis(co_array, [subscript(2), subscript(end), subscript(3:end-1)]);f = (f_-r*r_)/sqrt((1-r.^2)*(1-r_.^2));elsedisp('subscript参数错误,程序退出');return;endend
本文章仅为个人理解,如有不足请指出,如有建模高手,请多多指教(最近在参加数学建模哈哈)!
偏相关分析MATLAB代码实现相关推荐
- matlab的灰色关联,五种灰色关联度分析matlab代码
<五种灰色关联度分析matlab代码>由会员分享,可在线阅读,更多相关<五种灰色关联度分析matlab代码(3页珍藏版)>请在人人文库网上搜索. 1.灰色邓关联分析% p12- ...
- 【机器人学习】Delta机器人三维模型+正逆运动学分析+matlab代码
模型与代码下载地址 https://download.csdn.net/download/yjw0911/85003278 1.三维模型 正逆运动分析过程: Delta机器人三维模型+正逆运动学分析+ ...
- matlab中怎么灰色关联度,五种灰色关联度分析matlab代码
五种灰色关联度分析matlab代码 灰色邓氏关联度分析% P12 -- The Study on the Grey Relational Degree and Its Application func ...
- 【数学建模】灰色关联分析 + Matlab代码实现
文章目录 一.学习内容: 二.学习时间: 三.学习产出: 3.1 灰色关联分析基本思想 3.2 运用灰色关联分析的基本步骤 3.3 灰色关联分析代码实现(Matlab) 3.3.1 应用一:分析产业对 ...
- 【机器人学习】三自由机器人正逆运动学分析+matlab代码+直线轨迹与圆轨迹
下载链接:https://download.csdn.net/download/yjw0911/85111527 clc clear all close all M=300;l=150;L=50; t ...
- matlab灰色关联代码,灰色关联分析matlab代码
灰色关联分析 function output=grayrela(x0) %参考因子与比较因子共同存储在一个矩阵x0中,参考因子位于第一列 %斜率序列 for i=2:length(x0(:,1)) x ...
- 成矿远景叠加分析matlab代码,深水潭金成矿地质特征及成矿远景分析
1红旗沟深水潭金成矿的地质研究青海省都兰县红旗沟深水潭是青海省重要的金矿之一,红旗沟深水潭位于矿藏丰富的柴达木盆地的南源,红旗沟深水潭金成矿区的构造位置在昆中和昆北这两个较深的大断裂之间的构造带上,同 ...
- FASTICA独立成分分析matlab代码实现
%求数据的独立成分 %输入参数:Z建模数据(矩阵 n*m维,m为样本个数,n为变量个数) % perica提取信息含量(一个百分数) %输出参数:W分解矩阵(矩阵 n*q维,n为变量个数,q为独立成分 ...
- 回归预测分析MATLAB代码
%reg.m %回归分析法解 %---------------- clc %原始数据 %公路客运量(万人) glkyl=[12815 15543 19326 22864 26150 28468 308 ...
- matlab包络分析函数,数据包络分析 (MATLAB代码).doc
数据包络分析 (MATLAB代码) X=[] % 输入多指标输入矩阵X Y=[] % 输入多指标输出矩阵Y [n,m]=size(X) % n为输入的指标数,m为决策单元数 s=size(Y) % s ...
最新文章
- 基于深度学习的点云配准Benchmark
- (0099)iOS开发之Xcode编译工程报错问题汇总
- java多线程测试框架(含入参和返回值)
- echarts图表实例
- VMware 如何通过现有虚拟机克隆新的虚拟机
- phpstuday 修改网站访问目录
- SQLServer表内自关联级联删除
- java好学还是ui好学_java编程和ui设计哪个好学?
- qt c语言socket,c – Qt双向客户端服务器使用QTcpSocket和QTcpServer
- QCA9886降低功耗指令
- 易用宝项目记录day3-curd
- Discuz手机视频模板 AIUI7.3.0 商业版
- 7z001怎么解压在安卓手机上面_安卓手机用户换iPhone11怎么转移手机便签内容?...
- Cadence 应用注意事项--转载
- Ceph RBD:条带(stripe)详解
- 安防网络摄像头互联网直播录像设置永久储存却只能播放一半问题解决
- pychar调试报错:Cython extension speeds up Python debugging
- c语言中各类型所占字节
- 计算机系统漫游:贯穿计算机系统所有方面的重要概念
- python实现一个简单的项目建议书范文_建议收藏,22个Python迷你项目(附源码)