4/22/2007 CH NG VII STATA TRONG PHÂN TÍCH D CH T H C C B N Ths Bùi V n Tr ng Gi i thi u STATA m t cơng c tính tốn phân tích r t có hi u qu phân tích d ch t h c Trong khuôn kh ch ng ch gi i thi u s d ng STATA cho m t s phân tích c b n th ng dùng nghiên c u d ch t h c ti n hành nghiên c u v đ c m t nhiên c a b nh nh t l m c, t vong ho c tàn t t c a m t lo i b nh ho c hi n t ng s c kh e đó, nhà d ch t h c ti n hành nhi u lo i hình nghiên c u hay g p nghiên c u thu n t p nghiên c u b nh ch ng V i yêu c u đ n gi n c a nghiên c u này, nhà nghiên c u có th t o b ng c b n th c hi n m t s tính tốn liên quan Các b ng có th b ng ti p liên 2x2, b ng 2xK, ho c b ng phân t ng Khi phân tích có k t qu đ u “có – khơng”, STATA có th t o đ u b ng 2x2 Khi phân tích có k t qu th h ng, STATA có th t o b ng x K đ cho m c đích phân tích lo i Khi phân tích địi h i so sánh b ng, STATA có th t o b ng phân t ng Trong ch ng s gi i thi u m t s b ng c b n đ h c dùng đ phân tích đ c trình bày c nhà d ch t ti n cho vi c theo dõi, chúng tơi s trình bày l nh STATA theo lo i hình nghiên c u: nghiên nghiên c u thu n t p, b nh ch ng, b nh ch ng có ghép c p Sau chúng tơi s trình bày dùng l nh STATA phân tích phân t ng đ ki m sốt y u t nhi u (confounding factors) ho c xác đ nh y u t t ng tác v i y u t nguy c (interaction) ho c y u t làm thay đ i m i quan h ph i nhi m – b nh (effect modiffiction): g i chung y u t tác đ ng © Bùi V n Tr ng 4/22/2007 Phân tích nghiên c u thu n t p (cohort studies) Nghiên c u thu n t p (cohort studies) m t lo i nghiên c u quan sát, m t hay nhi u nhóm kh i đ u khơng có b nh, chia thành nhóm nh d a c s có ph i nhi m hay khơng ph i nhi m v i y u t nguy c Nghiên c u thu n t p có th nghiên c u thu n t p t ng lai (Prospective cohort studies) ho c nghiên c u thu n t p h i c u (Restrospective cohort studies) Các nhóm đ đ xác đ nh tr c nghiên c u đ ng h p m c m i c theo dõi kho ng th i gian c hai nhóm ph i nhi m không ph i nhi m v i y u t nguy c Phân tích phiên gi i k t qu c a nghiên c u thu n t p liên quan đ n vi c tính tốn t l m c m i m cm i đ nhóm ph i nhi m so sánh v i t l nhóm khơng ph i nhi m v i y u t nguy c Các s li u nghiên c u c trình bày b ng b ng 2x2 | Y u t | Có ph i nhi m Không | | T ng -+ + -B nh | a c | a+c Không b nh | b d | b+d -+ + -T ng Qua b ng, ta có th tính đ | a+b c+d | a+b+c+d c: - T l m c m i (t xu t nguy c ) nhóm ph i nhi m: Ie = a/(a+b) - T l m c m i (t xu t nguy c ) nhóm khơng ph i nhi m:Iue = c/(c+d) - Chênh l ch v nguy c (Risk Difference) RD = Ie - Iue= [a/(a+b) - c/(c+d)] - Nguy c t ng đ i RR (Risk Ratio, Relative Risk)= Ie / Iue = [a/(a+b)]/[c/(c+d)] ng đ i - Kho ng tin c y c a nguy c t - T xu t chênh l ch v b nh nhóm ph i nhi m nhóm khơng ph i nhi m (OR) = (a/b)/(c/d) = (a.d)/(b.c) © Bùi V n Tr ng 4/22/2007 - T l nguy c thu c (Attributable Risk Proportion) ARP = (RR-1)/RR T l m c m i c a nhóm ph i nhi m v i y u t nguy c s đ c so sánh v i t l m c m i c a nhóm khơng ph i nhi m v i y u t nguy c T s c a s so sánh này đ ng đ i (Relative c g i t s nguy c hay nguy c t Risk, RR) thông s d ch t h c đ c quan tr ng nh t nghiên c u thu n t p Dùng STATA cho phép có th tính tốn đ c t t c thông s v i kho ng tin c y c a chúng Ngoài b ng ti p liên đ so sánh m i liên quan gi a hai bi n phân lo i nên STATA cho phép tính tốn Chi bình ph ng giá tr p 2.1 Tính tốn tr c ti p Trong tr ng h p tính tốn tr c ti p ta dùng cú pháp l nh nh sau: csi #a #b #c #d [, csi_options] Gi i thích cú pháp l nh: csi: t vi t t t c a c m t cohort study immediate #a #b #c #d: s tr csi_options: có th csi_options ng h p đ m i ô b ng 2x2 c dùng mô t sau: Mô t or a k t qu t xu t chênh (OR) tb Tính tốn ki m nghi m d a kho ng tin c y exact Tính tốn p c a Fisher's exact level(#) t m c tin c y; m c đ nh m c đ tin c y 95 Ví d 1: Trong m t nghiên c u thu n t p theo dõi nguy c c a hút thu c đ i v i b nh m ch vành, ngu i ta ti n hành nghiên c u hai nhóm kh e m nh t 30 tu i tr lên: nhóm hút thu c khơng hút thu c K t qu sau m t 10 n m theo dõi © Bùi V n Tr ng 4/22/2007 cho th y s 4893 ng i hút thu c có 657 ng i xu t hi n b nh m ch vành 4236 không m c b nh nhóm khơng hút thu c 3480 ng i có 367 tr ng h p xu t hi n b nh m ch vành, 3480 không m c Dùng STATA đ tính tr c ti p cho k t qu nh sau: csi 657 367 4236 3480 | Exposed Unexposed | Total -+ + -Cases | 657 367 | 1024 Noncases | 4236 3480 | 7716 -+ + -Total | 4893 3847 | 8740 | | Risk | 1342735 095399 | 1171625 | | | Point estimate | [95% Conf Interval] | + -Risk difference | 0388744 | 0255539 052195 Risk ratio | 1.407493 | 1.247657 1.587805 Attr frac ex | 2895169 | 1984979 3701999 Attr frac pop | 1857545 | + chi2(1) = 31.47 Pr>chi2 = 0.0000 T b ng k t qu ta có : - Nguy c m c b nh 13,43%; nhóm ph i nhi m (hút thu c lá) 0,1343 (hay nhóm khơng ph i nhi m (khơng hút thu c) 0,0954 (hay 9,54%) - Chênh l nh v nguy c c a hai nhóm hút thu c không hút thu c là: 0,0389 (3,89%), 95% CI (2,56% - 5,22%) - Nguy c t ng đ i hay t xu t nguy c c a hai nhóm hút thu c không hút thu c (RR): 1,408, 95% CI (1,248 – 1,588) - T l nguy c m c b nh tim m ch hút thu c (Attributable Risk Proportion, ARP) 28,95% , 95% CI (19,85 - 37,02) - Chi bình ph ng = 31,47; giá tr p = 0,0000 Và có th k t lu n r ng: © Bùi V n Tr ng 4/22/2007 Nh ng ng i hút thu c có nguy c b b nh m ch vành cao g p 1,408 l n (95% CI: 1,248 – 1,588) so v i ng i không hút thu c (p= 0,000) Chênh l ch v t xu t nguy c 3,89% T l nguy c b nh m ch vành hút thu c 28,95%, 95% CI (19,85 - 37,02) 2.2 Tính tóan v i b s li u s n có Trong tr ng h p có b ng s li u, ta có th dùng cú pháp: cs var_case var_exposed [, cs_options] Gi i thích cú pháp l nh: cs: đ c vi t t t t c m t “cohort study” var_case = bi n b nh (1 = có b nh, 0= khơng có b nh) var_exposed = bi n ph i nhi m (1 = ph i nhi m; 0= không ph i nhi m) [, cs_options]: nh ph n tính tốn tr c ti p Ví d Trong b s li u nghiên c u v nguy c hút thu c v i b nh tim m ch (tên c a b s li u fuldata.dta) Câu h i nghiên c u hút thu c có làm t ng nguy c m c b nh tim m ch hay không? Dùng STATA, ta th c hi n b c nh sau: a M b s li u: use ful.dta b Tìm bi n ph i nhi m bi n b nh: - Bi n ph i nhi m (Bi n hút thu c): có tên cur_smok Trong bi n có gía tr : 1= hi n cịn hút thu c; 0= khơng hút thu c) - Bi n b nh (b nh tim m ch) có tên cvd, có hai giá tr : 1= có b nh; 0= khơng b nh © Bùi V n Tr ng 4/22/2007 Chú ý: Các bi n ph i nhi m bi n b nh bi n phân lo i có hai giá tr đ nh d ng N u gía tr khác nh t thi t ph i chuy n v đ nh d ng c Xác đinh ch s đánh giá: Dùng l nh STATA cho nghiên c u thu n t p cho k t qu : cs cvd cur_smok if form_smok!=1 | cur_smok | | Exposed Unexposed | Total -+ + -Cases | 354 367 | 721 Noncases | 2384 3113 | 5497 -+ + -Total | 2738 3480 | 6218 | | Risk | 1292915 1054598 | 1159537 | | | Point estimate | [95% Conf Interval] | + -Risk difference | 0238317 | 0076428 0400206 Risk ratio | 1.225979 | 1.06885 1.406207 Attr frac ex | 1843253 | 0644149 2888672 Attr frac pop | 0905009 | + chi2(1) = 8.49 Pr>chi2 = 0.0036 K t qu t phân tích cho th y : - T l nguy c m c b nh tim m ch c a nhóm hi n hút thu c 0,1293 (hay 12,93%), c a nhóm khơng hút thu c 0,1055 (hay 10,55%) - Chênh l nh v t xu t nguy c c a hai nhóm hút thu c khơng hút thu c là: 0,0238 (2,38%) - Nguy c t ng đ i (RR): 1,226 , 95% CI (1,069 – 1,406) - T l nguy c m c b nh tim m ch hút thu c (ARP) 18,43% - Khi bình ph ng = 8,49, giá tr p = 0,0036 Và ta có th k t lu n r ng : Nh ng ng i hút thu c có nguy c m c b nh tim m ch cao g p 1,226 l n (95% CI: 1,069 – 1,406) so v i nh ng ng i khơng hút © Bùi V n Tr ng 4/22/2007 thu c (p= 0,0036) Chênh l ch v t xu t nguy c 2,38% T l nguy c b nh tim m ch hút thu c 18,43% 2.3 Nghiên c u thu n t p tính theo th i gian- ng i i v i nhi u nghiên c u thu n t p không ph i lúc nhóm ph i nhi m khơng ph i nhi m có th theo dõi đ c toàn b cá th v i y u t ph i nhi m toàn b th i gian nghiên c u Vì đa s tr cá th đ c đ a vào nhóm nghiên c u (entry) c ng có nhi u tr nhi u th i m khác ng h p cá th có th kh i nhóm (withdrawal) tr nghiên c u k t thúc ho c không th theo dõi đ ng h p c cu c c (lost of follow-up) Do th i gian theo dõi cá th nhóm nghiên c u không gi ng Trong tr h p khơng th tính đ c s cá th m u s Vì v y tr S tr ng ng h p ng h p m c m i T xu t m c m i (Incidence Rate) = T ng s đ n v th i gian –ng ng i ta dùng m u s t ng s đ n v th i gian-ng i (person-times) i (person-times) đ c theo dõi c a cá th : Khi trình bày k t qu thu đ c s đ n v th i gian-ng i cá th ti p xúc ho c không ti p xúc v i y u t ph i nhi m ch không ph i t ng s cá th m i nhóm nghiên c u ng i th i gian có th là: n m-ng Trong tr n v th i gian-ng i, tháng-ng i th ng n m –ng i, tu n-ng ng h p s li u nghiên c u đ i (đ n v i) c trình bày b ng b ng 2x2 nh sau: © Bùi V n Tr ng 4/22/2007 Ph i nhi m | Có Khơng | T ng -+ + -Tr ng h p b nh | Th i gian-ng a i | b | a+b PTue PTe | PTe + PTe -+ + -T xu t m c m i| b/PTue a/PTe Qua b ng, ta có th tính đ | (a+b)/PTe + PTe) c: - T xu t m c m i nhóm ph i nhi m: IRe = a/PTe - T xu t nguy c nhóm khơng ph i nhi m: IRue = b/PTe - Chênh l ch v t xu t m c m i (Rate Difference) RD = IRe – IRue - T s t xu t m c m i c a hai nhóm RR (Rate Ratio)= IRe / IRue - Các kho ng tin c y c a ch s RD, RR - T l nguy c thu c (Attributable Rate Proportion) ARP = (RR-1)/RR S d ng STATA ta có th d dàng tính tốn ch s d ch t h c 2.3.1 Tính tốn tr c ti p Trong tr ng h p tính tốn tr c ti p ta dùng cú pháp l nh nh sau: iri #a #b #Na #Nb Gi i thích cú pháp l nh: iri: Vi t t t c a c m t “incidence rate immediate” #a: s tr ng h p m c m i nhóm ph i nhi m #b: s tr ng h p m c m i nhóm khơng ph i nhi m #Na: T ng s đ n v th i gian –ng iđ c theo dõi nhóm ph i nhi m #Nb: T ng s đ n v th i gian –ng iđ c theo dõi nhóm khơng ph i nhi m © Bùi V n Tr ng 4/22/2007 Ví d tìm hi u v nguy c c a hút thu c v i b nh tim m ch, ng ti n hành m t nghiên c u thu n t p t ng lai Sau m t th i gian theo dõi nhóm hút thu c khơng hút thu c cho k t qu xu t hi n b nh tim m ch s có 30 ng 36527,87 n m-ng iđ c theo dõi i ta hai nhóm hút thu c có 354 ng i i m c b nh m ch vành t ng s nhóm khơng hút thu c có 367 ng b nh m ch vành t ng s 48270,23 n m-ng im c i theo dõi Dùng STATA đ tính thông s d ch t nghiên c u thu n t p v i đ n v theo dõi n m - ng i cho k t qu sau: iri 354 367 36527.87 48270.23 | Exposed Unexposed | Total -+ + -Cases | 354 367 | 721 Person-time | 36527.87 48270.23 | 84798.1 -+ + -| | Incidence Rate | 0096912 007603 | 0085025 | | | Point estimate | [95% Conf Interval] | + -Inc rate diff | 0020882 | 0008137 0033627 Inc rate ratio | 1.274654 | 1.098363 1.479088 (exact) Attr frac ex | 2154732 | 0895545 3239077 (exact) Attr frac pop | 1057941 | + (midp) Pr(k>=354) = 0.0006 (exact) (midp) 2*Pr(k>=354) = 0.0012 (exact) K t qu cho th y : - T xu t m c m i nhóm ph i nhi m: IRe = 0,00969 - T xu t nguy c nhóm khơng ph i nhi m: IRue = 0,007603 - Chênh l ch v t xu t m c m i (Rate Difference) RD = 0,0020882 - T s t xu t m c m i c a hai nhóm RR (Rate Ratio)= IRe / IRue = 1,275 kho ng tin c y m c 95% = (1,098 – 1,479) © Bùi V n Tr ng 4/22/2007 - T l nguy c thu c (Attributable Rate Proportion) ARP = (RR-1)/RR = 0.2155 (hay 21,55%) - Giá tr p m t phía = 0,0006; p hai phía = 0,0012 K t lu n: Nh ng ng i hút thu c có nguy c m c b nh tim m ch cao g p 1,275 l n (95% CI: 1.098 - 1.479) so v i ng i không hút thu c (p= 0,0012 ) T l nguy c b nh m ch vành hút thu c 21,55% 2.3.2 Tính tốn b s li u cho s n Trong tr ng h p có s li u, cú pháp: ir case exposed person_time Gi i thích l nh: ir: vi t t t c a c m t “incidence rate” case: tr ng h p b nh ( = có b nh; 0= khơng có b nh) exposed: y u t ph i nhi m (1 = ph i nhi m; = không ph i nhi m) person-time: s đ n v th i gian-ng iđ c theo dõi c a cá th nghiên c u Phân tích nghiên c u b nh - ch ng (Case-Control studies) Nghiên c u b nh ch ng nghiên c u đ i t ng nghiên c u đ ch n c s có b nh (cases) khơng b nh- nhóm ch ng (controls) c lo i b nh mà d đ nh ti n hành nghiên c u Tùy thu c vào gi thuy t nhân qu mà nhà d ch t h c ch n đ ki m đ nh mà nghiên c u b nh ch ng có th thi t k m t nhóm b nh v i m t nhóm ch ng, m t nhóm b nh v i nhóm ch ng ho c nhóm ch ng Thơng th ng ng i ta © Bùi V n Tr ng 10 4/22/2007 th ng thi t k m t nhóm b nh v i m t nhóm ch ng Trong ph n phân tích STATA ch gi i h n thi t k m t nhóm b nh m t nhóm ch ng Phân tích nghiên c u b nh ch ng so sánh t n s ph i nhi m v i m t y u t nguy c gi a nhóm b nh v i nhóm ch ng s li u dich t h c đ tính tốn s k t h p này, c trình bày thành b ng 2x2 Ph i nhi m | Có Khơng | Total -+ + B nh | a c | a + c Ch ng | b d | b + d -+ + Total | Qua b ng, ta có th tính đ a+b c+d | a+b+c+d c - Chênh l ch v ph i nhi m nhóm b nh (Odds among cases): = [a/(a+c)]/[c/(a+c)] - Chênh lêch v ph i nhi m nhóm ch ng (Odds among controls): =[b/(b+d)/[d/(b+d)] -T s chênh OR (Odds Ratio): = Odds among cases/Odds among controls = [a/(a+c)]/[c/(a+c)]:[b/(b+d)/[d/(b+d)]=(a/c)/b/d)=(a.d)/(b.c) - Kho ng tin c y c a OR - Giá tr Chi bình ph ng, m c ý ngh a p 3.1 Tính tốn tr c ti p N u bi t tr c s li u c a b ng 2x2 ta có th dùng l nh sau tính tr c ti p nh sau: cci #a #b #c #d [, cci_options] Trong đó: © Bùi V n Tr ng 11 4/22/2007 cci: vi t t t c a c m t “ case control immediate” #a: s tr ng h p b nh có y u t #b: S tr ng h p b nh khơng có y u t #c: s tr ng h p ch ng có y u t #d: s tr ng h p ch ng khơng có y u t cci_options ph i nhi m ph i nhi m ph i nhi m ph i nhi m mô t cornfield S d ng c l ng Cornfield đ tính SE c a t s chênh tb Tính tốn d a ki m đ nh kho ng tin c y exact Tính tốn p xác c a Fisher's level(#) đ t m c đ tin c y; m c đ nh m c đ tin c y 95% Ví d Nghiên c u v nh h Hill (1953) nghiên c u v i 1357 tr tr ng c a thu c đ n b nh ung th ph i, R Doll va hai nhóm: nhóm b nh nhóm ch ng ng h p ung th ph i, có 1350 tr ng h p có hút thu c lá, ng h p b b nh không hút thu c Nhóm ch ng có 1296 tr 61 tr nhóm b nh ng h p hút thu c ng h p khơng hút thu c Tìm t s chênh c a nhóm b nh nhóm ch ng, chênh l ch có đáng tin c y khơng? Dùng STATA đ tính tốn tr c ti p cho k t qu nh sau: cci 1350 1296 61 Proportion | Exposed Unexposed | Total Exposed -+ + -Cases | 1350 | 1357 0.9948 Controls | 1296 61 | 1357 0.9550 -+ + -Total | 2646 68 | 2714 0.9749 | | | Point estimate | [95% Conf Interval] | + -Odds ratio | 9.077381 | 4.126815 23.59071 (exact) © Bùi V n Tr ng 12 4/22/2007 Attr frac ex | 8898361 | 7576824 9576104 (exact) Attr frac pop | 8852459 | + chi2(1) = 43.98 Pr>chi2 = 0.0000 T s chênh OR = 9,077, 95% CI ( 4,23 – 23,59) Khi bình ph ng = 43,98, giá tr p = 0,000 K t lu n: Nh ng ng i hút thu c có kh n ng b ung th ph i cao g p 9,077, 95% CI ( 4,23 – 23,59) l n so v i nh ng ng i không hút thu c lá, p = 0,000 3.2 Phân tích d a b s li u s n có i v i b s li u nghiên c u, đ tính ch s d ch t h c nghiên c u b nh ch ng ta dùng l nh nh sau: cc var_case var_exposed [, cc_options] Gi i thích cú pháp l nh: Cc: vi t t t c a c m t “case control” var_case: bi n b nh (1 = b nh, 0= ch ng) var_exposed: bi n ph i nhi m (1= có y u t ph i nhi m, 0= khơng có y u t ph i nhi m) cc_options Mơ T cornfield S d ng c l ng Cornfield đ tính SE c a t s chênh tb Tính tốn d a ki m đ nh kho ng tin c y exact Tính tốn p xác c a Fisher's level(#) đ t m c đ tin c y; m c đ nh m c đ tin c y 95% Ví d S d ng b s li u nghiên c u b nh ch ng đ tìm hi u v m i liên quan c a u ng r u v i ung th th c qu n B s li u s d ng k-tquan.dta.Dùng STATA, ta th c hi n b c nh sau: © Bùi V n Tr ng 13 4/22/2007 a M b s li u: use k-tquan.dta b Tìm bi n b nh/ch ng bi n ph i nhi m/không ph i nhi m: - Bi n ph i nhi m (u ng r tr : 1= có u ng r u): có tên alcohol Trong bi n có gía u; 0= khơng khơng u ng r u) - Bi n b nh (Ung th th c qu n) có tên case, có hai giá tr : 1= có b nh; 0= tr ng h p ch ng(khơng b nh) 3.3 Phân tích phiên gi i k t qu nghiên c u b nh ch ng Dùng l nh STATA cho nghiên c u thu n t p cho k t qu : cc case alcohol Proportion | Exposed Unexposed | Total Exposed -+ + -Cases | 96 104 | 200 0.4800 Controls | 109 666 | 775 0.1406 -+ + -Total | 205 770 | 975 0.2103 | | | Point estimate | [95% Conf Interval] | + -Odds ratio | 5.640085 | 3.937435 8.061794 (exact) Attr frac ex | 8226977 | 7460276 8759581 (exact) Attr frac pop | 3948949 | + chi2(1) = 110.26 Pr>chi2 = 0.0000 T s chênh OR = 5,64, 95% CI ( 3,94 - 8,06 ) Chi bình ph ng = 110,26 , giá tr p = 0,000 K t lu n: Nh ng ng i u ng r u có kh n ng b ung th th c qu n cao g p 5,64 l n, 95% CI ( 3,94 - 8,06 ) so v i nh ng ng i không hút thu c lá, p = 0,000 Phân tích b nh ch ng ghép c p Ghép c p m t n l c đ lo i b sai s ch n có th khơng th l a ch n vào nhóm M c tiêu đ ghép b nh ch ng v i y u t ngo i tr y u t © Bùi V n Tr ng 14 4/22/2007 nguy c N u có m t s khác bi t v s tr ng b nh tr ng h p nghiên c u ghép c p ph i nhi m v i y u t nguy c ch không ph i y u t ghép B ng x2 cho nghiên c u b nh ch ng ghép c p nh sau: Ch ng (controls) T ng Ph i nhi m Không ph i nhi m Ph i nhi m a b a+b Không ph i nhi m c d c+d a+c b+d a+b+c+d B nh (cases) T ng S d ng ki m đ nh chi bình ph ng c a McNemar cho nghiên c u b nh ch ng ghép c p Công th c: Chi bình ph ng c a McNemar = (b-c)2/(b+c) Trong STATA nghiên c u b nh ch ng ghép c p có th dùng l nh nh sau: mcci a b c d Gi i thích l nh: mcci: vi t t t c a c m t “matching case control intemediate” a: s c p có tr ng h p b nh ch ng ph i nhi m v i y u t nguy c b: s c p có tr ng h p c p b nh ph i nhi m v i y u t nguy c ch ng không ph i nhi m v i y u t nguy c c: s c p có tr ng h p b nh khơng ph i nhi m v i y u t nguy c ch ng có ph i nhi m v i y u t nguy c d: s c p có tr ng h p b nh ch ng không ph i nhi m v i y u t nguy c Ví d Trong m t nghiên c u b nh ch ng tìm hi u nguy c c a hút thu c th đ ng v i ung th ph i, nhà nghiên c u ti n hành tìm hi u nguy c ung th ph i 500 ng (case) đ i đàn ông đ c ch n đoán ung th ph i M i tr ng h p ch n đoán c ghép c p theo tu i, ch ng t c ngh nghi p v i m t tr ng h p © Bùi V n Tr ng 15 4/22/2007 ch ng (control) không b b nh ung th ph i K t qu v tình tr ng b nh ph i nhi m đ c th hi n b ng d i B nh có/khơng hút Ch ng có/khơng hút S c p thu c th đ ng thu c th đ ng Khơng Khơng 100 Khơng Có 33 Có Khơng 117 Có Khơng 250 S d ng l nh STATA tr c ti p cho nghiên c u ghép c p cho k t qu nh sau: mcci 250 117 33 100 | Controls Cases | Exposed | Unexposed | Total -+ + -Exposed | 250 117 | 367 Unexposed | 33 100 | 133 -+ + -Total | 283 McNemar's chi2(1) = 47.04 217 | 500 Prob > chi2 = 0.0000 Exact McNemar significance probability = 0.0000 Proportion with factor Cases 734 Controls 566 [95% Conf Interval] - difference 168 120305 215695 1.29682 1.203736 1.397101 rel diff .3870968 3004944 4736992 odds ratio 3.545455 2.391508 5.388252 ratio (exact) T k t qu ta có: Khi bình ph ng c a McNemar = 47,04 Giá tr p = 0,000 © Bùi V n Tr ng 16 4/22/2007 T l v i v iy ut : - B nh = (a+b)/n = 0,73 - Ch ng = (a+c)/n = 0,57 S khác nhau: B nh – Ch ng = (0,73- 0,57 ) = 0,17 T s : B nh/Ch ng = 1,29 ng đ i = (b-c)/(b+d) = 0,387 S khác t T s chênh (OR) = ad/bc = 3,55, 95% CI (2,39 – 5,39) i hút thu c th đ ng có kh bi ung Và ta có th k t lu n r ng nh ng ng i không hút thu c th đ ng (p th ph i cao g p 3,55 l n so v i nh ng ng =0,000, Khi bình ph ng= 47,04) Nhi u kh ng ch nhi u Trong nghiên c u d ch t h c ng i ta đ c p r t nhi u đ n y u t nhi u (confounding factors) Y u t nhi u y u t làm l ch m i quan h gi a b nh yêu t ph i nhi m Bi n nhi u m t bi n th ba có liên quan c nhi m b nh nghiên c u M t bi n đ y u t ph i c g i bi n nhi u n u: Nó m t bi n đ c l p, không ph i bi n trung gian c a bi n b nh Nó y u t liên quan đ n b nh v i y u t ph i nhi m Ki m soát kh ng ch nhi u r t đ c quan tâm nghiên c u d ch t h c y u t nhi u có th làm nh n đ nh sai m i quan h gi a b nh y u t ph i nhi m Chúng ta có th ki m soat nhi u vào th i m thi t k nghiên c u nh thu h p ph m vi nghiên c u, ghép c p ch n ng u nhiên Trong q trình phân tích có th ki m sốt nhi u b ng phân t ng, chu n hóa, ghép c p, ho c phân tích đa bi n Trong ch ng gi i thi u s d ng STATA phân tích phân t ng nghiên c u d ch t h c đ ki m soát y u t nhi u (ph n chu n hố, phân tích đa bi n s đ c gi i thi u m t d p khác) © Bùi V n Tr ng 17 4/22/2007 5.1 Phân tích phân t ng Trong b ng ti p liên 2x2 không xem xét đ n y u t nhi u (confounding) ho c y u t t ng tác (interaction) v i y u t ph i nhi m nghiên c u v i y u t nguy c khác Trong phân tích phân t ng, tìm hi u m i quan h ph i nhi m b nh có xem xét đ n y u t t nhi u ho c y u t ng tác ho c y u t làm thay đ i tác đ ng (effect modification) c a m i quan h gi a ph i nhi m b nh STATA cho phép th c hi n th c hi n phân tích phân t ng, qua có th xác đ nh đ c m i quan h c a y u t ph i nhi m b nh sau ki m soát (ho c hi u ch nh) cho y u t nhi u L nh đ t ng nh sau: - c khái quát nghiên c u thu n t p b nh ch ng có phân i v i nghiên c u thu n t p: cs case exposed, by (stavar) - i v i nghiên c u thu n t p tính theo đ n v th i gian-ng i: ir case, exposed person_year, by (stavar) - i v i nghiên c u b nh ch ng: cc case exposed, by (stavar) Gi i thích l nh: cc: vi t t t c m t “case control” , cs: vi t t t c m t “cohort study” ir: vi t t t c m t case: tr “ incidence rate” ng h p b nh (1= có b nh, = khơng có b nh) exposed: y u t ph i nhi m (1= ph i nhi m; 0= không ph i nhi m) person_year: th i gian-ng i theo dõi stavar: bi n phân t ng ây bi n nghi ng có th y u t gây nhi u ho c y u t làm thay đ i tác đ ng c a m i quan h gi a ph i nhi m b nh © Bùi V n Tr ng 18 4/22/2007 Ví d Trong m t nghiên ng r i ta mu n tìm hi u c u m i liên quan gi a u ng u ung th th c qu n (S d ng b s li u có tên k-tquan.dta) Trong q trình phân tích s li u ng i ta mu n tìm xem hút thu c có ph i y u t nhi u ho c y u t làm thay đ i tác đ ng c a m i quan h gi a u ng r u ung th th c qu n (n u có) hay khơng? Ta ti n hành nh sau: M b s li u: use k-tquan.dta" Tìm bi n b s li u: - case = tr ng h p b nh (1 = b ung th th c qu n, 0= không b b nh) - alcohol = u ng r u (1= có u ng r u; = khơng u ng r u) - tod4 = hút thu c (1= có hút thu c lá, = khơng hút thu c lá) Tìm hi u m i liên quan gi a u ng r cc u ung th th c qu n case alcohol Proportion | Exposed Unexposed | Total Exposed -+ + -Cases | 96 104 | 200 0.4800 Controls | 109 666 | 775 0.1406 -+ + -Total | 205 770 | 975 0.2103 | Point estimate | [95% Conf Interval] | + -Odds ratio | 5.640085 | 3.937435 8.061794 (exact) Attr frac ex | 8226977 | 7460276 8759581 (exact) Attr frac pop | 3948949 | + chi2(1) = 110.26 Pr>chi2 = 0.0000 Ta có: OR = 5,64, 95% CI ( 3.94 - 8.06 ) Chi bình ph ng = 110,26, giá tr p = 0,000 K t lu n: © Bùi V n Tr ng 19 4/22/2007 Nh ng ng i u ng r u có kh n ng b ung th th c qu n cao h n 5,64 l n, 95%CI ( 3,94 – 8,06) so v i nh ng ng i không u ng r u Nhà nghiên c u nghi ng hút thu c có th y u t nhi u ho c y u t làm thay đ i tác đ ng, ti n hành phân tích phân t ng theo tình tr ng hút thu c không hút thu lá, dùng l nh STATA cho k t qu sau: cc case alcohol, by ( tob4) Tob4 | OR [95% Conf Interval] M-H Weight -+ | 5.414411 3.678267 7.946703 10.17917 (exact) | 6.526786 2.080228 21.13344 1.365854 (exact) -+ Crude | 5.640085 3.937435 8.061794 (exact) M-H combined | 5.546013 3.920895 7.844703 Test of homogeneity (M-H) chi2(1) = 0.11 Pr>chi2 = 0.7384 Test that combined OR = 1: Mantel-Haenszel chi2(1) = 105.09 Pr>chi2 = 0.0000 Phiên gi i: - Tính OR t ng: Trong tr 5,41 Trong tr ng h p không hút thu c (tod4 =0 ) ta có OR = ng h p hút thu c (tod4 =1): OR = 6.53 - Xác đ nh xem y u t hút thu c có ph i y u t thay đ i tác đ ng ho c y u t t ng tác v i y u t ph i nhi m hay không ta dùng ki m đ nh tính thu n nh t c a OR qua t ng (Test of homogeneity of ORs) Qua b ng cho th y Chi bình ph t ng = 0,11, giá tr p = 0,7384 Do hút thu c không ph i y u t ng tác (interaction) v i y u t ph i nhi m ho c làm thay đ i tác đ ng (effect modification) c a m i quan h ph i nhi m b nh - V y hút thu c có th y u t nhi u Trong phân tích phân t ng cho th y: OR thơ (không xét đ n y u t hút thu c lá) = 5,64, 95% CI (3,94 - 8,06) OR c a Mantel-Haenszel hi u ch nh cho y u t nhi u (hút thu c lá)= 5,55, 95% CI © Bùi V n Tr ng 20 4/22/2007 (3,83 - 8,02) Chi bình ph ng c a Mantel-Haenszel = 105,09, giá tr p = 0,0000 K t lu n: Sau ki m soát/hi u ch nh cho y u t hút thu c cho th y nh ng ng u ng r i u có kh n ng b ung th th c qu n cao g p 5,55 l n (OR = 5.546013) so v i nh ng ng i nh ng ng i không u ng r u, v i chi bình ph ng = 105.09, p = 0,000 Ví d Trong nghiên c u v nguy c c a hút thu c v i b nh tim m ch ví d 2, phân t ng theo ch ng t c đ hi u ch nh cho y u t ch ng t c Dùng l nh STATA cho k t qu sau: cs cvd cur_smok, by (race) chung toc | RR [95% Conf Interval] M-H Weight -+ | 1.263619 1.090991 1.463562 138.6631 | 9580051 6458025 1.421137 22.23586 | 8571429 2340296 3.13932 2.333333 -+ Crude | 1.219202 1.06347 1.397739 M-H combined | 1.216177 1.06076 1.394365 Test of homogeneity (M-H) chi2(2) = 1.946 Pr>chi2 = 0.3779 - Ki m đ nh tính thu n nh t RR (Test of homogeneity) gi a t ng đ xem y u t ch ng t c có ph i y u t t ng tác v i y u t ph i nhi m (hút thu c lá) không? Qua b ng cho th y chi2 = 1,946 giá tr p = 0,3779 DO v y ch ng t c không ph i y u t t ng tác (interaction) ho c y u t hay đ i tác đ ng (effect modification) m i quan h gi a b nh tim m ch hút thu c c ng không ph i y u t t ng tác v i hút thu c đ làm thay đ i m i quan h - So sánh gi a RR thô (không xét đ n y u t ch ng t c) OR c a MantelHaenszel hi u ch nh cho y u t c ng t c xem có s khác bi t khơng Qua b ng cho © Bùi V n Tr ng 21 4/22/2007 th y khơng có s khac bi t (RR thô = 1,219; RR hi u ch nh c a Mantel-Haenszel = 1.,22) Do v y, y u t ch ng t c không ph i y u t nhi u Ta có th k t lu n d a vào RR thô (ho c RR hi u ch nh) K t lu n: Nh ng ng i hút thu c có nguy c b b nh tim m ch cao g p 1,22 l n, 95% CI (1,06 1,40 ) so v i nh ng ng i không hút thu c © Bùi V n Tr ng 22