Sample Size สำหรับ 2 Proportions เช่น Cross-Sectional, Cohort, Randomized Clinical Trial

STATA11 ใช้คำสั่ง sampsi เช่น sampsi 0.4 0.6 , STATA 14 มีคำสั่ง power twoprop 0.4 0.6 ,

STATA11
ใช้คำสั่ง sampsi เช่น sampsi 0.4 0.6
ถ้าไม่ระบุ r(#) ใน option sampsi จะใช้ค่า r(1) คือ n2/n1 = 1
ถ้าไม่ระบุ power ใน option sampsi จะใช้ค่า power = 0.9
ถ้าไม่ระบุ nocontin ใน option sampsi จะทำ continuity correction มาให้ด้วย


. sampsi 0.4 0.6
Estimated sample size for two-sample comparison of proportions
Test Ho: p1 = p2, where p1 is the proportion in population 1
and p2 is the proportion in population 2
Assumptions:
. . . . alpha = 0.0500 (two-sided)
. . . . power = 0.9000
. . . . p1 = 0.4000
. . . . p2 = 0.6000
. . . . n2/n1 = 1.00
Estimated required sample sizes:
. . . . n1 = 140
. . . . n2 = 140


คำสั่ง power twoprop p1 p2 ใน STATA 14
ระบุ power ตามต้องการ, ถ้าไม่ระบุุ power จะใช้ power = 0.8

การคำนวณจะไม่ใช้ continuity correction แต่ถ้าต้องการ continuity correction ให้ระบุ continu ในส่วน option


. power twoprop 0.4 0.6, power(0.9) continu
Performing iteration ...
Estimated sample sizes for a two-sample proportions test
Pearson's chi-squared test
Ho: p2 = p1 versus Ha: p2 != p1
Study parameters:
. . . . alpha = 0.0500
. . . . power = 0.9000
. . . . delta = 0.2000 (difference)
. . . . p1 = 0.4000
. . . . p2 = 0.6000
Estimated sample sizes:
. . . . N = 280
. . . . N per group = 140
.


ค่าที่ได้จะตรงกันกับ Openepi (Fleiss with Continuity Correction)
Exposed = 140
Unexposed = 140
http://www.openepi.com/SampleSize/SSCohort.htm






ถ้าจำนวน Unexposed เป็น 2 เท่าของ กลุ่ม Exposed
Openepi ใช้ Unexposed / Exposed ค่า (n1/n2) = 2

STATA 11 ใช้ option r(0.5) หมายถึง n2/n1 = 0.5
STATA 14 ใช้ option nratio(0.5) หมายถึง n2/n1 = 0.5


สูตรที่ใช้ในการคำนวณขนาดตัวอย่าง สำหรับ 2 Proportions
เช่น Cross-Sectional, Cohort, RCT

https://dl.dropboxusercontent.com/u/1999671/EPI_56/ss/s92.htm

ต้วอย่าง
alpha = 0.01, Za = 2.576
power = 0.9, Zb = 1.282
กลุ่มควบคุมป่วย 10% (P0 = 0.1)
กลุ่มทดลองป่วย 2.5% (P1 = 0.025)
RR = P1/P0 = 0.25

ตอบ
n0 = 308
n1 = 308


Source: ทวิป กิตยาภรณ์. Sample Size Estimation. บทที่ 9 วิจัยคลินิก. หน้า 73


sampsi 0.1 0.025 , alpha(0.01) power(0.9) nocont

(STATA 11 ยังไม่มีคำสั่ง power)


power twoprop 0.1 0.025 , alpha(0.01) power(0.9)


บันทึกนี้เขียนที่ GotoKnow โดย  ใน EPISTAT



ความเห็น (0)