Johansenメソッドを使用して共積分ベクトルを取得する


9

私はより良いヨハンセンの方法を理解しようとしているので、 3つのプロセスがある、Likelihood-Based-Inference-Cointegrated-Autoregressive-Econometricsという本によって与えられた例3.1を開発しました。

X1t=i=1tϵ1i+ϵ2t

X2t=αi=1tϵ1i+ϵ3t

X3t=ϵ4t

したがって、共積分ベクトルは[a、-1、0]と[0、0 1]になるはずですが、ヨハンセン法を実行すると、それらを取得できません。

私が試しているコードは次のとおりです:

import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
from statsmodels.tsa.stattools import adfuller
from statsmodels.tsa.johansen import coint_johansen

mu, sigma = 0, 1 # mean and standard deviation
n = 1000

s1 = np.random.normal(mu, sigma, n)
s2 = np.random.normal(mu, sigma, n)
s3 = np.random.normal(mu, sigma, n)

x_1t = np.cumsum(s1)+s2
x_2t = 7*np.cumsum(s1)+s3
x_3t = s3

#Creating Pandas object
todays_date = datetime.datetime.now().date()
index = pd.date_range(todays_date-datetime.timedelta(10), periods=n, freq='D')
y = pd.DataFrame(index=index, data={'col1': x_1t, 'col2': x_2t, 'col3':x_3t} )

p = 4
jres = coint_johansen(y, 0, p)

私はいくつかのp値を試してみましたが、共積分ベクトルを取得できません。何か間違っていることはわかっています。ありがとう。


参照さノートブックはここにある:github.com/mapsa/seminario-doc-2014/blob/master/...
AB3

回答:


6

答えを見つけました。誰かに役立つ場合は、次のノートブックを確認できます。

http://nbviewer.ipython.org/github/mapsa/seminario-doc-2014/blob/master/cointegration-example.ipynb


1
リンクが機能しなくなる可能性があるため、ここにその要旨をコピーしてみませんか?
gung-モニカの回復

coint_johansenメソッドが存在しないようですか?ここで何が欠けていますか?私のnoobの質問に対する私の謝罪を受け入れてください。
RAY

2
このリポジトリのcointブランチにありますgithub.com/josef-pkt/statsmodels
mapsa
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.