---

可下载资源

随着时间的推移，世界人口不断的增长，为了更好地把握世界人口的进展速度与规律。

模型一：logistic模型

logistic模型又称作阻滞增长模型，主要用来描述在环境资源有限制的情况下，人口数量的增长规律。由于一些因素的影响世界人口数量最终会达到一个饱和值。阻滞作用上体现在对的影响上，使得随着年份的增加而下降。若将表示为的函数，则它应是减函数。则有

由于bgistic回归模型就是基于二项分布族的广义线性模型，因此在R软件中，Logistic回归分析可以通过调用广义线性回归模型函数glm()来实现，其调用格式为

Log<一glm(formula,family=binomial,data)其中，formula为要拟合的模型，family=binomial说明分布为二项分布，data为可选择的数据框。

我们将1950年到2100年的人口数据进行录入，并调用glmnet包来进行拟合。

模型二：  AFRIMA模型

时间序列模型可分为段记忆模型和长记忆模型。一般的时间序列分析模型有自回归(AR)模型、滑动平均(MA)模型、自回归滑动平均(ARMA)模型、自回归整合滑动平均（ARIMA）模型等，这些模型主要是短记忆模型。

 arfi(Diut\[,2\]\[)#建立arfima模型

plot(Discnt\[,1\],Dscunt\[,2\])#原始数据

points(Dicount\[,1\]\[1:66\],f$fittedcol="red")#拟合数据

---

---

模型稳定性改进

对数据进行一阶差分使数据更加稳定。

模型三：  ARIMA模型 

 ARIMA模型定义

ARIMA模型全称为差分自回归移动平均模型。是由博克思和詹金斯于70年代初提出的一著名时间序列预测方法，博克思-詹金斯法。其中ARIMA(p，d，q)称为差分自回归移动平均模型，AR是自回归，p为自回归项；MA为移动平均，q为移动平均项数，d为时间序列成为平稳时所做的差分次数。

ARIMA模型的基本思想是：将预测对象随时间推移而形成的数据序列视为一个随机序列，用一定的数学模型来近似描述这个序列。这个模型一旦被识别后就可以从时间序列的过去值及现在值来预测未来值。

ARIMA模型的基本思想是：将预测对象随时间推移而形成的数据序列视为一个随机序列，用一定的数学模型来近似描述这个序列。这个模型一旦被识别后就可以从时间序列的过去值及现在值来预测未来值。

建模过程主要包括：

第一步：自回归过程

令Yt表示t时期的GDP。如果我们把Yt的模型写成

(Y\_t-δ)=α\_1 (Y_(t-1)-δ)+u_t

其中δ是Y的均值，而ut是具有零均值和恒定方差σ^2的不相关随机误差项(即ut是白噪音)，则成Yt遵循一个一阶自回归或AR(1)随机过程。

P阶自回归函数形式写成：

(Y\_t-δ)=α\_1 (Y_(t-1)-δ)+α\_2 (Y\_(t-2)-δ)+α\_3 (Y\_(t-3)-δ)+⋯+α\_p2 (Y\_(t-p)-δ)+u_t

模型中只有Y这一个变量，没有其他变量。可以理解成“让数据自己说话”。

第二步：移动平均过程

上述AR过程并非是产生Y的唯一可能机制。如果Y的模型描述成

Y\_t=μ+β\_0 u\_t+β\_1 u_(t-1)

其中μ是常数，u为白噪音(零均值、恒定方差、非自相关)随机误差项。t时期的Y等于一个常数加上现在和过去误差项的一个移动平均值。则称Y遵循一个一阶移动平均或MA(1)过程。

q阶移动平均可以写成：

Y\_t=μ+β\_0 u\_t+β\_1 u_(t-1)+β\_2 u\_(t-2)+⋯+β\_q u\_(t-q)

自回归于移动平均过程

如果Y兼有AR和MA的特性，则是ARMA过程。Y可以写成

Y\_t=θ+α\_1 Y_(t-1)+β\_0 u\_t+β\_1 u\_(t-1)

其中有一个自回归项和一个移动平均项，那么他就是一个ARMA(1，1)过程。Θ是常数项。

ARMA(p，q)过程中有p个自回归和q个移动平均项。

第三步：自回归求积移动平均过程

上面所做的都是基于数据是平稳的，但是很多时候时间数据是非平稳的，即是单整(单积)的，一般非平稳数据经过差分可以得到平稳数据。因此如果我们讲一个时间序列差分d次，变成平稳的，然后用AEMA(p，q)模型，则我们就说那个原始的时间序列是AEIMA(p，d，q)，即自回归求积移动平均时间序列。AEIMA(p，0，q)=AEMA(p，q)。

通常，ARIMA 模型建模步骤有4个阶段: 序列平稳性检验，模型初步识别，模型参数估计和模型诊断分析。

模型实现

步骤一：识别。找出适当的p、d、和q值。通过相关图和偏相关图可以解决。

步骤二：估计。估计模型周所含自回归和移动平均项的参数。有时可以用最小二乘法，有时候需要用非线性估计方法。(软件可以自动完成)

步骤三：诊断(检验)。看计算出来的残差是不是白噪音，是，则接受拟合；不是，则重新在做。

步骤四：预测。短期更为可靠。

具体来说：

首先，看一下有没有很明显的trend，需不需要differencing之后再建模。

从ACF和PACF的结果来看，序列没有很快地落入虚线范围之内，因此，序列不平稳。对序列进行差分。

画出ACF 和PACF，通过看图来决定用哪个模型（ARMA（p，q），ARIMA之类的）。

从差分后的数据结果来看，ACF在8阶后开始落入虚线范围，PACF在2阶后很快落入虚线范围，因此p=8，q=2，d=1。

xz1=automa(Dist\[1:66,2\],ic=c('bic'),trace=T)#自动查找最优的arima模型

 ARIMA(2,2,2)                    : -1058.701

 ARIMA(0,2,0)                    : -1026.504

 ARIMA(1,2,0)                    : -1049.834

 ARIMA(0,2,1)                    : -1048.83

 ARIMA(1,2,2)                    : -1062.532

 ARIMA(1,2,1)                    : -1050.673

 ARIMA(1,2,3)                    : -1058.723

 ARIMA(2,2,3)                    : Inf

 ARIMA(0,2,2)                    : -1060.99



 Best model: ARIMA(1,2,2)

从残差的ACF结果来看，序列很快稳定地落入虚线范围，模型稳定。

plot(Discnt\[,1\],Disnt\[,2\])

points(c(2016:2100) ,Diunt\[66,2\]+diffinv(xzrcast$mean)\[-1\],col="red")#预测值2015到2100年

从残差的ACF结果来看，序列很快稳定地落入虚线范围，模型稳定。

为了检验预测误差是均值为零的正态分布，我们可以画出预测误差的直方图，并覆盖上均值为零、标准方差的正态分布的曲线图到预测误差上。

points(mst$mids, myst$density, type="l", col="blue", lwd=2)

}

plotForecrrors(xzrrecast $residuals)

参考文献

【1】统计模拟及其R实现，肖枝洪，武汉大学出版社

【2】Logistic模型在人口预测中的应用，阎慧臻，大连工业大学学报，第27卷第4期