10. Kinh tế lượng không gian, sai số chuẩn Conley và mạng cặp đôi
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 không gian, sai số chuẩn Conley và mạng cặp đôi 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 (Cameron & Trivedi (2022); Baltagi (2021); Pesaran (2015); Greene (2019)). 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 (Anselin (1988); LeSage & Pace (2009); Elhorst (2014); Fotheringham, Yang & Kang (2017); Conley (1999); Bramoullé, Djebbari & Fortin (2009)).
Khung biểu diễn toán học tổng quát cho các mô hình thuộc chương 10 được xác lập tại phương trình (1):
$$Y = \rho W Y + X\beta + W X \theta + u, \quad u = \lambda W u + \varepsilon \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 (Anselin (1988); LeSage & Pace (2009); Elhorst (2014); Fotheringham, Yang & Kang (2017); Conley (1999); Bramoullé, Djebbari & Fortin (2009); Cameron & Trivedi (2022); Baltagi (2021); Pesaran (2015); Greene (2019)).
2. Cấu trúc các nhóm phương pháp trong chương
Chương 10 được tổ chức thành 4 nhóm bài giảng theo thứ tự từ nền tảng đến nâng cao:
3. Danh mục các mô hình kinh tế lượng
Bảng dưới đây đối chiếu toàn bộ 10 bài giảng thuộc chương 10 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:
| 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 |
|---|---|---|---|---|
| 10.1.1 | Mô hình trễ không gian của biến giải thích (SLX) | Stata spregress · spatialreg · spreg |
Bảng phân rã Direct, Indirect (Spillover $\theta$) và Total Effects kèm $p$-value, kiểm định Wald đồng thời $H_0: \theta = 0$, kiểm định Moran’s $I$ & LM-Error trên phần dư. | Lan tỏa đầu tư hạ tầng giữa các tỉnh giáp ranh Việt Nam |
| 10.1.2 | Mô hình tự hồi quy trễ không gian của biến phụ thuộc (SAR) | base Stata sp · spatialreg · spreg / scipy |
Kiểm định tự tương quan không gian toàn cục Moran’s I, kiểm định nhân tử Lagrange Anselin LM-Lag và Robust LM-Lag, điều kiện ổn định không gian |ρ| < 1 / λ_max(W) và phân rã tác động Trực tiếp / Gián tiếp / Tổng của LeSage–Pace | Tương tác tăng trưởng GRDP 63 tỉnh thành |
| 10.1.3 | Mô hình tự tương quan không gian trong sai số (SEM) | base Stata sp · spatialreg · spreg / scipy |
Kiểm định Moran’s I phần dư OLS, kiểm định Anselin LM-Error và Robust LM-Error, kiểm định Wald/LR cho H_0: λ = 0 và kiểm định giả thuyết nhân tử chung (Common Factor test) so với SDM | Cú sốc thời tiết/môi trường lan tỏa không gian |
| 10.1.4 | Mô hình không gian Durbin và phân rã tác động LeSage–Pace (SDM) | base Stata sp · spatialreg · spreg / scipy |
Phân rã ma trận đạo hàm riêng (I – ρ W)^(-1)(I β_k + W θ_k) thành Tác động Trực tiếp (Direct), Tác động Lan tỏa Gián tiếp (Indirect Spillover) và Tổng tác động (Total) kèm sai số chuẩn mô phỏng MCMC/Delta, và kiểm định LR rút gọn về SAR (H_0: θ = 0) hoặc SEM (H_0: θ + ρ β = 0) | Tác động trực tiếp và tràn của FDI cấp tỉnh |
| 10.1.5 | Mô hình kết hợp tự hồi quy trễ và sai số không gian (SARAR) | Stata spregress · spatialreg · spreg |
Kiểm định ý nghĩa riêng rẽ và đồng thời của $(\rho, \lambda)$, phân rã LeSage–Pace Direct/Indirect/Total impacts có hiệu chỉnh vết ma trận $(I-\rho W)^{-1}$, Anselin–Kelejian Moran dư. | Mô hình hóa đồng thời tương tác giá đất và nhiễu vùng |
| 10.2.1 | Mô hình dữ liệu bảng không gian hiệu ứng cố định | base Stata sp · splm / spatialreg · linearmodels / scipy |
Kiểm định Hausman không gian của Mutl–Pfaffermayr (Spatial FE vs Spatial RE), kiểm định LR cho hiệu ứng cố định không gian và thời gian, và phân rã tác động trực tiếp / gián tiếp ngắn hạn và dài hạn | Bảng không gian 63 tỉnh thành giai đoạn 2010–2025 |
| 10.2.2 | Mô hình dữ liệu bảng không gian hiệu ứng ngẫu nhiên | base Stata sp · splm / spatialreg · linearmodels / scipy |
Kiểm định Hausman không gian (nếu p > 0.05 thì chọn Spatial RE), kiểm định LM không gian cho hiệu ứng ngẫu nhiên của Baltagi–Song–Koh và phân rã tác động biên không gian LeSage–Pace | Kiểm định Hausman không gian FE vs RE |
| 10.3.1 | Hồi quy trọng số địa lý và đa thang đo (GWR, MGWR) | gwr / Mata adaptive Gaussian/bisquare kernel · GWmodel · mgwr |
Băng thông tối ưu riêng $bw_k$ cho từng biến, $\text{AIC}_c$ so sánh MGWR vs GWR vs Global OLS, kiểm định Monte-Carlo / Leung $F_3$ về sự biến thiên không gian thực sự của từng $\beta_k(u_i, v_i)$. | Biến thiên hệ số theo tọa độ địa lý các tỉnh/huyện Việt Nam |
| 10.3.2 | Sai số chuẩn bền vững không gian–thời gian Conley và mạng cặp đôi | acreg · fixest · Conley Haversine Bartlett Spatial-Temporal HAC & Aronow–Samii–Assenova Dyadic Sandwich Estimator |
Đối chiếu tỷ số độ phóng đại sai số chuẩn $\text{SE}_{\text{Conley}} / \text{SE}_{\text{HC1}}$ (lượng hóa mức độ tương quan không gian ngầm), số cặp láng giềng trung bình trong bán kính $d_{\text{cutoff}}$, kiểm định nhạy cảm qua đa bán kính {50, 100, 250} km. | Hồi quy dữ liệu tọa độ GPS và thương mại cặp đôi |
| 10.4.1 | Mô hình hiệu ứng đồng đẳng trên mạng lưới Bramoullé | base Stata · AER / spatialreg · linearmodels |
Kiểm định điều kiện nhận dạng độc lập tuyến tính của các ma trận I, G, G^2 (Bramoullé–Djebbari–Fortin condition), thống kê F giai đoạn 1 cho biến công cụ G^2 X và kiểm định Hansen J | Lan tỏa quyết định tài chính trong mạng lưới chuỗi cung ứng |
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 10, nhà nghiên cứu cần thực hiện đầy đủ danh mục kiểm tra học thuật sau (Cameron & Trivedi (2022); Baltagi (2021); Pesaran (2015); Greene (2019)):
- Câu hỏi nhận dạng trung tâm: Làm thế nào để lượng hóa hiệu ứng lan tỏa địa lý (Direct vs Indirect Spillovers), biến thiên hệ số đa thang đo (MGWR) và sai số chuẩn không gian Conley HAC / mạng cặp đôi?
- Hệ thống kiểm định chẩn đoán bắt buộc: Global Moran’s I & Local LISA, Anselin LM-Lag & Robust LM-Error, LeSage–Pace Direct/Indirect/Total Impacts, AICc Bandwidth Selection, Conley Distance Cutoff Sensitivity.
- 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+ (Anselin (1988); LeSage & Pace (2009); Elhorst (2014); Fotheringham, Yang & Kang (2017); Conley (1999); Bramoullé, Djebbari & Fortin (2009)).
5. Tài liệu tham khảo
- Anselin, L. (1988). Spatial Econometrics: Methods and Models. Kluwer Academic Publishers.
- LeSage, J., & Pace, R. K. (2009). Introduction to Spatial Econometrics. Chapman & Hall/CRC.
- Elhorst, J. P. (2014). Spatial Econometrics: From Cross-Sectional Data to Spatial Panels. Springer.
- Fotheringham, Yang, & Kang (2017). Econometric specification, identification, and asymptotic inference.
- Conley (1999). Econometric specification, identification, and asymptotic inference.
- Bramoullé, Y., Djebbari, H., & Fortin, B. (2009). Identification of peer effects through social networks. Journal of Econometrics, 150(1), 41–55.
- Cameron, A. C., & Trivedi, P. K. (2022). Microeconometrics Using Stata (2nd ed., Vol. 2, Chương 26: Spatial Regression Models). Stata Press.
- Baltagi, B. H. (2021). Econometric Analysis of Panel Data (6th ed., Chương 13: Spatial Panel Data Models). Springer.
- Pesaran, M. H. (2015). Time Series and Panel Data Econometrics (Chương 30: Spatial Econometrics and Network Interdependence). Oxford University Press.
- Greene, W. H. (2019). Econometric Analysis (8th ed., Chương 11: Spatial Autoregression and Robust Covariance Estimation). Pearson.
↑ Tài liệu gốc: KINH TẾ LƯỢNG ỨNG DỤNG TOÀN DIỆN TỪ A – Z
← Chương trước: 9. Suy luận nhân quả và đánh giá tác động chính sách
→ Nhóm đầu tiên trong chương: 10.1. Hồi quy không gian cắt ngang
→ Bài học đầu tiên trong chương: 10.1.1. Mô hình trễ không gian của biến giải thích (SLX)
→ Chương tiếp theo: 11. Phân tích sống sót, thời gian sự kiện và rủi ro cạnh tranh
Bài liên quan
- 10.1. Hồi quy không gian cắt ngang
- 10.1.1. Mô hình trễ không gian của biến giải thích (SLX)
- 10.1.2. Mô hình tự hồi quy trễ không gian của biến phụ thuộc (SAR)
- 10.1.3. Mô hình tự tương quan không gian trong sai số (SEM)
- 10.1.4. Mô hình không gian Durbin và phân rã tác động LeSage–Pace (SDM)
- 10.1.5. Mô hình kết hợp tự hồi quy trễ và sai số không gian (SARAR)
- 10.2. Dữ liệu bảng không gian
- 10.3. Hồi quy địa lý đa thang đo và sai số chuẩn Conley
- 10.4. Kinh tế lượng mạng lưới