5. Kinh tế lượng dữ liệu bảng tĩnh, động và phụ thuộc chéo
TÓM TẮT CHƯƠNG 5:
- Quy mô chương: Chương 5. Kinh tế lượng dữ liệu bảng tĩnh, động và phụ thuộc chéo bao gồm 6 nhóm phương pháp và 28 bài giảng chuyên sâu.
- Bài toán nhận dạng cốt lõi: “Làm sao để kiểm soát đại lượng không quan sát được theo cá thể, xử lý nội sinh động (GMM), phụ thuộc chéo thế hệ 2 (ARDL) và hiệu ứng ngưỡng trên dữ liệu bảng?”
- Hệ thống kiểm định bắt buộc: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR bậc 2 & Hansen J, Bootstrap Threshold F-test.
- Nguồn tham chiếu học thuật: Tổng hợp từ các công trình gốc của Arellano & Bond (1991); Blundell & Bond (1998); Pesaran (2006, 2007); Chudik & Pesaran (2015); Hansen (1999); Machado & Santos Silva (2019) và hệ thống giáo trình chuẩn quốc tế của Baltagi (2021); Wooldridge (2010); Pesaran (2015); Verbeek (2021).
1. Đặt vấn đề và bối cảnh nghiên cứu
Trong nghiên cứu thực nghiệm kinh tế, tài chính và chính sách công, việc lựa chọn đúng cấu trúc mô hình thuộc chương Kinh tế lượng dữ liệu bảng tĩnh, động và phụ thuộc chéo quyết định trực tiếp tính không chệch, tính nhất quán và hiệu lực suy luận thống kê của các ước lượng theo Baltagi (2021); Wooldridge (2010); Pesaran (2015); Verbeek (2021). Khi cấu trúc dữ liệu thực tế vi phạm các giả định cổ điển, nhà nghiên cứu cần áp dụng chính xác các phương pháp chuyên biệt được trình bày trong chương này theo Arellano & Bond (1991); Blundell & Bond (1998); Pesaran (2006, 2007); Chudik & Pesaran (2015); Hansen (1999); Machado & Santos Silva (2019).
Khung biểu diễn toán học tổng quát cho các mô hình thuộc chương 5 được xác lập tại phương trình 1:
$$Y_{it} = \rho Y_{i,t-1} + X_{it}^\top \beta + \mu_i + \lambda_t + \varepsilon_{it} \tag*{[1]}$$
Từ phương trình 1, mỗi nhóm phương pháp con trong chương sẽ cụ thể hóa điều kiện nhận dạng, cấu trúc ma trận hiệp phương sai của sai số và thuật toán tối ưu hóa phù hợp với từng thiết kế dữ liệu tại Việt Nam theo Arellano & Bond (1991); Blundell & Bond (1998); Pesaran (2006, 2007); Chudik & Pesaran (2015); Hansen (1999); Machado & Santos Silva (2019); Baltagi (2021); Wooldridge (2010); Pesaran (2015); Verbeek (2021).
2. Cấu trúc các nhóm phương pháp trong chương
Chương 5 được tổ chức thành 6 nhóm bài giảng theo thứ tự từ nền tảng đến nâng cao:
Bảng 5.1Cấu trúc các nhóm phương pháp trực thuộc chương 5
3. Danh mục các mô hình kinh tế lượng
Bảng dưới đây đối chiếu toàn bộ 28 bài giảng thuộc chương 5 cùng bộ thư viện thực thi trên Stata 18, R 4.x, Python 3.11+ và các kiểm định chẩn đoán bắt buộc:
Bảng 5.2Danh mục các mô hình kinh tế lượng và thư viện thực thi thuộc chương 5
| Mã mục | Tên bài giảng | Thư viện thực thi Stata · R · Python | Kiểm định chẩn đoán bắt buộc | Dữ liệu thực nghiệm Việt Nam |
|---|---|---|---|---|
| 5.1.1 | Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm | base Stata · plm · linearmodels |
Kiểm định nhân tử Lagrange Breusch–Pagan cho hiệu ứng ngẫu nhiên (LM), kiểm định $F$ cho hiệu ứng cố định và sai số chuẩn phân cụm theo đơn vị | Bảng đầu tư & giá trị doanh nghiệp niêm yết HOSE |
| 5.1.2 | Mô hình hiệu ứng cố định và phép biến đổi nội nhóm (FE) | base Stata · plm / fixest · linearmodels |
Kiểm định $F$ cho sự tồn tại của hiệu ứng cố định cá thể, kiểm định Hausman / Mundlak đối chiếu (FE) với (RE), kiểm định phương sai thay đổi Modified Wald và tự tương quan Wooldridge | Tăng trưởng GRDP và đầu tư công 63 tỉnh thành GSO |
| 5.1.3 | Mô hình hiệu ứng ngẫu nhiên và kiểm định Hausman (RE) | base Stata · plm · linearmodels |
Kiểm định Breusch–Pagan (LM) đối chiếu (RE) và Pooled (OLS), kiểm định Hausman và kiểm định Hausman vững phân cụm Mundlak–Wooldridge | Hiệu quả hoạt động ngân hàng thương mại Việt Nam |
| 5.1.4 | Ước lượng trung bình giữa các đơn vị cắt ngang (BE) | base Stata · plm · linearmodels |
Kiểm định phương sai thay đổi White trên hồi quy trung bình nhóm và so sánh sự chênh lệch hệ số Between (BE) với Within (FE) để phát hiện nội sinh cấp đơn vị | Khác biệt cấu trúc dài hạn giữa các tỉnh thành |
| 5.1.5 | Ước lượng sai phân bậc nhất khử hiệu ứng cố định (FD) | base Stata · plm · linearmodels |
Kiểm định tự tương quan bậc nhất của phần dư sai phân $\Delta \varepsilon_{it}$, dưới $H_0$ không tự tương quan ở mức thì $\text{Corr}(\Delta \varepsilon_{it}, \Delta \varepsilon_{i,t-1}) = -0.5$, và kiểm định ngoại sinh chặt | Bảng hộ gia đình VHLSS qua các kỳ điều tra |
| 5.1.6 | Mô hình hiệu ứng ngẫu nhiên tương quan Mundlak (CRE) | Stata 18 xtreg · plm · linearmodels |
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. | Giữ biến cố định vùng miền/ngành trong mô hình FE |
| 5.1.7 | Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth | Stata asreg / xtreg · plm · linearmodels |
Hiệu chỉnh sai số chuẩn Newey–West (HAC) cho chuỗi thời gian hệ số $\hat{\beta}_t$, thống kê $t$ Fama–MacBeth và $R^2$ cắt ngang trung bình | Định giá chéo lợi suất cổ phiếu Việt Nam |
| 5.1.8 | Sai số chuẩn hiệu chỉnh dữ liệu bảng Beck–Katz (PCSE) | base Stata xtpcse · pcse · linearmodels |
Hệ số tự tương quan $\hat{\rho}$ chung hoặc riêng theo đơn vị, Wald $\chi^2$ toàn mô hình, Breusch–Pagan LM kiểm định tương quan đồng thời. | Bảng tài khóa – tăng trưởng ASEAN — $T \ge N$ — |
| 5.1.9 | Sai số chuẩn bền vững không gian–thời gian Driscoll–Kraay | xtscc · plm · linearmodels |
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$. | Bảng 63 tỉnh thành Việt Nam có tương quan chéo |
| 5.2.1 | Ước lượng GMM sai phân dữ liệu bảng động Arellano–Bond | base Stata / xtabond2 · plm · pydynpd / linearmodels |
Kiểm định tự tương quan bậc 2 của phần dư sai phân Arellano–Bond AR bậc 2 yêu cầu $p \gt 0.05$, kiểm định Hansen $J$ cho tính hợp lệ của biến công cụ và quy tắc số công cụ nhỏ hơn hoặc bằng số nhóm $N$ | Động học cấu trúc vốn doanh nghiệp Việt Nam |
| 5.2.2 | Ước lượng GMM hệ thống dữ liệu bảng động Blundell–Bond | base Stata / xtabond2 · plm · pydynpd / linearmodels |
Kiểm định Arellano–Bond AR bậc 1 và AR bậc 2, kiểm định Hansen $J$ toàn hệ thống, kiểm định Difference-in-Hansen cho tính hợp lệ của các công cụ phương trình mức và hiệu chỉnh sai số chuẩn hữu hạn mẫu Windmeijer (2005) | Tăng trưởng kinh tế cấp tỉnh và tín dụng động |
| 5.2.3 | Ước lượng biến giả hiệu chỉnh chệch Nickell (LSDVC) | Stata xtlsdvc · plm · statsmodels / numpy |
Đánh giá độ hội tụ từ bộ ước lượng khởi tạo Arellano–Bond hoặc Blundell–Bond và khoảng tin cậy Bootstrap tham số | Bảng 25 ngân hàng thương mại Việt Nam — $N$ nhỏ — |
| 5.3.1 | Ước lượng biến công cụ hai giai đoạn dữ liệu bảng | base Stata · plm · linearmodels |
Thống kê $F$ giai đoạn 1 Kleibergen–Paap rk Wald $F$ trên dữ liệu đã khử trung bình cá thể, kiểm định nội sinh Davidson–MacKinnon / Hausman và kiểm định Hansen $J$ | Nội sinh FDI và năng suất lao động cấp tỉnh |
| 5.3.2 | Ước lượng phương pháp moment tổng quát dữ liệu bảng (GMM) | Stata xtivreg2 · plm / fixest · linearmodels |
Kiểm định nhận dạng yếu Kleibergen–Paap rk (LM) và Wald $F$, kiểm định Hansen $J$ overidentification và kiểm định $C$ difference-in-Sargan cho tính ngoại sinh của tập con công cụ | Quản trị công ty và hiệu quả tài chính под phương sai thay đổi |
| 5.3.3 | Ước lượng biến công cụ nội tại Hausman–Taylor | base Stata · plm · linearmodels / statsmodels |
Kiểm định nhận dạng quá mức Sargan–Hansen yêu cầu số biến ngoại sinh biến thiên theo thời gian $k_1 \ge g_2$ số biến nội sinh bất biến theo thời gian, và kiểm định Hausman đối chiếu với (FE) | Thương mại và khoảng cách địa lý nội sinh |
| 5.4.1 | Hồi quy OLS hiệu chỉnh toàn phần trên bảng đồng tích hợp (FMOLS) | Stata xtcointreg · plm / cointReg · statsmodels / numpy |
Kiểm định đồng tích hợp bảng Pedroni / Kao trước ước lượng, kiểm tra độ nhạy đối với cửa sổ hạt nhân Bartlett/Parzen và so sánh hệ số Pooled vs Group-Mean (FMOLS) | Tiêu thụ năng lượng, GDP và phát thải $\text{CO}_2$ ASEAN |
| 5.4.2 | Hồi quy OLS động dẫn–trễ trên bảng đồng tích hợp (DOLS) | Stata xtcointreg · plm / cointReg · statsmodels / numpy |
Lựa chọn độ dài dẫn–trễ tối ưu theo (AIC) và (BIC), sai số chuẩn dài hạn (HAC) và kiểm định tính dừng của phần dư đồng tích hợp | Độ mở thương mại và tăng trưởng dài hạn |
| 5.4.3 | Ước lượng trung bình nhóm gộp dữ liệu bảng (PMG) | xtpmg · plm · statsmodels / scipy |
Hệ số hiệu chỉnh sai số trung bình $\bar{\phi} \lt 0$ và có ý nghĩa thống kê, kiểm định Hausman so sánh PMG vs MG, kiểm định giả thuyết đồng nhất dài hạn. | Truyền dẫn lãi suất dài hạn các ngân hàng |
| 5.4.4 | Ước lượng trung bình nhóm hệ số không đồng nhất (MG) | xtmg · plm · statsmodels |
Kiểm định tính không đồng nhất hệ số dốc Swamy / Pesaran–Yamagata $\tilde{\Delta}$, kiểm định Pesaran CD trên phần dư, Hausman MG vs PMG/DFE. | Độ co giãn thu nhập theo từng quốc gia |
| 5.4.5 | Mô hình hiệu chỉnh sai số hiệu ứng cố định động (DFE) | xtpmg · fixest · linearmodels |
Tốc độ hội tụ $\phi \in (-2, 0)$, Delta-method $z$-test cho hệ số dài hạn $\theta$, kiểm định Hausman tam giác MG vs PMG vs DFE. | Đầu tư và tiết kiệm vùng kinh tế |
| 5.5.1 | Kiểm định nghiệm đơn vị dữ liệu bảng thế hệ 2 Pesaran (CIPS) | xtcips · plm · statsmodels |
Thống kê CIPS và CIPS, giá trị tới hạn chính xác theo cặp kích thước mẫu $N$ và $T$ ở mức $1\%, 5\%, 10\%$ — không chặn / có chặn / có xu thế , phân bố $t_i$ CADF. | Kiểm định tính dừng vĩ mô các nền kinh tế mới nổi |
| 5.5.2 | Kiểm định đồng tích hợp bảng hiệu chỉnh sai số Westerlund | xtwest · plm · linearmodels |
Bộ 4 thống kê $G_\tau, G_\alpha, P_\tau, P_\alpha$, giá trị $Z$-stat tiệm cận chuẩn và robust bootstrap $p$-value dưới phụ thuộc chéo. | Đồng tích hợp tài chính – tăng trưởng dưới cú sốc chung |
| 5.5.3 | Ước lượng trung bình nhóm hiệu ứng tương quan chung (CCEMG) | xtdcce2 · plm · statsmodels |
Kiểm định Pesaran CD trước và sau lọc CCE — kỳ vọng phần dư hết phụ thuộc chéo , bậc mũ phụ thuộc chéo $\hat{\alpha}_{\text{CD}}$, kiểm định $\tilde{\Delta}$ slope homogeneity. | Tăng trưởng & chuyển dịch năng lượng xanh |
| 5.5.4 | Ước lượng gộp hiệu ứng tương quan chung (CCEP) | xtdcce2 · plm · linearmodels |
Kiểm định Pesaran CD phần dư, Pesaran–Yamagata $\hat{\Delta}_{\text{adj}}$ kiểm chứng điều kiện gộp hệ số, $R^2_{\text{MG/CCEP}}$. | Bảng doanh nghiệp/ngân hàng chịu cú sốc vĩ mô toàn cầu |
| 5.5.5 | Mô hình ARDL bảng tăng cường trung bình chéo (CS-ARDL) | xtdcce2 · plm · statsmodels |
Kiểm định CD phần dư sau gắn độ trễ trung bình chéo, hệ số hội tụ ECM $\hat{\phi} \lt 0$, kiểm định nghiệm đơn vị CIPS trên phần dư dài hạn. | Nợ công và tăng trưởng dài hạn dưới phụ thuộc chéo |
| 5.6.1 | Mô hình hồi quy ngưỡng dữ liệu bảng tĩnh Hansen | xthreg · pdR · linearmodels |
Kiểm định Bootstrap $F_1$ tồn tại hiệu ứng ngưỡng, vì $\gamma$ không định danh dưới $H_0$, khoảng tin cậy tỷ số hợp lý $LR_1(\gamma) \le c(\alpha)$ cho điểm ngưỡng $\hat{\gamma}$. | Ngưỡng đòn bẩy nợ tối ưu của doanh nghiệp |
| 5.6.2 | Mô hình ngưỡng dữ liệu bảng động có nội sinh Seo–Shin | xthenreg · gmm · linearmodels |
Kiểm định Sup-Wald bootstrap về tính tuyến tính / tồn tại bước nhảy ngưỡng, kiểm định Hansen $J$ tính hợp lệ biến công cụ trễ, khoảng tin cậy $95\%$ cho ngưỡng $\hat{\gamma}$. | Ngưỡng lạm phát / nợ xấu trong bảng động |
| 5.6.3 | Hồi quy phân vị dữ liệu bảng theo phương pháp moment (MM-QR) | xtqreg · quantreg · statsmodels |
Hệ số Location — $\hat{\beta}$ , hệ số Scale — $\hat{\gamma}$ , phân vị chuẩn hóa $\hat{q}(\tau)$, kiểm định Wald tính bất biến theo phân vị — $\gamma = 0$ . | Tác động FDI & đổi mới công nghệ trên các phân vị phát thải |
4. Quy tắc chọn mô hình và kiểm định
Trước khi công bố kết quả ước lượng thuộc chương 5, nhà nghiên cứu cần thực hiện đầy đủ danh mục kiểm tra học thuật sau theo Baltagi (2021); Wooldridge (2010); Pesaran (2015); Verbeek (2021):
- Câu hỏi nhận dạng trung tâm: Làm sao để kiểm soát đại lượng không quan sát được theo cá thể, xử lý nội sinh động (GMM), phụ thuộc chéo thế hệ 2 (ARDL) và hiệu ứng ngưỡng trên dữ liệu bảng?
- Hệ thống kiểm định chẩn đoán bắt buộc: Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR bậc 2 & Hansen J, Bootstrap Threshold F-test.
- Suy luận thống kê vững: Luôn kiểm tra cấu trúc phương sai thay đổi, tự tương quan hoặc phụ thuộc cụm để lựa chọn ma trận hiệp phương sai sai số chuẩn phù hợp trên Stata 18, R 4.x và Python 3.11+ theo Arellano & Bond (1991); Blundell & Bond (1998); Pesaran (2006, 2007); Chudik & Pesaran (2015); Hansen (1999); Machado & Santos Silva (2019).
LƯU Ý HỌC THUẬT: Khi trình bày bảng kết quả ước lượng thuộc chương 5 trong khóa luận, luận văn, luận án hoặc bài báo khoa học, cần báo cáo đầy đủ giá trị thống kê và $p$-value của bộ kiểm định gồm Hausman & Mundlak CRE Test, Pesaran CD Test, Pesaran CIPS Unit Root, Westerlund Cointegration, Arellano–Bond AR bậc 2 & Hansen J, Bootstrap Threshold F-test ngay dưới chân bảng để minh chứng mô hình thỏa mãn điều kiện nhận dạng theo Baltagi (2021); Wooldridge (2010); Pesaran (2015); Verbeek (2021).
5. Tài liệu tham khảo
- Arellano, M, & Bond, S.(1991). Some tests of specification for panel data: Monte Carlo evidence and an application to employment equations. The Review of Economic Studies, 58, số 2, 277–297.
- Blundell, R, & Bond, S.(1998). Initial conditions and moment restrictions in dynamic panel data models. Journal of Econometrics, 87, số 1, 115–143.
- Pesaran(2006, 2007). Econometric specification, identification, and asymptotic inference.
- Chudik, & Pesaran(2015). Econometric specification, identification, and asymptotic inference.
- Hansen(1999). Econometric specification, identification, and asymptotic inference.
- Machado, & Santos Silva(2019). Econometric specification, identification, and asymptotic inference.
- Baltagi, B. H.(2021). Econometric Analysis of Panel Data, 6th ed, Chương 2–12: Static, Dynamic and Cross-Sectionally Dependent Panels. Springer.
- Wooldridge, J. M.(2010). Econometric Analysis of Cross Section and Panel Data, 2nd ed, Chương 10–11. MIT Press.
- Pesaran, M. H.(2015). Time Series and Panel Data Econometrics, Chương 26–29: Large Heterogeneous and Cross-Sectionally Dependent Panels. Oxford University Press.
- Verbeek, M.(2021). Panel Methods for Finance: A Guide to Empirical Practice. De Gruyter.
ĐIỀU HƯỚNG HỌC TẬP:
↑ Tài liệu gốc: KINH TẾ LƯỢNG ỨNG DỤNG TOÀN DIỆN TỪ A – Z
← Chương trước: 4. Mô hình đếm, lực hấp dẫn thương mại và biến phụ thuộc giới hạn
→ Nhóm đầu tiên trong chương: 5.1. Dữ liệu bảng tĩnh
→ Bài học đầu tiên trong chương: 5.1.1. Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm
→ Chương tiếp theo: 6. Chuỗi thời gian đơn biến, phi tuyến, điểm gãy và bong bóng
↑ Tài liệu gốc: KINH TẾ LƯỢNG ỨNG DỤNG TOÀN DIỆN TỪ A – Z
← Chương trước: 4. Mô hình đếm, lực hấp dẫn thương mại và biến phụ thuộc giới hạn
→ Nhóm đầu tiên trong chương: 5.1. Dữ liệu bảng tĩnh
→ Bài học đầu tiên trong chương: 5.1.1. Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm
→ Chương tiếp theo: 6. Chuỗi thời gian đơn biến, phi tuyến, điểm gãy và bong bóng
Bài liên quan
- 5.1. Dữ liệu bảng tĩnh
- 5.1.1. Hồi quy gộp dữ liệu bảng và sai số chuẩn phân cụm
- 5.1.2. Mô hình hiệu ứng cố định và phép biến đổi nội nhóm (FE)
- 5.1.3. Mô hình hiệu ứng ngẫu nhiên và kiểm định Hausman (RE)
- 5.1.4. Ước lượng trung bình giữa các đơn vị cắt ngang (BE)
- 5.1.5. Ước lượng sai phân bậc nhất khử hiệu ứng cố định (FD)
- 5.1.6. Mô hình hiệu ứng ngẫu nhiên tương quan Mundlak (CRE)
- 5.1.7. Hồi quy hai bước theo kỳ cắt ngang Fama–MacBeth
- 5.1.8. Sai số chuẩn hiệu chỉnh dữ liệu bảng Beck–Katz (PCSE)
- 5.1.9. Sai số chuẩn bền vững không gian–thời gian Driscoll–Kraay
- 5.2. Dữ liệu bảng động
- 5.3. Biến công cụ dữ liệu bảng
- 5.4. Đồng tích hợp và ARDL dữ liệu bảng
- 5.4.1. Hồi quy OLS hiệu chỉnh toàn phần trên bảng đồng tích hợp (FMOLS)
- 5.4.2. Hồi quy OLS động dẫn–trễ trên bảng đồng tích hợp (DOLS)
- 5.4.3. Ước lượng trung bình nhóm gộp dữ liệu bảng (PMG)
- 5.4.4. Ước lượng trung bình nhóm hệ số không đồng nhất (MG)
- 5.4.5. Mô hình hiệu chỉnh sai số hiệu ứng cố định động (DFE)
- 5.5. Dữ liệu bảng phụ thuộc chéo thế hệ 2
- 5.5.1. Kiểm định nghiệm đơn vị dữ liệu bảng thế hệ 2 Pesaran (CIPS)
- 5.5.2. Kiểm định đồng tích hợp bảng hiệu chỉnh sai số Westerlund
- 5.5.3. Ước lượng trung bình nhóm hiệu ứng tương quan chung (CCEMG)
- 5.5.4. Ước lượng gộp hiệu ứng tương quan chung (CCEP)
- 5.5.5. Mô hình ARDL bảng tăng cường trung bình chéo (CS-ARDL)
- 5.6. Ngưỡng và phân vị dữ liệu bảng