function?dyzbhg（xy）
%多元逐步回歸分析
%作者：唐世星
%2006.11.20
%xy為待輸入的原始數(shù)據(jù)，按照先x后y按列排列的數(shù)組
%如：x1?x2?x3?x4?y等等

%clc;%clear?all;
%計(jì)算離差陣R（mm）
[nm]=size（xy）;

%F1=0;F2=0;
%disp（‘均值為:‘）
xy_aver=mean（xy）%求均值
for?i=1:m
????for?j=1:i
????????R（ij）=0;
????????for?k=1:n
????????????R（ij）=R（ij）+（xy（ki）-xy_aver（i））*（xy（kj）-xy_aver（j））;
????????end
????????R（ji）=R（ij）;
????end
????SR（i）=sqrt（R（ii））;%計(jì)算對(duì)角線元素的平方根
end
%disp（‘************?Deviation?Matrix?&?Value?of?SR?（離差陣R&SR）?***********‘）?%輸出離差陣R及SR
%[R??SR‘]?
%計(jì)算相關(guān)系數(shù)R（mm）
for?i=1:m
????for?j=1:i
????????R（ij）=R（ij）/（SR（i）*SR（j））;
????????R（ji）=R（ij）;
????end
end
%disp（‘**********?Correlation?Coefficient?Matrix?（相關(guān)系數(shù)陣R）?**********‘）%輸出相關(guān)系數(shù)陣R
%R
flag=1;%是否重復(fù)進(jìn)行逐步回歸的標(biāo)志
while（flag）
????disp（‘********?Stepwise?Regression?Analysis?Start?*************‘）
????F1=input（‘剔除門坎值:F1=‘）;
????F2=input（‘引入門坎值:F2=‘）;
????S=0;%計(jì)算步數(shù)
????L=0;%引入方程的自變量個(gè)數(shù)
????FQ=n-1;%殘差平方和的自由度????
????disp（‘**************?Discriminant?Value?o

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????305??2006-12-12?12:03??xiaoqu.m

?????文件???????7261??2006-12-12?12:02??dyzbhg.m

-----------?---------??----------?-----??----

?????????????????7566????????????????????2