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)参见《现代气候统计诊断与