Hàm điều khiển (control function) là một phương pháp ước lượng kinh tế lượng và thống kê được sử dụng để khắc phục hiện tượng biến giải thích nội sinh (endogeneity) bằng cách mô hình hóa và kiểm soát thành phần sai số chưa quan sát được thông qua các biến phụ trợ hoặc phần dư ước lượng từ phương trình rút gọn. Thay vì tìm cách thay thế biến giải thích nội sinh bằng giá trị dự báo tuyến tính như trong phương pháp biến công cụ truyền thống, phương pháp hàm điều khiển trực tiếp bổ sung một hàm của phần dư vào phương trình hồi quy chính để làm sạch tương quan giữa biến giải thích và sai số ngẫu nhiên.
Bản chất lý thuyết và cơ chế kiểm soát tính nội sinh
Trong phân tích kinh tế lượng và suy diễn nhân quả, vấn đề biến nội sinh nảy sinh khi một hoặc nhiều biến giải thích có tương quan với sai số ngẫu nhiên của mô hình cấu trúc. Nguyên nhân của hiện tượng này thường bắt nguồn từ biến bị bỏ sót chưa quan sát được, sai số đo lường trong dữ liệu, hoặc tính đồng thời giữa nguyên nhân và kết quả. Nếu ước lượng bằng phương pháp bình phương tối thiểu thông thường (Ordinary Least Squares, OLS), các hệ số hồi quy thu được sẽ bị chệch và không nhất quán, dẫn đến các kết luận sai lệch về mặt chính sách và thực nghiệm.
Phương pháp tiếp cận biến công cụ truyền thống, điển hình là ước lượng bình phương tối thiểu hai bước (Two-Stage Least Squares, 2SLS), giải quyết vấn đề này bằng cách phân rã biến nội sinh thành hai thành phần: một thành phần ngoại sinh được dự báo từ các biến công cụ và một thành phần nội sinh nằm trong phần dư. Trong mô hình tuyến tính chuẩn, ước lượng 2SLS sử dụng giá trị khớp của biến nội sinh từ bước thứ nhất để thay thế cho chính biến đó ở bước thứ hai.
Ngược lại, phương pháp hàm điều khiển tiếp cận bài toán từ một góc nhìn đối ngẫu nhưng mang lại sự linh hoạt vượt trội trong các cấu trúc phi tuyến. Thay vì loại bỏ phần nội sinh để chỉ giữ lại phần dự báo ngoại sinh, phương pháp này giữ nguyên biến giải thích ban đầu trong phương trình cấu trúc, đồng thời đưa thêm một hàm toán học của phần dư ước lượng từ phương trình rút gọn vào mô hình. Hàm toán học bổ sung này được gọi là hàm điều khiển. Vai trò cốt lõi của hàm điều khiển là điều kiện hóa thành phần sai số cấu trúc theo phần dư của phương trình rút gọn, từ đó triệt tiêu hoàn toàn sự tương quan giữa biến giải thích nội sinh và phần sai số còn lại.
Trong bối cảnh mô hình hồi quy tuyến tính chuẩn, sự tương đương về mặt số học giữa ước lượng hàm điều khiển và ước lượng bình phương tối thiểu hai bước 2SLS là một đồng nhất thức đại số OLS thuần tuý xuất phát từ định lý Frisch-Waugh-Lovell. Tính chất trùng khớp này nghiệm đúng trên mọi mẫu dữ liệu thực nghiệm mà không đòi hỏi bất kỳ giả định nào về dạng phân phối của sai số ngẫu nhiên. Tuy nhiên, sức mạnh thực sự và giá trị ứng dụng độc đáo của phương pháp hàm điều khiển nằm ở khả năng mở rộng tự nhiên sang các lớp mô hình phi tuyến phức tạp, mô hình lựa chọn rời rạc, mô hình dữ liệu bảng có hiệu ứng cố định không quan sát được và các mô hình kinh tế lượng phi tham số.
Mô hình chọn mẫu hai bước Heckman và tỷ số nghịch đảo Mills
Nguồn gốc lịch sử và ứng dụng kinh điển nhất của tư duy hàm điều khiển bắt nguồn từ công trình mang tính bước ngoặt của nhà kinh tế học đoạt giải Nobel James Heckman vào năm 1979 công bố trên tạp chí Econometrica. Trong nghiên cứu này, Heckman (1979) đã chứng minh rằng hiện tượng sai lệch chọn mẫu (sample selection bias) thực chất có thể được nhìn nhận và xử lý như một dạng sai số định dạng mô hình do bỏ sót biến giải thích.
Xét bài toán ước lượng phương trình tiền lương của người lao động, trong đó tiền lương chỉ được quan sát đối với những cá nhân thực sự tham gia vào lực lượng lao động. Quyết định tham gia thị trường lao động là một quá trình tự lựa chọn tối ưu của cá nhân, dẫn đến việc nhóm người đi làm có thể sở hữu những đặc điểm năng lực hoặc động lực chưa quan sát được vượt trội hơn so với nhóm không đi làm. Nếu chỉ hồi quy OLS đơn thuần trên mẫu những người có việc làm, sai số của phương trình tiền lương sẽ không có kỳ vọng bằng không do điều kiện chọn mẫu.
Heckman (1979) đã đề xuất quy trình ước lượng hai bước mẫu mực đặt nền móng cho phương pháp hàm điều khiển:
- Bước thứ nhất: Ước lượng một mô hình probit về xác suất tham gia thị trường lao động của toàn bộ mẫu quan sát dựa trên các biến ngoại sinh và ít nhất một biến công cụ ngoại trừ ảnh hưởng đến quyết định đi làm nhưng không tác động trực tiếp lên mức lương. Từ các hệ số ước lượng probit, tính toán giá trị tỷ số nghịch đảo Mills (inverse Mills ratio) cho từng quan sát.
- Bước thứ hai: Đưa tỷ số nghịch đảo Mills đã tính toán vào phương trình hồi quy tiền lương như một biến giải thích bổ sung và tiến hành ước lượng bằng phương pháp OLS trên mẫu quan sát được chọn.
Trong cấu trúc này, tỷ số nghịch đảo Mills chính là một hàm điều khiển tham số cụ thể. Bằng cách đưa hàm này vào phương trình hồi quy, kỳ vọng có điều kiện của sai số cấu trúc được kiểm soát hoàn toàn, biến phần sai số còn lại thành một đại lượng có kỳ vọng có điều kiện bằng không nhưng mang đặc tính phương sai thay đổi (heteroskedasticity) phụ thuộc vào các biến giải thích. Hệ số hồi quy đi kèm tỷ số nghịch đảo Mills cung cấp một kiểm định thống kê trực tiếp về sự hiện diện của sai lệch chọn mẫu: nếu hệ số này khác không có ý nghĩa thống kê, giả thuyết về tính ngoại sinh của quá trình chọn mẫu bị bác bỏ.
Mô hình phi tuyến và ước lượng probit đồng thời Rivers-Vuong
Khi biến phụ thuộc là biến nhị phân hoặc biến định tính rời rạc, việc áp dụng phương pháp biến công cụ tuyến tính thông thường (như 2SLS tuyến tính cho mô hình xác suất tuyến tính) thường bộc lộ nhiều khiếm khuyết nghiêm trọng, chẳng hạn như dự báo xác suất nằm ngoài khoảng không một hoặc vi phạm giả định phương sai sai số đồng nhất. Trong bối cảnh này, phương pháp hàm điều khiển cung cấp một giải pháp ước lượng thanh lịch và nhất quán.
Nghiên cứu kinh điển của Rivers và Vuong (1988) công bố trên Journal of Econometrics đã phát triển phương pháp ước lượng thông tin giới hạn hai bước (hai bước hợp lý cực đại có điều kiện) cho mô hình probit chứa biến giải thích nội sinh liên tục. Quy trình ước lượng Rivers-Vuong vận hành theo các bước chặt chẽ:
- Bước thứ nhất (Phương trình rút gọn): Hồi quy tuyến tính biến giải thích nội sinh trên toàn bộ tập hợp các biến ngoại sinh và các biến công cụ ngoại trừ bằng phương pháp OLS. Từ kết quả hồi quy, trích xuất chuỗi phần dư ước lượng đại diện cho thành phần biến thiên chưa được giải thích của biến nội sinh.
- Bước thứ hai (Phương trình cấu trúc probit): Ước lượng mô hình probit tiêu chuẩn của biến phụ thuộc nhị phân với tập biến giải thích bao gồm các biến độc lập ban đầu, biến nội sinh liên tục và phần dư ước lượng được trích xuất từ bước thứ nhất.
Theo chứng minh của Rivers và Vuong (1988), phần dư đưa vào phương trình probit ở bước thứ hai đóng vai trò là một hàm điều khiển tuyến tính. Khi kiểm soát phần dư này, biến nội sinh ban đầu trở nên ngoại sinh có điều kiện đối với phần sai số ngẫu nhiên còn lại của mô hình probit. Một ưu điểm thực nghiệm vô cùng to lớn của tiếp cận Rivers và Vuong (1988) là kiểm định tính ngoại sinh trở nên cực kỳ đơn giản: chỉ cần thực hiện kiểm định giả thuyết thống kê chuẩn (chẳng hạn kiểm định t hoặc kiểm định tỷ số khả dĩ) cho hệ số của phần dư trong phương trình probit bước thứ hai. Nếu hệ số này không có ý nghĩa thống kê, nhà nghiên cứu có thể an tâm kết luận về tính ngoại sinh mà không cần lo ngại về sai lệch nội sinh.
Mở rộng phi tham số và mô hình tam giác không cộng tính
Trong các ứng dụng kinh tế lượng hiện đại, các giả định phân phối tham số ngặt nghèo (như giả định sai số phân phối chuẩn đa biến) thường bị hoài nghi vì nếu phân phối thực tế bị sai lệch định dạng, các ước lượng tham số của hàm điều khiển có thể mất đi tính nhất quán. Để vượt qua rào cản này, lý thuyết kinh tế lượng vi mô đã chứng kiến bước phát triển vượt bậc hướng tới các phương pháp tiếp cận bán tham số và phi tham số hoàn toàn.
Đóng góp nền tảng của Imbens và Newey (2009) trên tạp chí Econometrica đã mở rộng phương pháp hàm điều khiển cho lớp mô hình phương trình đồng thời dạng tam giác không có tính cộng tính (nonseparable triangular simultaneous equations models). Trong mô hình này, biến phụ thuộc và biến nội sinh liên hệ với nhau thông qua một hàm cấu trúc tổng quát phi tuyến và không thể tách rời sai số ngẫu nhiên theo phép cộng đơn thuần.
Imbens và Newey (2009) đã chứng minh rằng một hàm điều khiển phi tham số hợp lệ có thể được xây dựng trực tiếp từ hàm phân phối tích lũy có điều kiện (conditional cumulative distribution function) của biến giải thích nội sinh khi biết giá trị của các biến công cụ. Cụ thể, biến điều khiển được định nghĩa là giá trị xác suất tích lũy mà biến nội sinh đạt được trong phân phối có điều kiện của nó ứng với các biến công cụ. Khi điều kiện hóa theo biến điều khiển này, tính độc lập có điều kiện giữa biến nội sinh và sai số cấu trúc được thiết lập đầy đủ.
Định lý nhận dạng phi tham số của Imbens và Newey (2009) đòi hỏi hai điều kiện quan trọng: phương trình dạng rút gọn ở bước thứ nhất phải đơn điệu ngặt theo một sai số vô hướng duy nhất và biến giải thích nội sinh phải là một biến liên tục. Khi các điều kiện này được thỏa mãn, phương pháp hàm điều khiển phi tham số cho phép nhận dạng và ước lượng nhất quán các đại lượng kinh tế lượng quan trọng như hàm cấu trúc trung bình (average structural function) và hàm cấu trúc phân vị (quantile structural function) mà không cần áp đặt bất kỳ giả định tham số nào lên hàm hồi quy hay phân phối sai số.
Hồi quy phân vị và ứng dụng kinh tế lượng thực nghiệm
Bên cạnh việc ước lượng tác động trung bình, các nhà kinh tế học thực nghiệm ngày càng quan tâm đến tác động phân phối không đồng nhất của các can thiệp chính sách lên các phân vị khác nhau của biến kết quả. Nghiên cứu của nhà kinh tế lượng Sokbae Lee vào năm 2007 trên Journal of Econometrics đã tiên phong thiết lập phương pháp tiếp cận hàm điều khiển cho các mô hình hồi quy phân vị có biến giải thích nội sinh.
Theo Lee (2007), việc xử lý biến nội sinh trong mô hình phân vị thông qua biến công cụ truyền thống gặp phải thách thức toán học lớn do tính chất phi tuyến của toán tử phân vị. Bằng cách ước lượng phần dư từ phương trình dạng rút gọn ở bước thứ nhất và đưa biến điều khiển này vào phương trình hồi quy phân vị ở bước thứ hai theo dạng phi tham số chuỗi, Lee (2007) đã chứng minh rằng phương pháp hàm điều khiển cho phép ước lượng nhất quán hiệu ứng can thiệp phân vị (quantile treatment effects) trên toàn bộ dải phân vị của tổng thể phân tích mà không cần áp đặt ràng buộc dạng hàm lên tương quan giữa các sai số.
Trong bài báo tổng quan toàn diện công bố trên Journal of Human Resources, Wooldridge (2015) đã hệ thống hóa các phương pháp hàm điều khiển trong kinh tế lượng ứng dụng hiện đại. Wooldridge (2015) nhấn mạnh rằng phương pháp hàm điều khiển đặc biệt vượt trội trong các tình huống sau:
- Mô hình biến phụ thuộc dạng tỷ lệ (fractional response models) khi có biến giải thích nội sinh.
- Mô hình đếm (count data models) như Poisson hay nhị thức âm có biến nội sinh liên tục.
- Mô hình có hệ số ngẫu nhiên và tính không đồng nhất chưa quan sát tương tác với biến nội sinh.
- Các kiểm định tính nội sinh mạnh mẽ, có khả năng kháng lại hiện tượng phương sai thay đổi mà không đòi hỏi các thuật toán tối ưu hóa phi tuyến phức tạp.
So sánh kỹ thuật giữa biến công cụ chuẩn và hàm điều khiển
Mặc dù cùng chia sẻ mục tiêu triệt tiêu chệch do biến nội sinh gây ra thông qua việc sử dụng các biến công cụ ngoại trừ, phương pháp biến công cụ truyền thống (IV/2SLS) và phương pháp hàm điều khiển (CF) sở hữu những đặc tính kỹ thuật và phạm vi ứng dụng rất khác biệt.
Bảng dưới đây tổng hợp các điểm tương đồng và dị biệt cốt lõi giữa hai phương pháp tiếp cận kinh tế lượng này:
| Tiêu chí so sánh | Biến công cụ truyền thống (IV / 2SLS) | Phương pháp hàm điều khiển (Control Function) |
|---|---|---|
| Cơ chế xử lý biến nội sinh | Thay thế biến nội sinh bằng giá trị dự báo từ phương trình cấp một | Giữ nguyên biến nội sinh và bổ sung hàm của phần dư phương trình cấp một |
| Hành vi trong mô hình tuyến tính | Cho ước lượng hệ số góc nhất quán chuẩn tắc | Cho ước lượng hệ số góc hoàn toàn trùng khớp đại số theo định lý Frisch-Waugh-Lovell |
| Khả năng thích ứng mô hình phi tuyến | Khó mở rộng, dễ rơi vào bẫy ước lượng chệch (cấm kỵ 2SLS ngây thơ) | Rất tự nhiên và hiệu quả cho các mô hình probit, tobit, hồi quy phân vị |
| Yêu cầu giả định phân phối | Chỉ đòi hỏi tính trực giao của biến công cụ và điều kiện hạng | Thường đòi hỏi giả định về cấu trúc sai số hoặc tính đơn điệu vô hướng |
| Kiểm định tính nội sinh | Cần tính toán các kiểm định Hausman phức tạp qua ma trận hiệp phương sai | Kiểm định trực tiếp bằng hệ số t của biến phần dư trong phương trình chính |
| Hiệu chỉnh sai số chuẩn | Được các phần mềm thống kê tự động chuẩn hóa | Thống kê t kiểm định ngoại sinh hợp lệ chuẩn tắc; suy diễn hệ số cấu trúc cần sandwich hoặc bootstrap |
Theo phân tích phương pháp luận chi tiết của Wooldridge (2015), một ưu thế thực nghiệm vượt bậc của phương pháp hàm điều khiển nằm ở kiểm định tính ngoại sinh: dưới giả thuyết vô hiệu biến giải thích là ngoại sinh, kiểm định t thông thường đối với hệ số của phần dư ở bước thứ hai là hoàn toàn chuẩn xác và hợp lệ tiệm cận mà không cần bất kỳ sự hiệu chỉnh sai số chuẩn nào. Tuy nhiên, khi giả thuyết ngoại sinh bị bác bỏ và nhà nghiên cứu muốn suy diễn thống kê trên các hệ số cấu trúc khác, việc phần dư là một biến số được ước lượng sinh ra từ bước thứ nhất (generated regressor) sẽ khiến sai số chuẩn OLS hoặc probit thông thường bị chệch xuống. Trong tình huống đó, việc áp dụng công thức điều chỉnh ma trận thông tin dạng sandwich hoặc sử dụng kỹ thuật lấy mẫu lặp lại bootstrap theo cụm là bắt buộc để đảm bảo khoảng tin cậy và kiểm định thống kê chuẩn xác.
Thách thức của biến giải thích nội sinh rời rạc và tính không đồng nhất chưa quan sát
Một trong những giới hạn phương pháp luận quan trọng nhất của phương pháp hàm điều khiển trong kinh tế lượng hiện đại là việc xử lý các biến giải thích nội sinh có tính chất rời rạc (chẳng hạn như biến nhị phân chỉ thị can thiệp chính sách hoặc biến đếm nguyên). Như Wooldridge (2015) đã chỉ ra, khi biến nội sinh là biến rời rạc, phần dư từ phương trình dạng rút gọn cấp một không thể liên tục và không thể đóng vai trò như một hàm điều khiển hoàn hảo để trực giao hóa hoàn toàn sai số cấu trúc trong các mô hình phi tuyến không cộng tính.
Trong các tình huống này, việc áp dụng máy móc phương pháp hàm điều khiển tiêu chuẩn có thể dẫn đến các ước lượng không nhất quán nếu dạng phân phối thực tế của sai số bị sai lệch định dạng. Để khắc phục hạn chế này, các nhà kinh tế lượng thường áp dụng các giải pháp xấp xỉ tham số linh hoạt hoặc kết hợp phương pháp hàm điều khiển với ước lượng hiệu ứng can thiệp cục bộ (Local Average Treatment Effects, LATE) thông qua biến công cụ truyền thống. Sự kết hợp hài hòa giữa biến công cụ tuyến tính và hàm điều khiển phi tuyến mang lại cho các nhà nghiên cứu một bộ công cụ thực nghiệm đa năng và đáng tin cậy để giải mã các mối quan hệ nhân quả phức tạp trong kinh tế học và khoa học xã hội hiện đại.