1、对于多元线性回归方法系数的求解,通常也可以使用MAT雉搽妤粲LAB中自带的函数regress来求。虽然可以调用函数,直接求解多元线性回归系数,但是求解的过程还需要劲忧商偌了解一下。首先需要创建一个函数:function dy_beta=dyxxgxs(X,y);X,y是输入的两个参数,自变量和因变量,如图所示
2、然后检测输入的因变量和自变量是否正确,如图所示~=是不等于的意思,如果不等于则报错。if nargin~=2 erro('需要输入因变量和自变量')end
3、接下来,需要求X自变量和y因变量的行列数,使用[m n]=size();如图所示[n q]=size(X);[n1 q1]=size(y);%求变量的行列数
4、然后再检验自变量和因变量的行是否相等,检测自变量y是否为列向量,如图所示if 艘绒庳焰n~=n1, error('输入参数y的行数,必须等于输入参数X的行数.');endif collhs~=1, error('输入参数y应该是一个列向量');end
5、需要对自变量X进行添加一列1的数据,将数据写成矩阵形式,需要一列1的数据才能得到正确的结果。B=[ones(n,1),X];
6、最后就是求多元线性回归模型的系数,即公式dy_beta=B\y;即可得到系数参数,如图所示