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)
ไม่มีความเห็น