全回归程序框图

读此段有困难时可略去，因其与回归应用操作无涉；或待读过附录三后再读。

本法将表1-1全部16个模型类型的回归，编入一个程序，仅需给出原始观测数据，建模自动进行。运行规制分两种，一是单点某一模型类型运行，要求给定模型号，最后输出所选模型的回归结果；二是全部运行16个模型类型并自动择优，此时不给定模型号，仅 “回车”（相当于输入“0”），程序运行中不输出各模型的中间结果，只输出最优模型的运行结果。程序运行过程见框图1-1。

图1-1中主要符号定义：N——问题维数（变量数）；M0——初始回归元数；（X，Y ）——观测数据，X为自变量；R——与回归中变量筛选有关的回归轮次数；P——模型类型编号（顺序同模型表）。以上各量，在程序初始运行前，皆由向程序输入的数据块（详后）中自动取得。T ——运行方式标记变量，手选模型标T=1，自动选择模型标T=0（T值由程序根据选定的开机方法自动确定，不用回归者特别给定）。F与Fo分别表示回归中算得的变量显著性程度的实际值与界限值；当当前模型所有变量的F值皆大于Fo，即变量都显著时，该模型类型的回归结束。T=1时，程序只做一个模型的运行；若T=0，自第一个模型类型做起，一个模型做完，程序记取该模型的残差s（也称“标准偏差” ），并转向下运行P=P+1模型，直到全部16个模型类型都做完，选s最小的模型作最优模型，此后程序令T=1，自动转入手选方式，再按选定模型做一次带中间输出的回归，即得所有16个模型类型中最优回归方程的全部信息。框图中有“*”号标记者为“最优”之意。

图1-1　全回归法程序框图

本法实施需使用奔-1以上微机。程序运行仅需输入观测数据，免除了一切不必要的中间显示，以缩短计算机运行时间。所得模型为过程规律的客观反映，无悬念，不会带入人为因素，无论多少人做回归，只要使用相同观测数据和全回归法程序，都得相同结果；也最大限度地减轻了建模者的劳动量，只要懂得简单的过程模型化常识，会取数据，会开机调用程序，每一位需求者都“会做回归分析”。特别令人放心的是，所得模型为当前条件下的最优的不可改进的模型。回归建模只应运行一遍（次），更多遍运行无益（因得相同结果）。