2.2.5　程序设计

2.2.5　程序设计

2.2.5.1　计算流程图

多元线性回归计算的程序流程如图2.6 所示。其中高斯—约当消去法的计算流程很多文献都有介绍，这里没有详细列出。

2.2.5.2　程序代码

（1）多元线性回归子程序。输入原始资料X_Y（m＋1，n），其中前m个为因子，第m＋1 个为预报对象，因子数m，系列长度n；输出回归系数B（0 tom）、总离差Syy、回归离差和U。在计算过程中调用高斯—约当消去法子程序GAUSSJ求解回归系数。其程序代码如下：

（2）高斯—约当消去法子程序。该子程序在这里是用来求解多元回归系数，子程序输入正规方程组的元数m，数组A（1 tom，1 tom）输入时按列存放正规方程组系数矩阵A，输出时存放A的逆阵；数组B （1tom）输入时按列存放正规方程组右侧常数项矩阵，输出时存放回归系数。其程序代码如下：

（3）F 检验子程序。F 检验子程序需要调用ShiSuanF过程。程序输入因子数m、样本容量n、信度r0 以及总离差Syy、回归离差和U。其VisualBasic程序代码如下：

函数过程ShiSuanF需要调用GAMMLN、BETACF和BETAI三个过程（三个过程的程序代码见2.1.4.2）。程序输入组间离差平方和的自由度f1、组内离差平方和的自由度f2以及信度r，输出为F 分布表中的上分位点F （α）值f3。其VisualBasic程序代码如下：

（4）t 检验子程序。t 检验子程序需要调用ShiSuanT过程（见本章2.1.4.2）。程序输入因子数m、样本容量n、信度r0以及总离差Syy、回归离差和U。其Visual Basic程序代码如下：

计算出回归系数后，应用多元回归模型式（2.34），输入各预报因子值就能计算出预报对象的预测值。(https://www.daowen.com)