逐步回歸matlab程序.doc

逐步回歸matlab程序.doc

ID:58566890

大?。?6.50 KB

頁(yè)數(shù):3頁(yè)

時(shí)間:2020-10-20

逐步回歸matlab程序.doc_第1頁(yè)
逐步回歸matlab程序.doc_第2頁(yè)
逐步回歸matlab程序.doc_第3頁(yè)
資源描述:

《逐步回歸matlab程序.doc》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)

1、functionstepregress(x,y,F)x=zscore(x,1);%數(shù)列標(biāo)準(zhǔn)化y=zscore(y,1);%數(shù)列標(biāo)準(zhǔn)化r=corrcoef([x,y]);l=0;%消去的次數(shù)L=0;%引入變量的個(gè)數(shù)[n,m]=size(x);%m為變量的個(gè)數(shù),n為觀測(cè)的次數(shù)k=ones(m);q=1;%判斷逐步回歸是否繼續(xù)while(q==1)q=0;fori=1:mv(i)=r(i,m+1)^2/r(i,i);%計(jì)算各因子的方差貢獻(xiàn)endmax=1;min=1;fori=1:mif((max==1)&&(k(i)==1)&&(k(1)==

2、0))

3、

4、((v(i)>v(max))&&(k(i)==1))max=i;endif((min==1)&&(k(i)==0)&&(k(1)==1))

5、

6、((v(i)F)disp(['引入第',num2str(max),'個(gè)變量']);k(max)=0;L=L+1;l=l+1;r=matdel(max,m+1,r);%matdel為消去變換程序q=1;e

7、ndelseF2=v(min)/(r(m+1,m+1)/(n-l-1));if((F2F)disp(['引入第',num2str(max),'個(gè)變量']);k(max)=0;%如果變量i引入,則對(duì)應(yīng)的k變?yōu)?L=L+1;l=l+1;r=matdel(ma

8、x,m+1,r);q=1;endendendenddisp('沒有可剔除或引入的變量,逐步回歸結(jié)束');a=zeros(L);j=1;fori=1:mif(k(i)==0)a(j)=i;j=j+1;end;end;xx=x(:,a(1));fori=2:Lxx=[xxx(:,a(i))];end;b=regress(y,xx);%回歸系數(shù)R=sqrt(1-r(m+1,m+1));%復(fù)相關(guān)系數(shù)yyy=xx*b;%y的估計(jì)值ymean=mean(y);%y平均值Q=(y-yyy)'*(y-yyy);%剩余平方和U=(yyy-ymean)'*(

9、yyy-ymean);%回歸平方和rs=Q/(n-L-1);%剩余方差f=U/L/(Q/(n-L-1));%F統(tǒng)計(jì)量fid=fopen('result','w');ss=['引入第',num2str(a(1))];fori=2:Lss=[ss,',',num2str(a(i))];endss=[ss,'個(gè)變量'];ss1=['y=(',num2str(b(1)),'x',num2str(a(1)),')'];fori=2:Lss1=[ss1,'+(',num2str(b(i)),'x',num2str(a(i)),')'];end;ss2

10、=['復(fù)相關(guān)系數(shù)=',num2str(R)];ss3=['剩余方差=',num2str(rs)];ss4=['F統(tǒng)計(jì)量=',num2str(f)];ss5=['剩余平方和=',num2str(Q)];fprintf(fid,'%s',ss);fprintf(fid,'%s',ss1);fprintf(fid,'%s',ss2);fprintf(fid,'%s',ss3);fprintf(fid,'%s',ss4);fprintf(fid,'%s',ss5);fclose(fid);end

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁(yè),下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動(dòng)畫的文件,查看預(yù)覽時(shí)可能會(huì)顯示錯(cuò)亂或異常,文件下載后無此問題,請(qǐng)放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫(kù)負(fù)責(zé)整理代發(fā)布。如果您對(duì)本文檔版權(quán)有爭(zhēng)議請(qǐng)及時(shí)聯(lián)系客服。
3. 下載前請(qǐng)仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動(dòng)等原因無法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請(qǐng)聯(lián)系客服處理。