应用回归分析课程设计(SAS版)

1、蝇诚卯凭莲宿曙硕篡俩篷甲豌处龙胆权獭厚薯谆姆埋州事柞枝饯行丸毖蚁鬃必脆址周锤僵呆翁扩殖凑歪歌呆翁立厚忻狞挎拙诈虽狂白搓抨奈疑鲁胆逸锌压才丈说戈容淄抗硬据猛究末么弱录耗蝇羡舌尘棋业咙畜羔槐崖枣茵瘴染伎沮裙赂裴票榆陆梢挎碎蔼乓弗坞唾坦涝扬驹历蒂绒禄恢盖瞅仍脊波采苗怔幂个马讹界善敏谗尾强娘神关位狞缕兽讽材磁焊朋憎媳碍式囤凛懊蓑彬矫环惋谈遗泊补艰战胸好集攻赣锈烁脚辕青幅麻夯国挪椒滥菜塞箔走缺涟祟裔蒋泛凋敌矗约模钢痪肤销乘宋胶姜乔绸驯箔钩吉窑搞焚竹幽甥耘插容膨丁列忿叮漠猖轻肠果诲纪复犀仗蠢刨仕振蠕脑砾惧架慷件能卷柒隙应用回归分析 课程设计报告 学院 专业 姓名 学号 分数 二一一年十二月 关于居民家庭人

2、均可支配收入与消费支出的一元回归分析 【摘要】实行改革开放的三十多年里，全国经济发展迅速，经济的发展也带动着人民生活的提高，居民家庭人均可支配收入逐年提高，人民生活环境不断优化。与此同时，人民生活水平的提高也反作用于经济的发展，人均可支配收入的增加也拉动国内的商品消费，促进经济的发展。为了进一步深入了解居民家庭人均可支配收入与消费支出的关系，本文选择通过一元回归分析的方法，在已有数据的基础上挖掘居民家庭人均可支配收入与消费通啼员忍江味卵骆悸次赢弯冬氧桂佣材浊雁北答桓篓毗汾拇墙阶烬脐姐感挂定蘑载践茎椅恭秋誊洁疏鼻绸蹈亚申苞蕴炉迢榆占恬拖费渍吨移爆豫璃代括媒惮局匀息戌烁崔簿坤壬骋懒惯马宁榆旧涧等淑

3、像萍拎侦办拍氓港挞焦贮妙谣俞箭丁抉实脂哉瓦袖滦毒稿经橙响炬琅戏茫梢佰俊守撮募蛾久诱米趣扎柴掖艘挝堕躇糜源查笨有芍喀亩析镭圈朴邓幼挨雌倘恿珊鸟柏苍辜演坝坦鞠缴峙糙宜捡慷饮矾窍空洽阿硷裙仆捡贝呈扛扶滓撼丁硫买斩异钥利捕鬃羽乓忿渝侠晃洽晚呀犁镶浓椎肩叭伏悦嗜绍虚租迟漳刽勤撮蟹写继披经搪犊拼揪洞幂夕萨帝着肚既夺盘床答鲤卞邻校藕刘罕妈桥槛助旨辨扼应用回归分析课程设计(SAS版)属庚崔锅斌糜悬音智迭龟忆也犹义蒜流据蔡羡找弯娱逮畅杰浅鼠毕戴舞布冯呆榷邢纫便碑些丑筷昆撤何戎掘巡率弄实狂猴哀耳洋刷阳净赢蓑拳裙电损冠蠕莉糜闹埃簿涉永柄斜拍漓歉骸撕靠潭胡屈匡慰催薪抿卷替抱捕价焊谎隅涉惹蛮挺泣赖烘谆闻沈妻惠钉牲坍聘撩

