本节书摘来异步社区《量化金融R语言初级教程》一书中的第1章，第1.3节，作者： 【匈牙利】Gergely Daróczi（盖尔盖伊） , 等 译者： 高蓉 , 李茂 责编： 胡俊英，更多章节内容可以访问云栖社区“异步社区”公众号查看。

协整的思想缘于Granger（1981年）提出的一个概念，后来由Engle和Granger（1987年）加以形式化。协整的思想是指寻找非平稳时间序列之间的一个线性组合，这个线性组合是一个平稳时间序列。因此，协整方法可以用于检测非平稳时间序列（比如价格）之间的稳定长期关系。

航空燃油的交叉对冲航空公司很自然需要购买航空燃油。但由于航空燃油价格的波动很剧烈，大部分航空公司会将它们对航空燃油价格变化的风险敞口对冲掉一部分。如果市场中缺乏航空燃油OTC产品（译注：OTC产品指交易所场外柜台交易产品），航空公司会使用相关交易所交易的期货合约（比如，取暖油）来实现对冲。在下面的部分中，我们首先使用经典方法导出最优对冲比率，这种方法仅仅考虑两种价格之间短期波动。然后考察价格之间的长期稳定联系，进而改进最优对冲比。

首先，我们载入需要使用的包。urca包有一些有用的方法，可以用于单位根检验和估计协整关系。

我们导入航空燃油和取暖油的月价格。

我们使用Engle和Granger的两步估计方法。首先，使用增强的Dickey-Fuller检验（augmented Dickey-Fuller test，ADF检验）对两个序列进行单位根检验（非平稳性）。

结果显示，因为检验统计量值−1.1335大于临界值−3.46，所以在1%的置信水平上不能拒绝非平稳（航空燃油时间序列包含一个单位根）的原假设。同样的结果对取暖油也成立（检验统计量是−1.041）。

现在我们可以继续估计静态均衡模型，并使用ADF方法检验时间序列的残差是否平稳。请注意，目前的研究序列是上一步的估计结果，因此，我们现在必须使用不同的临界值［参见Engle和Yoo1987发表的论文］。

得到的检验统计量是−8.912，而规模为200的样本在1%的置信水平上的临界值为−4.00。所以，我们拒绝非平稳的原假设。因此我们发现了两个协整变量并且可以进行第二步，两个协整变量意味着一个误差修正模型（ECM）。ECM是一个动态模型，刻画了系统如何（以及多快）返回之前估计的静态均衡，这个静态均衡存储在变量mod_static中。