Сравнение пропорций с двумя выборками, оценка размера выборки: R против Stata

Сравнение пропорций с двумя выборками, оценка размера выборки: R против Stata

Я получил разные результаты для размеров выборки, а именно:

В R

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

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

В стате

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] ). Это можно изменить с помощью nocontinuityопции, например,

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

что дает размер выборки 161 на группу. Использование коррекции непрерывности дает более консервативный тест (т. Е. Больший размер выборки) и, очевидно, имеет меньшее значение при увеличении размера выборки.

Фрэнк Харрелл, в документации для bpower(часть его пакета Hmisc ), указывает, что формула без исправления непрерывности является довольно точной, тем самым обеспечивая некоторое оправдание для отказа от исправления.