利用Matlab进行线性回归分析

9利用Matlab进行线性回归分析回归分析是处理两个及两个以上变量间线性依存关系的统计方法。可以通过软件Matlab实现。1.利用Matlab软件实现在Matlab中，可以直接调用命令实现回归分析，（1）[b,bint,r,rint,stats]=regress(y,x)，其中b是回归方程中的参数估计值，bint是b的置信区间，r和rint分别表示残差及残差对应的置信区间。stats包含三个数字，分别是相关系数，F统计量及对应的概率p值。（2）recplot(r,rint)作残差分析图。（3）rstool(x,y)一种交互式方式的句柄命令。以下通过具体的例子来说明。例现有多个样本的因变量和自变量的数据，下面我们利用Matlab，通过回归分析建立两者之间的回归方程。%一元回归分析x=[1097128415021394130315551917205121112286231120032435262529483,553372];%自变量序列数据y=[698872988807738102513161539156117651762196019022013244627362825];%因变量序列数据X=[ones(size(x')),x'],pause[b,bint,r,rint,stats]=regress(y',X,0.05),pause%调用一元回归分析函数rcoplot(r,rint)%画出在置信度区间下误差分布。10%多元回归分析%输入各种自变量数据x1=[5.52.58332.98946.55.55653.58647.57]';x2=[3155675038713056427360445039557040506259]';x3=[108127812125851112610106111199]';x4=[86916151781041671264414681311]';%输入因变量数据y=[79.3200.1163.1200.1146.0177.730.9291.9160339.4159.686.3237.5107.2155201.4100.2135.8223.3195]';X=[ones(size(x1)),x1,x2,x3,x4];[b,bint,r,rint,stats]=regress(y,X)%回归分析Q=r'*rsigma=Q/18rcoplot(r,rint);%逐步回归X1=[x1,x2,x3,x4];stepwise(X1,y,[1,2,3])%逐步回归%X2=[ones(size(x1)),x2,x3];%X3=[ones(size(x1)),x1,x2,x3];%X4=[ones(size(x1)),x2,x3,x4];%[b1,b1int,r1,r1int,stats1]=regress(y,X2)%[b2,b2int,r2,r2int,stats2]=regress(y,X3);%[b3,b3int,r3,r3int,stats3]=regress(y,X4);