Nhóm 5.1: Dữ liệu bảng tĩnh (static)
static) là phân nhánh chuyên sâu thuộc Chương 05: Kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, quy tụ 9 phương pháp ước lượng (từ B056 đến B064). 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 5.1 (Dữ liệu bảng tĩnh)
Khi triển khai nghiên cứu định lượng trong miền Kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, chúng ta thường đứng trước câu hỏi: giữa 9 mô hình cùng thuộc nhóm phương pháp con (econometric sub-group, static), đâ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 5.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):
$$Y_{it} = rho Y_{i,t-1} + X_{it}^top beta + mu_i + lambda_t + varepsilon_{it} tag{1}$$
Sự khác biệt giữa các mô hình từ B056 (pooled_ols) đến B064 (driscoll_kraay_fe) 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 9 Mô hình thuộc Nhóm 5.1
Bảng đối chiếu dưới đây giúp bạn nhìn tổng thể toàn bộ 9 bài học trong Nhóm 5.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 |
|---|---|---|---|---|---|
| B056 | pooled_ols |
Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm (Pooled OLS & Clustered SE) Pooled OLS |
L1 | Stata: regress y x1 x2, vce(cluster id)R: plm::plm(y ~ x1 + x2, data = pdata, model = 'pooling')Python: PooledOLS.from_formula('y ~ 1 + x1 + x2', data=pdata).fit(cov_type='clustered', cluster_entity=True) |
Phương pháp Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm (Pooled OLS & Clustered SE) (Pooled OLS) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B057 | panel_fe |
Mô hình hiệu ứng cố định dữ liệu bảng và phép biến đổi Within (Panel Fixed Effects — FE) Panel Fixed Effects |
L1 | Stata: xtreg y x1 x2, fe vce(cluster id)R: plm::plm(y ~ x1 + x2, data = pdata, model = 'within')Python: PanelOLS.from_formula('y ~ 1 + x1 + x2 + EntityEffects', data=pdata).fit(cov_type='clustered', cluster_entity=True) |
Phương pháp Mô hình hiệu ứng cố định dữ liệu bảng và phép biến đổi Within (Panel Fixed Effects — FE) (Panel Fixed Effects) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B058 | panel_re |
Mô hình hiệu ứng ngẫu nhiên GLS và kiểm định Hausman (Panel Random Effects — RE) Panel Random Effects |
L1 | Stata: xtreg y x1 x2, re vce(cluster id)R: plm::plm(y ~ x1 + x2, data = pdata, model = 'random')Python: RandomEffects.from_formula('y ~ 1 + x1 + x2', data=pdata).fit(cov_type='clustered', cluster_entity=True) |
Phương pháp Mô hình hiệu ứng ngẫu nhiên GLS và kiểm định Hausman (Panel Random Effects — RE) (Panel Random Effects) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B059 | between_ols |
Ước lượng trung bình giữa các đơn vị cắt ngang (Between Effects Estimator — BE) Between OLS |
L1 | Stata: xtreg y x1 x2, beR: plm::plm(y ~ x1 + x2, data = pdata, model = 'between')Python: BetweenOLS.from_formula('y ~ 1 + x1 + x2', data=pdata).fit() |
Phương pháp Ước lượng trung bình giữa các đơn vị cắt ngang (Between Effects Estimator — BE) (Between OLS) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B060 | first_difference |
Ước lượng sai phân bậc nhất khử hiệu ứng cố định (First-Difference Estimator — FD) First Difference |
L1 | Stata: regress D.y D.(x1 x2), vce(cluster id)R: plm::plm(y ~ x1 + x2, data = pdata, model = 'fd')Python: FirstDifferenceOLS.from_formula('y ~ x1 + x2', data=pdata).fit(cov_type='clustered', cluster_entity=True) |
Phương pháp Ước lượng sai phân bậc nhất khử hiệu ứng cố định (First-Difference Estimator — FD) (First Difference) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B061 | mundlak_cre |
Mô hình hiệu ứng ngẫu nhiên tương quan Mundlak–Wooldridge (Correlated Random Effects — CRE) Mundlak / Wooldridge Correlated Random Effects (CRE) |
L2 | Stata: Stata 18 xtreg, cre / Mundlak (xtreg y x1 x2 z_fixed, cre vce(cluster id))R: plm / fixest (plm(y ~ x1 + x2 + mean_x1 + mean_x2 + z_fixed, model = "random"))Python: linearmodels (RandomEffects bổ sung $bar{X}_i$ group means + Cluster-robust Wald test trên $gamma_{bar{X}} = 0$) |
Đưa trung bình theo thời gian của từng biến thay đổi $bar{X}_i$ vào mô hình Random Effects, cho phép ước lượng chính xác hệ số within giống hệt Fixed Effects **đồng thời** giữ lại và ước lượng được các biến bất biến theo thời gian $Z_i$. |
| B062 | fama_macbeth |
Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth và hiệu chỉnh Newey–West Fama-MacBeth |
L2 | Stata: asreg y x1 x2, fmb newey(3)R: plm::pmg(y ~ x1 + x2, data = pdata, index = c('year', 'id'))Python: FamaMacBeth.from_formula('y ~ 1 + x1 + x2', data=pdata).fit(cov_type='kernel') |
Phương pháp Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth và hiệu chỉnh Newey–West (Fama-MacBeth) trong họ mô hình kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến, 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. |
| B063 | pcse_panel |
Sai số chuẩn hiệu chỉnh dữ liệu bảng Beck–Katz (Panel-Corrected Standard Errors — PCSE) Beck–Katz Panel-Corrected Standard Errors (PCSE) |
L2 | Stata: xtpcse (xtpcse y x1 x2, correlation(ar1) hetonly)R: pcse (pcse::pcse(ols_fit, groupN = df$entity, groupT = df$time))Python: Prais–Winsten $text{AR}(1)$ + Beck–Katz contemporaneous $hat{Sigma}_{Ntimes N} otimes I_T$ sandwich covariance |
Ước lượng dữ liệu bảng chuỗi thời gian – chéo (TSCS) xử lý đồng thời phương sai sai số khác nhau theo bảng (panel heteroskedasticity), tương quan đồng thời giữa các đơn vị và tự tương quan $text{AR}(1)$. |
| B064 | driscoll_kraay_fe |
Ước lượng bảng với sai số chuẩn bền vững không gian–thời gian Driscoll–Kraay (xtscc) Driscoll–Kraay Spatial-Temporal SCC Panel |
L2 | Stata: xtscc (ssc install xtscc: xtscc y x1 x2, fe lag(m))R: plm (lmtest::coeftest(plm_fe, vcov = plm::vcovSCC(plm_fe, maxlag = m)))Python: linearmodels (PanelOLS(..., entity_effects=True).fit(cov_type="kernel", kernel="bartlett", bandwidth=m)) |
Ước lượng bảng hiệu ứng cố định/gộp với ma trận hiệp phương sai Newey–West trên chuỗi trung bình chéo của điều kiện moment, bền vững trước cả tự tương quan, phương sai thay đổi và phụ thuộc không gian/chéo tổng quát. |
3. Hướng dẫn Chọn nhanh Mô hình Phù hợp trong Nhóm 5.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 B056 — Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm (Pooled OLS & Clustered SE) (
pooled_ols) khi: Áp dụng khi nghiên cứu Bảng đầu tư & giá trị doanh nghiệp niêm yết HOSE và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B057 — Mô hình hiệu ứng cố định dữ liệu bảng và phép biến đổi Within (Panel Fixed Effects — FE) (
panel_fe) khi: Áp dụng khi nghiên cứu Tăng trưởng GRDP và đầu tư công 63 tỉnh thành GSO và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B058 — Mô hình hiệu ứng ngẫu nhiên GLS và kiểm định Hausman (Panel Random Effects — RE) (
panel_re) khi: Áp dụng khi nghiên cứu Hiệu quả hoạt động ngân hàng thương mại Việt Nam và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B059 — Ước lượng trung bình giữa các đơn vị cắt ngang (Between Effects Estimator — BE) (
between_ols) khi: Áp dụng khi nghiên cứu Khác biệt cấu trúc dài hạn giữa các tỉnh thành và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B060 — Ước lượng sai phân bậc nhất khử hiệu ứng cố định (First-Difference Estimator — FD) (
first_difference) khi: Áp dụng khi nghiên cứu Bảng hộ gia đình VHLSS qua các kỳ điều tra và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B061 — Mô hình hiệu ứng ngẫu nhiên tương quan Mundlak–Wooldridge (Correlated Random Effects — CRE) (
mundlak_cre) khi: Khi nhà nghiên cứu cần ước lượng cả tác động của biến cố định theo thực thể (khoảng cách địa lý, thể chế gốc, giới tính, ngành) mà không muốn chịu chệch tương quan $text{Cov}(alpha_i, X_{it}) neq 0$ của RE chuẩn. (Kiểm định kèm theo: Kiểm định Wald cụm đa chiều $H_0: gamma_{bar{X}} = 0$ (kiểm định Mundlak bền vững dưới phương sai thay đổi — vượt trội Hausman cổ điển khi clustered SE).). - Chọn B062 — Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth và hiệu chỉnh Newey–West (
fama_macbeth) khi: Áp dụng khi nghiên cứu Định giá chéo lợi suất cổ phiếu Việt Nam và các bài toán thuộc nhóm dữ liệu bảng tĩnh. (Kiểm định kèm theo: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR(2) & Hansen J, Bootstrap Threshold F-test). - Chọn B063 — Sai số chuẩn hiệu chỉnh dữ liệu bảng Beck–Katz (Panel-Corrected Standard Errors — PCSE) (
pcse_panel) khi: Bảng kinh tế chính trị, tài khóa công, vĩ mô so sánh có số đơn vị $N$ vừa phải ($N in [10, 50]$) và chuỗi thời gian dài hơn hoặc tương đương ($T ge N$, nơi FGLS Parks thổi phồng độ tin cậy giả). (Kiểm định kèm theo: Hệ số tự tương quan $hat{rho}$ chung hoặc riêng theo đơn vị (psar1), Wald $chi^2$ toàn mô hình, Breusch–Pagan LM kiểm định tương quan đồng thời.). - Chọn B064 — Ước lượng bảng với sai số chuẩn bền vững không gian–thời gian Driscoll–Kraay (xtscc) (
driscoll_kraay_fe) khi: Rất phù hợp cho bảng tỉnh thành Việt Nam (63/34 tỉnh) hoặc bảng ngành/quốc gia có tương quan không gian – thời gian phức tạp mà không cần biết trước ma trận trọng số không gian $W$. (Kiểm định kèm theo: Kiểm định Pesaran CD, Wooldridge panel autocorrelation, within-$R^2$, độ rộng băng thông Bartlett tự động $m = lfloor 4(T/100)^{2/9} rfloor$.).
↑ Quay lại Tổng quan Chương 05: Chương 05 — Kinh tế lượng dữ liệu bảng: tĩnh, động, phụ thuộc chéo và phi tuyến
→ Đi tới Bài học đầu tiên của Nhóm 5.1: Bài B056 — Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm (Pooled OLS & Clustered SE) (
pooled_ols)→ 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 B056: Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm (Pooled OLS & Clustered SE) (Pooled OLS — pooled_ols)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 B057: Mô hình hiệu ứng cố định dữ liệu bảng và phép biến đổi Within (Panel Fixed Effects — FE) (Panel Fixed Effects — panel_fe)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 B058: Mô hình hiệu ứng ngẫu nhiên GLS và kiểm định Hausman (Panel Random Effects — RE) (Panel Random Effects — panel_re)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 B059: Ước lượng trung bình giữa các đơn vị cắt ngang (Between Effects Estimator — BE) (Between OLS — between_ols)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 B060: Ước lượng sai phân bậc nhất khử hiệu ứng cố định (First-Difference Estimator — FD) (First Difference — first_difference)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 B061: Mô hình hiệu ứng ngẫu nhiên tương quan Mundlak–Wooldridge (Correlated Random Effects — CRE) (Mundlak / Wooldridge Correlated Random Effects (CRE) — mundlak_cre)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 B062: Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth và hiệu chỉnh Newey–West (Fama-MacBeth — fama_macbeth)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 B063: Sai số chuẩn hiệu chỉnh dữ liệu bảng Beck–Katz (Panel-Corrected Standard Errors — PCSE) (Beck–Katz Panel-Corrected Standard Errors (PCSE) — pcse_panel)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 B064: Ước lượng bảng với sai số chuẩn bền vững không gian–thời gian Driscoll–Kraay (xtscc) (Driscoll–Kraay Spatial-Temporal SCC Panel — driscoll_kraay_fe)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)