Câu lệnh Stata tính tương quan giữa hai biến định lượng

5 0 0
Câu lệnh Stata tính tương quan giữa hai biến định lượng

Đang tải... (xem toàn văn)

Thông tin tài liệu

Muoán bieát caùc bieán coù moái quan heä vôùi nhau nhö theá naøo, chuùng ta phaûi tính heä soá töông quan (kyù hieäu laø r) Heä soá töông quan naøy coù giaù trò töø –1 ñeán 1. Keát quaû heä soá töông quan seõ coù yù nghóa khaùc nhau nhö sau r=0 :khoâng coù moái töông quan giöõa hai bieán quan taâm r 0 : coù moái töông quan thuaän 0 < |r| < 0.2 : coù moái töông quan raát yeáu 0.2< |r| < 0.5 : coù moái töông quan yeáu 0.5 < |r| < 0.7 : coù moái töông quan trung bình 0.7 < |r| < 0.9 : coù moái töông quan chaët cheõ > 0.9 : coù töông quan raát chaët cheõ. Nhöng treân thöïc teá 0.5 < |r| < 0.7 ñaõ ñöôïc goïi laø töông quan khaù chaët cheõ. Chuùng ta coù theå tính ñöôïc heä soá töông quan baèng leänh corr hay pwcorr vaø chuùng ta coù theå xem ñöôïc caùc moái töông quan naøy baèng ño àthò khi goõ leänh graph matrix ñeå veõ ma traän ñoà thò chaám ñieåm. Cuù phaùp

