实验七-多元回归模型

《实验七-多元回归模型》由会员分享，可在线阅读，更多相关《实验七-多元回归模型（18页珍藏版）》请在装配图网上搜索。

1、实验七 多元回归模型（2学时）一、实验目的和要求1. 熟练掌握多元线性回归模型的建立方法，掌握并能检验所建立回归方程的显著性与方程系数的显著性，能根据实际问题作预测与控制；2.掌握平方和分解公式，会编程求总离差平方和TSS、回归平方和RSS、残差平方和ESS、复相关系数平方等统计量； 3.会根据实际问题对建立多元非线性回归模型，掌握多元线性回归的regress命令格式二、实验内容1.多元线性回归模型（1）多元线性回归模型多元线性回归模型待定常数，回归系数，.矩阵表示对进行次独立观测，得组数据 则有 ，其中 相互独立，且采用矩阵记号 -观测向量- 设计矩阵-待估回归参数向量 -随机误差向量 多元

2、线性回归模型（2）参数估计及性质-的最小二乘估计-随机误差项方差的无偏估计-回归方程给出，可由的观测值和经验回归方程求得的预测值%求回归参数命令（3）复相关系数及相关性检验总离差平方和分解 总离差残差平方和（Total Sum of Squares） 残差平方和（Error Sum of Squares） 回归平方和（Regression Sum of squares）复相关系数平方，回归愈越显著%求复相关系数平方命令TSS=sum(y-mean(y).2) %计算总离差平方和，y是因变量Y数据RSS=sum(y1-mean(y).2) %计算回归平方和ESS=sum(y-y1).2) %计算

3、残差平方和R2=RSS/ESS; %计算样本决定系数R2=RSS/TSS（4）回归方程的显著性检验检验假设： 统计量 给出显著性水平，检验值，当拒绝，认为与线性回归显著；否则线性关系不显著%回归方程显著性检验命令F=(n-p-1)*SSR/SSE %计算的F统计量,n是样本容量F1=finv(0.95,p,n-p-1) %查F统计量0.05的分位数F2=finv(0.99,p,n-p-1) %查F统计量0.01的分位数p=1-fcdf(F,p,n-p-1) %求检验P值，F是上面计算结果（5）回归系数的统计推断检验假设 统计量 检验值当，拒绝,认为与线性回归显著；否则不显著%回归系数显著性的t

