Nhóm 9.4: Tác động điều trị & Nhân quả ML (treatment_effects)
treatment_effects) là phân nhánh chuyên sâu thuộc Chương 09: Suy luận nhân quả và đánh giá tác động chính sách, quy tụ 6 phương pháp ước lượng (từ B138 đến B160). 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 9.4 (Tác động điều trị & Nhân quả ML)
Khi triển khai nghiên cứu định lượng trong miền Suy luận nhân quả và đánh giá tác động chính sách, 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, treatment_effects), đâ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 9.4 đề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):
$$text{ATT}(g, t) = text{E}[Y_{it}(g) – Y_{it}(infty) mid G_{ig} = 1] tag{1}$$
Sự khác biệt giữa các mô hình từ B138 (ipw) đến B160 (causal_forest) 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 9.4
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 9.4 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 |
|---|---|---|---|---|---|
| B138 | ipw |
Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) IPW |
L1 | Stata: teffects ipw (y) (treat x1 x2), ateR: survey / WeightIt IPW Horvitz–Thompson EstimatorPython: Horvitz–Thompson & Hajek Normalized IPW Treatment Estimator |
Phương pháp Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) (IPW) trong họ mô hình suy luận nhân quả và đánh giá tác động chính sách, 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. |
| B139 | aipw |
Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) AIPW |
L2 | Stata: teffects aipw (y x1 x2) (treat x1 x2), ateR: AIPW / Doubly-Robust Robins–Rotnitzky–Zhao EstimatorPython: Doubly-Robust Augmented IPW (AIPW) Estimator |
Phương pháp Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) (AIPW) trong họ mô hình suy luận nhân quả và đánh giá tác động chính sách, 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. |
| B150 | synthetic_control |
Phương pháp đối chứng tổng hợp Abadie–Diamond–Hainmueller (Synthetic Control Method — SCM) Synthetic Control |
L2 | Stata: synth y x1 x2, trunit(1) trperiod(2018) nestedR: Synth::synth(dataprep.res)Python: Abadie–Diamond–Hainmueller Constrained Simplex Donor Weight Solver |
Phương pháp Phương pháp đối chứng tổng hợp Abadie–Diamond–Hainmueller (Synthetic Control Method — SCM) (Synthetic Control) trong họ mô hình suy luận nhân quả và đánh giá tác động chính sách, 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. |
| 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) Ridge-Augmented & Matrix-Completion Synthetic Control (ASCM / FECT) |
L3 | Stata: synth + Ridge debiasing / fect (ssc install synth; Ridge outcome bias-correction + conformal inference)R: Synth + glmnet / fect (Synth::synth + glmnet Ridge bias-correction & fect::fect(..., method = "mc") trên CRAN active; tương thích augsynth)Python: pysyncon (pysyncon.AugSynth + Chernozhukov–Wüthrich–Zhu conformal $p$-value & in-space placebo) |
Khắc phục chệch khớp không hoàn hảo trước can thiệp (imperfect pre-treatment fit) của Synthetic Control cổ điển khi đơn vị xử lý nằm gần hoặc ngoài bao lồi (convex hull) của nhóm đối chứng thông qua hiệu chỉnh chệch mô hình kết quả Ridge hoặc hoàn thiện ma trận nhân tố ẩn. |
| B159 | double_ml |
Học máy khử chệch trực giao Neyman hai bước Chernozhukov (Double / Debiased Machine Learning — DML) Double ML |
L3 | Stata: xporegress y d, controls(x1-x10) xfolds(5)R: DoubleML::DoubleMLPLR$new(data, ml_l, ml_m)$fit()Python: DoubleML Cross-Fitted Neyman-Orthogonal Partialling-Out PLR |
Phương pháp Học máy khử chệch trực giao Neyman hai bước Chernozhukov (Double / Debiased Machine Learning — DML) (Double ML) trong họ mô hình suy luận nhân quả và đánh giá tác động chính sách, 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. |
| 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) Causal Forest |
L4 | Stata: Causal Forest CATE via Python/R GRF integrationR: grf::causal_forest(X, Y, W) -> grf::average_treatment_effect(cf)Python: Orthogonal R-Learner / Causal Forest CATE(x) Estimator |
Phương pháp 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) (Causal Forest) trong họ mô hình suy luận nhân quả và đánh giá tác động chính sách, 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 9.4
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 B138 — Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) (
ipw) khi: Áp dụng khi nghiên cứu Đánh giá tác động bảo hiểm y tế hộ gia đình và các bài toán thuộc nhóm tác động điều trị & nhân quả ml. (Kiểm định kèm theo: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check). - Chọn B139 — Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) (
aipw) khi: Áp dụng khi nghiên cứu Tác động hỗ trợ lãi suất lên đầu tư doanh nghiệp và các bài toán thuộc nhóm tác động điều trị & nhân quả ml. (Kiểm định kèm theo: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check). - Chọn B150 — Phương pháp đối chứng tổng hợp Abadie–Diamond–Hainmueller (Synthetic Control Method — SCM) (
synthetic_control) khi: Áp dụng khi nghiên cứu Đánh giá cơ chế đặc thù TP.HCM / Hải Phòng và các bài toán thuộc nhóm tác động điều trị & nhân quả ml. (Kiểm định kèm theo: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check). - Chọn B151 — Đố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) (
augmented_synthetic_control) khi: Nghiên cứu tình huống so sánh (1 hoặc một vài tỉnh/quốc gia thực thi cơ chế đặc thù như TP.HCM, Đà Nẵng, Hải Phòng) với $T_0 ge 8$ kỳ trước chính sách. (Kiểm định kèm theo: Độ giảm chệch L2 imbalance trước/sau Ridge augmentation, Conformal inference $p$-value, kiểm định hoán vị không gian (in-space placebo ratio) & leave-one-donor-out.). - Chọn B159 — Học máy khử chệch trực giao Neyman hai bước Chernozhukov (Double / Debiased Machine Learning — DML) (
double_ml) khi: Áp dụng khi nghiên cứu Nhận dạng nhân quả phi tuyến với biến kiểm soát chiều cao và các bài toán thuộc nhóm tác động điều trị & nhân quả ml. (Kiểm định kèm theo: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check). - Chọn B160 — 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) (
causal_forest) khi: Áp dụng khi nghiên cứu Cá nhân hóa chính sách hỗ trợ doanh nghiệp và hộ gia đình và các bài toán thuộc nhóm tác động điều trị & nhân quả ml. (Kiểm định kèm theo: Joint Pre-Trend Wald Test, Goodman-Bacon Negative Weight Diagnostic, Rambachan–Roth HonestDiD Breakdown M*, Cattaneo rddensity Manipulation Test, SMD Balance Check).
↑ Quay lại Tổng quan Chương 09: Chương 09 — Suy luận nhân quả và đánh giá tác động chính sách
→ Đi tới Bài học đầu tiên của Nhóm 9.4: Bài B138 — Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) (
ipw)→ 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 B138: Trọng số nghịch đảo xác suất điều trị Horvitz–Thompson (Inverse Probability Weighting — IPW) (IPW — ipw)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 B139: Ước lượng bền vững kép Robins–Rotnitzky–Zhao (Augmented IPW — Doubly-Robust AIPW) (AIPW — aipw)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 B150: Phương pháp đối chứng tổng hợp Abadie–Diamond–Hainmueller (Synthetic Control Method — SCM) (Synthetic Control — synthetic_control)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 B151: Đố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) (Ridge-Augmented & Matrix-Completion Synthetic Control (ASCM / FECT) — augmented_synthetic_control)Cấp độ 3: Nâng cao (Advanced)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B159: Học máy khử chệch trực giao Neyman hai bước Chernozhukov (Double / Debiased Machine Learning — DML) (Double ML — double_ml)Cấp độ 3: Nâng cao (Advanced)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)
- Bài B160: 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) (Causal Forest — causal_forest)Cấp độ 4: Chuyên sâu (Specialized)Dữ liệu Việt Nam (EcoData)Thực hành Trực tiếp (EcoLab)Tri-Engine Verified (Stata · R · Python)