regress多元回归拟合代码

admin

n=10;m=5;alpha=0.05;
x1=[10521667313.33,13880396970.82,12685660578.41,14048932040.99,15010921209.65,16103786351.41,16334811436.83,17277037490.47,20351939978.58,21610405402.53];

x2=[1041709112.47,1440816247.17,1268359415.28,1263787426.53,1368801995.54,1484218837.50,1490854303.85,1520329814.80,1798177056.32,1857104874.30];
x3=[0.0145,0.2215,0.0975,0.01472,0.0851,0.0254,0.0461,0.0788,0.0415,-0.0332];

x4=[3.7127,2.5508,1.8585,1.4840,1.6380,2.2226,4.4028,2.94,3.1862,4.6272];
x5=[37021382625.531,32045534525.480,23974111617.647,16287217039.656,22618981489.215,35640964099.679,116369639594.802,40158146232.970,35071362434.280,26814763200.536];
y=[9.37,15.18,13.02,6.56,6.07,8.99,20.36,14.39,7.93,4.97];
x=[ones(n,1),x1',x2',x3',x4',x5'];


[b,bint,r,rint,s]=regress(y',x);
s2=sum(r.^2)/(n-m-1);
%b,bint,s,s2;
rcoplot(r,rint);
title('残差图的绘制');
xlabel('数据 ');ylabel('残差');