4、检验命令T=b1/sqrt(SSE/(n-2)*sqrt(sum(x-mean(x).2) %t统计量观测值to, x是自变量，b1是X的回归系数T1=tinv(0.975,n-p-1) %t统计量0.05的分位数T2=tinv(0.995,n-p-1) %t统计量0.01的分位数p=2-2*tcdf(T,n-p-1) %t检验的p值（6）预测及统计推断 因变量的点估计和区间估计给出，的预测值的置信区间 4.多元线性回归建模的基本步骤(1)对问题进行直观分析，选择因变量与解释变量，作出因变量与各解释变量散点图，初步设定多元线性回归模型参数个数；(2) 多元回归建模命令输入因变量与自变量的观测数

5、据(y,X), 计算参数的估计regeress，调用格式有以下三种：（1）b = regress(Y,X)（2）b,bint,r,rint,stats = regress(Y,X)（3）b,bint,r,rint,stats = regress(Y,X,alpha) 输入参数：因变量观测向量;矩阵，第一列元素全为1，第j列是自变量Xj观测向量，对一元线性回归，取p=1即可;alpha为显著性水平.输出参数：向量b-回归系数估计值bint-回归系数的(1-alpha)置信区间；向量r-残差列向量；rint-模型的残差的(1- a)的置信区间；stats-用于检验回归模型的统计量，有4个分量值：第

6、一个是复相关系数平方，第二个是F统计量值，第三个是与统计量F对应的概率P，当P n,p=size(A); %读取A的行数和列数X=A(:,1:3) %X1-X3构成的数据矩阵 y=A(:,4)mmdl=LinearModel.fit(X,y) %调用LinearModel类的fit方法求解模型 结果输出：结果输出 mmdl = Linear regression model: y 1 + x1 + x2 + x3Estimated Coefficients:参数估计表 Estimate SE tStat pValue 参数估计 残差平方 t统计量 检验P值 (Intercept) 162.06

7、 346.15 0.46818 0.64686 x1 7.2739 1.3523 5.3787 9.7273e-05 x2 13.957 3.167 4.4071 0.00059659 x3 -4.3996 19.76 -0.22265 0.82702 以下是关于一些统计量的结果：Number of observations: 18, Error degrees of freedom: 14Root Mean Squared Error: 100R-squared: 0.957, Adjusted R-Squared 0.948F-statistic vs. constant model: 1

8、05, p-value = 7.75e-10 其他结果，点击工作区mmdl,显示很多结果。包括给出3个平方和SSE、SST、SSR及MSE，残差Residuals，因变量拟合Fittde等诸多结果。说明：1）,修正后的复相关系数平方 均认为和的线性关系显著.2）F统计量观测值,及检验概率均拒绝原假设，认为和的线性关系显著.3）从参数估计表看出, 认为, 认为认为与关系显著，与关系不显著，可进一步考虑是否剔除.此外，该类结果还给出了一系列统计量值，如Fitted-因变量拟合结果，Residual残差结果，平方和SSESSTSSR等。例3.2.3现代服务业是社会分工不断深化的产物，随着经济的发展，

9、科学技术的进步，现代服务业的发展受到多种因素和条件的影响。不仅受到经济总体发展水平的影响，还受到第二产业、就业、投入等因素的影响，从这几个主要方面出发，利用江苏省统计年鉴的有关数据，通过建立多元线性回归模型对1990-2008年各种因素对现代服务业的影响进行回归分析。假如构建如下江苏省服务业增长模型：，Y代表江苏省服务业的增加值(单位：亿元)，反映了江苏省服务业发展的总体水平。x1x4表示影响江苏省服务业发展的四种主要因素和影响，其中x1代表江苏省人均GDP(单位：元)，说明江苏省总体经济发展水平对服务业的影响；x2代表江苏省第二产业的增加值(单位：亿元)，说明了工业发展对服务业的影响，体现了

10、生产性服务业的需求规模；x3表示江苏省服务业的就业人数(单位：万人)； x4表示江苏省服务业资本形成总额(单位：亿元)，主要体现服务业投资的经济效应。表3.11 江苏省关于服务业发展及各影响因素相关数据年份服务业增加值Y省人均GDPX1第二产业增加值X2服务业就业人数X3服务业资本形成总额X4198937.76203870.24589.74252.01199028.13210935.53623.19275.82199193.582353101.33640.95330.711992160.623106325.34706.39439.321993286.584321478.79786.37620.

11、971994277.125801588.72855.97858.911995387.117319528.49920.451102.711996367.168471358.86975.661293.431997291.779371337.741025.221370.211998280.0110049228.241102.311624.741999227.6110695280.051151.681773.372000329.1611765515.741192.021903.372001385.4412882471.571263.772131.872002437.0214396697.031341.

12、862189.782003601.39168301182.621407.632686.572004704.72202231650.881443.373362.1920051291.11245601917.051542.463930.5620061360.09288141895.81625.064628.5920071769.28339282055.561713.335287.91解：（1）建立回归模型，,n=19,p=4回归分析程序%输入各影响因素的数据x0=203870.24589.74252.01210935.53623.19275.822353101.33640.95330.713106

13、325.34706.39439.324321478.79786.37620.975801588.72855.97858.917319528.49920.451102.718471358.86975.661293.439371337.741025.221370.2110049228.241102.311624.7410695280.051151.681773.3711765515.741192.021903.3712882471.571263.772131.8714396697.031341.862189.78168301182.621407.632686.57202231650.881443.

14、373362.19245601917.051542.463930.56288141895.81625.064628.59339282055.561713.335287.91;y=37.76,28.13,93.58,160.62,286.58,277.12,387.11,367.16,291.77,280.01,227.61,329.16,385.44,437.02,601.39,704.72,1291.11,1360.09,1769.28; %Y服务业增加值列向量n,p=size(x0); %矩阵的行数即样本容量n，列数即回归变量个数x=ones(n,1),x0; % 构造设计矩阵,第一列为1

15、，后4列为X1-X4db,dbint,dr,drint,dstats=regress(y,x); % 调用多元回归分析命令，输出回归参数db，回归参数置信区间dbinr，残差dr,残差置信区间drint,及一些统计量dstats程序运行结果：1）回归参数的估计n=19, p=4输出：db =345.2493 0.1672 0.1962 -0.7012 -0.6537得到回归参数估计为所以，服务业增加值Y对4个自变量的线性回归方程为回归参数的区间估计结果dbint = 22.8409 667.6576 的置信区间 0.0731 0.2613 的置信区间为0.0731, 0.2631 0.0201

16、 0.3722 的置信区间 -1.1650 -0.2374 的置信区间 -1.2867 -0.0207 的置信区间2）残差输出结果：dr = -83.8119 -59.4867 0.5889 14.6534 82.1633 7.9945 80.5396 64.6092 -72.1976 44.5876 5.7768 -4.5670 73.2709 -79.9409 -46.9710 -136.1633 113.7938 -10.20225.3624即残差向量残差的置信区间结果drint = -221.9639 54.3402 即的置信区间-221.9639, 54.3402 -205.9396

17、 86.9662 -151.4255 152.6034 -140.5198 169.8267 -62.4916 226.8182 -144.3484 160.3374 -72.5509 233.6301 -92.9528 222.1712 -214.9770 70.5819 -103.0712 192.2464 -137.4586 149.0121 -162.2339 153.0999 -64.6853 211.2271 -195.3209 35.4392 -191.1233 97.1813 -228.7041 -43.6224 -12.7351 240.3227 -149.3628 128.

18、9584 -85.7809 96.5057 的置信区间-85.7809, 96.50573）输出 统计量值dstats = 1.0e+003 * %下列数据要100.00010 0.1727 0.0000 5.7926可决系数R2=0.00010接近0， F统计量观测值,检验P值， Fa,或者由F0的p值为p0=00.05，所以拒绝H0，接受备择假设H1，说明总体回归系数bi不全为零，即表明模型的线性关系在95%的置信水平下显著成立.2）回归系数的t 检验%t检验程序S=MSE*inv(x*x); % 计算回归参数的协方差矩阵T0=db./sqrt(diag(S); % 每个回归参数的T统计量

19、 Ta=tinv(0.975,n-p-1); % t分布的分位数pp=tpdf(T0,n-p-1); % 每个回归参数的T统计量对应的概率t检验程序结果：S =1.0e+04 * 2.2597 -0.0003 -0.0004 -0.0032 0.0024 -0.0003 0.0000 0.0000 0.0000 -0.0000 -0.0004 0.0000 0.0000 0.0000 -0.0000 -0.0032 0.0000 0.0000 0.0000 -0.00000.0024 -0.0000 -0.0000 -0.0000 0.0000说明：,对角元素开方为，t统计量值 T0 = 2.

20、2967 3.8122 2.3899 -3.2424 -2.2149Ta =2.1448 T检验P值 pp = 0.0356 0.0019 0.0301 0.0059 0.0412t检验：统计量,对角元为，表3.13 参数估计表变量值标准差t值p值常数项省人均GDP第二产业增加值服务业就业人数服务业资本形成总额345.24930.16720.1962-0.7012-0.6537150.3220.0440.0820.2160.2952.29673.81222.3899-3.2424-2.21480.03560.00190.03010.00590.0412给a=0.05，表3.13 最后一列概率均小于0.05，拒绝H0，接受备择假设H1,即回归系数bi(i=0,1,2,3,4)显著不为零.3）复相关系数及相关性检验%可决系数检验 程序：R2=RSS/TSS; % 计算样本决定系数程序结果R2 =0.9801该方程的拟合优度判定系数调整后的拟合优度判定系数由此说明该多元线性回归方程的拟合程度比较理想.