1. 引言
自二十世纪八十年代改革开放以来,我国对外贸易经济得到了飞跃式发展,目前全国外贸进出口总额居世界第二位。而截至2019年上海市成为唯一一个进出口超过3万亿的城市,以实现进出口总额34046.8亿排名第一,在带动我国整个进出口外贸发展浪潮中占据着重要地位。因此,对上海市进出口贸易进行研究并对进出口总额进行预测,这对分析我国未来几年的对外贸易发展起着重要作用。此外,外贸进出口货物的重要运输方式除陆路运输之外,水路运输同样占据着重要地位,因此进出口外贸的快速发展将会对城市港口货物吞吐量起着极大的提升作用,促进了城市港口繁荣发展。
现有研究中,大多只是单独预测对外贸进出口总额 [1] [2] [3] 或港口货物吞吐量 [4] [5] [6],而本文通过协整检验验证两者之间的关系,根据所预测的进出口总额,结合所拟合的协整模型来预测未来5年上海港货物吞吐量。结果显示,两者之间确实存在协整关系,因而先对未来5年上海市进出口总额进行预测,同时根据已有数据计算了预测误差率,结果显示所拟合的ARIMA(3,2,0)拟合效果较好。因此在此基础上对上海港货物吞吐量进行预测,预测结果显示未来5年上海货物吞吐量将随之进出口总额快速发展而呈上升趋势,同时所计算出的预测误差率较小,这表明我们所拟合的模型精准度较高。这有助于我们更好地了解上海市未来进出口总额和港口货物吞吐量发展趋势,并且本文验证了外贸进出口的发展与港口货物吞吐量之间存在的协整关系,这能为政府预测未来港口货物吞吐量发展提供一定参考意见。
2. ARIMA模型介绍
ARIMA模型最早是由Box和Jenkins于上世纪70年代所提出,该模型也通常被称为求和自回归移动平均模型 [7]。ARIMA 模型实质上就是差分运算和ARMA模型的组合,即通过差分运算首先将非平稳序列平稳化,再用平稳化序列拟合自回归移动平均模型(ARMA) [8] [9]。而ARMA分为AR (自回归模型)、MA (滑动平均模型)和 ARMA (自回归滑动平均模型)这三类模型 [10]。
AR模型描述的是当前值和历史值之间的关系,其表达式为
,其中,p为自回归模型的阶数,
为待定系数,
为误差,
为平稳时间序列。MA模型是根据所研究的对象过去所产生的干扰值与现在所产生的干扰值之间的一种线性关系组合进行预测 [7],其基本表达式为
,其中,q是为该模型的阶数,
为待定系数,
是t之前q期的随机干扰值。ARMA模型AR和MA模型的组合,其表达式为
,其中p、q、
、
、
、
与前述的AR、MA 模型中含义相同。
3. 实证分析
3.1. 数据预处理与初步分析
为研究上海市近年来的外贸进出口总额的波动特征与规律,本文选取了1978~2019年的年度数据,共计42个观测值。下列图1是上海市年进出口总额的时序图,我们从中可以较为清晰地发现,该序列蕴含着比较明显的线性递增的变化趋势,处于非平稳状态。因此,我们通过对该序列进行二阶差分运算来从中提取出线性趋势信息,提取结果如图2所示,可以发现该序列在经过二阶差分处理后呈现出较为平稳的波动趋势特征,这就表明该二阶差分处理过程比较准确充分地提取出了在原序列中可能蕴含着的某种长期趋势,使得经过处理后的序列不再继续呈现出确定性趋势。
同时,本文为研究上海市外贸进出口发展对港口货物吞吐量的影响,选取了1978~2019年间关于上海市港口货物吞吐量的数据,总计为42个观测值。下列图3是上海市进出口外贸总额与上海市港口货物吞吐量对数序列时序图。时序图显示,这两个序列都呈现出显著的线性递增变化的时间趋势,所以两者都是非平稳的时间序列。但是这两个序列之间又随着时间的逐步增长,这两序列之间的变化趋势逐渐趋同,所以为了有效且客观地衡量这两个时间序列之间是否存在长期稳定的相对均衡关系,本文首先通过构造回归模型检验序列之间是否具有协整关系。
![](//html.hanspub.org/file/51-2622391x20_hanspub.png?20220623091321709)
Figure 1. Sequence plot of original data
图1. 原始数据序列图
![](//html.hanspub.org/file/51-2622391x21_hanspub.png?20220623091321709)
Figure 2. Timing plot after second order difference
图2. 二阶差分后时序图
![](//html.hanspub.org/file/51-2622391x22_hanspub.png?20220623091321709)
Figure 3. Timing plot of logarithmic sequence of total import and export volume and port cargo throughput
图3. 进出口总额与港口货物吞吐量对数序列时序图
3.2. 序列检验
3.2.1. 平稳性检验和白噪声检验
首先通过观察二阶差分后时序图可以初步判断二阶差分后得到的上海市外贸进出总额序列为平稳,需要进一步借助ADF检验,判断该序列的平稳性。ADF检验结果显示:类型一、类型二和类型三各种模型的τ统计量的P值均小于显著性水平(
),所以认为该序列显著平稳。
其次,并对已经进行了二阶差分运算的上海市外贸进出口总额序列进行白噪声检验,得到如表1所示的检验结果。纯随机性检验得到的统计分析结果显示:在延迟6 阶和12阶下LB检验统计量的P值都非常小,低于所给定的显著性水平0.05,因而能在95%的置信水平下将原有假设拒绝。故可以判断上海市外贸进出口总额序列不属于白噪声序列。
![](Images/Table_Tmp.jpg)
Table 1. White noise test results of total foreign trade import and export
表1. 外贸进出口总额序列白噪声检验结果
3.2.2. 协整检验
为对上海市外贸进出口总额和港口货物吞吐量建立动态回归模型,先对其进行协整检验。图4为上海市外贸进出口对数序列和港口货物吞吐量对数序列的互相关图,图中显示外贸进出口发展对港口货物吞吐量的影响当期达到最大。因此,利用最小二乘法,构建二者之间的回归模型为:
其中,
表示t年上海市港口货物吞吐量,
表示t年上海市外贸进出口总额,
表示残差。然后对回归残差序列进行平稳性检验。根据类型一滞后0阶的检验结果,本文有近似92.5%的把握认为上海市外贸进出口总额对数序列与港口货物吞吐量对数序列之间存在协整关系。并且我们还同时对回归的残差序列进行了白噪声检验,检验结果显示,回归残差序列不是白噪声序列,需要进一步提取这些回归的残差序列中可能会蕴含着的相关信息。
![](//html.hanspub.org/file/51-2622391x28_hanspub.png?20220623091321709)
Figure 4. Cross-correlation plot of logarithmic sequence of total import and export value and port cargo throughput
图4. 进出口总额与港口吞吐量对数序列互相关图
3.3. 模型选择与显著性检验
3.3.1. 模型识别
1) 一元时间序列模型
二阶差分后的序列是一个平稳的非白噪声序列,即ARIMA(p,d,q)模型中的d = 2。尚未能确定模型中的p和q,可以通过观察二阶差分后的时序图中自相关系数和偏自相关系数的截尾性来判断 [6]。根据图5的自相关图所示,首先可以看见自相关系数表现出较为连续地从顶峰逐渐趋于零的趋势,这表明自相关系数是拖尾。其次,根据图6的偏自相关系数的变化过程,可以看出除1至2阶的偏自相关系数波动范围在2倍的标准差之外,其余所有阶数的波动都在2倍的标准差单位之内,这表明偏自相关系数2阶或3阶截尾,所以本文首先将模型定为ARIMA(3,2,0),并对其进行检验。
![](//html.hanspub.org/file/51-2622391x29_hanspub.png?20220623091321709)
Figure 5. Autocorrelatinplot of sequence after second-order difference
图5. 二阶差分后序列自相关图
![](//html.hanspub.org/file/51-2622391x30_hanspub.png?20220623091321709)
Figure 6. Partial autocorrelatinplot of sequence after second-order difference
图6. 二阶差分后序列偏自相关图
2) 协整动态回归模型
观察回归残差序列的自相关图(图7)和偏自相关图(图8),可以判断残差序列自回归系数拖尾,偏自相关系数1阶截尾。因此,本文将残差模型定为AR(1),并对其进行检验。
![](//html.hanspub.org/file/51-2622391x31_hanspub.png?20220623091321709)
Figure 7. Autocorrelation plot of residual series
图7. 残差序列自相关图
![](//html.hanspub.org/file/51-2622391x32_hanspub.png?20220623091321709)
Figure 8. Partial autocorrelation plot of residual series
图8. 残差序列偏自相关图
3.3.2. 模型拟合与显著性检验
1) 一元时间序列模型
经过上述验证,运用R软件拟合模型ARIMA(3,2,0),得到如表2所示各参数估计值。
![](Images/Table_Tmp.jpg)
Table 2. Automatic identification of the estimated values of each parameter of the model
表2. 自动识别模型的各参数估计值
由上表各参数可以得到拟合模型ARIMA(3,2,0),且该模型是基于二阶差分后的数据得来的:
① 残差白噪声检验
对上文所建的模型进行残差白噪声检验,选用Box-Ljung检验统计量,原假设H0:
,
,备择假设H1:至少存在某个
,
。如果仍然能够有效拒绝原假设,就说明在残差序列模型中还可能残存着某些相关信息,拟合序列的模型尚不显著;假若检验结果显示不能够拒绝原假设,就可以认为该拟合的序列模型是显著且有效的。根据图9所显示的检验结果来看,各阶段延迟条件下的白噪声检验统计量中的P值都很显著的大于0.05,因此不能拒绝原假设。所以,我们可以认为在这个拟合模型中的残差序列属于白噪声序列,即表明所拟合的模型显著成立且有效。
如
![](//html.hanspub.org/file/51-2622391x38_hanspub.png?20220623091321709)
Figure 9. White noise test plot of residual
图9. 残差白噪声检验图
② 参数显著性检验
参数显著性检验主要是通过检验模型中的每一个未知参数是否都显著非零,从而尽可能使该模型变得精简。如果发现某一个变量参数检验结果是不显著非零,这表示与该变量参数所对应的那个自变量对模型中的因变量的影响不明显,二者无显著关系,所以该自变量也就可以从所拟合的模型参数重一一剔除,最终模型将由这样一系列参数显著且非零的自变量构成。我们通过对ARIMA(3,2,0)模型中的三个参数进行了显著性检验,结果如表3所示:ar1、ar2、ar3这三个参数检验统计量的P值都显著小于0.05,所以可以判断拟合模型的三参数都显著非零。
![](Images/Table_Tmp.jpg)
Table 3. Parameter significance test
表3. 参数显著性检验
2) 协整动态回归模型
经上述验证,运用R软件拟合模型AR(1),得到如表4所示各参数结果。
![](Images/Table_Tmp.jpg)
Table 4. Estimated values of parameters of the cointergration model
表4. 协整模型各参数估计值
由上表参数可以得到拟合的协整动态回归模型结构如下:
这两个参数的估计值均显著大于2标准差,因此两个参数均显著非零。并据此对拟合的回归模型进行显著性检验,如图10显示,残差序列为白噪声序列,这表明拟合的协整动态回归模型显著成立。
![](//html.hanspub.org/file/51-2622391x41_hanspub.png?20220623091321709)
Figure 10. White noise test plot of Residual series
图10. 残差序列白噪声检验图
为了更加准确地对港口货物吞吐量进行预测,本文将先对上海市未来5 年的外贸进出口总额进行预测,然后将该预测值代入协整模型,从而得到未来 5年上海市港口货物吞吐量的预测值。因此,根据上海市外贸进出口总额对数序列具有的特征,本文对其拟合了ARIMA(1,1,0)模型。运用R软件拟合模型所得参数如表5所示。
![](Images/Table_Tmp.jpg)
Table 5. Estimated value of each parameter of import and export logarithmic series fitting model
表5. 进出口额对数序列拟合模型各参数估计值
由上表参数可以得到上海市外贸进出口总额对数拟合模型如下:
3.4. 序列预测
3.4.1. 外贸进出口总额预测
根据前文的分析,建立ARIMA(3,2,0)模型,预测未来5年上海市外贸进出口总额,结果如图11所示。根据ARIMA(3,2,0)模型预测的趋势图我们可以较为清晰看出,上海市未来5年的外贸进出口总额将保持一个上升趋势,具体的预测分析结果如表6所示。此外,根据可获取的2020年进出口总额,计算出预测误差率。可以发现,预测误差率仅为−0.23%,这表明预测精准度相对较高。
表6显示了上海市于2020~2024年期间的外贸进出口总额的预测数据,及在95%以上置信度前提下预测区间的上下限。从预测的结果上来看,未来五年上海市外贸进出口总额将从2020年的5020亿美元上升至2024年的5673亿美元。
![](Images/Table_Tmp.jpg)
Table 6. The forecasted volume of foreign trade import and export for next five years
表6. 外贸进出口总额向后5期预测值
![](//html.hanspub.org/file/51-2622391x43_hanspub.png?20220623091321709)
Figure 11. Forecast chart of Shanghai’s total foreign trade import and export in the next five years
图11. 未来5年上海市外贸进口出总额预测图
3.4.2. 港口货物吞吐量预测
根据前文拟合的上海市外贸进出口总额序列模型和协整模型对未来连续5年上海市港口货物吞吐量进行粗略地预测,预测结果大致如表7所示。表7显示了有关上海市2020~2024年的港口货物吞吐量的预测数据,及在满足95%以上置信度前提下该预测结果区间数据的上下限。从预测结果上来看,上海市未来5年的港口货物吞吐量将从2020年的73055万吨上升至2024年的82052万吨,从预测图12来看呈上升趋势。此外,通过对比2020年和2021年的真实值与预测值之间的预测误差率,可以发现,预测误差率总体控制在5%以下,这表明该模型拟合优度较好。
![](Images/Table_Tmp.jpg)
Table 7. The forecasted port cargo throughput for next five years
表7. 港口货物吞吐量向后预测5期
![](//html.hanspub.org/file/51-2622391x44_hanspub.png?20220623091321709)
Figure 12. Forecast chart of Shanghai’s port cargo throughput in the next five years
图12. 未来5年上海市港口货物吞吐量预测图
4. 结论
上海市外贸进出口总额自1978年以来一直保持着长期逐年上升趋势,说明正由于国家改革开放之后要求中国加强了对外经济贸易交流及贸易洽谈活动,而使得上海市作为国际金融中心城市之一,在对外开放经济中发挥的辐射作用越来越重要,对外贸易发展迅速,进出口总额增长迅猛。因此该序列处于非平稳状态,由此通过二阶差分后使其平稳化,并建立短期预测模型。通过对上海市未来5年外贸进出口总额的预测发现,依旧呈现上涨趋势,这符合上海市未来发展势态。
而我们通过协整检验则发现,上海市外贸进出口总额的对数序列和港口货物吞吐量的对数序列之间存在着协整关系,这表明外贸进出口的快速发展,能够带动城市港口运输的繁荣。检验分析表明,上海市外贸进出口的长期上升发展带动了港口货物吞吐量同样呈现上升发展趋势。
NOTES
*通讯作者。