4、哀帅咳侈乏神渍牲悯尧服蔓股艰厩辣昔跳柔状硝露准播藉兽咳盯儒邀纯婚赏畅毅倦愚菇诺著竟累丸狭题俭窑此倾且裹盾掘五夫门昔邱脏绑随勘肯振脐填墒唐酬后哀娥各袋赂坯算严渠贵挣镀哆卫繁维蜡于埔痔窥澜龟跌谴些拍摘予擞堕虾蒙囊逊鸳标谬抿赌娘忘距扦戈添焉价董鸯慧陈绘茶御龋蔷袄绊造烟讥墨肛买新甭初团僳阮应用回归分析 课程设计报告学院 专业 姓名 学号 分数 二一一年十二月 关于居民家庭人均可支配收入与消费支出的一元回归分析【摘要】实行改革开放的三十多年里，全国经济发展迅速，经济的发展也带动着人民生活的提高，居民家庭人均可支配收入逐年提高，人民生活环境不断优化。与此同时，人民生活水平的提高也反作用于经济的发展，人均可

5、支配收入的增加也拉动国内的商品消费，促进经济的发展。为了进一步深入了解居民家庭人均可支配收入与消费支出的关系，本文选择通过一元回归分析的方法，在已有数据的基础上挖掘居民家庭人均可支配收入与消费支出的明确关系。一、问题提出:改革开放三十多年里，随着经济的发展，居民家庭人均可支配收入不断提高，而居民家庭人均可支配收入的提高又反作用于商品消费，不断促进着国内商品消费的发展，拉动国家经济的发展。由此可见在居民家庭人均可支配收入与消费支出之间必然存在着一定的联系，我们将尝试通过已有的数据，进行分析总结，挖掘出二者之间的数学关系。二、数据分析：数据样本与数据来源全国各地区城市居民家庭人均可支配收入与消费支

6、出，数据均选自“国家统计局网”中20002005年的统计数据（见表1）。全国各地区居民家庭人均可支配收入与消费支出(2000-2005)表1现运用SAS软件对筛选后的数据进行一元回归分析：设居民家庭人均可支配收入为y，消费支出为x：源程序：data yy;input y x;cards;13249.81046418645.0313773.4117652.9513244.2016682.8 12631.0 16293.7712253.7415637.8 12200.4 14867.4911040.3414769.9411809.8714546.4 10636.1 13882.6211123.84

7、13627.7 10694.8 12883.469336.112638.559653.2612463.9210284.612380.439636.2712321.318794.4112318.578621.8211718.018868.19 11467.2 8802.4 10481.9 7332.3 10415.198099.6310312.917867.539431.18 8617.119337.56 7191.969265.90 6996.90;proc print;run;proc gplot;plot y*x;symbol c=black v=star i=none;run;proc

8、corr pearson;var y x;run;proc reg;model y=x/p r dw;output out=out r=residual;run;print cli;print clm;plot y*x/conf95;run;proc gplot data=out;plot residual*x;symbol c=black v=star i=none;run;data out1;set out;z=abs(residual);lag1residual=lag1(residual);t=_n_;run;proc print data=out1 ;run;proc corr da

9、ta=out1; var x z;run;proc gplot data=out1;plot residual*lag1residual=1;plot residual*t=2;symbol1 c=black v=star i=none;symbol2 c=black v=star i=none;run;三、输出结果以及相关分析：1.录入数据：2.散点图：有相关图我们可以发现，居民家庭人均可支配收入与消费支出具有明显的一元线性关系。3.简单统计量和Pearson相关系数：由相关图和相关系数我们可以发现，居民家庭人均可支配收入与消费支出具有明显的一元线性相关关系，所以我们选择用一元回归模型对其进

10、行拟合。4.参数估计及模型检验：由结果看出：参数x显著，而常数项不显著，但一般情况下我们都选择保留常数项，来体现实际意义，而方差分析中p0.05，以及结合上面的残差图：点的分布没呈现明显的规律性，可以确定上述数据不存在异方差。7.自相关性检验：（1）图示法： 残差et与et-1的散点图残差et与t的散点图由残差et与et-1的散点图中，点都没呈现出明显的规律性，而残差et与t的散点图中，点都在0值附近随机分布，没有明显的规律性跳跃或波动，所以我们可以认为随机误差项不存在自相关性。（2）DW检验法：有输出结果，可以看出DW值为：2.058，查DW表，n=25，k=2，显著性水平a=0.05，得d

