matlab求偏相关系数,计算皮尔逊相关系数t检验——Matlab
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
说明：正文贴出matlab计算皮尔逊相关系数及t检验(95%)的子函数，以便查阅；
附件中包含相关系数、t检验(99%，95%，90%)M文件，自行取用；
如有疏漏，欢迎批评指正~
Part 1：% 本函数实现了皮尔逊相关系数的计算操作
% 皮尔逊相关系数计算公式(P18)及无偏相关系数矫正公式(P19)参见《现代气候统计诊断与预测技术》-魏凤英
%
% 输入：
% X：输入的数值序列
% Y：输入的数值序列
%
% 输出：
% coeff：两个输入数值序列X，Y的相关系数
%
if length(X) ~= length(Y)
error('两个数值数列的维数不相等');
return;
end
fenzi = sum((X-sum(X)/length(X)).*(Y-sum(Y)/length(Y)));
fenmu = sqrt(sum((X-sum(X)/length(X)).^2))*sqrt(sum((Y-sum(Y)/length(Y)).^2));
coef = fenzi / fenmu;
if length(X)>30
coeff = coef;
elseif length(X)>4 && length(X)<30
coeff = coef*(1+(1-coef^2)/(2*(length(X)-4))); %计算无偏相关系数加以矫正
else
coeff = coef;
fprintf('数据长度小于5 %8.4f\n',coeff);
end
end %函数myPearson结束
Part 2：
function ttest = myttest005(X , Y)
% 本函数实现了皮尔逊相关系数的t检验(α=0.05)操作
% t检验公式(P31)及t分布表(P288)参见《现代气候统计诊断与预测技术》-魏凤英
%
% 输入：
% X：输入的皮尔逊相关系数r
% Y：输入的样本量n
%
% 输出：
% ttest：通过统计检验的相关系数
%
if Y <= 2
error('计算相关系数的样本量不够');
return;
end
if X == NaN
ttest = X;
return;
end
Z = Y-2;
tvalue = sqrt(Z)*abs(X)/sqrt(1-X.^2);%用abs(X)可做正or负相关的t检验
if Z == 1
tnorm = 12.71;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 2
tnorm = 4.30;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 3
tnorm = 3.18;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 4
tnorm = 2.78;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 5
tnorm = 2.57;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 6
tnorm = 2.45;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 7
tnorm = 2.37;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 8
tnorm = 2.31;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 9
tnorm = 2.26;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 10
tnorm = 2.23;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 11
tnorm = 2.20;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 12
tnorm = 2.18;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 13
tnorm = 2.16;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 14
tnorm = 2.15;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 15
tnorm = 2.13;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 16
tnorm = 2.12;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 17
tnorm = 2.11;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 18
tnorm = 2.10;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 19
tnorm = 2.09;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 20
tnorm = 2.09;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 21
tnorm = 2.08;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 22
tnorm = 2.07;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 23
tnorm = 2.07;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 24
tnorm = 2.06;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 25
tnorm = 2.06;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 26
tnorm = 2.06;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 27
tnorm = 2.05;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 28
tnorm = 2.05;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 29
tnorm = 2.04;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z == 30
tnorm = 2.04;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z > 30 && Z < 40
tnorm = 2.03;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z >= 40 && Z < 50
tnorm = 2.02;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z >= 50 && Z < 60
tnorm = 2.01;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z >= 60 && Z < 80
tnorm = 2.00;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z >= 80 && Z < 100
tnorm = 1.99;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
elseif Z >= 100 && Z <= 120
tnorm = 1.98;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
else
tnorm = 1.96;
if tvalue > tnorm
ttest = X;
else
ttest = NaN;
end
end
end%函数myttest005结束
matlab求偏相关系数,计算皮尔逊相关系数t检验——Matlab相关推荐
- python怎么计算相关系数_Python三种方法计算皮尔逊相关系数
1.皮尔逊相关系数在统计学中,皮尔逊相关系数( Pearson correlation coefficient),又称皮尔逊积矩相关系数(Pearson product-moment correlat ...
- js计算皮尔逊相关系数
没找到网上有什么现成的,所以就自己写了个,挺简单的,记录一下. /*** 计算皮尔逊相关系数* @list1:数字列表1* @list2:数字列表2* @return number * */ func ...
- Python计算皮尔逊相关系数
** Python计算皮尔逊相关系数 ** 皮尔逊相关系数是衡量随机变量X与Y相关程度的一种方法,相关系数的取值范围是[-1,1].相关系数的绝对值越大,则表明X与Y相关度越高.当X与Y线性相关时,相 ...
- 推荐系统(2)--计算皮尔逊相关系数
计算皮尔逊相关系数 代码如下: #-*- coding:utf-8 -*- '''皮尔逊相关系数 ''' import json import numpy as np# 计算user1 和 user2 ...
- 在matlab中可以计算式子的最大值吗,matlab求最大值
matlab遗传算法求最大值问题例题完整代码_IT认证_资格考试/认证_教育专区.matlab遗传算法求最大值问题例题完整代码 本文件的目的是减少您打字的烦恼 %主程序:用...... 11MATLA ...
- pandas计算皮尔逊相关系数
前置说明: 本文中使用的数据在文末提供下载,本文使用IH和IC主力合约的数据为例 目录 1. 概念 2. 计算过程 1. 概念 皮尔逊(pearson)相关系数:用于判断数据是否线性相关的方法 皮尔逊 ...
- 皮尔逊相关系数和检验P值
看两者是否算相关要看两方面:显著水平以及相关系数 (1)显著水平,就是P值,这是首要的,因为如果不显著,相关系数再高也没用,可能只是因为偶然因素引起的,那么多少才算显著,一般p值小于0.05就是显著了 ...
- Python--如何计算皮尔逊相关系数（Pearson correlation coefficient）
皮尔逊系数 在统计学中,皮尔逊相关系数( Pearson correlation coefficient),又称皮尔逊积矩相关系数(Pearson product-moment correlation ...
- Python三种方法计算皮尔逊相关系数以及实现给定数据集，返回数据集中每个特征和标签的相关系数
特征预处理完之后,我们需要选择有意义的特征作为输入机器学习的算法和模型进行训练 在统计学中,皮尔逊相关系数( Pearson correlation coefficient),又称皮尔逊积矩相关系数( ...
- matlab求偏微分方法解析解,偏微分数值解(2,MATLAB求解方法)学案.ppt
这部分主要讨论如何用MATLAB实现对偏微分方程的数值仿真求解．MATLAB的偏微分方程工具箱(PDE Toolbox)的出现,为偏微分方程的求解以及定性研究提供了捷径．主要步骤为: 2.1 用偏微分 ...
最新文章
- 悖论对计算机科学影响,科学网—基于对角线引理和维特根斯坦思想对于悖论的分析 - 庄朝晖的博文...
- mysql 被关闭了_mysql被关闭简析
- 8个高效的Python爬虫框架分享
- 自定义控件android特效,Android自定义控件eBook实现翻书效果实例详解
- C语言那点事——如何从零学好C语言？
- C/C++高级算法之绘制曼德布洛特集
- html计时器组件,vue 计时器组件的实现代码
- 水面反光如何拍摄_拍摄水面反射的创意
- 【JAVA基础篇】IO流
- PocoClassGenerator：RDBMS所有表/视图生成Dapper POCO类代码
- linux汇编指令输出到屏幕,Linux 汇编语言（GNU GAS汇编）开发指南
- matlab gui输入数据库,从数据库值填充Matlab GUI列表框
- 如何利用即时通讯工具进行营销
- N沟道与P沟道增强型MOS管电压、原理、导通条件！
- Amazon亚马逊 Vendor Central Label详解
- 手机号不用后，这操作必须要做，不然银行卡里的钱就别人花了
- 【项目整理】安卓应用商店评论监控平台
- Path环境变量是什么？有什么用？怎么工作？JDK,JRE,JVM,集成开发工具是什么？有什么用？怎么工作？JDK安装包文件结构？有必要设置Path环境变量？Eclipse和IDEA开发工具的运作
- 基层社会治理综合管理平台智慧街道Java商业源码
- 《计算机科学与探索》期刊 从投稿到以为录用，再到拒稿过程