Töông quan giöõa hai hay nhieàu bieán ñònh löôïng Muoán bieát caùc bieán coù moái quan heä vôùi nhau nhö theá naøo, chuùng ta phaûi tính heä soá töông quan (kyù hieäu laø r) Heä soá töông quan naøy coù giaù trò töø –1 ñeán 1 Keát quaû heä soá töông quan seõ coù yù nghóa khaùc nhau nhö sau r=0 :khoâng coù moái töông quan giöõa hai bieán quan taâm r 0 : coù moái töông quan thuaän 0 < |r| < 0.2 : coù moái töông quan raát yeáu 0.2< |r| < 0.5 : coù moái töông quan yeáu 0.5 < |r| < 0.7 : coù moái töông quan trung bình 0.7 < |r| < 0.9 : coù moái töông quan chaët cheõ > 0.9 : coù töông quan raát chaët cheõ Nhöng treân thöïc teá 0.5 < |r| < 0.7 ñaõ ñöôïc goïi laø töông quan khaù chaët cheõ Chuùng ta coù theå tính ñöôïc heä soá töông quan baèng leänh corr hay pwcorr vaø chuùng ta coù theå xem ñöôïc caùc moái töông quan naøy baèng ño àthò khi goõ leänh graph matrix ñeå veõ ma traän ñoà thò chaám ñieåm Cuù phaùp Tính heä soá töông quan giöõa hai bieán hay corr varlist hoaëc ma traän caùc heä soá töông quan cuûa pwcprr varlist , sig moät nhoùm caùc bieán ñònh löôïng star(5) Tính heä soá töông quan cho töøng caëp bieán ñònh löôïng, caùc option coù yù nghóa: sig : theå hieän möùc yù nghóa thoáng keâ cuûa heä soá töông quan star(5) : ñaùnh daáu sao vaøo nhöõng heä soá coù yù nghóa graph matrix varlist…., thoáng keâ options Caùc option söû duïng trong caâu leänh naøy laø: half : chæ veõ moät beân ma traän jitter( ) :loaïi boû ngaãu nhieân nhöõng chaám töôïng tröng cho soá löôïng raát nhoû traùnh vieäc nhöõng ñieåm naøy truøng laép msymbol( ) : bieåu töôïng cuûa daáu chaám (hình troøn daïng loã, troøn…) msize( …) : kích thöôùc cuûa chaám diagonal : daùn nhaõn cho caùc bieán corr weight iq age (obs=100) | weight iq age -+ - weight | 1.0000 iq | -0.2920 1.0000 age | 0.4131 -0.4363 1.0000 Dieãn giaûi keát quaû:  Heä soá töông quan naøy ñöôïc tính treân soá quan saùt chung (obs=100) traû lôøi ñuû ba caâu hoûi veà weight, iq, vaø age Soá  Giöõa iq vaø weight coù töông quan nghòch yeáu  Giöõa weight vaø age coù töông quan thuaän yeáu  Giöõa iq vaø age coù töông quan nghòch yeáu pwcorr weight iq age, obs sig star(5) | weight iq age -+ - weight | 1.0000 | | 106 | iq | -0.3006* 1.0000 | 0.0025 | 99 109 | age | 0.4152* -0.4366* 1.0000 | 0.0000 0.0000 | 106 109 117 | Dieãn giaûi keát quaû  Leänh pwcorr duøng ñeå tính heä soá töông quan cho töøng caëp bieán vì vaäy daân soá ñeå tính seõ khaùc nhau tuøy vaøo töøng soá quan saùt traû lôøi ñuû cho caëp bieán ñoù  Coù 99 quan saùt traû lôøi ñuû hai caâu hoûi veà iq vaø weight Giöõa bieán iq vaø bieán weight coù töông quan nghòch yeáu coù yù nghóa thoáng keâ (p=0,00025)  Coù 106 quan saùt traû lôøi ñuû hai caâu hoûi veà age vaø weight Giöõa bieán age vaø bieán weight coù töông quan thuaän yeáu coù yù nghóa thoáng keâ (p=0,000)  Coù 109 quan saùt traû lôøi ñuû hai caâu hoûi veà iq vaø age Giöõa bieán iq vaø bieán age coù töông quan nghòch yeáu coù yù nghóa thoáng keâ (p=0,00) graph matrix weight iq age, half jitter(1) msymbol(circle_hollow) msize(small) Dieãn giaûi keát quaû:  OÂ ñoà thò chaám ñieåm giöõa bieán weight vaø bieán iq ñi xuoáng chöùng toû töông quan nghòch Caùc chaám ñieåm rôøi raïc khoâng thaønh ñöôøng thaúng chöùng toû hai bieán töông quan yeáu  OÂ ñoà thò chaám ñieåm cuûa bieán age vaø bieán weight ñi leân chöùng toû chuùng töông quan thuaän Tuy nhieân caùc chaám ñieåm cuõng rôøi raïc neân chuùng ta cuõng keát luaän raèng chuùng töông quan yeáu  OÂ ñoà thò chaám ñieåm cuûa bieán iq vaø bieán age cuõng coù hình daïng töông töï nhö giöõa bieán weight vaø bieán iq Töông quan giöõa moät bieán ñònh löôïng vaø moät bieán ñònh tính Trong tröôøng hôïp moät hoaëc hai bieán ñònh löôïng coù phaân boá khoâng chuaån hay moät bieán laø bieán ñònh tính, chuùng ta khoâng theå söû duïng leänh corr hay pwcorr ñeå tìm heä soá töông quan maø chuùng ta phaûi söû duïng leänh Spearman vaø leänh ktau ñeå tính heä soá töông quan spearman vaø heä soá töông quan ktau Chuùng ta cuõng coù theå söû duïng ñoà thò ñeå moâ taû moái töông quan baèng leänh twoway Cuù phaùp spearman var1 var2 Tính heä soá töông quan spearman ktau var1 var2 Tính heä soá tau-b cuûa kendall twoway (scatter bpt1 bpt2 if Caùc option ñöôïc söû duïng trong bñl==#, option) (scatter bpt1 leänh naøy: bpt2 if bñ;==##, option) , msymbol( ) : bieåu options töôïng cuûa chaám (troøn, troøn loã, daáu x ) mcolor (…) : maøu cuûa chaám (ñen, traéng ) jitter( ) : Ví duï: spearman weight life Number of obs = 105 Spearman's rho = 0.0622 Test of Ho: weight and life are independent Prob > |t| = 0.5282 ktau weight life Number of obs = 105 Kendall's tau-a = 0.0359 Kendall's tau-b = 0.0513 Kendall's score = 196 SE of score = 308.818 (corrected for ties) Test of Ho: weight and life are independent Prob > |z| = 0.5278 (continuity corrected) Dieãn giaûi keát quaû  Keát quaû cho thaáy moái töông quan Spearman(0.0622) vaø Kendall(0.0513) giöõa hai bieán weight vaø bieán life ñeàu raát yeáu  Heä soá ktau-b baèng 0,0513 vôùi p=0,5278 chöùng toû giöõa hai bieán weight vaø life khoâng coù moái töông quan twoway (scatter weight age if life==1, msymbol(circle) mcolor(black) jitter(1)) (sca if life==2, msymbol(x) mcolor(black) jitter(1)), legend > tter weight age (order(1 "no > " 2 "yes")) Diễn giải kết quả Ñoà thò chaám cho thaáy ôû caû nhoùm coù yù ñònh töï töï vaø khoâng coù yù ñònh töï töï coù söï taêng troïng nhieàu ôû nhoùm phuï nöõ coù ñoä tuoåi cao

Ngày đăng: 19/03/2024, 16:57

Tài liệu cùng người dùng

Tài liệu liên quan