coint

语法

coint(Y0, Y1, [trend="c"], [method="aeg"], [maxLag], [autoLag="aic"])

详情

检验单变量方程中是否存在协整性。结果以字典的形式返回,包括:

  • tStat:浮点数标量,表示残差的单位根检验的 t 统计量
  • pValue:浮点数标量,表示 MacKinnon 近似 p 值
  • criticalValues:字典,表示基于回归曲线,在1%、5%和10%水平上的检验统计量的临界值

参数

Y0 数值向量,元素不可为空值,表示协整系统中的第一个元素。

Y1 数值向量或数值矩阵,元素不可为空值,且元素个数与 Y0 相同。表示协整系统中的其余元素。

trend 字符串,指定协整方程的回归中使用的趋势项。取值有如下选择:

  • "c":默认值,表示只使用常数。
  • "ct":使用常数和趋势。
  • "ctt":使用常数、线性趋势和二次趋势。
  • "n":不使用常数和趋势。

method 字符串,指定协整检验所使用的方法,支持取值为 “aeg”,代表增广 Engle-Granger 两步协整性检验。

maxLag 非负整数,作为调用 adfuller 时传递的参数,指定 adfuller 检验中使用的最大滞后期。

autoLag 字符串,作为调用 adfuller 时传递的参数,指定在 0~maxLag 中自动确定滞后期长度时使用的方法。取值有如下选择:

  • "aic":默认值,表示使用 Akaike Information Criterion 来确定滞后期数值。
  • "bic":表示使用 Bayesian information criterion 来确定滞后期数值。
  • "tstat":将滞后期初始值设为 maxLag,然后逐步减 1,直到上一个滞后期数值的 t 统计量在 5% 显著性水平上显著。
  • "max":将滞后期数值设置为 maxLag。

例子

Y0 = 234 267 289 301 312 323 334 345 356;
Y1 = 267 289 301 312 323 334 345 356 367;
coint(Y0, Y1);

输出为字典:

tValue->-1.498236972489574
pValue->0.761867238199341
criticalValues->[-5.789286875000001,-4.206501875,-3.6171]