


下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、偏最小二乘回歸matlab程序代碼偏最小二乘回歸matlab程序代碼單因變量function y=pls(pz)row,col=size(pz);aver=mean(pz);stdcov=std(pz); %求均值和標(biāo)準(zhǔn)差rr=corrcoef(pz); %求相關(guān)系數(shù)矩陣%data=zscore(pz); %數(shù)據(jù)標(biāo)準(zhǔn)化stdarr = ( pz - aver(ones(row,1),:) )./ stdcov( ones(row,1),:); % 標(biāo)準(zhǔn)化數(shù)據(jù)結(jié)果與zscore()一致x0=pz(:,1:col-1);y0=pz(:,end); %提取原始的自
2、變量、因變量數(shù)據(jù)e0=stdarr(:,1:col-1);f0=stdarr(:,end); %提取標(biāo)準(zhǔn)化后的自變量、因變量數(shù)據(jù)num=size(e0,1);%求樣本點(diǎn)的個(gè)數(shù)temp=eye(col-1);%對(duì)角陣for i=1:col-1%以下計(jì)算 w,w*和 t 的得分向量, w(:,i)= ( e0'* f0 )/ norm( e0'*f0 ); t(:,i)=e0*w(:,i) %計(jì)算成分 ti 的得分 &
3、#160; alpha(:,i)=e0'*t(:,i)/(t(:,i)'*t(:,i) %計(jì)算 alpha_i ,其中(t(:,i)'*t(:,i)等價(jià)于norm(t(:,i)2 e=e0-t(:,i)*alpha(:,i)' %計(jì)算殘差矩陣 e0=e; %計(jì)算w*矩陣 if i=1
4、0; w_star(:,i)=w(:,i); else for j=1:i-1 temp=temp*(eye(col-1)-w(:,j)*alpha(:,j)');
5、0; end w_star(:,i)=temp*w(:,i); end %以下計(jì)算 ss(i)的值 beta=t(:,1:i),ones(num,1)f0 %求回歸方程的系數(shù) beta(end,:)=;
6、160; %刪除回歸分析的常數(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
7、=t1(j,:);she_f=f1(j,:); %把舍去的第 j個(gè)樣本點(diǎn)保存起來(lái) t1(j,:)=;f1(j,:)=; %刪除第j個(gè)觀測(cè)值 beta1=t1,ones(num-1,1)f1; %求回歸分析的系數(shù) beta1(end,:)=;
8、 %刪除回歸分析的常數(shù)項(xiàng) cancha=she_f-she_t*beta1; %求殘差向量 press_i(j)=sum(cancha.2); end press(i)=sum(press_i)
9、; if i>1 q_h2(i)=1-press(i)/ss(i-1) else q_h2(1)=1 end if q_h2(i)< fprintf('提出的成分個(gè)數(shù) r=%d',i); &
10、#160; r=i; break endendbeta_z=t,ones(num,1)f0; %求標(biāo)準(zhǔn)化y關(guān)于主成分得分向量t的回歸系數(shù)beta_z(end,:)=; %刪除常數(shù)項(xiàng)xishu=w_star*beta_z; %求標(biāo)準(zhǔn)化y關(guān)于x的回歸系數(shù), 且是針對(duì)標(biāo)準(zhǔn)數(shù)據(jù)的回歸系數(shù),每一列是一個(gè)回歸方
11、程mu_x=aver(1:col-1);mu_y=aver(end);sig_x=stdcov(1:col-1);sig_y=stdcov(end);ch0=mu_y-mu_x./sig_x*sig_y*xishu; %計(jì)算原始數(shù)據(jù)的回歸方程的常數(shù)項(xiàng)xish=xishu'./sig_x*sig_y; %計(jì)算原始數(shù)據(jù)的回歸方程的系數(shù),每一列是一個(gè)回歸方程rc=corrcoef(x0*xish'+ch0,y0)sol=ch0;xish' %顯示回歸方程的系數(shù),每一列是一個(gè)方程,每一列的第一個(gè)
12、數(shù)是常數(shù)項(xiàng)多因變量function y=pls(pz,xnum,ynum)row,col=size(pz);aver=mean(pz);stdcov=std(pz); %求均值和標(biāo)準(zhǔn)差rr=corrcoef(pz); %求相關(guān)系數(shù)矩陣data=zscore(pz); %數(shù)據(jù)標(biāo)準(zhǔn)化stdarr = ( pz - aver(ones(row,1),:) )./ stdcov( ones(row,1),:); % 標(biāo)準(zhǔn)化自變量n=xnum;m=ynum; %n 是自變量的個(gè)數(shù),m是因變量的個(gè)數(shù)x0=pz(:,1:n);y0=pz(:,n+1
13、:end); %提取原始的自變量、因變量數(shù)據(jù)e0=data(:,1:n);f0=data(:,n+1:end); %提取標(biāo)準(zhǔn)化后的自變量、因變量數(shù)據(jù)num=size(e0,1);%求樣本點(diǎn)的個(gè)數(shù)temp=eye(n);%對(duì)角陣for i=1:n%以下計(jì)算 w,w*和 t 的得分向量, matrix=e0'*f0*f0'*e0; vec,val=eig(matrix) %求特征值和特征向量 val=diag(val); %提出對(duì)角線元素&
14、#160; val,ind=sort(val,'descend'); w(:,i)=vec(:,ind(1) %提出最大特征值對(duì)應(yīng)的特征向量 t(:,i)=e0*w(:,i) %計(jì)算成分 ti 的得分 alpha(:,i)=e0'*t(:,i)/(t(:,i)'*t(:,i) %計(jì)算 alpha_i ,其中(t(:,i)'*t(:,i)等價(jià)于
15、norm(t(:,i)2 e=e0-t(:,i)*alpha(:,i)' %計(jì)算殘差矩陣 e0=e; %計(jì)算w*矩陣 if i=1 w_star(:,i)=w(:,i);
16、else for j=1:i-1 temp=temp*(eye(n)-w(:,j)*alpha(:,j)'); end
17、60; w_star(:,i)=temp*w(:,i); end %以下計(jì)算 ss(i)的值 beta=t(:,1:i),ones(num,1)f0 %求回歸方程的系數(shù) beta(end,:)=; %刪除回歸分析的常數(shù)項(xiàng) cancha=f0-t(:,1:i)*beta; %求殘差矩陣
18、60; 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,:)=;f
19、1(j,:)=; %刪除第j個(gè)觀測(cè)值 beta1=t1,ones(num-1,1)f1; %求回歸分析的系數(shù) beta1(end,:)=; %刪除回歸分析的常數(shù)項(xiàng)
20、60; cancha=she_f-she_t*beta1; %求殘差向量 press_i(j)=sum(cancha.2); end press(i)=sum(press_i) if i>1 q_h2(i)=1-press(i)/ss(i-1) else q_h2(1)=1 end if q_h2(i)< fprintf('提出的成
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 項(xiàng)目的可行性分析與評(píng)估計(jì)劃
- 2025年電子廚房秤項(xiàng)目建議書(shū)
- 2025年微球載體材料合作協(xié)議書(shū)
- 醫(yī)院銷售合同
- 電商交易平臺(tái)的商品展示與服務(wù)免責(zé)協(xié)議書(shū)
- Oxadixyl-Standard-生命科學(xué)試劑-MCE
- Dimethenamide-P-Standard-生命科學(xué)試劑-MCE
- 4-Aminonicotinic-acid-生命科學(xué)試劑-MCE
- 2-3-Isopropylideneguanosine-生命科學(xué)試劑-MCE
- 幼兒繪本小藍(lán)和小黃讀后感
- 信息技術(shù)基礎(chǔ)ppt課件(完整版)
- 幼兒繪本故事:奇怪的雨傘店
- 弘揚(yáng)與傳承中華傳統(tǒng)文化課件(共16張PPT)
- 鋼琴基礎(chǔ)教程教案
- 電子課件-《飯店服務(wù)心理(第四版)》-A11-2549
- 糖基轉(zhuǎn)移酶和糖苷酶課件(PPT 111頁(yè))
- (語(yǔ)文A版)四年級(jí)語(yǔ)文下冊(cè)課件跳水 (2)
- 自來(lái)水業(yè)務(wù)辦理授權(quán)委托書(shū)
- 電子商務(wù)客戶服務(wù)ppt課件匯總(完整版)
- 2022危險(xiǎn)化學(xué)品及石化化工項(xiàng)目必須進(jìn)園區(qū)的規(guī)定文件匯總
- 【單元設(shè)計(jì)】第七章《萬(wàn)有引力與宇宙航行》單元教學(xué)設(shè)計(jì)及教材分析課件高一物理人教版(2019)必修第二冊(cè)
評(píng)論
0/150
提交評(píng)論