pcorrxy.m

来自「wang xiao ping 版遗传算法」· M 代码 · 共 33 行

function [partial]=pcorrxy(x,var1,var2)
%
% PCORRXY
%
% Returns partial correlation between var1 and var2
% contained in a matrix x
%
% All matrices are standardised prior to calculations.
%
% Called: [partial]=pcorrxy(x,var1,var2)
%
% Calls: standm(x)
%	 reduce(x,var1,var2)
%

stdx=standm(x);
xx=reduce(stdx,var1,var2);

temp=inv(xx'*xx)*xx';

par1=temp*stdx(:,var1);
par2=temp*stdx(:,var2);

res1=stdx(:,var1)-xx*par1;
res2=stdx(:,var2)-xx*par2;

stdres1=standvec(res1);
stdres2=standvec(res2);

[nvals,ncols]=size(stdres1);

partial=(stdres1'*stdres2)/(nvals-1);