时间序列ARIMA期末论文-ARIMA模型在总人口预测中的应用

=1\*ROMANIARIMA模型在总人口预测中的应用【摘要】人口发展与社会经济的发展是密不可分的,研究我国总人口的发展,对我国人口数进行分析和预测,有利于及时控制人口的增长调节人口平衡,利于政府及时了解发展趋势并做出反应对策使我国人口发展步入健康的轨道。本文利用时间序列建模原理和思路，并结合R3.2.1软件对1962年——2014年我国年底总人口数据做分析和预测。找到对原始数据有着较好的拟合度和较高的预测精度的模型。利用此模型可对我国年底总人口进行合理的预测。【关键词】ARIMA建模总人口人口预测目录TOC\o"1-3"\u一、引言 3TOC\o"1-3"\u1.1研究背景 31.2研究现状 4二、模型建立 5TOC\o"1-3"\u2.1模型识别 52.2模型的参数估计 82.3模型的诊断 102.模型的预测 12模型的优缺点及推广 133.1模型的优缺点 133.2模型的推广 13结束语14【参考文献】15附录16引言1.1研究背景我国是世界上人口最多的国家，自1980年开始，年末中国大陆总人口就已经超过了10亿，并一直保持约占世界总人口的五分之一，亚洲人口的三分之一。中国人口的发展同中国社会的发展一样经过了漫长而曲折的道路。在世纪的进程中，目前我国进入了一个全新的时代，要想在21世纪——这个充满竞争与挑战的时代中变的富强、屹立于世界民族之林，实现我们的中国梦，这全取决于人。能否顺利解决人口现状等问题，是我国乃自世界共同面临的问题，由于地球的资源是有限的，它不可能无限制的容纳人口，当人口过多，会由于经济跟不上，工作岗位欠缺，医疗等水平不足，从而导致整个社会处于一种动荡之中；然而如果人口过少，又会由于人员不足，导致各方面人力资源不足，无法正常完成各项必须社会活动，这也会极大地限制一个国家的发展，因此，对人口的研究是具有相当的意义的。我国由于幅员广阔，民族众多，各民族发展水平不一，同时作为世界第一人口大国，我国的耕地面积却相对不足，因此我国每年都需要从国外大量进口粮食，由于过分依赖于进口这对我国的发展影响巨大，为此甚至有国外反华势力叫嚣只要断绝给中国供粮，三五年之内中国必定大乱。当然那只是敌对势力的一厢情愿与恶意诋毁，但我们自己却必须认识到在由于人口的问题而导致的一系列问题，关于人口问题我国必须重视,并根据其趋势做出反应对策。因此，认真分析我国当前人口现状，从中发现其变化的趋势，并对未来总人口进行短期预测，及时采取必要的政治及经济措施来解决人口发展问题，对树立未来的发展目标很有必要。总之，人口是构成社会的主体，在我国社会主义现代化建设中，人口问题始终是极为重要的问题，而人口问题的本质是发展问题。人口发展与社会经济的发展也是密不可分的。基于此，我们利用时间序列中的ARMA模型对我国人口进行预测，对人口的控制起到指导作用，有利于政府采取必要的政治及经济措施来进行调控。所以，对其进行分析和测试是非常有意义的工作。1.2研究现状在对人口问题的研究上，国内外学者做了相当多的工作。在国内程等利用自限模型对我国的人口增长进行了预测，认为中国在2010年-2019年人口数依次会缓慢增加，2016年突破14亿大关，且未来15年人口净增加量不会超过1亿；蒋慧基于多元统计模型对广西人口增长进行了分析，得出了人口增长的综合因子，并提出了稳定人口增长的建议；丁明等运用相空间重构神经网络模型对我国人口增长进行预测，很好的解决了非线性的问题，为我国人口增长预测提供了一种新的方法；王保等，利用Logistic模型进行人口预测，并检验了2005年—2007年的数据误差，取得了理想的效果。在国外，\o"SearchforStocks,B.J."\o"SearchforRosenF"Rosen利用Malthusian模型对人口进行研究，也取得了不错的成绩。本文基于时间序列在研究时间相关问题上的优势，以1949年-2014年的年末总人口数据，利用时间序列知识建模，找到适合人口增长的过程的时间序列模型（模型识别），然后利用参数估计估计出模型的参数（参数估计），再对模型进行诊断，判断模型的好坏（模型诊断），最后利用已经建立的模型对未来的给定的时间进行预测（预测）。\o"SearchforStocks,B.J."\o"SearchforStocks,B.J."模型建立2.1模型识别\o"SearchforStocks,B.J."首先需要对数据进行预先的处理，观察其时间序列图是否为平稳序列，可以用R软件完成（具体程序见于附录1），可得到时间序列图为：\o"SearchforStocks,B.J."图1：人口时间序列图由图一可以看出，年底总人口数随着时间增加的同时也在逐年上涨，有着明显的上升趋势。因此可得出这列数据是不平稳的、方差也是不平稳的结论。并且数据大致是呈线性变化的，因此可以考虑做差分变\o"SearchforStocks,B.J."换。先对数据做一阶差分变换后再观察序列是否平稳，可运用R软件编程得其变换后的序列及变换后序列的自相关图和偏自相关图(具体程序见附录2)：图二：一阶差分序列图及其ACF、PACF图由图二可知，对原始数据在进行一阶差分之后的时间序列图显示并不算平稳，并且一阶差分后的ACF图仍呈现出近似直线下降趋势，因此可以考虑再做一次差分运算。通过编程可以得到（具体程序见附录3）：图三：对数二阶差分序列图及其ACF、PACF图再观察其二阶差分后的时间序列图，基本上趋于平稳，而其自相关图(ACF)和偏自相关图（PACF）呈现出明显的拖尾形式，且PACF图在滞后6阶比较显著，ACF在滞后1,3,6阶相对较显著，由此可以认为原序列基本上可以用ARIMA(1,2,1)，ARIMA(1,2,2)，ARIMA(1,2,3)，进行拟合。对ARIMA(1,2,1)模型序列满足：Yt+Yt-2-2Yt-1=∅Yt-1+Y由=1\*GB3①式可以得到：Yt=2+∅Yt-1+1-2∅Yt-2-∅Yt-3同理可以得到ARIMA(1,2,2)模型：Yt=2+∅Yt-1+1-2∅YARIMA(1,2,3)的模型为：Yt=2+∅Yt-1+1-2∅Y2.2模型的参数估计根据2.1节，已经找到几个可能用于拟合的模型，分别为ARIMA(1,2,1)，ARIMA(1,2,2)，ARIMA(1,2,3)，那么接下来就应该对其进行参数估计，在对模型进行参数估计时有多种方法可以选择，这里选择用最小二乘估计，通过R编程序可以得到如下结果（具体程序见于附录4）：Call:arima(x=x,order=c(1,2,1),method="CSS")Coefficients:ar1ma1-0.47940.5150s.e.0.10330.1901sigma^2estimatedas21752:partloglikelihood=-327.05Call:arima(x=x,order=c(1,2,2),method="CSS")Coefficients:ar1ma1ma2-0.52160.8442-0.3147s.e.0.00090.09120.1014sigma^2estimatedas16939:partloglikelihood=-320.67Call:arima(x=x,order=c(1,2,3),method="CSS")Coefficients:ar1ma1ma2ma3-0.52030.7646-0.11110.3405s.e.0.00080.13270.14350.1278sigma^2estimatedas15191:partloglikelihood=-317.89由上面的输出可以确定在=2\*GB3②式中的系数，并且由估计值的误差项se判断该系数是否显著（通过判定在0是否在区间[∅-2se,∅+2se]内，若在则系数不显著，如不在则显著）。从而代入=2\*GB3②可以确定ARIMA(1,2,1)模型为：Yt=1.5206Yt-1-0.0412Yt-2-0.4794Y类似代入=3\*GB3③式得ARIMA(1,2,2)表达式为：Yt=1.4784Yt-1+0.0432Yt-2-0.5216代入=4\*GB3④式得ARIMA(1,2,3)模型的表达式：Yt=1.4797Yt-1+0.0406Yt-2-0.52032.3模型的诊断在对模型完成了识别和参数估计之后，需要对模型进行诊断，诊断模型是否具有对原时间序列数据的很好的拟合效果。主要进行拟合模型的残差分析和分析过度参数化；对模型进行过度参数分析主要看在进行差分时是否出现过度差分的状况，而对残差进行分析主要需要做以下几个方面的工作：=1\*GB2⑴，检验残差是否是随机的，一个模型如果能很好的拟合，那么拟合后的残差基本上是随机的，残差是应该围绕在某条平行于x=0这一条直线上下波动的，并且波动的幅度不会很大，这可以用做残差序列图观察得到；=2\*GB2⑵，检验残差是否呈正态性，一个模型如果能很好的拟合，那么其残差应该是呈正态性的，这里用残差QQ图和S-W正态性检验（原假设为：H0:数据是呈正态性的）进行；=3\*GB2⑶，判断残差之间是否是相互独立的，一个模型如果能很好的拟合，那么其残差之间相对是比较独立的，这里主要用残差的自相关序列图和L-B检验（原假设为：H0:原数据的残差之间是不相关的）进行。诊断ARIMA(1,2,1)模型，用R软件编程序（具体程序见于附录5）输出为：图4：ARIMA(1,2,1)模型残差的序列图、ACF图和QQ图Shapiro-WilknormalitytestW=0.93417,p-value=0.00592Box-LjungtestX-squared=23.846,df=23,p-value=0.4124由图4中的残差序列图可以看出残差基本基本上是围绕x=0这条直线上下波动的，因此可以认为ARIMA(1,2,1)模型拟合满足残差是随机的条件；又由图4中的QQ可以看出残差基本上是集中在一条直线上的，由S-W检验的的p=0.00592<0.05，检验也可以认为残差是非正态的；又在残差的自相关图中，只有滞后二阶是是显著的，因此可以认为残差之间基本上也是不相关的，特别由B-L检验的p=0.4121>0.05，因此没有充分的理由拒绝原假设，应该认为残差是相互独立的。到此就已经对模型行进行了诊断，由于残差正态性不足，说明用ARIMA(1,2,1)模型拟合原数据不是十分合适。同理可对ARIMA(1,2,2)和ARIMA(1,2,3),进行诊断，这里图形和数据检验就不再一一呈现在论文中。在诊断中发现模型ARIMA(1,2,2)用于拟合是满足条件的，而ARIMA(1,2,3)模型也由于残差正态性不足而不适合用于拟合合。因此在对人口模型进行预测时采用ARIMA(1,2,1)模型。并且在进行参数冗余分析时发现，对于该时间序列用ARIMA(1,2,2)模型拟合后，对ARIMA(p,d,q)中p,d,q任意一个变小都不能再满足条件，因此用ARIMA(1,2,2)模型是适合且简化的，因此后文选用ARIMA(1,2,2)模型进行预测。2.4模型预测由2.1,2.2,2.3三节已经完成了对模型的识别、参数估计和模型诊断，现在就需要运用该模型对原时间序列趋势进行预测。在实际情况中，过去的状态已经是即成的事实，我们关心的更多的是未来是什么状态，从而就可以判断在目前的情况未来的情况，进而可以进行相应的措施予以应对。这里假设对未来5年进行预测，编程[7]得到（具体程序见附录6TimeSeries:Start=2015End=2019Frequency=1[1]137487.2138208.4138936.0139666.2140397.4TimeSeries:Start=2015End=2019Frequency=1[1]162.5725305.4282518.2338796.22561129.1609由上输出，可知由模型ARIMA(3,2,6)预测的2015年-2019年年末总人口数分别为：137487.2,138208.4,138936.0,139666.2,140397.4（单位：万人）。并且还可以对人口数进行区间估计，在此就不列出。\o"SearchforStocks,B.J."模型优缺点及推广3.1模型优缺点优点：ARIMA模型在对许多时间序列都适用，并且在建模过程中有多种方法可以选择，并且在模型的诊断中能过对模型的好坏进行评价，对预测起到很好的作用。缺点：ARIMA模型对数据进行预测时，对短期预测效果还是比较好的，但随着时间的延长，其预测误差就比较大了。3.2模型的推广ARIMA模型对大多以时间为变量的趋势都有比较好的拟合效果，可广泛运用于气温、股票、降水量、商品价格等的研究之中。结束语人口问题是人类社会伴随始终的问题，在新中国成立之初，由于对人口问题的认识不足，导致我国大量的人员过剩，从而影响到后面几十年的发展，给教育、医疗社会带来了巨大的压力，即便后来在意识到问题的严重性后，在计划生育政策下，人口得到了有效的控制，然而由于人口基数过大，人们传统观念的根深蒂固，我国人口还是处于不断增长的时期，最近几年又由于伴随着人口老龄化情况加剧，社会经济压力又有了很大的压力，这对我国的经济在一定程度上有着不良影响，慢慢政府在对计划生育问题上又有了一定放松，在一定条件下允许第二胎。这一系列的重大举措反应了当代人对人口问题的重视。本文对未来几年利用模型进行了预测，对政策是一定的补充说明。论文结束之际，在此向各位老师表示感谢，对在相关领域做出贡献的学者表示感谢。参考文献：[1]程华，高孝成，幺焕民《我国人口增长预测问题的自限模型》，[J],哈尔滨师范大学自然科学学报，2009年，第25卷第3期：34-36[2]蒋慧萍，《基于多元统计模型的广西人口增长分析》，[J],企业科技与发展，2010年，第10期：190-192[3]丁明磊，杨晓娜，曹连海《相位空间重构神经网络模型在我国人口增长预测中的应用》，[J],华北水利水电学院学报，2008年6月，第29卷第3期：95-98[4]王保学，蔡果兰，《Logistic模型的参数估计及人口预测》[J]，北京工商大学学报（自然科学版），2009年11月，第27卷第6期：75-79[5]\o"SearchforRosenF"RosenF,《Theprincipleofpopulationaspoliticaltheory:Godwin'sOfpopulationandtheMalthusiancontroversy》,[J],\o"SearchforJournalOfTheHistoryOfIdeas"JournalOfTheHistoryOfIdeas[JHistIdeas]1970Jan-Mar;Vol.31(1),pp.33-48.[6]潘红宇等译，《时间序列分析及应用》[M]，机械工业出版社，2011年1月[7]薛毅，陈丽萍，《统计建模与R软件（下册）》[M]，清华大学出版社附录1、做时间序列的序列图plot(x)2、做序列的一阶差分时间按序列图及处理后的ACF图,PACF图程序：mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)x<-ts(mydata,frequency=1,start=1962)par(mfrow=c(3,1))plot(x)plot(diff(x,difference=1))#做差分后的时间序列图acf(diff(x,difference=1),ci.type='ma')#做差分后的自相关图pacf(diff(x,difference=1))#做差分后的偏自相关图3、做序列的二阶差分时间按序列图及处理后的ACF图,PACF图程序：mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)x<-ts(mydata,frequency=1,start=1962)par(mfrow=c(3,1))plot(diff(x,difference=2))#做差分后的时间序列图acf(diff(x,difference=2),ci.type='ma')#做差分后的自相关图pacf(diff(x,difference=2))#做差分后的偏自相关图4、估计模型中的参数程序：mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)x<-ts(mydata,frequency=1,start=1962)arima(x,order=c(1,2,1),method="CSS")#CSS代表条件最小二乘估计arima(x,order=c(1,2,2),method="CSS")#CSS代表条件最小二乘估计arima(x,order=c(1,2,3),method="CSS")#CSS代表条件最小二乘估计5、对模型的残差进行检验程序：library(TSA)mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)x<-ts(mydata,frequency=1,start=1962)m<-matrix(c(1:3,2),2,2)layout(m);layout.show(3);model=arima(x,order=c(1,2,1))#当需要改诊断其他模型时就改变此处的p,d,q值即可plot(rstandard(model,infl=iflSR),ylab='StandardizedResiduals',type='o');abline(h=0);#画残差图qqnorm(residuals(model))#画残差的QQ图qqline(residuals(model))#画残差图QQ图shapiro.test(residuals(model))#对残差做正态性检验acf(residuals(model))#对残差做自相关图LB.test(model,lag=25)#求残差Ljung-Box统计量6、对未来5年年末人口数的预测程序：library(TSA)mydata<-read.table("C:/Users/Administrator/Desktop/shuju.txt",header=T)x<-ts(mydata,frequency=1,start=1962)model=arima(x,order=c(1,2,2))#确定模型predict(model,c(5))#对未来五年在model模型下的预测基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统基于51单片机的嵌入式Internet系统的设计与实现单片机监测系统在挤压机上的应用MSP430单片机在智能水表系统上的研究与应用基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用单片机在高楼恒压供水系统中的应用HYPERLINK"/deta