PHASE 02 · VAR STABILITY
VAR 안정성
동반행렬 고유값, 스펙트럼 반경, 수준 대 차분. ρ가 1보다 작아야 Granger의 전제입니다.
ρ(C) = max |λ|
VAR 안정성 검증
동반행렬의 모든 고유값이 단위원 안에 있어야 충격반응이 수렴하고 Granger F가 유효합니다. 수준과 Δp를 같이 봅니다.
ρ Δp
0.966
ρ 수준
1.022
시차
4
T
92
DW
1.82
판정
단위근 근접
러·우 휴전 합의 시점 · 9월 또는 10월 인상 · ρ가 1에 가깝습니다. 검정은 유지하되 충격반응은 오래 갑니다.
VAR 충분VARMA(3,1) ΔAIC=0.23
적합성 방법PIPELINE
8단계 검증
01 · Align
정렬
두 선두 확률을 같은 날짜로 맞춥니다. VAR은 균형 패널이 필요합니다.
공통 날짜 97
02 · Difference
차분
수준 VAR은 거의 단위근입니다. Δp로 정상 근방에 올린 뒤 안정성을 봅니다.
Δp 96
03 · Lag
시차
시스템 AIC = log|Σ| + 2k²p/T. 방정식마다 다른 p를 쓰지 않습니다.
AIC p=4
04 · Fit
추정
방정식별 OLS. 가우스 가정 아래 VAR은 SUR과 같습니다.
T=92 · log|Σ|=-22.49
05 · Companion
동반행렬
C는 kp×kp. 첫 블록 행이 A₁…Aₚ, 아래는 이동항등입니다.
C는 8×8
06 · Eigenvalues
고유값
det(I − A₁z − ⋯ − Aₚzᵖ)=0의 역수가 C의 λ. QR로 구합니다.
고유값 8개
07 · Radius
반경
ρ(C)=max|λ|. Lütkepohl: 모든 |λ|<1 일 때만 안정.
ρ=0.966 · 수준 1.022
08 · Verdict
판정
폭발이면 충격반응이 발산하고 Granger F는 유효하지 않습니다.
단위근 근접
고유값
|λ| ≥ 1 이면 단위원 밖. 복소쌍은 진동, 실수근은 지수 감쇠·발산.
| # | Re | Im | |λ| | 위치 |
|---|---|---|---|---|
| 1 | 0.966 | 0.000 | 0.966 | 안 |
| 2 | 0.226 | 0.748 | 0.781 | 안 |
| 3 | 0.226 | -0.748 | 0.781 | 안 |
| 4 | -0.265 | 0.720 | 0.767 | 안 |
| 5 | -0.265 | -0.720 | 0.767 | 안 |
| 6 | -0.690 | 0.000 | 0.690 | 안 |
| 7 | -0.593 | 0.309 | 0.669 | 안 |
| 8 | -0.593 | -0.309 | 0.669 | 안 |
쌍 원장
ρ가 큰 쌍이 위로. 행을 누르면 단위원이 바뀝니다.
| 상대 | ρ Δp | ρ 수준 | p | DW | Granger | 판정 |
|---|---|---|---|---|---|---|
| 0.966 | 1.022 | 4 | 1.82 | 유효 | 단위근 근접 | |
| 0.953 | 1.008 | 3 | 1.57 | 유효 | 단위근 근접 | |
| 0.926 | 0.994 | 3 | 1.87 | 유효 | 안정 | |
| 0.831 | 0.991 | 3 | 1.83 | 유효 | 안정 | |
| 0.826 | 1.005 | 4 | 1.66 | 유효 | 안정 | |
| 0.779 | 0.996 | 4 | 1.82 | 유효 | 안정 | |
| 0.768 | 0.997 | 4 | 1.92 | 유효 | 안정 | |
| 0.752 | 1.001 | 3 | 1.76 | 유효 | 안정 | |
| 0.736 | 0.986 | 3 | 1.85 | 유효 | 안정 | |
| 0.707 | 0.997 | 3 | 1.77 | 유효 | 안정 | |
| 0.687 | 0.931 | 1 | 1.76 | 유효 | 안정 | |
| 0.619 | 0.985 | 2 | 1.80 | 유효 | 안정 |