Chương 09: Suy luận nhân quả và đánh giá tác động chính sách (Causal Inference & Modern Policy Evaluation)
causal_policy), chúng ta sẽ cùng nhau làm chủ 25 mô hình kinh tế lượng được phân bổ có hệ thống vào 5 nhóm phương pháp chuyên sâu. Trọng tâm cốt lõi của chương là giải quyết câu hỏi: “Làm sao để nhận dạng tác động nhân quả tin cậy bằng Staggered DID thế hệ mới, Đối chứng tổng hợp (SCM/SDID), Hồi quy gián đoạn (rdrobust/RKD), MTE, Bunching và Học máy nhân quả?” thông qua thực hành đối chứng 100% trên bộ ba công cụ Stata 18 MP, R 4.x và Python 3.11+ với dữ liệu thực tế từ EcoData và phòng thí nghiệm EcoLab.1. Đặt vấn đề & Bối cảnh Nghiên cứu của Chương 09
Trong hành trình thực hiện khóa luận tốt nghiệp, luận văn thạc sĩ, luận án tiến sĩ hay các công bố quốc tế chuẩn ISI/Scopus giai đoạn 2021–2026, việc lựa chọn đúng họ mô hình kinh tế lượng (econometric model family, causal_policy) quyết định trực tiếp tính hợp lệ nội tại của toàn bộ nghiên cứu. Rất nhiều bản thảo bị tạp chí từ chối ngay vòng đầu tiên không phải vì thiếu dữ liệu, mà vì nhà nghiên cứu áp dụng ước lượng cổ điển khi cấu trúc dữ liệu đòi hỏi các kỹ thuật chuyên biệt của Suy luận nhân quả và đánh giá tác động chính sách.
Ở Chương 09 này, chúng ta sẽ tiếp cận từng phương pháp từ trực giác kinh tế gần gũi nhất, đi qua cấu trúc đại số và xác suất vừa đủ để hiểu rõ bản chất nhận dạng, rồi trực tiếp thực thi trên cả ba nền tảng phần mềm phổ biến nhất hiện nay. Như biểu diễn tổng quát trong công thức (1), cấu trúc toán học nền tảng của họ mô hình causal_policy được đặt trên hệ tiên đề sau:
$$text{ATT}(g, t) = text{E}[Y_{it}(g) – Y_{it}(infty) mid G_{ig} = 1] tag{1}$$
Từ công thức (1), chúng ta có thể thấy rõ cách thức mà từng nhóm mô hình con trong chương mở rộng hoặc nới lỏng các giả định khắt khe về kỳ vọng có điều kiện, ma trận hiệp phương sai của nhiễu và cơ chế nội sinh trong dữ liệu kinh tế Việt Nam.
2. Bản đồ Cấu trúc 5 Nhóm Phương pháp Con thuộc Chương 09
Để giúp bạn dễ dàng định vị phương pháp phù hợp nhất với đề tài và cấu trúc dữ liệu của mình, bảng dưới đây tổng hợp toàn bộ 5 nhóm con và 25 bài học trong Chương 09:
3. Danh mục Chi tiết 25 Mô hình Kinh tế lượng & Bộ Công cụ Tri-Engine
Mỗi bài học dưới đây đều được thiết kế theo bố cục 8 phần chuẩn sư phạm: từ giới thiệu phương pháp, tình huống dữ liệu thực tế tại Việt Nam, diễn giải toán học, gói lệnh Stata 18 / R 4.x / Python 3.11+, đến thực hành từng bước và kiến giải chính sách:
| Mã Bài | Khóa Mô hình | Tên Bài học Chuyên sâu | Cấp độ | Gói Thư viện Tri-Engine (Stata · R · Python) | Dữ liệu Thực nghiệm Việt Nam (EcoData) |
|---|---|---|---|---|---|
| B136 | propensity_score_matching |
Ghép cặp điểm xu hướng Rosenbaum–Rubin (Propensity Score Matching — PSM) | L1 | base Stata · MatchIt / stats · statsmodels / scikit-learn |
Đánh giá tác động chương trình tín dụng ưu đãi hộ VHLSS |
| B137 | entropy_balancing |
Cân bằng trọng số Entropy chính xác moment Hainmueller (Entropy Balancing ebal) | L2 | Stata ebalance · ebal / stats · scipy / statsmodels |
Tác động đào tạo nghề lên thu nhập lao động Việt Nam |
| B138 | ipw |
Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) | L1 | base Stata · stats · statsmodels |
Đánh giá tác động bảo hiểm y tế hộ gia đình |
| B139 | aipw |
Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) | L2 | base Stata · stats · statsmodels |
Tác động hỗ trợ lãi suất lên đầu tư doanh nghiệp |
| B140 | twfe_did |
Sai biệt trong sai biệt hai chiều hiệu ứng cố định cổ điển (Two-Way Fixed Effects DID — TWFE) | L1 | base Stata · fixest / plm · linearmodels |
Tác động nâng lương tối thiểu vùng tại Việt Nam |
| B141 | event_study_did |
Thiết kế nghiên cứu sự kiện động dẫn–trễ và kiểm định tiền xu hướng (Dynamic Event-Study DID) | L2 | Stata reghdfe / xtreg · fixest · linearmodels / statsmodels |
Động học trước và sau khi tỉnh có cao tốc kết nối |
| B142 | bacon_decomposition |
Định lý phân rã trọng số $2times 2$ và so sánh bị cấm của Goodman-Bacon (bacondecomp) | L3 | bacondecomp · bacondecomp · Exact Goodman-Bacon Theorem 1 closed-form $2times 2$ DID variance-weight & within/between decomposition solver |
Chẩn đoán chệch TWFE khi chính sách triển khai lệch pha |
| B143 | callaway_santanna_did |
Ước lượng DID theo nhóm–thời gian không đồng nhất Callaway–Sant’Anna (csdid / did) | L3 | Stata csdid · did · statsmodels / numpy |
Lộ trình chuyển đổi số lệch pha các tỉnh thành |
| B144 | sun_abraham_did |
Nghiên cứu sự kiện trọng số tương tác đoàn hệ của Sun–Abraham (eventstudyinteract / sunab) | L3 | eventstudyinteract · fixest · pyfixest / Interaction-WeightedFE |
Tác động niêm yết lệch năm lên quản trị công ty |
| B145 | borusyak_imputation_did |
Ước lượng DID suy diễn hiệu quả cao Borusyak–Jaravel–Spiess và Gardner (did_imputation / did2s) | L3 | did_imputation / did2s · didimputation / did2s · pyfixest |
Đánh giá cải cách hành chính cấp tỉnh với phương sai tối thiểu |
| B146 | dechaisemartin_did |
Ước lượng DID động cho xử lý bật–tắt đảo chiều De Chaisemartin–D’Haultfœuille (did_multiplegt_dyn) | L4 | did_multiplegt_dyn · DIDmultiplegtDYN · Cohort-baseline switching DID_M & dynamic DID_l solver với cluster-level influence-function SE |
Tác động các đợt giãn thuế / kiểm soát dịch bật–tắt |
| B147 | wooldridge_etwfe |
Mô hình hiệu ứng cố định hai chiều mở rộng bão hòa Mundlak của Wooldridge (ETWFE / jwdid) | L3 | Stata 18 xthdidregress twfe / jwdid · etwfe + marginaleffects · Saturated Mundlak Cohort $times$ Post-Period FE estimator |
Staggered DID cho biến đếm và tỷ trọng |
| B148 | lp_did |
Phương pháp chiếu cục bộ kết hợp mẫu đối chứng sạch Dube–Girardi–Jordà–Taylor (LP-DID) | L4 | lpdid · fixest clean-control LP loop · pyfixest |
Đánh giá chính sách có phụ thuộc động học trễ |
| B149 | honest_did_sensitivity |
Phân tích độ nhạy vi phạm giả định xu hướng song song Rambachan–Roth (HonestDiD) | L3 | honestdid · HonestDiD · ARP / Fixed-Length Confidence Interval |
Kiểm định độ bền kết luận chính sách |
| B150 | synthetic_control |
Phương pháp đối chứng tổng hợp Abadie–Diamond–Hainmueller (Synthetic Control Method — SCM) | L2 | Stata synth · Synth · scipy / numpy |
Đánh giá cơ chế đặc thù TP.HCM / Hải Phòng |
| B151 | augmented_synthetic_control |
Đối chứng tổng hợp tăng cường hiệu chỉnh chệch Ridge và hoàn thiện ma trận (ASCM / FECT) | L3 | synth + Ridge debiasing / fect · Synth + glmnet / fect · pysyncon |
Đối chứng tổng hợp khi khớp tiền xử lý chưa hoàn hảo |
| B152 | synthetic_did |
Sai biệt trong sai biệt tổng hợp Arkhangelsky–Athey–Hirshberg–Imbens–Wager (SDID) | L3 | sdid · fect / Synth + Frank–Wolfe SDID · pysyncon / Frank–Wolfe regularized unit $hat{omega}$ |
Đánh giá chính sách khu vực kinh tế trọng điểm |
| B153 | sharp_rdd |
Thiết kế hồi quy gián đoạn sắc nhọn tại ngưỡng cắt (Sharp Regression Discontinuity — Sharp RDD) | L1 | Stata rdrobust · rdrobust / stats · statsmodels |
Ngưỡng điểm chuẩn học bổng / ngưỡng quy mô ưu đãi thuế |
| B154 | fuzzy_rdd |
Thiết kế hồi quy gián đoạn mờ sử dụng ngưỡng cắt làm biến công cụ cục bộ (Fuzzy RDD) | L2 | Stata rdrobust · rdrobust / AER · linearmodels / statsmodels |
Tác động chương trình hỗ trợ hộ cận nghèo quanh chuẩn nghèo |
| B155 | rdrobust_sharp_fuzzy |
Hồi quy gián đoạn đa thức cục bộ hiệu chỉnh chệch bền vững Calonico–Cattaneo–Titiunik (rdrobust & rddensity) | L3 | rdrobust + rddensity · rdrobust + rddensity · rdrobust + rddensity |
Chuẩn vàng ước lượng RDD và kiểm định thao túng mật độ |
| B156 | regression_kink_rkd |
Thiết kế gãy khúc hồi quy tại điểm đổi độ dốc công thức chính sách Card–Lee–Pei–Weber (RKD) | L4 | rdrobust, deriv · rdrobust · rdrobust |
Tác động trần trợ cấp bảo hiểm thất nghiệp |
| B157 | bunching_estimator |
Ước lượng tích tụ hành vi tại điểm gãy khúc và bước nhảy thuế suất Saez & Kleven–Waseem (Bunching) | L4 | bunchit / polynomial counterfactual density · bunching · Saez |
Phản ứng của doanh nghiệp tại ngưỡng thuế SME |
| B158 | marginal_treatment_effects |
Hiệu ứng điều trị biên và biến công cụ cục bộ Heckman–Vytlacil (MTE & PRTE) | L4 | mtefe / margte · ivmte / sampleSelection · Semiparametric Polynomial & Normal Local-IV MTE |
Lợi suất giáo dục đại học theo kháng cự lựa chọn ẩn |
| B159 | double_ml |
Học máy khử chệch trực giao Neyman hai bước Chernozhukov (Double / Debiased Machine Learning — DML) | L3 | Stata 18 lasso · DoubleML / glmnet · scikit-learn / statsmodels |
Nhận dạng nhân quả phi tuyến với biến kiểm soát chiều cao |
| B160 | causal_forest |
Rừng ngẫu nhiên nhân quả ước lượng tác động điều trị không đồng nhất cá thể Wager–Athey (Causal Forest) | L4 | Stata Python/R bridge · grf / randomForest · scikit-learn / econml |
Cá nhân hóa chính sách hỗ trợ doanh nghiệp và hộ gia đình |
4. Cây Quyết định Lựa chọn Mô hình & Bộ Kiểm định Chẩn đoán Bắt buộc
Khi đứng trước bộ dữ liệu thực tế của đề tài, chúng ta cần tuân thủ quy trình kiểm định tiền ước lượng (pre-estimation) và hậu ước lượng (post-estimation) nghiêm ngặt. Hãy kiểm tra các tiêu chí sau trước khi báo cáo kết quả của Chương 09:
- Câu hỏi nhận dạng trung tâm: Làm sao để nhận dạng tác động nhân quả tin cậy bằng Staggered DID thế hệ mới, Đối chứng tổng hợp (SCM/SDID), Hồi quy gián đoạn (rdrobust/RKD), MTE, Bunching và Học máy nhân quả?
- Hệ thống kiểm định bắt buộc: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check
- Nguyên tắc đối chứng Tri-Engine: Luôn kiểm tra độ hội tụ của hàm hợp lý / điều kiện moment và đối chiếu sai số chuẩn vững (robust / clustered standard errors) giữa Stata, R và Python.
- Tích hợp dữ liệu thực nghiệm Việt Nam: Tận dụng trực tiếp các bộ dữ liệu chuẩn hóa từ EcoData (VHLSS, GSO 63 tỉnh thành, Hải quan, Doanh nghiệp niêm yết HOSE/HNX) và chạy kiểm chứng một chạm trên EcoLab.
5. Bộ Prompt Copy-Ready Sinh Infographic Tổng quan Chương 09
Chúng ta có thể sao chép trực tiếp đoạn prompt chuẩn hóa dưới đây vào Gemini NotebookLM hoặc Gemini Nano Banana Pro để tạo ngay bản đồ Infographic 4K tóm tắt toàn bộ cây quyết định của Chương 09 phục vụ slide bảo vệ khóa luận, luận văn hoặc bài giảng trên lớp:
🎨 Copy-Ready Prompt: Gemini NotebookLM / Nano Banana Pro — Infographic Chương 09
[GEMINI NOTEBOOKLM / NANO BANANA PRO — INFOGRAPHIC BLUEPRINT CHƯƠNG 09]
Vai trò: Chuyên gia Thiết kế Đồ họa Thông tin Học thuật (Academic Information Designer) & Giảng viên Kinh tế lượng 18 năm kinh nghiệm.
Chủ đề: Bản đồ Tư duy & Cây Quyết định Lựa chọn Mô hình cho Chương 09: Suy luận nhân quả và đánh giá tác động chính sách (Causal Inference & Modern Policy Evaluation — Họ `causal_policy`).
Bố cục 4 tầng (Tỷ lệ 16:9 chuẩn 4K 3840x2160, nền trắng ngà #F8FAFC, bảng màu xanh Navy #0F172A, xanh cổ vịt #0D9488, vàng đồng #D97706):
1. Tầng Tiêu đề (Top Banner): 'CHƯƠNG 09: SUY LUẬN NHÂN QUẢ VÀ ĐÁNH GIÁ TÁC ĐỘNG CHÍNH SÁCH' — Tổng hợp 5 nhóm con & 25 mô hình Tri-Engine (Stata 18 · R 4.x · Python 3.11+).
2. Tầng Cây Quyết định (Central Decision Flowchart): Sơ đồ rẽ nhánh giúp nhà nghiên cứu chọn đúng mô hình trong 5 nhóm con: (1) Nhóm 9.1 Ghép cặp & Cân bằng mẫu (2 bài: B136, B137); (2) Nhóm 9.4 Tác động điều trị & Nhân quả ML (6 bài: B138, B139, B150, B151, B159, B160); (3) Nhóm 9.2 DID cổ điển & Nghiên cứu sự kiện (3 bài: B140, B141, B143); (4) Nhóm 9.5 DID cổ điển & Nghiên cứu sự kiện (8 bài: B142, B144, B145, B146, B147, B148, B149, B152); (5) Nhóm 9.3 Hồi quy gián đoạn, Gãy khúc, Bunching & MTE (6 bài: B153, B154, B155, B156, B157, B158).
3. Tầng Công thức & Kiểm định Trụ cột (Core Identification & Diagnostics): Câu hỏi trọng tâm (Làm sao để nhận dạng tác động nhân quả tin cậy bằng Staggered DID thế hệ mới, Đối chứng tổng hợp (SCM/SDID), Hồi quy gián đoạn (rdrobust/RKD), MTE, Bunching và Học máy nhân quả?), công thức lõi (text{ATT}(g, t) = text{E}[Y_{it}(g) - Y_{it}(infty) mid G_{ig} = 1]), kiểm định bắt buộc (Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check) và các mô hình tiêu biểu (B136 Ghép cặp điểm xu hướng Rosenbaum–Rubin (Propensity Score Matching — PSM) (propensity_score_matching), B137 Cân bằng trọng số Entropy chính xác moment Hainmueller (Entropy Balancing ebal) (entropy_balancing), B138 Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) (ipw), B139 Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) (aipw), B140 Sai biệt trong sai biệt hai chiều hiệu ứng cố định cổ điển (Two-Way Fixed Effects DID — TWFE) (twfe_did), B141 Thiết kế nghiên cứu sự kiện động dẫn–trễ và kiểm định tiền xu hướng (Dynamic Event-Study DID) (event_study_did)).
4. Tầng Chân trang (Footer): Hệ sinh thái Dữ liệu Kinh tế Việt Nam EcoData (ecodata.io.vn) · Phòng thí nghiệm Kinh tế lượng EcoLab (ecolab.io.vn) · Cổng Học thuật VietLOD (vietlod.tnsai.vn).6. Tài liệu Tham khảo Học thuật Trụ cột của Chương 09
- Tài liệu gốc & Giáo trình nền tảng: Callaway & Sant’Anna (2021); Sun & Abraham (2021); Borusyak, Jaravel & Spiess (2024); Arkhangelsky et al. (2021); Calonico, Cattaneo & Titiunik (2014); Chernozhukov et al. (2018)
- Wooldridge, J. M. (2010). Econometric Analysis of Cross Section and Panel Data (2nd ed.). MIT Press.
- Cameron, A. C., & Trivedi, P. K. (2022). Microeconometrics Using Stata (2nd ed., Vols. 1–2). Stata Press.
- Hansen, B. E. (2022). Econometrics. Princeton University Press.
← Chương trước: Chương 08: Độ biến động tài chính, định giá tài sản và rủi ro hệ thống
→ Bắt đầu Bài học đầu tiên của Chương 09: Bài B136 — Ghép cặp điểm xu hướng Rosenbaum–Rubin (Propensity Score Matching — PSM) (
propensity_score_matching)→ Chương tiếp theo: Chương 10: Kinh tế lượng không gian, Conley HAC và mạng cặp đôi
Hướng dẫn chọn mô hình và so sánh 2 phương pháp ước lượng thuộc Nhóm 9.1 (Ghép cặp & Cân bằng mẫu) trên
Hướng dẫn chọn mô hình và so sánh 6 phương pháp ước lượng thuộc Nhóm 9.4 (Tác động điều trị & Nhân quả
Hướng dẫn chọn mô hình và so sánh 3 phương pháp ước lượng thuộc Nhóm 9.2 (DID cổ điển & Nghiên cứu sự
Hướng dẫn chọn mô hình và so sánh 8 phương pháp ước lượng thuộc Nhóm 9.5 (DID cổ điển & Nghiên cứu sự
Hướng dẫn chọn mô hình và so sánh 6 phương pháp ước lượng thuộc Nhóm 9.3 (Hồi quy gián đoạn, Gãy khúc, Bunching