逐步回归matlab程序

~
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); %计算各因子的方差贡献
end
max=1;
min=1;
>
for i=1:m
if((max==1)&&(k(i)==1)&&(k(1)==0))||((v(i)>v(max))&&(k(i)==1)) max=i;
end
if((min==1)&&(k(i)==0)&&(k(1)==1))||((v(i)<v(min))&&(k(i)==0)) min=i;
end
end
《
if(l<3)&&(L+1<=m)
F1=v(max)/((r(m+1,m+1)-v(max))/(n-l-2));
if(F1>F)
disp( [ '引入第', num2str(max), '个变量']);
k(max)=0;
L=L+1;
l=l+1;
￥
r=matdel(max,m+1,r); %matdel为消去变换程序
q=1;
end
else
F2=v(min)/(r(m+1,m+1)/(n-l-1));
if((F2<F)&&(k(min)==0))
disp( [ '剔除第', num2str(min), '个变量']);
k(min)=1;
…
L=L-1;
l=l+1;
r=matdel(min,m+1,r);
q=1;
else
F1=v(max)/((r(m+1,m+1)-v(max))/(n-l-2));
if(F1>F)
disp( [ '引入第', num2str(max), '个变量']); /
k(max)=0; %如果变量i引入，则对应的k变为0 L=L+1;
l=l+1;
r=matdel(max,m+1,r);
q=1;
end
end
end
…
end
disp('没有可剔除或引入的变量，逐步回归结束');
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:L
xx=[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)'*(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))];
end
ss=[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,'%s\n',ss);
fprintf(fid,'%s\n',ss1);
fprintf(fid,'%s\n',ss2);
<
fprintf(fid,'%s\n',ss3);
fprintf(fid,'%s\n',ss4);
fprintf(fid,'%s',ss5);
fclose(fid);
end。