PHASE 02 · VAR STABILITY
VAR 안정성
동반행렬 고유값, 스펙트럼 반경, 수준 대 차분. ρ가 1보다 작아야 Granger의 전제입니다.
ρ(C) = max |λ|
VAR 안정성 검증
동반행렬의 모든 고유값이 단위원 안에 있어야 충격반응이 수렴하고 Granger F가 유효합니다. 수준과 Δp를 같이 봅니다.
ρ Δp
0.969
ρ 수준
1.040
시차
4
T
92
DW
1.80
판정
단위근 근접
연내 일본은행 추가 인상 · 9월 또는 10월 인상 · ρ가 1에 가깝습니다. 검정은 유지하되 충격반응은 오래 갑니다.
VAR 충분VARMA(3,1) ΔAIC=0.18
적합성 방법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.23
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.969 · 수준 1.040
08 · Verdict
판정
폭발이면 충격반응이 발산하고 Granger F는 유효하지 않습니다.
단위근 근접
고유값
|λ| ≥ 1 이면 단위원 밖. 복소쌍은 진동, 실수근은 지수 감쇠·발산.
| # | Re | Im | |λ| | 위치 |
|---|---|---|---|---|
| 1 | 0.969 | 0.000 | 0.969 | 안 |
| 2 | -0.280 | 0.709 | 0.762 | 안 |
| 3 | -0.280 | -0.709 | 0.762 | 안 |
| 4 | -0.713 | 0.000 | 0.713 | 안 |
| 5 | -0.374 | 0.473 | 0.603 | 안 |
| 6 | -0.374 | -0.473 | 0.603 | 안 |
| 7 | 0.022 | 0.329 | 0.330 | 안 |
| 8 | 0.022 | -0.329 | 0.330 | 안 |
쌍 원장
ρ가 큰 쌍이 위로. 행을 누르면 단위원이 바뀝니다.
| 상대 | ρ Δp | ρ 수준 | p | DW | Granger | 판정 |
|---|---|---|---|---|---|---|
| 0.969 | 1.040 | 4 | 1.80 | 유효 | 단위근 근접 | |
| 0.915 | 0.986 | 2 | 1.87 | 유효 | 안정 | |
| 0.788 | 0.993 | 4 | 1.90 | 유효 | 안정 | |
| 0.779 | 0.920 | 2 | 1.57 | 유효 | 안정 | |
| 0.746 | 0.994 | 4 | 1.81 | 유효 | 안정 | |
| 0.738 | 0.996 | 3 | 1.80 | 유효 | 안정 | |
| 0.733 | 0.999 | 2 | 1.48 | 유효 | 안정 | |
| 0.689 | 0.987 | 2 | 1.62 | 유효 | 안정 | |
| 0.640 | 0.994 | 3 | 1.83 | 유효 | 안정 | |
| 0.619 | 0.985 | 2 | 1.80 | 유효 | 안정 | |
| 0.564 | 0.995 | 2 | 1.86 | 유효 | 안정 | |
| 0.544 | 1.000 | 2 | 1.68 | 유효 | 안정 |