起来的模型。(3 )先观察观察模型残差 Matlab 代码如下: n=20; m=2; X=[ones(n,1),x1',x2']; [b,bint,r,rint,s]=regress(y',X); b,bint,s rcoplot(r,rint) 结果是: 如图所示,应该剔除第 8 组和第 20 组数据。在执行如下 Matlab 代码: y=[11.2 13.4 40.7 5.3 24.8 12.7 20.9 8.7 9.6 14.5 26.9 15.7 36.2 18.1 28.9 14.9 25.8 21.7 ]; x1=[16.5 20.5 26.3 16.5 19.2 16.5 20.2 17.2 14.3 18.1 23.1 19.1 24.7 18.6 24.9 17.9 22.4 20.2 ]; x2=[6.2 6.4 9.3 5.3 7.3 5.9 6.4 4.9 6.4 6.0 7.4 5.8 8.6 6.5 8.4 6.7 8.6 8.4 ]; x3=[587 643 635 692 1248 643 1964 713 749 7895 762 2793 741 625 854 716 921 595 ]; n=18; m=2; X =[ones(n,1),x1',x2']; [b,bint,r,rint,s]=regress(y',X); b,bint,s rcoplot(r,rint) stepwise(x,y') % 进行逐步回归 pause ; n=18; X =[ones(n,1) x1' x2']; % 由前面的逐步回归可以得到包含 2 个变量 x1,x2 时s 最小[b,bi,r,ri,s]=regress(y',x); s2=sum(r.^2)/(n-3); b,bi,s,s2 rcoplot(r,ri) % 残差分析剔除之后结果如下:
全文预览
