Двопробне порівняння пропорцій, оцінка розміру вибірки: R проти Stata

Двопробне порівняння пропорцій, оцінка розміру вибірки: R проти Stata

Я отримав різні результати щодо розмірів вибірки:

В Р

power.prop.test(p1 = 0.70, p2 = 0.85, power = 0.90, sig.level = 0.05)

Результат: $n = 160.7777$ (так 161) для кожної групи.

У штаті

sampsi 0.70 0.85, power(0.90) alpha(0.05)

Результат: для кожної групи.$n = 174$

Чому різниця? Дякую.

До речі, я провів той самий розрахунок розміру вибірки в SAS JMP , результат: (майже такий же, як результат R).$n = 160$

Різниця пов’язана з тим, що sampsiкоманда Stata (застаріла як Stata 13 і замінена на power) використовує корекцію безперервності за замовчуванням, тоді як R - power.prop.test()ні (детальніше про формулу, яку використовує Stata, див. У [PSS] power twoproportions ). Це можна змінити за допомогою nocontinuityпараметра, наприклад,

sampsi 0.70 0.85, power(0.90) alpha(0.05) nocontinuity

що дає розмір вибірки 161 на групу. Використання корекції безперервності дає більш консервативний тест (тобто, більший розмір вибірки), і, очевидно, має значення менше, оскільки розмір вибірки збільшується.

Френк Харрелл у документації до bpower(частини свого пакету Hmisc ) вказує, що формула без корекції безперервності є досить точною, тим самим надаючи певне обґрунтування відмови від виправлення.