逐步回归matlab程序(共3页).doc

《逐步回归matlab程序(共3页).doc》由会员分享，可在线阅读，更多相关《逐步回归matlab程序(共3页).doc（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上function stepregress(x,y,F)x=zscore(x,1); %数列标准化y=zscore(y,1); %数列标准化r=corrcoef(x,y);l=0; %消去的次数L=0; %引入变量的个数n,m=size(x); %m为变量的个数，n为观测的次数k=ones(m);q=1; %判断逐步回归是否继续while(q=1)q=0; for i=1:m v(i)=r(i,m+1)2/r(i,i); %计算各因子的方差贡献endmax=1;min=1;for i=1:m if(max=1)&(k(i)=1)&(k(1)=0)|(v(i)v(max)

2、&(k(i)=1) max=i; end if(min=1)&(k(i)=0)&(k(1)=1)|(v(i)v(min)&(k(i)=0) min=i; endendif(l3)&(L+1F) disp( 引入第, num2str(max), 个变量); k(max)=0; L=L+1; l=l+1; r=matdel(max,m+1,r); %matdel为消去变换程序 q=1; endelse F2=v(min)/(r(m+1,m+1)/(n-l-1); if(F2F) disp( 引入第, num2str(max), 个变量); k(max)=0; %如果变量i引入，则对应的k变为0 L

3、=L+1; l=l+1; r=matdel(max,m+1,r); q=1; end endendenddisp(没有可剔除或引入的变量，逐步回归结束);a=zeros(L);j=1;for i=1:m if (k(i)=0) a(j)=i; j=j+1; end;end;xx=x(:,a(1);for i=2:Lxx=xx x(:,a(i);end;b=regress(y,xx); %回归系数R=sqrt(1-r(m+1,m+1); %复相关系数yyy=xx*b; %y的估计值ymean=mean(y); %y平均值Q=(y-yyy)*(y-yyy); %剩余平方和U=(yyy-ymean)

4、*(yyy-ymean); %回归平方和rs=Q/(n-L-1); %剩余方差f=U/L/(Q/(n-L-1); %F统计量fid=fopen(result,w);ss=引入第,num2str(a(1);for i=2:L ss=ss,num2str(a(i);endss=ss,个变量;ss1=y=（,num2str(b(1),x,num2str(a(1),);for i=2:L ss1=ss1,+(,num2str(b(i),x,num2str(a(i),);end;ss2=复相关系数=,num2str(R);ss3=剩余方差=,num2str(rs);ss4=F统计量=,num2str(f);ss5=剩余平方和=,num2str(Q);fprintf(fid,%sn,ss);fprintf(fid,%sn,ss1);fprintf(fid,%sn,ss2);fprintf(fid,%sn,ss3);fprintf(fid,%sn,ss4);fprintf(fid,%s,ss5);fclose(fid);end 专心-专注-专业