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

10

Сравнение пропорций с двумя выборками, оценка размера выборки: 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

dwstu
источник

Ответы:

7

Различие связано с тем, что 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 ), указывает, что формула без исправления непрерывности является довольно точной, тем самым обеспечивая некоторое оправдание для отказа от исправления.

Фил Шумм
источник
2
Отличный ответ. Кажется, что причина заключается не в разнице между двумя методами в моем посте, а в том, что один из этих методов использует коррекцию непрерывности, а другой - нет.
Майкл М
1
Спасибо. При наличии только двух пропорций (то есть таблицы 2х2) не имеет значения, указали ли вы альтернативу в виде двух пропорций или одной пропорции и отношения шансов. А поскольку точный тест Фишера является консервативным для биномиальной задачи с двумя выборками, оценки мощности, основанные на этом, ближе к оценкам формулы с исправлением непрерывности.
Фил Шумм
1
bsamsize(0.70, 0.85, alpha=0.05, power=0.90)N1знак равноN2знак равно160.7777