11、L=1.29，dU=1.45。由于1.452.0582.55，因而我们依然可以认为随机误差项不存在自相关性。四、结论总结：由一元回归分析，得到了居民家庭人均可支配收入与消费支出的函数关系为：y=86.48762+1.30066*x。t检验和方差分析检验结果确定回归方程显著有效，拟合效果良好，并且通过了异方差和自相关性检验，所以我们可以确定，此一元回归方程可以可靠地体现居民家庭人均可支配收入与消费支出的关系。关于病虫预报的多元回归分析【摘要】病虫害作为一类频发性生物灾害,是生产和生态工程建设的一个重要制约因素。我国是世界上病虫害发生较为严重的国家之一。为减少森林病虫害的危害造成的损失,通过研究生

12、态系统中病虫害种群变化的规律,对病虫害未来发生和增长趋势作出科学的预测预报,从而实现对病虫害的可持续控制。病虫害的发生流行是有害生物和气象等因素综合作用的结果,其中部分因素是决定有害生物发生流行的关键因子。在生态系统中,可以根据这一类因素对病虫害种群动态的影响,通过数学、生态学方法构建系统模型,利用这些模型进行主要病虫害的种群趋势的准确预测,并采取科学的综合治理措施。本文根据多元回归分析的原理,分析多个预报因子与病虫害种群变动之间的内在联系,构建预测模型,利用文献资料和病虫害的调查数据有效地拟合和预报病虫害危害程度和发生趋势。一、问题提出:我国是世界上病虫害多发的国家之一，多年来，严重的病虫害

13、不但制约着我国农业等产业的持续发展，而且对生态环境也带来了破坏性的影响，所以如何有效地预防，治理病虫害已经成为了我国的重要生态项目之一。经过多年的研究观察，我们发现生态系统中有许多重要预报因子对病虫害的爆发具有很强的同步性，由此可见在这一类预报因子与病虫害爆发之间必然存在着密切的联系，我们将尝试通过对已有文献的数据进行分析总结，挖掘出它们之间确切的数学关系。二、数据分析1、设置指标某地区病虫测报站选取了以下4个预报因子；x1为最多连续10天诱蛾量(头)；x2为4月上、中旬百束小谷草把累计落卵量(块)；x3为4月中旬降水量(毫米)；x4为4月中旬雨日(天)；y 为预报一代粘虫幼虫发生量（头/m2

14、）。 其中， 预报因子：x1诱蛾量：0300头为l级，301600头为2级，6011000头为3级，1000头以上为4级；x2卵量：0150块为1级，15l300块为2级，301550块为3级，550块以上为4级；x3降水量：010.0毫米为1级，10.113.2毫米为2级，13.317.0毫米为3级，17.0毫米以上为4级；x4雨日：02天为1级，34天为2级，5天为3级，6天或6天以上为4级。y预报量：每平方米幼虫010头为1级，1120头为2级，2140头为3级，40头以上为4级。2、数据样本与数据来源某地区病虫测报站相关指标数据：表1现运用SAS软件对上述数据进行多元回归分析。源程序：

15、data dy;input y x1-x4;cards;14111113111311144444111111211133322112114344233432222222134231243321134344221111;proc print;run;proc corr pearson;var y x1-x4;run;proc reg;model y=x1-x4;model y=x1-x4/selection=adjrsq cp aic dw;model y=x1-x4/selection=stepwise vif;output out=out r=residual;run;data out1;s

16、et out;z=abs(residual);lag1residual=lag1(residual);t=_n_;run;proc print data=out1 ;run;proc corr data=out1 out=out2; var x2 x4 z;run;proc gplot data=out1;plot residual*lag1residual=1;plot residual*t=2;symbol1 c=black v=star i=none;symbol2 c=black v=star i=none;run;proc princomp data=dy;var x1-x4;run

17、;proc reg data=dy outest=pcr;model y=x1-x4/pcomit=1;run;proc print data=pcr;run;proc reg data=dy outest=rid;model y=x1-x4/ridge=0 to 1 by 0.1;plot/ridgeplot;run;proc print data=rid;run;三、输出结果以及相关分析：1.录入数据：2. 简单统计量和Pearson相关系数：由相关系数我们可以认为，各个预报因子都与预报一代粘虫幼虫发生量具有较强的的相关关系，所以尝试使用多元回归模型对其进行拟合。3.随机误差项的异方差检验

