Thực hành Stata - Phần 1: Làm quen với Stata

PXTThực hành Stata*Thực hành StataPhần 1: Làm quen với StataPXTThực hành Stata*LÀM QUEN VỚI STATA 1. Cửa sổ kết quả (ctrl 1) Stata results: Hiển thị kết quả các câu lệnh vừa thực hiện. PXTThực hành Stata*LÀM QUEN VỚI STATA2. Cửa sổ viết câu lệnh (ctrl 4) Stata Command: Cho phép người thực hành viết câu lệnh thực hành stata.PXTThực hành Stata*3. Cửa sổ hiển thị câu lệnh đã thực hiện (ctrl 5) Stata Review: hiển thị các câu lệnh mà người thực hành đã thực hiện.LÀM QUEN VỚI STATAPXTThực hành Stata*4. Cửa sổ hiển thị tên các biến (ctrl 6) Stata variables: hiển thị tên các biến.LÀM QUEN VỚI STATAPXTThực hành Stata* Lệnh đặt bộ nhớ: Set mem 3200 [m/k]TẠO BIẾN, LƯU NHỚ VÀ LẬP DATA Lệnh nhập dữ liệu: Input [Tên biến1] [Tên biến2] [Tên biến3] [Giá trị biến1] [giá trị biến2] [giá trị biến3] endPXTThực hành Stata*SỬA CHỮA DỮ LIỆU Sửa chữa một giá trị quan sát nào đó của biến bị nhầm lẫn: Replace [tên biến có giá trị quan sát bị sai]=[giá trị đúng] in [STT quan sát mà biến bị sai].(Lệnh này chỉ sửa chữa được một giá trị của biến.)Ví dụ: replace X2=170 in 1Lệnh này sẽ thay thế giá trị cũ của biến X2 tại quansát số 1 bằng giá trị mới (170).PXTThực hành Stata*SỬA CHỮA DỮ LIỆU Lệnh đổi tên biến:Rename [tên biến cũ] [tên biến mới]Rename X1 Hoa Lệnh hiển thị thông tin về dữ liệu đã sửa chữa:List SỬA CHỮA DỮ LIỆUSửa chữa dữ liệu trên cửa sổ editor: Lệnh edit: editChú ý: sau khi sửachữa dữ liệu phải lưu dữ liệu.PXTThực hành Stata*PXTThực hành Stata*SỬA CHỮA DỮ LIỆU Lệnh xóa biến, xóa quan sát Drop [tên các biến cần xóa]Drop in [số thứ tự của quan sát cần xóa]PXTThực hành Stata*Drop điều kiệnVí dụ: drop if Y>48 Lệnh trên sẽ xóa các quan sát mà tại đó biến Y có giá trị lớn hơn 48.SỬA CHỮA DỮ LIỆUPXTThực hành Stata*HỒI QUYCấu trúc câu lệnh:Reg [tên biến phụ thuộc] [danh sách các biến độc lập]Ví dụ: reg Y X1 X2PXTThực hành Stata*Hồi quy có điều kiện:Cấu trúc câu lệnh:Reg [biến phụ thuộc] [danh sách biến độc lập] [in/if phạm vi quan sát/biểu thức điều kiện]Ví dụ:Reg Y X1 X2 in 1/10 HỒI QUYPXTThực hành Stata*Reg Y X1 X2 if Y>150Hồi quy với độ tin cậy tùy chọn:Reg Y X1 X2, level(độ tin cậy)Ví dụ:Reg Y X1 X2, level(98) HỒI QUYPXTThực hành Stata*Predict ymu, xbPredict ei, resid Hoặc có thể lấy giá trị của ei bằng lệnh tạo biến:gen ei=Y-ymuHỒI QUYLấy giá trị của Y^, ei:PXTThực hành Stata*○Ví dụ: gen D=1 if Y>150 HỒI QUY VỚI BIẾN GIẢLệnh tạo biến giả (hai câu lệnh)Gen D=1 if biểu thức điều kiện○Replace D=0 if Y150)PXTThực hành Stata* HỒI QUY VỚI BIẾN GIẢLệnh hồi quy với giả:Reg [biến phụ thuộc] [các biến số lượng biến giả] Ví dụ: reg Y X1 X2 D Chú ý: Hồi quy với biến giả cũng có thể thực hiện lệnh hồi quy có điều kiện (in/if, level) PXTThực hành Stata* Lệnh hồi quy với biến giả và biến tương tác: Reg Y X1 X2 D DX2 Chú ý: Hồi quy với biến giả cũng có thể thực hiện lệnh hồi quy có điều kiện (in/if, level)HỒI QUY VỚI BIẾN GIẢ Tạo biến tương tác: Gen DX2=D*X2PXTThực hành Stata*ĐA CỘNG TUYẾN1.Dựa vào hàm hồi qui phụ:Phát hiện dựa vào R2 cao và tỷ số t thấp;Thực hiện hồi quy phụ;Kiểm định sự phù hợp của hàm hồi quy phụ;Khắc phục đơn giản nhất là bỏ biến phụ thuộc của HHQ phụ.ĐA CỘNG TUYẾN2. Dựa vào nhân tử phóng đại phương saiVif + Nếu vif ≥10 thìHHQ có ĐCT; + Lệnh vif phải thực hiện sau lệnh HQ và Có tác dụng đối vớilệnh HQ gần nhất.PXTThực hành Stata*KIỂM ĐỊNH GIẢ THIẾT1. Kiểm định giá trị trung bình của biếnTtest [tên biến]=[giá trị so sánh, level()]Ttest [tên biến1]=[tên biến2, level()]PXTThực hành Stata*KIỂM ĐỊNH GIẢ THIẾTPXTThực hành Stata*Ví dụ: 1. ttest Y=145, level(90) 2. ttest x1=x2, level(90)Chú ý: khi không chỉ ra độ tin cậy thì Stata sẽ mặc định độ tin cậy là 95%.KIỂM ĐỊNH GIẢ THIẾTPXTThực hành Stata*2. Kiểm định hệ số test [tến biến]=[giá trị so sánh] test [tến biến1]=[tên biến2]Chú ý: kiểm định hệ số của biến phải thực hiện sau lệnh hồi quy và có tác dụng với lệnh hồi quy gần nhất. KIỂM ĐỊNH GIẢ THIẾTPXTThực hành Stata*Ví dụ 1: test X2=0.5(Kiểm định hệ số của biến X2=0.5) Ví dụ 2: test X2=X3(Kiểm định hs của biến X2=hs của biến X3) PHƯƠNG SAI SAI SỐ THAY ĐỔIHettest [danh sách các biến độc lập]Ví dụ: hettest X2 X3 DBpagan [danh sách các biến độc lập]Ví dụ: bpagan X2 X3 DwhitetstPXTThực hành Stata*TỰ TƯƠNG QUANBước 1: tạo biến số thứ tự: gen stt=_nBước 2: tạo biến thời gian: tsset sttBước 3: kiểm định DW: dwstatPXTThực hành Stata*KIỂM ĐỊNH VỀ DẠNG HÀMKiểm định thừa biến (=kiểm định hệ số của biến độc lập bằng 0)Kiểm định thiếu biến (kiểm định Ramsey): ovtestPXTThực hành Stata*