本文通过构建VECM模型和TVECM模型对上海铜期货价格和现货价格之间的协整关系进行了分析。实证研究结果表明:在VECM (2)模型中,铜的现货市场与期货市场之间存在着现货向期货的单项协整关系,非均衡误差对铜的现货对数化价格影响更加显著;在滞后2阶的情况下,铜现货对数化价格与期货对数化价格互为格兰杰因果,二者之间存在双向引导关系。在TVECM模型中引入投资者的交易成本后,门限估计值分别为0.032、0.117。通过进一步分析后发现在短期时间内当铜的现货对数化价格与期货对数化价格之间的基差大于0.117时,可以通过铜的期货对数化价格判断现货的价格走势。最后根据研究结论,提出了相应的对策建议。 This paper analyses the co-integration relationship between Shanghai copper futures price and spot price by constructing VECM model and TVECM model. The results of the empirical study show that in the VECM (2) model, there is a spot-to-futures monotonic co-integration relationship between the spot market and the futures market of copper, and the imbalance error has a more significant impact on the spot logarithmised price of copper. At lag order 2, the spot logarithmised price of copper and the futures logarithmised price are Granger-causal to each other, and there is a two-way bootstrapping relationship between them. After the introduction of in-vestors’ transaction costs in the TVECM model, the threshold estimates are 0.032, 0.117. Through further analysis, it is found that in the short-term time when the basis difference between the spot logarithmic price of copper and the futures logarithmic price is greater than 0.117, the price trend of the spot can be judged by the futures logarithmic price of copper. Finally, based on the conclusions of the study, corresponding countermeasures are proposed.
本文通过构建VECM模型和TVECM模型对上海铜期货价格和现货价格之间的协整关系进行了分析。实证研究结果表明:在VECM (2)模型中,铜的现货市场与期货市场之间存在着现货向期货的单项协整关系,非均衡误差对铜的现货对数化价格影响更加显著;在滞后2阶的情况下,铜现货对数化价格与期货对数化价格互为格兰杰因果,二者之间存在双向引导关系。在TVECM模型中引入投资者的交易成本后,门限估计值分别为0.032、0.117。通过进一步分析后发现在短期时间内当铜的现货对数化价格与期货对数化价格之间的基差大于0.117时,可以通过铜的期货对数化价格判断现货的价格走势。最后根据研究结论,提出了相应的对策建议。
VECM模型,TVECM模型,协整检验,铜期现货价格
—Based on the VECM and TVECM Models
Chongbiao Tang, Yongmei Ding*, Jiaxiong Yu
College of Science, Wuhan University of Science and Technology, Wuhan Hubei
Received: Nov. 13th, 2023; accepted: Dec. 19th, 2023; published: Dec. 28th, 2023
This paper analyses the co-integration relationship between Shanghai copper futures price and spot price by constructing VECM model and TVECM model. The results of the empirical study show that in the VECM (2) model, there is a spot-to-futures monotonic co-integration relationship between the spot market and the futures market of copper, and the imbalance error has a more significant impact on the spot logarithmised price of copper. At lag order 2, the spot logarithmised price of copper and the futures logarithmised price are Granger-causal to each other, and there is a two-way bootstrapping relationship between them. After the introduction of investors’ transaction costs in the TVECM model, the threshold estimates are 0.032, 0.117. Through further analysis, it is found that in the short-term time when the basis difference between the spot logarithmic price of copper and the futures logarithmic price is greater than 0.117, the price trend of the spot can be judged by the futures logarithmic price of copper. Finally, based on the conclusions of the study, corresponding countermeasures are proposed.
Keywords:VECM Model, TVECM Model, Co-Integration Test, Futures and Spot Prices of Copper
Copyright © 2023 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
铜是一种重要的基础工业金属,它广泛应用于建筑、电力、机械制造、交通运输和家电等多个重要行业中,铜对中国这个电力消费大国至关重要 [
在国内,使用向量误差修正模型(Vector Error Correction Model,以下简称VECM)和门限向量误差修正模型(Threshold Vector Error Correction Model,以下简称TVECM)研究协整关系的相关研究已经逐渐增多,从不同角度对协整关系进行了探讨。其中康赞亮等(2006)基于协整理论建立了VECM模型,研究了外部投资、国际贸易投资与中国经济增长之间存在的长期均衡关系,通过实证分析表明中国经济增长与外部投资之间存在着显著的双向因果关系,但影响程度存在差异 [
为了获取铜期货和现货价格协整关系,本文将VECM模型和TVECM模型结合使用,目的是为了解这两种价格的波动、行为等方面是否存在关联性,以及不同变量之间的长期均衡关系是否存在。VECM模型是一种多元回归模型,用于研究不同变量之间的因果关系。在铜期货和现货价格协整分析中,使用VECM模型以期望有效地分析铜期货价格与现货价格的长期均衡关系,推算铜期货价格对现货价格的影响,以及两者之间的波动幅度。TVECM模型针对现有时间序列数据的非线性特征,引入了阈值,将数据分为不同的区域,然后在每个区域中估计参数。在铜期货和现货价格协整分析中,使用TVECM模型可以更好地探究变量之间非线性关系的变化和门限,以了此解铜的期货和现货的价格波动的周期性和频率,预测价格变化方向,以及识别市场变化趋势,并为投资者提供稳健的投资决策依据。
为了研究铜的期货价格与现货价格的长期均衡关系,本文以沪铜主力的日收盘价作为铜的期货价格,(数据来源于东方财富网),以上海阴极铜的日结算价作为铜的现货价格(数据来源于万得资讯),数据的时间跨度为2019年9月25日至2022年11月30日。由于期货交易并不是与现货交易同步,铜的期货价格与现货价格并不是日日对应,所以需要进一步筛选相互匹配的数据,经过筛选、去除缺失值后共得到445对有效数据。
图1. 铜的现货价格与期货价格趋势对比图
本文研究变量的命名与数据分布特征见表1。
变量 | 符号表示 | 样本数 | 均值 | 方差 | 最大值 | 最小值 |
---|---|---|---|---|---|---|
铜的现货价格 | s | 445 | 50822.49 | 5256.33 | 59,483 | 35,135 |
铜的期货价格 | f | 445 | 45108.49 | 5056.81 | 52,800 | 34,390 |
表1. 铜的现期货价格分布特征表
通过表1可以发现铜的期货价格与现货价格差异较大,并且铜的现货价格均值要高于铜的期货价格。为了更好地直观展示铜的期货价格与现货价格之间的价格趋势,利用R进行了数据可视化,如图1所示。
通过图1可以发现,铜的期货价格与现货价格波动较大,从整体趋势来看,铜的现货价格高于铜的现货价格。为了降低数据的波动和可能存在的异方差性,并提升后续建模的效果,本文对铜的期货价格和现货价格均进行对数化处理,分别记为lnf、lns,后续的建模均使用对数化后的数据。
在进行时序建模分析时,首先需要对时序数据进行平稳性检验。平稳性是指时间序列的统计特性在时间上保持不变,即时序图形中的均值、方差和自相关系数等不随时间而变化。如果时间序列数据不平稳,可能会导致误差项存在自相关性,模型无法准确描述历史数据和预测未来数据。因此平稳性检验是一项必要的步骤,可以帮助我们确定是否需要进行平稳性转化。常用的平稳性检验方法有绘制时间序列的时序图、自相关图和偏自相关图和单位根检验,本文选用单位根检验中的ADF检验方法对变量进行单位根检验。
ADF单位根检验共有三个模型,分别为:
模型1不包含趋势项和截距项: Δ X t = δ X t − 1 + ∑ i = 1 m β i Δ X t − 1 + ε t (1)
模型2包含截距项: Δ X t = α + δ X t − 1 + ∑ i = 1 m β i Δ X t − 1 + ε t (2)
模型3包含截距项和趋势项: Δ X t = α + θ t + δ X t − 1 + ∑ i = 1 m β i Δ X t − 1 + ε t (3)
在上面三个模型表达式中,其中 α 表示截距项, θ t 表示趋势项, ε t 表示独立同分布的扰动项, t 表示时间变量。
ADF检验的假设为:
H 0 : δ = 0 ,即存在一单位根,说明序列平稳
H 1 : δ ≠ 0 ,即不存在一单位根,说明序列不平稳
对铜的现货对数价格 ln s 、期货对数价格 ln f 和差分后的现货对数价格 Δ ln s 、差分后的期货对数价格 Δ ln f ,分别按模型3、模型2、模型1依次做ADF检验,何时检验拒绝原假设,即原序列不存在单位根,为平稳序列时,停止单位根检验,最佳的滞后阶数以赤池信息(AIC)准则为准,检验结果见表2。
通过表2可以看出,铜的现货对数价格 ln s 、期货对数价格 ln f ,经过ADF检验可以发现P值均大于0.05,故不能拒绝原假设,即序列 ln s 、 ln f 均为非平稳序列。而一阶差分序列 Δ ln s 、 Δ ln f 的ADF检验P值均为0.002,小于显著性水平0.05,故拒绝原假设,认为一阶差分序列 Δ ln s 、 Δ ln f 为平稳序列。通过ADF单位根检验,可以发现铜的现货对数价格 ln s 与期货对数价格 ln f 是一阶单整序列,可以继续进行协整检验。
序列 | 变量 | 检验模型 | ADF检验t统计量 | P值 | 结论 |
---|---|---|---|---|---|
原序列 | ln s | 模型3 | −1.448 | 0.148 | 不平稳 |
ln s | 模型2 | 0.144 | 0.886 | 不平稳 | |
ln s | 模型1 | −1.705 | 0.089 | 不平稳 | |
ln f | 模型3 | −1.729 | 0.085 | 不平稳 | |
ln f | 模型2 | −1.669 | 0.096 | 不平稳 | |
ln f | 模型1 | −0.422 | 0.673 | 不平稳 | |
一阶差分序列 | Δ ln s | 模型3 | −16.035 | 0.002 | 平稳 |
Δ ln f | 模型3 | −14.816 | 0.002 | 平稳 |
表2. ADF单位根检验结果表
建立VECM模型并确定最优的滞后阶数,本文根据施瓦茨信息(SC)准则,确定VECM模型的最佳滞后阶数为2阶,故建立VECM (2)模型,表达式如式(4)、式(5)所示。
Δ ln s t = α 1 + δ 1 E C T t − 1 + ∑ i = 1 2 θ 1 , i Δ ln s t − i + ∑ i = 1 m γ 1 , i Δ ln f t − i + e 1 , t (4)
Δ ln f t = α 2 + δ 2 E C T t − 1 + ∑ i = 1 2 θ 2 , i Δ ln s t − i + ∑ i = 1 2 γ 2 , i Δ ln f t − i + e 2 , t (5)
其中式(4)、(5)中, E C T t − 1 为非均衡误差项,并且 E C T t − 1 = φ s t − 1 + τ f t − 1 , θ 、 γ 为模型待估系数, Δ 表示序列一阶差分,其他符号含义同上。
为了进一步了解铜的现货对数价格 ln s 与期货对数价格 ln f 之间是否存在长期均衡关系,就需要进行协整检验,常用的协整检验方法有Engle-Granger检验和Johanson检验。本文选取Engle-Granger协整检验对 ln s 、 ln f 序列进行协整检验。协整检验结果见表3。
参数估计值 | 标准误 | P值 | |
---|---|---|---|
δ 1 | −0.023 | 0.010 | 0.023 |
δ 2 | −0.006 | 0.009 | 0.466 |
表3. 铜的现货对数价格lns与期货对数价格lnf的Engle-Granger协整检验结果
通过表3可以发现,在5%显著性水平下, δ 1 的参数估计值通过了显著性检验,说明铜的现货对数价格 ln s 与期货对数价格 ln f 之间存在着现货对数价格向期货对数价格的单项协整关系。在其他条件不变的情况下,当铜的现货对数价格相对于期货对数价格超过均衡关系1%时,那么下一期铜的现货对数价格将会平均下降0.023%。这进一步说明非均衡误差对铜的现货对数化价格影响更加显著,当短期时间内铜的期货对数化价格与现货对数化价格超过或者低于其均衡关系时,可以通过铜的期货对数化价格来研判现货对数化价格的走势。
VECM模型的参数估计结果如表4所示。
待估参数 | 参数估计值 | 标准误 | T统计量 | P值 |
---|---|---|---|---|
α 1 | −0.001 | 0.001 | −1.724 | 0.085 |
θ 1 , 1 | −0.191 | 0.048 | −3.988 | 0.000 |
γ 1 , 1 | 0.146 | 0.057 | −2.530 | 0.017 |
θ 1 , 2 | 0.021 | 0.048 | 0.418 | 0.676 |
γ 1 , 2 | −0.071 | 0.062 | −1.184 | 0.002 |
α 2 | −0.001 | 0.001 | −0.651 | 0.515 |
θ 2 , 1 | 0.010 | 0.040 | 0.257 | 0.797 |
γ 2 , 1 | 0.014 | 0.049 | 0.292 | 0.771 |
θ 2 , 2 | −0.176 | 0.041 | −4.340 | 0.000 |
γ 2 , 2 | −0.025 | 0.050 | −0.503 | 0.615 |
表4. VECM模型的参数估计系数表
通过表4可以看出,待估参数 α 1 、 α 2 、 θ 1 , 2 、 θ 2 , 1 、 γ 2 , 1 、 γ 2 , 2 对应t统计量的P值均大于显著性水平0.05,故参数估计不显著,而其他待估参数均显著不为0。将表3和表4中参数估计值代入式(4)、式(5)可得VECM模型:
Δ ln s t = − 0.023 E C T t − 1 − 0.191 Δ ln s t − 1 + 0.01 Δ ln s t − 2 + 0.146 Δ ln f t − 1 − 0.071 Δ ln f t − 2 − 0.001 (6)
Δ ln f t = − 0.006 E C T t − 1 + 0.01 Δ ln s t − 1 − 0.176 Δ ln s t − 2 + 0.014 Δ ln f t − 1 − 0.025 Δ ln f t − 1 − 0.001 (7)
通过VECM模型的参数估计结果可以发现,当铜的期现货对数化价格在偏离均衡态时,向量误差模型对铜的现货对数化价格和期货对数化价格具有反向调节作用。铜的期现货对数化价格的1期滞后系数、铜的期货对数化价格的2期滞后系数对当期铜现货价格影响显著,而铜的现货对数化价格的2期滞后系数对当期铜现货价格影响不显著。分析表明,铜的期货对数化价格能够显著影响铜期货对数化价格,而在短期时间内铜现货对数化价格会影响铜期货对数化价格。
利用Engle-Granger可以求出协整向量,其结果如表5所示。
通过表5可知,非均衡误差 E C T t − 1 = ln s t − 1 − 1.011 ln f t − 1 ,即非均衡误差为铜的现货对数化价格和期货对数化价格的基差。
估计方法 | s t − 1 ( φ 的估计值) | f t − 1 ( τ 的估计值) |
---|---|---|
Engle-Granger | 1.000 | −1.011 |
表5. Engle-Granger的协整向量
在此之前,已经通过协整检验得出了铜现货对数化价格与期货对数化价格之间存在着均衡协整关系,并建立了VECM模型,为了更好地检验铜现货对数化价格与期货对数化价格之间是否存在着因果关系,本文采用格兰杰因果检验,检验结果如表6所示。
检验原假设 | F统计量 | P值 |
---|---|---|
铜现货对数化价格不是期货对数化价格的格兰杰原因 | 14.856 | 0.000 |
铜期货对数化价格不是现货对数化价格的格兰杰原因 | 55.369 | 0.000 |
表6. 格兰杰因果关系检验结果表
通过表6可以看出,在5%的显著性水平下,铜现货对数化价格不是期货对数化价格的格兰杰原因的原假设和铜期货对数化价格不是现货对数化价格的格兰杰原因均被拒绝,说明在滞后2阶的情况下,铜现货对数化价格与期货对数化价格互为格兰杰因果,二者之间存在双向引导关系。
为了确保所构建模型的有效性和精度,本文利用ARIMA模型来VECM (2)模型的残差时间序列性质,如果残差序列是平稳的,则说明VECM (2)模型的残差满足时间序列模型的基本假设,无自相关性和异方差性。如果残差序列不平稳,则需要重新调整VECM模型的参数。通过借助R中的auto.arima()函数拟合残差 e 1 , t 、 e 2 , t ,拟合结果见表7。
残差序列 | 拟合模型 | σ 2 | 对数似然比 | AIC |
---|---|---|---|---|
e 1 , t | ARIMA(0,0,0) | 0.000 | 1169.520 | −2337.040 |
e 2 , t | ARIMA(0,0,0) | 0.000 | 1247.420 | −2492.830 |
表7. VECM (2)模型残差检验表
通过表7可以看出,两个残差 e 1 , t 、 e 2 , t 拟合的最佳模型均为ARIMA(0,0,0),所以这两个残差序列均是平稳序列,说明不需要重新调整VECM(2)模型的参数,所构建的模型有效。
为了判断构建的VECM (2)模型是否合适加入时间趋势项,本文加入时间趋势项后进行VECM (2)模型估计,估计结果见表8。
时间趋势项 | 参数估计值 | 标准误差 | T统计量值 | P值 |
---|---|---|---|---|
现货时间趋势项 | − 3.64 × 10 − 6 | 6.61 × 10 − 6 | −0.549 | 0.583 |
期货时间趋势项 | 6.12 × 10 − 6 | 5.53 × 10 − 6 | 1.105 | 0.269 |
表8. 加入时间趋势项的VECM (2)估计结果
通过表8可以发现,在5%的显著性水平下,铜的现货对数化价格和期货对数化价格均没有显著的时间趋势,由此可以判断VECM (2)模型并不适合加入时间趋势项。
前面所构建的VECM (2)模型,没有考虑到在实际中的交易成本对模型的影响,对投资者而言,当铜的现货对数化价格与期货对数化价格之差的绝对值大于交易成本时,投资者才会进行套利。在考虑实际交易成本后,VECM (2)模型对铜的期货对数化价格与现货对数化价格之间关系的可解释性将会减弱。考虑交易成本后,为了使铜的期货对数化价格与现货对数化价格之间关系研究更加科学,本文构建TVECM模型,模型表达式如式(8)、(9)所示。
Δ ln s t = α 1 + δ 1 E C T t − 1 + δ 3 D t − 1 E C T t − 1 + ∑ i = 1 3 θ 1 , i Δ ln s t − i + ∑ i = 1 3 γ 1 , i Δ ln f t − i + e 1 , t (8)
Δ ln f t = α 2 + δ 2 E C T t − 1 + δ 4 D t − 1 E C T t − 1 + ∑ i = 1 3 θ 2 , i Δ ln s t − i + ∑ i = 1 3 γ 2 , i Δ ln f t − i + e 2 , t (9)
在表达式(8)、(9)中, δ 3 、 δ 4 表示投资者存在套利时的修正误差系数,而 δ 1 、 δ 2 表示投资者不存在套利时的误差修正系数, D t − 1 为虚拟变量,表达式如下:
D t − 1 : { 1 , 表 示 存 在 套 利 机 会 0 , 表 示 不 存 在 套 利 机 会 (10)
在TVECM模型中设定两个门限值,门限值估计结果如表9所示。
门限值划分区间 | 现货 | P值 | 期货 | P值 |
---|---|---|---|---|
E C T t − 1 ≥ 0.117 | −0.151 | 0.008 | 0.057 | 0.244 |
0.032 < E C T t − 1 < 0.117 | −0.139 | 0.012 | 0.019 | 0.776 |
E C T t − 1 ≤ 0.032 | −0.108 | 0.599 | −0.026 | 0.138 |
表9. TVECM模型的门限值估计结果
通过表9可以发现,在5%的显著性水平下,当 E C T t − 1 ≥ 0.117 时,铜的期货与现货之间存在着现货向期货的单项协整,铜的现货对数化价格相比于期货对数化价格超出均衡关系1%,那么下一期的现货对数化价格将会平均下降0.151%。这表明在短期时间内当铜的现货对数化价格与期货对数化价格之间的基差大于0.117时,可以通过铜的期货对数化价格判断现货的价格走势。
本文通过构建VECM模型和TVECM模型对上海铜期货价格和现货价格之间的协整关系进行了分析,通过对实证研究结果,得出了以下结论:
第一,在VECM (2)模型中,铜的现货市场与期货市场之间存在着现货向期货的单项协整关系。在其他条件不变的情况下,当铜的现货对数价格相对于期货对数价格超过均衡关系1%时,那么下一期铜的现货对数价格将会平均下降0.023%。进一步说明非均衡误差对铜的现货对数化价格影响更加显著,当短期时间内铜的期货对数化价格与现货对数化价格超过或者低于其均衡关系时,可以通过铜的期货对数化价格来研判现货对数化价格的走势。
第二,在滞后2阶的情况下,铜现货对数化价格与期货对数化价格互为格兰杰因果,二者之间存在双向引导关系。
第三,在TVECM模型中引入投资者的交易成本后,门限估计值分别为0.032、0.117,通过进一步分析后发现在短期时间内当铜的现货对数化价格与期货对数化价格之间的基差大于0.117时,可以通过铜的期货对数化价格判断现货的价格走势。
基于上述结论,本文提出如下对策建议:
第一,加强市场监管。铜期货市场的监管需要得到加强,以防止价格操纵、内幕交易等不法行为的发生。需要建立完善的市场监管机构,加强对铜期货市场的监管,及时发现和打击不良行为。
第二,增强市场透明度。加强铜期货与现货市场之间的信息交流和透明度。双方应加强信息披露,尤其是现货市场价格的公开透明,以便交易者及时调整交易策略,降低市场风险。
第三,推进现货市场化。加快进口关税下调、资源税改革、内外贸资源平衡配置等举措,提高中国现货市场化程度,让市场经济机制更好地调节供求关系,减少现货价格波动。
第四,加强供应链管理。实施合理的供应链管理,以降低铜价格波动的风险。例如,在采矿、冶炼、加工和销售等环节,加强企业间合作和沟通,同时从全球范围内进行生产和销售,以减少价格波动的风险。
第五,加强铜产业政策的制定和调整。制定和实施有针对性的政策,以适应市场需求的变化。政策可以包括相关税费的政策调整、资源配置优化、环保政策的落实等。
湖北省教育厅科学研究计划指导性项目(编号:B2022001)。
唐崇彪,丁咏梅,余佳雄. 上海铜期货价格与现货价格的协整关系研究——基于VECM和TVECM模型Study on the Co-Integration Relationship between Shanghai Copper Futures Price and Spot Price—Based on the VECM and TVECM Models[J]. 运筹与模糊学, 2023, 13(06): 7294-7302. https://doi.org/10.12677/ORF.2023.136716
https://doi.org/10.3390/su9101768