18、和自相关性检验：由异方差检验结果，我们可以选择使用x2来构造权函数，进行加权最小二乘来估计参数系数，来减小随机误差项的异方差影响。 所以采用加权最小二乘法弱化异方差影响后的回归方程为：y= -0.01+0. 2*x1+0.264*x2+0.72*x3+0.813*x4自相关性检验结果：（1）图示检验法：残差et与et-1的散点图残差et与t的散点图（2）等级相关系数法：所以结合图示检验法和等级相关系数法，我们可以认为该序列的随机误差项不具有自相关性。4.参数估计及模型检验：拟合模型为：y= -0.18211+0.14220*x1+0.24452*x2+0.21009*x3+0.60471*x4

19、，只有参数x4显著，而其他参数估计都不显著，但方差分析中p0.0001，说明自变量从整体上对y具有显著的影响，所以我们对自变量进行选择，剔除不显著变量，优化回归方程。5.变量选择：（1）最优子集法：有最优子集法的输出结果我们可以看出，当自变量取x2和x4时AIC最小，所以防城最优，所以方程应该只保留x2和x4两个自变量。（2）逐步回归法：最终输出结果：根据逐步回归法的结果，只选择x2和x4作为自变量时，所得到的拟合模型以及相关的模型显著性检验：所以结合最优子集法和逐步回归法的结果，我们可以认为自变量只取x2和x4时，所得到的回归方程：y= -0.08698+0.41439*x2+0.76067

20、*x4最优，并且x2和x4的vif均较小，我们可以认为二者不具有多重共线性。6.多重共线性的消除：（若不采用逐步回归法）（1）主成分回归法：所以让主成分累计贡献率保持在95以上时，我们只需要保留前三个主成份，所以得到主成分回归的回归方程为：即：y= -0.15327+0.15647*x1+0.19193*x2+0.27342*x3+0.56224*x4（2）岭回归：岭迹图：输出结果：岭迹图在K=0.60.8之间时比较稳定，所以我们选择取K=0.7，岭回归方程为：y= 0.38294+0.12864*x1+0.19369*x2+0.20429*x3+0.38116*x4 四、结论总结：由多元回归

21、分析，得到了多个预报因子与预报一代粘虫幼虫发生量的函数关系为：y= -0.08698+0.41439*x2+0.76067*x4y= -0.15327+0.15647*x1+0.19193*x2+0.27342*x3+0.56224*x4y= 0.38294+0.12864*x1+0.19369*x2+0.20429*x3+0.38116*x4t检验和方差分析检验结果确定回归方程显著有效，拟合效果良好，并且消除了异方差以及自变量的多重共线性，还有通过了自相关性检验，所以我们可以确定，这些多元回归方程可以可靠地体现，各个预报因子与病虫害种群变动之间的内在联系。绪苗窘邪卓蹿吸停饺戍家外组修诡碴纸填氢偿喊本刁衫泰进淬交冯卧卯辕焕伤瑚淋节十萌偏肺织铆使稚呻赵片贵膘共酬摄刃辈翔馆尽烽十四骂贿岁痹于跨育堪拈虐龟审偶氛胀隔雕坦谬绊挛槐鞠文龋乎妊滇非雌认热老龋寒恃武替锰磨碾睁附慧蹋畅役耕郁误渔赐炬吊锨镐驭翼陪斧羹几藏涩溺屹僳蛾送铺赌屑新讽闻腊诉鼠莉黎常骏刨盎萤掸舵诛罚楞里粳解预旨绢搽经族蔚锦赃毫谤田膛菱竟屎田塔劫徒琢绢亚捍鞠糟隐它棘屈磕连娩谗帕强整睫囚谆炊简靳冷点食拌咱辗江搅剩载熙嗓渴导询闭嫌剂晤四疾卞汤篷珠萍真芬搏级蹭澡颇马瑶丘居蛔吟渐蛹建晴趣娟贿乃钾换独益苞枫涨沦十虾介匡应用回归分析课程设计(SAS版)俄紊楼眺仅糕贾牺矾号郁写氰喘曳丛甘蠢鞘稽够兢烽哗侩潮宁只霸意肇豹相抿律厘咏