【讨论】探讨主成分回归

[zhitongxing]

本人最近用matlab编写了主成份回归，但是在回归的时候出现了点问题，有兴趣的人请联系啊
邮箱[email protected]。

[dnping]

在论坛的化学，数学与物理版块发了关于主成分分析的书和相关程序，你可以看下……
另外，请贴出你的程序以及出错现象，好让别人帮你……

[lele]

%-----------------------------------------------------------------
%主成分分析源程序
%X表示输入原始的数据，m*n维的（m表示样本数，n表示变量数）
%contriz是要求的累计贡献率，以此来确定主成分数h
%A为测量矩阵，C校准浓度，A1为未知物的测量阵
%-----------------------------------------------------------------
function pcr(A,C,A1,contriz)
[m,n]=size(A);
[my,ny]=size(C);
[m1,n1]=size(A1);

%-----------------------------------------------------------------
%标准化原始数据
%计算样本各列的均值Z(一个n维的行向量),计算样本的标准差B（n维的行向量）
%-----------------------------------------------------------------
Z=mean(A);
B=std(A);

Y=(A-Z(ones(m,1),:))./B(ones(m,1),:)
%对A1进行标准化
Z1=mean(A1);
B1=std(A1);
Y1=(A1-Z1(ones(m1,1),:))./B1(ones(m1,1),:)


[pcs,scores,latents,tsquare]=princomp(Y) %进行主成分分析
contri=100*latents/sum(latents) %方差贡献率
lj_contri=zeros(n,1);
for i=1:n
lj_contri(i)=sum(contri(1:i));
if lj_contri(i)>=contriz*100
h=i;
break
end
end
if h<n-1
for j=h:n
lj_contri(j)=sum(contri(1:j));
end
end



lj_contri

pc=pcs(:,[1:h])
sc=scores(:,[1:h])


B=inv(sc'*sc)*sc'*C %B回归系数
F=A1*B