1、偏最小二乘法代碼clcclearload aa.txt ;%原始數(shù)據(jù)存放在純文本文件aa.txt中disp(期望)mu=mean(aa)disp(和標(biāo)準(zhǔn)差 )sig=std(aa) %求均值和標(biāo)準(zhǔn)差disp(相關(guān)系數(shù)矩陣)rr=corrcoef(aa) %求相關(guān)系數(shù)矩陣disp(數(shù)據(jù)標(biāo)準(zhǔn)化)data=zscore(aa) % 數(shù)據(jù)標(biāo)準(zhǔn)化n=3;m=3; %n是自變量的個(gè)數(shù),m是因變量的個(gè)數(shù)x0=aa(:,1:n);y0=aa(:,n+1:end);e0=data(:,1:n);f0=data(:,n+1:end);num=size(e0,1);%求樣本點(diǎn)的個(gè)數(shù)chg=eye(n); %w到w*

2、變換矩陣的初始化for i=1:n%以下計(jì)算w, w*和t的得分向量，matrix=e0*f0*f0*e0;vec,val=eig(matrix);%求特征值和特征向量val=diag(val); %提出對(duì)角線(xiàn)元素val,ind=sort(val,descend);w(:,i)=vec(:,ind(1) %提出最大特征值對(duì)應(yīng)的特征向量w_star(:,i)=chg*w(:,i); %計(jì)算 w*的取值t(:,i)=e0*w(:,i); %計(jì)算成分ti的得分alpha=e0*t(:,i)/(t(:,i)*t(:,i); %計(jì)算 alpha_ichg=chg*(eye(n)-w(:,i)*alpha

3、); %計(jì)算 w 到 w*的變換矩陣e=e0-t(:,i)*alpha; % 計(jì)算殘差矩陣e0=e;%以下計(jì)算ss(i)的值beta=t(:,1:i),ones(num,1)f0; % 求回歸方程的系數(shù)beta(end,:)=; %刪除回歸分析的常數(shù)項(xiàng)cancha=f0-t(:,1:i)*beta; % 求殘差矩陣ss(i)=sum(sum(cancha.八2); %求誤差平方和%以下計(jì)算press(i)for j=1:num t1=t(:,1:i);f1=f0;she_t=t1(j,:);she_f=f1(j,:); %把舍去的第j個(gè)樣本點(diǎn)保存起來(lái)t1(j,:)=;f1(j,:)=; %刪除

4、第 j 個(gè)觀測(cè)值beta1=t1,ones(num-1,1)f1; % 求回歸分析的系數(shù)beta1(end,:)=; %刪除回歸分析的常數(shù)項(xiàng)cancha=she_f-she_t*beta1; % 求殘差向量press_i(j)=sum(cancha.A2);endpress(i)=sum(press_i);if i1Q_h2(i)=1-press(i)/ss(i-1);elseQ_h2(1)=1;endif Q_h2(i)0.0975fprintf(提出的成分個(gè)數(shù) r=%d,i);r=i;breakendenddisp(Y關(guān)于t的回歸系數(shù))beta_z=t(:,1:r),ones(num,1)

5、f0 %求 Y 關(guān)于 t 的回歸系數(shù)beta_z(end,:)=; % 刪除常數(shù)項(xiàng)disp(求Y關(guān)于X的回歸系數(shù)，且是針對(duì)標(biāo)準(zhǔn)數(shù)據(jù)的回歸系數(shù)) xishu=w_star(:,1:r)*beta_z %求Y關(guān)于X的回歸系數(shù)，且是針對(duì)標(biāo)準(zhǔn)數(shù)據(jù)的回歸系數(shù)，%每一列是一個(gè)回歸方程mu_x=mu(1:n);mu_y=mu(n+1:end);sig_x=sig(1:n);sig_y=sig(n+1:end);for i=1:mch0(i)=mu_y(i)-mu_x./sig_x*sig_y(i)*xishu(:,i); %計(jì)算原始數(shù)據(jù)的回歸方程的常數(shù)項(xiàng)endfor i=1:mxish(:,i)=xishu

6、(:,i)./sig_x*sig_y(i); %計(jì)算原始數(shù)據(jù)的回歸方程的系數(shù)，每一列是一個(gè)回歸方程endsol=ch0;xish %顯示回歸方程的系數(shù)，每一列是一個(gè)方程，每一列的第一 個(gè)數(shù)是常數(shù)項(xiàng)disp(顯示回歸方程的系數(shù)，每一列是一個(gè)方程，每一列的第一個(gè)數(shù)是 常數(shù)項(xiàng));save mydata x0 y0 num xishu ch0 xishload mydatanumch0=repmat(ch0,num,1);yhat=ch0+x0*xish; % 計(jì)算 y 的預(yù)測(cè)值y1max=max(yhat);y2max=max(y0);ymax=max(y1max;y2max)cancha=yhat-y0; % 計(jì)算殘差subplot(2,2,1)plot(0:ymax (1) ,0:ymax(1),yhat(:,1),y0(:,1),*)subplot(2,2,2)plot(