Nhóm 6.1: Nghiệm đơn vị & Đồng tích hợp (unit_root_cointegration)
Tóm tắt Nhóm phương pháp 6.1 (AEO Quick Answer): Nhóm 6.1 — Nghiệm đơn vị & Đồng tích hợp (mã nhóm
unit_root_cointegration) là phân nhánh chuyên sâu thuộc Chương 06: Chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng, quy tụ 6 phương pháp ước lượng (từ B084 đến B089). Phần tổng quan này giúp chúng ta so sánh trực diện điều kiện áp dụng, cấu trúc toán học và cú pháp thực thi trên Stata 18 MP, R 4.x và Python 3.11+ để chọn đúng mô hình cho bài toán nghiên cứu.1. Bản chất Học thuật & Vị trí của Nhóm 6.1 (Nghiệm đơn vị & Đồng tích hợp)
Khi triển khai nghiên cứu định lượng trong miền Chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng, chúng ta thường đứng trước câu hỏi: giữa 6 mô hình cùng thuộc nhóm phương pháp con (econometric sub-group, unit_root_cointegration), đâu là công cụ khớp nhất với quá trình sinh dữ liệu (DGP) của đề tài?
Mọi phương pháp trong Nhóm 6.1 đều xuất phát từ khung biểu diễn đại số cơ sở của chương như trình bày tại công thức (1):
$$Delta Y_t = alpha_0 + rho Y_{t-1} + theta^+ X_{t-1}^+ + theta^- X_{t-1}^- + sum_{j=1}^{p-1} gamma_j Delta Y_{t-j} + varepsilon_t tag{1}$$
Sự khác biệt giữa các mô hình từ B084 (adf) đến B089 (johansen) nằm ở cách chúng ta đặt giả định lên phân phối của sai số ngẫu nhiên, cấu trúc nội sinh và phương pháp tối ưu hóa hàm mục tiêu từ công thức (1).
2. Bảng So sánh Trực diện 6 Mô hình thuộc Nhóm 6.1
Bảng đối chiếu dưới đây giúp bạn nhìn tổng thể toàn bộ 6 bài học trong Nhóm 6.1 kèm câu lệnh cốt lõi trên cả ba phần mềm:
| Mã Bài | model_key |
Tên Mô hình (VI / EN) | Cấp độ | Lệnh Thực thi Chuẩn (Stata 18 · R 4.x · Python 3.11+) | Ý nghĩa Kinh tế lượng Trọng tâm |
|---|---|---|---|---|---|
| B084 | adf |
Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) ADF Test |
L1 | Stata: dfuller y, lags(2) trendR: tseries::adf.test(y) / urca::ur.df(y, type = 'trend', lags = 2)Python: adfuller(y, maxlag=4, regression='ct', autolag='AIC') |
Phương pháp Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) (ADF Test) trong họ mô hình chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng, giúp ước lượng chính xác tham số cấu trúc và kiểm định giả thuyết kinh tế lượng. |
| B085 | phillips_perron |
Kiểm định nghiệm đơn vị phi tham số Phillips–Perron (PP Unit-Root Test) Phillips–Perron Nonparametric Unit-Root Test |
L1 | Stata: pperron (pperron y, trend lags(m))R: urca (urca::ur.pp(y, type = "Z-tau", model = "trend", lags = "short"))Python: arch (arch.unitroot.PhillipsPerron(y, trend="c", test_type="tau")) |
Kiểm định nghiệm đơn vị hiệu chỉnh phi tham số Newey–West trực tiếp trên thống kê $Z_alpha, Z_t$ mà không cần thêm nhiều biến trễ $Delta Y_{t-j}$, bền vững trước cả tự tương quan chuỗi và phương sai thay đổi tổng quát. |
| B086 | df_gls |
Kiểm định nghiệm đơn vị khử xu thế cục bộ lực kiểm định cao Elliott–Rothenberg–Stock (DF-GLS) Elliott–Rothenberg–Stock DF-GLS Unit-Root Test |
L2 | Stata: dfgls (dfgls y, maxlag(k))R: urca (urca::ur.ers(y, type = "DF-GLS", model = "trend", lag.max = k))Python: arch (arch.unitroot.DFGLS(y, trend="ct", method="maic")) |
Khử xu thế cục bộ bằng GLS gần biên nghiệm đơn vị ($bar{c} = -7$ hoặc $-13.5$) trước khi chạy hồi quy Dickey–Fuller, mang lại lực kiểm định (statistical power) cao vượt trội so với ADF và PP trên mẫu hữu hạn. |
| B087 | kpss |
Kiểm định tính dừng Kwiatkowski–Phillips–Schmidt–Shin đối chứng với ADF (KPSS Test) KPSS Test |
L1 | Stata: kpss y, maxlag(4)R: tseries::kpss.test(y, null = 'Trend')Python: kpss(y, regression='ct', nlags='auto') |
Phương pháp Kiểm định tính dừng Kwiatkowski–Phillips–Schmidt–Shin đối chứng với ADF (KPSS Test) (KPSS Test) trong họ mô hình chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng, giúp ước lượng chính xác tham số cấu trúc và kiểm định giả thuyết kinh tế lượng. |
| B088 | zivot_andrews |
Kiểm định nghiệm đơn vị có điểm gãy cấu trúc nội sinh Zivot–Andrews Zivot–Andrews Endogenous Structural-Break Unit-Root Test |
L2 | Stata: zandrews (ssc install zandrews: zandrews y, break(both) lagmethod(AIC))R: urca (urca::ur.za(y, model = "both", lag = k))Python: arch (arch.unitroot.ZivotAndrews(y, trend="ct", method="aic")) |
Kiểm định nghiệm đơn vị cho phép 1 điểm gãy cấu trúc nội sinh ở mức chặn (Model A), độ dốc xu thế (Model B) hoặc cả hai (Model C), khắc phục lỗi chấp nhận sai nghiệm đơn vị của ADF khi chuỗi có cú sốc gãy bậc. |
| B089 | johansen |
Kiểm định đồng tích hợp đa biến theo vết và trị riêng cực đại Johansen Johansen Test |
L2 | Stata: vecrank y x1 x2, trend(constant) lags(2)R: urca::ca.jo(df[, c('y', 'x1', 'x2')], type = 'trace', K = 2)Python: coint_johansen(df[['y', 'x1', 'x2']], det_order=0, k_ar_diff=1) |
Phương pháp Kiểm định đồng tích hợp đa biến theo vết và trị riêng cực đại Johansen (Johansen Test) trong họ mô hình chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng, giúp ước lượng chính xác tham số cấu trúc và kiểm định giả thuyết kinh tế lượng. |
3. Hướng dẫn Chọn nhanh Mô hình Phù hợp trong Nhóm 6.1
Dựa trên kinh nghiệm phản biện và hướng dẫn luận văn, luận án thực tế, chúng ta có thể áp dụng quy tắc chọn nhanh sau đây:
- Chọn B084 — Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) (
adf) khi: Áp dụng khi nghiên cứu Kiểm định tính dừng chuỗi CPI, GDP, VN-Index Việt Nam và các bài toán thuộc nhóm nghiệm đơn vị & đồng tích hợp. (Kiểm định kèm theo: ADF/PP/KPSS/DF-GLS/Zivot–Andrews Unit Root, Ljung–Box Q, Pesaran PSS Bounds F & t Test, Bai–Perron Sup-F, PSY GSADF Wild Bootstrap). - Chọn B085 — Kiểm định nghiệm đơn vị phi tham số Phillips–Perron (PP Unit-Root Test) (
phillips_perron) khi: Đối chứng bắt buộc cặp đôi với kiểm định ADF trong bảng kiểm định tính dừng chuẩn xuất bản cho chuỗi thời gian kinh tế và tài chính. (Kiểm định kèm theo: Thống kê $Z(tau)$ và $Z(alpha)$, MacKinnon approximate $p$-value, băng thông Bartlett Newey–West, giá trị tới hạn $1%, 5%, 10%$.). - Chọn B086 — Kiểm định nghiệm đơn vị khử xu thế cục bộ lực kiểm định cao Elliott–Rothenberg–Stock (DF-GLS) (
df_gls) khi: Chuỗi vĩ mô có nghiệm gần 1 (lãi suất thực, tỷ giá thực, tỷ lệ thất nghiệp) với cỡ mẫu trung bình ($T in [30, 150]$) nơi ADF/PP có lực kiểm định rất thấp. (Kiểm định kèm theo: Thống kê $tau_{text{DF-GLS}}$, chọn độ trễ Ng–Perron Modified AIC (MAIC — chống méo kích thước kiểm định), bảng giá trị tới hạn ERS.). - Chọn B087 — Kiểm định tính dừng Kwiatkowski–Phillips–Schmidt–Shin đối chứng với ADF (KPSS Test) (
kpss) khi: Áp dụng khi nghiên cứu Đối chứng cặp đôi ADF–KPSS cho chuỗi xuất khẩu và các bài toán thuộc nhóm nghiệm đơn vị & đồng tích hợp. (Kiểm định kèm theo: ADF/PP/KPSS/DF-GLS/Zivot–Andrews Unit Root, Ljung–Box Q, Pesaran PSS Bounds F & t Test, Bai–Perron Sup-F, PSY GSADF Wild Bootstrap). - Chọn B088 — Kiểm định nghiệm đơn vị có điểm gãy cấu trúc nội sinh Zivot–Andrews (
zivot_andrews) khi: Mọi chuỗi vĩ mô (GDP, CPI, tỷ giá, lãi suất, thương mại) có bước nhảy chế độ rõ rệt (như đổi mới, khủng hoảng 2008, COVID-19 2020). (Kiểm định kèm theo: Thống kê $min t_{text{ZA}}$, thời điểm gãy nội sinh $hat{T}_B$, giá trị tới hạn Zivot–Andrews ($-5.57, -5.08, -4.82$ cho Model C), độ trễ tối ưu AIC/BIC.). - Chọn B089 — Kiểm định đồng tích hợp đa biến theo vết và trị riêng cực đại Johansen (
johansen) khi: Áp dụng khi nghiên cứu Cầu tiền, lãi suất và thu nhập tại Việt Nam và các bài toán thuộc nhóm nghiệm đơn vị & đồng tích hợp. (Kiểm định kèm theo: ADF/PP/KPSS/DF-GLS/Zivot–Andrews Unit Root, Ljung–Box Q, Pesaran PSS Bounds F & t Test, Bai–Perron Sup-F, PSY GSADF Wild Bootstrap).
Điều hướng Học tập trong Nhóm 6.1:
↑ Quay lại Tổng quan Chương 06: Chương 06 — Chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng
→ Đi tới Bài học đầu tiên của Nhóm 6.1: Bài B084 — Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) (
→ Thực hành trực tiếp trên EcoLab: Mở Phòng thí nghiệm Kinh tế lượng EcoLab (Tri-Engine)
↑ Quay lại Tổng quan Chương 06: Chương 06 — Chuỗi thời gian đơn biến, phi tuyến, điểm gãy cấu trúc và bong bóng
→ Đi tới Bài học đầu tiên của Nhóm 6.1: Bài B084 — Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) (
adf)→ Thực hành trực tiếp trên EcoLab: Mở Phòng thí nghiệm Kinh tế lượng EcoLab (Tri-Engine)
Bài liên quan
- Bài B084: Kiểm định nghiệm đơn vị Dickey–Fuller tăng cường (Augmented Dickey–Fuller — ADF Test) (ADF Test — adf)Cấp độ 1: Nền tảng (Foundation)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B085: Kiểm định nghiệm đơn vị phi tham số Phillips–Perron (PP Unit-Root Test) (Phillips–Perron Nonparametric Unit-Root Test — phillips_perron)Cấp độ 1: Nền tảng (Foundation)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B086: Kiểm định nghiệm đơn vị khử xu thế cục bộ lực kiểm định cao Elliott–Rothenberg–Stock (DF-GLS) (Elliott–Rothenberg–Stock DF-GLS Unit-Root Test — df_gls)Cấp độ 2: Trung cấp (Intermediate)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B087: Kiểm định tính dừng Kwiatkowski–Phillips–Schmidt–Shin đối chứng với ADF (KPSS Test) (KPSS Test — kpss)Cấp độ 1: Nền tảng (Foundation)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B088: Kiểm định nghiệm đơn vị có điểm gãy cấu trúc nội sinh Zivot–Andrews (Zivot–Andrews Endogenous Structural-Break Unit-Root Test — zivot_andrews)Cấp độ 2: Trung cấp (Intermediate)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B089: Kiểm định đồng tích hợp đa biến theo vết và trị riêng cực đại Johansen (Johansen Test — johansen)Cấp độ 2: Trung cấp (Intermediate)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)