多元逐步回归的matlab代码

本程序利用matlab实现多元逐步回归的基本原理

function dyzbhg（xy）
%多元逐步回归分析
%作者：唐世星
%2006.11.20
%xy为待输入的原始数据，按照先x后y按列排列的数组
%如：x1 x2 x3 x4 y等等

%clc;%clear all;
%计算离差阵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））;%计算对角线元素的平方根
end
%disp（‘************ Deviation Matrix & Value of SR （离差阵R&SR） ***********‘） %输出离差阵R及SR
%[R  SR‘] 
%计算相关系数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 （相关系数阵R） **********‘）%输出相关系数阵R
%R
flag=1;%是否重复进行逐步回归的标志
while（flag）
    disp（‘******** Stepwise Regression Analysis Start *************‘）
    F1=input（‘剔除门坎值:F1=‘）;
    F2=input（‘引入门坎值:F2=‘）;
    S=0;%计算步数
    L=0;%引入方程的自变量个数
    FQ=n-1;%残差平方和的自由度    
    disp（‘************** Discriminant Value o

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

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

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

----------- ---------  ---------- -----  ----

                 7566                    2