Hồi quy tuyến tính STATA là việc ước lượng mô hình bình phương nhỏ nhất thông thường (OLS) bằng lệnh regress, sau đó kiểm tra các giả định như đa cộng tuyến, phương sai sai số thay đổi, dạng hàm và quan sát ảnh hưởng trước khi diễn giải hệ số. Bỏ qua phần kiểm tra, con số trong bảng có thể đúng về tính toán nhưng sai về suy luận.
Bài viết tổng hợp năm nhóm nội dung mà Luận Văn 360 thường đối chiếu khi rà soát chương kết quả: đặc tả mô hình, ước lượng, kiểm định giả định, điểm ảnh hưởng và cách báo cáo. Mỗi nhóm đều gắn với lệnh cụ thể, cách đọc đầu ra và hướng xử lý khi giả định bị vi phạm.
Tóm tắt nhanh: mô hình phải xuất phát từ lý thuyết và giả thuyết, không phải từ việc thử biến đến khi có ý nghĩa thống kê.
Sau regress, tối thiểu cần xem VIF, kiểm định phương sai thay đổi và kiểm định dạng hàm.
Báo cáo hệ số kèm sai số chuẩn, mức ý nghĩa, cỡ mẫu và loại sai số chuẩn đã dùng.
Hồi quy tuyến tính STATA cần kiểm tra những gì?
Theo các giáo trình kinh tế lượng phổ biến như Wooldridge (Introductory Econometrics), ước lượng OLS cho kết quả không chệch và hiệu quả khi thỏa các giả định Gauss–Markov. Trên thực tế, giả định hiếm khi thỏa trọn vẹn, nên việc kiểm tra nhằm biết mức độ vi phạm và chọn cách khắc phục phù hợp. Bảng dưới đây tóm tắt các kiểm tra hay dùng nhất.
| Vấn đề | Lệnh STATA sau regress | Giả thuyết H0 hoặc cách đọc | Hướng xử lý thường gặp |
|---|---|---|---|
| Đa cộng tuyến | estat vif | VIF cao cho thấy biến trùng thông tin | Xem lại biến, gộp hoặc bỏ có lý do |
| Phương sai sai số thay đổi | estat hettest, estat imtest | H0: phương sai sai số không đổi | Dùng sai số chuẩn vững |
| Sai dạng hàm | estat ovtest, linktest | H0: mô hình không bỏ sót dạng phi tuyến | Thêm bình phương, lấy log, tương tác |
| Phần dư không chuẩn | qnorm, swilk | H0: phần dư phân phối chuẩn | Xem cỡ mẫu, biến đổi biến phụ thuộc |
| Quan sát ảnh hưởng | predict kèm cooksd, lvr2plot | Điểm có đòn bẩy và phần dư lớn | Kiểm tra nguồn, phân tích độ nhạy |
Xác định biến và đặc tả mô hình
Trước khi gõ lệnh, hãy viết phương trình hồi quy ra giấy. Mỗi biến cần một lý do từ lý thuyết hoặc nghiên cứu trước, giúp bạn tránh thêm bớt biến tùy tiện sau khi nhìn thấy giá trị p. Bản đặc tả nên trả lời được các ý sau:
- Biến phụ thuộc đo lường khái niệm gì, theo đơn vị nào.
- Biến giải thích chính gắn với giả thuyết nào và dấu kỳ vọng ra sao.
- Biến kiểm soát nào được nghiên cứu trước sử dụng phổ biến.
- Mẫu ước lượng gồm những quan sát nào, giai đoạn nào.
Biến kiểm soát và đơn vị đo
Biến kiểm soát giúp tách ảnh hưởng của biến chính khỏi các yếu tố gây nhiễu, ví dụ quy mô và tuổi doanh nghiệp khi nghiên cứu tác động của đòn bẩy đến hiệu quả. Đơn vị đo ảnh hưởng trực tiếp đến cách đọc hệ số. Biến tiền tệ có độ lệch lớn thường được lấy logarit để hệ số mang nghĩa phần trăm gần đúng.
Biến phân loại và tương tác
STATA hỗ trợ biến nhân tố (factor variables), nên bạn không cần tự tạo từng biến giả. Tiền tố i. trước biến ngành sẽ tự sinh các biến giả và bỏ một nhóm làm gốc. Ký hiệu ## tạo cả hiệu ứng chính lẫn tương tác, còn c. đánh dấu biến liên tục khi tham gia tương tác. Cách này cũng giúp lệnh margins tính hiệu ứng biên chính xác.
Ước lượng bằng lệnh regress
Cú pháp cơ bản là regress theo sau bởi biến phụ thuộc rồi danh sách biến độc lập. Nên lưu mỗi mô hình bằng estimates store để so sánh và xuất bảng sau này. Nếu ước lượng nhiều đặc tả, hãy chạy chúng trên cùng một mẫu để kết quả so sánh được với nhau.
Đọc bảng kết quả regress
Phần đầu bảng cho số quan sát, thống kê F và giá trị p của nó, R bình phương và R bình phương hiệu chỉnh. Phần dưới liệt kê hệ số, sai số chuẩn, thống kê t, giá trị p và khoảng tin cậy 95%. Hãy đọc dấu và độ lớn hệ số trước, rồi mới xem mức ý nghĩa, vì một hệ số có ý nghĩa thống kê vẫn có thể rất nhỏ về mặt kinh tế.
Khi nào dùng sai số chuẩn vững?
Tùy chọn vce(robust) cho sai số chuẩn vững với phương sai thay đổi mà không làm thay đổi hệ số ước lượng. Khi quan sát được nhóm theo đơn vị, chẳng hạn nhiều nhân viên trong cùng một công ty, vce(cluster) theo mã nhóm sẽ phù hợp hơn. Hãy ghi rõ loại sai số chuẩn dưới bảng kết quả.
Kiểm tra giả định sau ước lượng
Các lệnh postestimation phải chạy ngay sau regress vì chúng dùng kết quả đang lưu trong bộ nhớ. Một số kiểm định như estat hettest được xây dựng cho mô hình ước lượng với sai số chuẩn thông thường, nên hãy chạy chúng trước khi chuyển sang phiên bản dùng sai số chuẩn vững.
Đa cộng tuyến và phương sai thay đổi
Lệnh estat vif in hệ số phóng đại phương sai cho từng biến. Hair và cộng sự (Multivariate Data Analysis) thường dẫn ngưỡng tham khảo VIF bằng 10, một số tác giả dùng ngưỡng chặt hơn là 5; hãy đối chiếu yêu cầu của giảng viên. Với phương sai thay đổi, giá trị p của estat hettest nhỏ hơn 0,05 cho thấy nên bác bỏ giả thuyết phương sai không đổi.
Dạng hàm và phân phối phần dư
Kiểm định Ramsey RESET qua estat ovtest gợi ý mô hình có thể bỏ sót dạng phi tuyến khi bác bỏ H0. Để xem phần dư, dùng predict với tùy chọn residuals, rồi vẽ rvfplot hoặc qnorm. Với mẫu lớn, sai lệch nhẹ khỏi phân phối chuẩn thường ít ảnh hưởng đến suy luận, nhưng vẫn nên mô tả trong luận văn.
Khi dữ liệu sắp xếp theo thời gian
Nếu quan sát là chuỗi thời gian, cần khai báo bằng tsset trước, sau đó dùng estat dwatson hoặc estat bgodfrey để kiểm tra tự tương quan. Với dữ liệu theo đơn vị và năm, kiểm định này không còn phù hợp; hãy chuyển sang công cụ của dữ liệu bảng trong STATA.
Điểm ảnh hưởng và quan sát bất thường
Một vài quan sát có thể kéo đường hồi quy về phía chúng. Lệnh predict với tùy chọn cooksd lưu khoảng cách Cook cho mỗi dòng; một quy tắc kinh nghiệm hay dùng là xem kỹ các điểm lớn hơn 4 chia cho cỡ mẫu. Biểu đồ lvr2plot đặt đòn bẩy cạnh phần dư bình phương, giúp nhận ra điểm vừa xa vừa lệch.
Kinh nghiệm từ Luận Văn 360: nhiều bản thảo chúng tôi rà soát loại bỏ điểm ảnh hưởng chỉ vì sau khi loại, hệ số chính trở nên có ý nghĩa.
Cách an toàn hơn là kiểm tra nguồn gốc của điểm đó, rồi báo cáo kết quả hồi quy tuyến tính STATA cả khi giữ và khi loại. Nếu kết luận đổi chiều, đó là phát hiện cần thảo luận, không phải điều nên che đi.
Nếu bạn cần một lượt rà soát độc lập cho đặc tả mô hình và các kiểm định, dịch vụ hỗ trợ phân tích dữ liệu bằng STATA của Luận Văn 360 sẽ góp ý cụ thể để bạn tự hoàn thiện. Mọi quyết định cuối cùng về mô hình vẫn thuộc về bạn và giảng viên hướng dẫn.
Báo cáo kết quả hồi quy tuyến tính STATA
Bảng kết quả nên đặt các mô hình cạnh nhau theo cột, ví dụ mô hình cơ sở rồi mô hình có thêm biến kiểm soát. Lệnh etable từ STATA 17 hoặc lệnh cộng đồng esttab giúp xuất bảng thống nhất. Nội dung tối thiểu gồm:
- Hệ số và sai số chuẩn hoặc thống kê t đặt trong ngoặc, kèm ký hiệu mức ý nghĩa.
- Số quan sát, R bình phương hiệu chỉnh và thống kê F của từng mô hình.
- Ghi chú về loại sai số chuẩn và các biến giả đã đưa vào.
- Kết quả tóm tắt các kiểm định giả định, có thể đặt ở bảng riêng hoặc phụ lục.
Khi diễn giải, hãy nói bằng đơn vị thật. Ví dụ minh họa: “Giữ các yếu tố khác không đổi, đòn bẩy tăng 1 điểm phần trăm gắn với ROA giảm khoảng 0,05 điểm phần trăm.” Phần giải thích ý nghĩa kiểm định có thể tham khảo thêm sổ tay thống kê của NIST.
Lỗi thường gặp
Diễn giải hệ số như quan hệ nhân quả
Hồi quy OLS trên dữ liệu quan sát cho biết mối liên hệ có điều kiện, không tự chứng minh nhân quả. Nếu nghi ngờ biến nội sinh, hãy nêu trong phần hạn chế hoặc cân nhắc phương pháp như biến công cụ với lệnh ivregress, sau khi trao đổi với người hướng dẫn.
Chạy kiểm định nhưng không xử lý kết quả
Không ít bản thảo báo cáo estat hettest có ý nghĩa nhưng vẫn trình bày sai số chuẩn thông thường. Mỗi kiểm định cần đi kèm một quyết định: giữ nguyên vì không vi phạm, hoặc điều chỉnh và nêu cách điều chỉnh. Hội đồng thường hỏi đúng điểm này.
Kết luận
Hồi quy tuyến tính STATA đáng tin cậy khi đi đủ năm nhóm việc: đặc tả mô hình từ lý thuyết, ước lượng bằng regress với loại sai số chuẩn phù hợp, kiểm tra giả định, soát điểm ảnh hưởng và báo cáo minh bạch. Hãy ghi toàn bộ lệnh vào do-file để kết quả có thể chạy lại.
Trước khi ước lượng, nên hoàn tất bảng thống kê mô tả trong STATA trên cùng mẫu. Nếu bạn dùng R song song, bài hồi quy tuyến tính trong R giúp đối chiếu kết quả. Để được góp ý trên chính tệp của mình, hãy tìm hiểu dịch vụ phân tích STATA của Luận Văn 360.
