资源描述:
《刘杰20120512013第9章:多元回归分析》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、《数学建模与数学实验》第9章:多元回归分析实验目的:1.直观了解多元回归分析的基本内容;2.掌握用数学软件求解多元回归分析问题。实验内容:1.多元回归分析的基本理论;2.用数学软件求解多元回归分析问题;3.课堂实验练习。一、多元线性回归(一)回归模型称如下模型为多元线性回归模型(高斯-马尔可夫线性模型)其中,,,。称为回归平面方程。%Y必须是列向量多元线性回归分析的主要任务是:(1)用样本值对未知参数和作点估计;(2)对回归系数作假设检验;(3)在处对y作预测与控制,并对y作区间估计。(二)模型参数估计1.
2、回归系数的最小二乘估计作离差平方和:最小二乘法就是选择的估计,使得《数学建模与数学实验》称为的最小二乘估计(LS估计)。2.的无偏估计残差平方和:,则,称为剩余方差,它是的无偏估计。(三)检验、预测与控制1.回归方程的显著性检验对回归方程的显著性检验,可提出如下原假设:若假设被拒绝,则回归显著,认为与之间确有线性关系,所求的线性回归方程有意义;否则回归不显著。2.预测(1)点预测:根据回归方程,给定自变量,预测.(2)区间预测:给定自变量,预测的置信区间.(四)回归分析的相关命令1.确定回归系数的点估计:b
3、=regress(Y,X)2.求回归系数的点估计和区间估计、并检验回归模型:《数学建模与数学实验》残差的点估计回归系数的区间估计残差的区间估计用于检验回归模型的四个数:相关系数r2、F值、接受原假设的概率p、残差方差的估计。回归系数的点估计显著性水平[b,bint,r,rint,stats]=regress(Y,X,alpha)(1)相关系数越接近1,说明回归方程越显著;(2)时拒绝,值越大,说明回归方程越显著;(3)时,拒绝,说明回归模型成立。3.画出残差及置信区间:rcoplot(r,rint)%筛选不
4、符合的数据4.方程回归系数的检验:stepwise(X,Y)【例9-1】火柴消费与各因素之间的关系分析火柴公司的火柴销售量与各方面因素有很大联系,根据往年的销售情况,收集到了以下的一些数据:年份火柴销售量y/万件煤气、液化气用户x1/万户卷烟销售量x2/万箱蚊香销售量x3/十万盒打火石销售量X4/万粒197117.8427.4321.4311.0925.78197218.2729.9524.9614.4828.16197320.2933.5328.3716.9724.26197422.6137.3142.5
5、720.1630.18197526.7141.1645.1626.3917.08197631.1945.7352.4627.047.39197730.5050.5945.323.083.88197829.6358.8246.8024.4610.53197929.6965.2851.1133.8220.09198029.2571.2553.2933.5721.22198131.0573.3755.3639.5912.63198232.2876.6854.0048.4911.17试分析火柴的销售量与各个因素之间
6、的关系。1.回归分析及其检验《数学建模与数学实验》运行结果如下:《数学建模与数学实验》从上述分析结果可得:回归方程:相关系数:,说明线性程度较好;检验统计量:,说明回归方程显著;接受原假设的概率:,说明回归方程显著。1.残差分析,作残差图从上图可以看出,第9个点和第11个点为奇异点,可以舍弃后重新求解回归方程。2.对系数进行逐步检验《数学建模与数学实验》从上表可以看出系数对影响不显著,可以认为蚊香的销售量与火柴的销售量基本无关,删除该变量,重新计算如下:经过重新计算,可得回归方程如下:《数学建模与数学实验》
7、二、一元多项式回归(一)回归1.确定多项式系数格式:[p,s]=polyfit(x,y,m)功能:利用m次多项式对数据x,y进行曲线拟合,p为多项式系数,s为一个矩阵,用来估计预测误差。2.一元多项式回归命令格式:polytool(x,y,m)功能:此命令产生一个交互式的画面,画面中有拟合曲线和y的置信区间。(二)预测和预测误差估计格式1:Y=polyval(p,x)功能:求回归多项式p在x处的预测值Y。格式2:[Y,delta]=polyconf(p,x,s,alpha)功能:求回归多项式p在x处的预测值
8、Y,以及预测值Y的置信区间Y±delta.【例9-2】观测物体降落的距离s与时间t的关系,得到数据如下表:t(s)1/302/303/304/305/306/307/30s(cm)11.8615.6720.6026.6933.7141.9351.13t(s)8/309/3010/3011/3012/3013/3014/30s(cm)61.4972.9085.4499.08113.77129.54146.481.作