tbd matlab,如何用matlab实现偏相关分析

matlab偏最小二乘分析(2009-11-20 16:43:07)转载▼

标签：杂谈

单因变量function y=pls(pz)

[row,col]=size(pz);

aver=mean(pz);

stdcov=std(pz); %求均值和标准差

rr=corrcoef(pz); %求相关系数矩阵

%data=zscore(pz); %数据标准化

stdarr = ( pz - aver(ones(row,1),:) )./ stdcov( ones(row,1),:); % 标准化数据结果与zscore()一致

x0=pz(:,1:col-1);y0=pz(:,end); %提取原始的自变量、因变量数据

e0=stdarr(:,1:col-1);f0=stdarr(:,end); %提取标准化后的自变量、因变量数据

num=size(e0,1);%求样本点的个数

temp=eye(col-1);%对角阵

for i=1:col-1

%以下计算 w，w*和 t 的得分向量，

w(:,i)= ( e0'* f0 )/ norm( e0'*f0 );

t(:,i)=e0*w(:,i) %计算成分 ti 的得分

alpha(:,i)=e0'*t(:,i)/(t(:,i)'*t(:,i)) %计算 alpha_i ,其中(t(:,i)'*t(:,i))等价于norm(t(:,i))^2

e=e0-t(:,i)*alpha(:,i)' %计算残差矩阵

e0=e;

%计算w*矩阵

if i