function?step_wise（）
%多元逐步回歸本程序在MATLAB2008a下運行通過.
clc;clear?all;
disp（‘----------姓名:滕傳亮學號:s08030482-----------‘）;
Fin=input（‘F進=‘）;%F進
Fout=input（‘F出=‘）;%F出
%XY=[x1?x2?x3?...?xn?y]，先x后y按列排列的數組
XY=[137261911.5;1511403419.8;218291713.7;1912153321.6;2711132722.3;3210211519.1;
	178181611.7;2610352319.4;146141810.6;2813213425.5;199132918.7;1210193819.3;
	238251715.6;2811333224.7;219181915.3;3514243429.8;166191410.2;2410322619.8;
	2211393825.3;10717209.7;188342214.8;2911282120.7;1811163219.6;1610153420.3;
	187231411.1;2311292920.7;2513414028.9;329121518.3;3611371821.5;319251417.7;
	2913143828.3;1810113521.6];
[mn]=size（XY）;
n=n-1;%自變量個數
f=0;%引入變量個數
s=0;%計算步數
mXY=mean（XY）;%求各列平均值
L=XY‘*XY-m*mXY‘*mXY;%計算離差陣L（n+1n+1）
%計算相關系數陣R（n+1n+1）
R=diag（1./sqrt（diag（L）））;
R=R*L*R;
In=[];%已引入的所有變量下標索引數組
Out=1:n;%未引入的所有變量下標索引數組
disp（‘開始逐步回歸計算:‘）;
while?1
%計算引入變量的偏回歸平方和，和未引入變量的引入貢獻
P=R（1:nend）.^2./diag（R（1:n1:n））;
if?s>1
????[pminimin]=min（P（In））;%尋找已引入自變量方差貢獻的最小值和索引
????Fx=（m-f-1）*pmin/R（endend）;%計算待剔除變量的F檢驗值
????if?Fx<=Fout
????????s=s+1;????????
????????disp（sprintf（‘=====================第%d步=====================‘s））;
????????disp（sprintf（‘變量x%d被剔除F%d=%f<=F出=%f‘In（imin）In（imin）FxFout））;????????
????????R=T（RIn（imin））;%消去變換
????????Out（end+1）=In（imin）;%將剔除的變量下標加入未引入變量索引數組
????????In（imin）=[];%將剔除的變量下標從引入變量索引數組中刪除????????????
????????f=f-1;%引入變量個數減一
????????continue;
????end
end
[pmaximax]=max（P（Out））;%尋找未引入變量方差貢獻最大值和索引