Hướng dẫn sử dụng Stata 8.0 UMP

Hướng dẫn sử dụng Stata 8.0 Khởi động Stata 8.0 • Vào Start – Programs – Medistat – Stata • Màn hình Stata bao gồm: – thanh: • Thanh tiêu đề • Thanh menu • Thanh cơng cụ: – cửa sổ: • Lệnh (command); Kết (results); xem lại (review) biến số (variable) Thanh tiêu đề Thanh menu Thanh công cụ Xem lại Kết Biến số Lệnh • Sắp xếp cửa sổ cho tiện việc sử dụng lưu cửa sổ với menu Prefs :: Save Windows Preferences • Nhật kí kết (lưu giữ kết quả): log cần thực trước thực phân tích thống kê với Stata • Quan sát (observation) biến (variable) • Tên biến (name) nhãn biến (label) • Format: string, general date • Biến định tính biến định lượng Bài tập • Nhập số liệu người (3 nam nữ) vào Stata: Giới – Nam – Nam – Nữ – Nữ – Nam Đường huyết 80 90 100 110 120 • Hãy tính tần suất giới tính, giá trị đường huyết trung bình, giá trị đường huyết trung bình phân tầng theo giới • Stata phân biệt chữ hoa chữ thường; cần lưu ý để gõ lệnh Stata (thường chữ thường) • Tên biến (name of variable) – Bắt đầu chữ – Các kí tự chữ cái, chữ số hay gạch (_) Khơng cho phép kí tự đặc biệt dấu trừ (-) hay dấu trắng( ) • Lưu ý: Stata khơng hiển thị font tiếng Việt (unicode, ABC, TCVN3, v.v) • Biến số: đặc tính hay đại lượng đối tượng – Đặc tính: biến số định tính – Đại lượng: biến số định lượng • Thống kê mơ tả: sum (định lượng) tab1 (định tính) • Để mơ tả phân tầng sử dụng lệnh bysort Statistics :: Summaries, tables, & tests :: Tables :: Two-way Tables with measures of association Biến số phụ thuộc theo hàng Biến số độc lập nằm theo cột Tính % hàng: NC bệnh chứng 1=case, | 0=no, 1=yes 0=control | khong co | Total -+ + -chung | 534 474 | 1008 | 52.98 47.02 | 100.00 -+ + -benh | 210 42 | 252 | 83.33 16.67 | 100.00 -+ + -Total | 744 516 | 1260 | 59.05 40.95 | 100.00 Pearson chi2(1) = 76.8300 Pr = 0.000 Statistics :: Observational/Epi Analysis :: Tables for epidemiologists :: Case-control odds ratio Biến số phụ thuộc theo hàng Biến số độc lập nằm theo cột Thực tập phân tích khử yếu tố gây nhiễu Contains data from C:\DATA\CHILUMBA.DTA obs: 1,260 Leprosy and BCG vaccination vars: 29 Jan 1996 17:02 size: 18,900 (99.8% of memory free) storage display value variable name type format label variable label id float %9.0g Identity number caco byte %8.0g 1=case, 0=control agegrp byte %8.0g Age group sex byte %8.0g 0=male, 1=female bcgscar byte %8.0g 0=no, 1=yes school byte %8.0g Duration of schooling mbcont byte %8.0g Contact with MB case pbcont byte %8.0g Contact with PB case - ` cc caco bcgscar Proportion | Exposed Unexposed | Total Exposed -+ + -Cases | 42 210 | 252 0.1667 Controls | 474 534 | 1008 0.4702 -+ + -Total | 516 744 | 1260 0.4095 | | | Point estimate | [95% Conf Interval] | + -Odds ratio | 2253165 | 1543594 3234627 Prev frac ex | 7746835 | 6765373 8456406 Prev frac pop | 3642857 | + chi2(1) = 76.83 Pr>chi2 = 0.0000 Biến số phụ thuộc theo hàng Biến số độc lập nằm theo cột Biến số phụ thuộc theo hàng Biến số độc lập nằm theo cột Biến số phân tầng Sử dụng phân tầng Mantel-Haenszel Trọng số Mantel-Haenzel OR tầng cc caco bcgscar, by( agegrp ) Age group | OR [95% Conf Interval] M-H Weight -+ 1/14 | 468 2004388 1.054927 10.03344 15/24 | 4457382 1821365 1.049273 9.334783 25/34 | 2826511 0942287 8728589 6.333333 35-44 | 1784882 0677976 425671 15.72251 45+ | 095122 0022768 6113332 6.507937 -+ Crude | 2253165 1543594 3234627 M-H combined | 293582 1968806 4377801 Test of homogeneity (M-H) chi2(4) = 5.12 Pr>chi2 = 0.2756 OR tổng hợp Test that combined OR = 1: Mantel-Haenszel chi2(1) = Pr>chi2 = Ảnh hưởng BCG lên bệnh phong có khác nhóm tuổi hay khơng Có ảnh hưởng BCG lên bệnh phong hay không (sau loại trừ bỏ yếu tố gây nhiễu) yếu tố nguy cơ) 38.42 0.0000 (exact) (exact) (exact) (exact) (exact) (exact) logistic caco bcgscar agegrp sex school mbcont pbcont Logistic regression Log likelihood = -529.16522 Number of obs LR chi2(6) Prob > chi2 Pseudo R2 = = = = 1209 126.84 0.0000 0.1070 -caco | Odds Ratio Std Err z P>|z| [95% Conf Interval] -+ -bcgscar | 3383164 0683668 -5.36 0.000 2276736 5027285 agegrp | 1.258598 0668384 4.33 0.000 1.134184 1.396659 sex | 1.029841 1682939 0.18 0.857 7475992 1.418637 school | 6599573 0764139 -3.59 0.000 525967 8280816 mbcont | 2.903683 1.512648 2.05 0.041 1.04598 8.060745 pbcont | 1.795166 3998049 2.63 0.009 1.160197 2.777649 Thực phân tích bảng 2x2 phân tầng Kiểm định tính đồng OR (RR) theo tầng p-value10% No Không nhiễu, không tương tác, báo cáo OR (RR) thô Có tương tác Báo cáo OR (RR) tầng ORMH − ORcrude 0.294 − 0.225 = = 23% ORMH 0.294 Yes Biến phân tầng biến gây nhiễu Báo cáo OR (RR) kết hợp MH Thực bảng x phân tầng Xem kết kiểm định tính đồng chất có ý nghĩa TK? Khơng đồng chất Báo cáo OR nhóm (tương tác) OR nhóm = ORKH ORKH đánh giá tốt mối liên quan ORthô ≈ ORKH Báo cáo ORThô (Không gây nhiễu) Khơng tương đương Báo cáo ORKH (Có gây nhiễu) ...Khởi động Stata 8.0 • Vào Start – Programs – Medistat – Stata • Màn hình Stata bao gồm: – thanh: • Thanh tiêu đề • Thanh menu • Thanh công... tích cho biến số định lượng – Thực lệnh thống kê khác v.v Thốt khỏi chương trình Stata Vào lại chương trình Stata; sử dụng menu File::log::View để xem kết phân tích thực bước Bài tập • • Mở tập... vào Stata: Giới – Nam – Nam – Nữ – Nữ – Nam Đường huyết 80 90 100 110 120 • Hãy tính tần suất giới tính, giá trị đường huyết trung bình, giá trị đường huyết trung bình phân tầng theo giới • Stata