Перегляд вихідного коду сторінки Статистична потужність

'''Потужність''' або '''чутливість''' перевірки подвійної [[Перевірка статистичних гіпотез|статистичної гіпотези]]&nbsp;— це ймовірність того, що тест правильно відкидає [[Нульова гіпотеза|нульову гіпотезу]] (Н<sub>0</sub>) коли альтернативна гіпотеза (Н<sub>1</sub>) істина. Це можна еквівалентно розглядати як ймовірність прийняття альтернативної гіпотези (Н<sub>1</sub>) Коли це істина, то тест може виявити ефект, якщо ефект дійсно існує. Тобто
: <math> \text{power} = \mathbb P\big( \text{reject} H_0 \big| H_1 \text{ is true} \big) </math>
Потужність тесту іноді, менш формально, належить до ймовірності відхилення нульової, якщо це не є істиною, хоча це не формальне визначення, як зазначено вище. Потужність в загальному випадку є функцією можливих розподілів, часто визначається параметром, згідно з альтернативною гіпотезою. Якщо потужність збільшується, то зменшується ймовірність [[Похибки першого і другого роду|помилки II типу]] (хибно негативна), яку також називають [[Похибки першого і другого роду|похибкою першого та другого роду]] (β), оскільки потужність дорівнює 1−β, знову ж таки, згідно з альтернативною гіпотезою. Аналогічне поняття [[Похибки першого і другого роду|помилки першого роду]] також називають «хибнопозитивною» або на рівні тесту при нульовій гіпотезі.

Аналіз потужності може бути використаний для розрахунку мінімального [[Визначення розмірів вибірки|обсягу вибірки]], необхідного, щоб можна з достатнім ступенем ймовірності визначити вплив заданого розміру. Наприклад: «скільки разів треба кинути монету, щоб зробити висновок, що це фальсифіковано ?»<ref>{{Cite web |url=http://www.statisticsdonewrong.com/power.html |title=Statistical power and underpowered statistics |website=statisticsdonewrong.com |author=Alex Reinhart |language=en |accessdate=10 грудня 2016 |archive-date=20 січня 2017 |archive-url=https://web.archive.org/web/20170120101701/https://www.statisticsdonewrong.com/power.html }}</ref> Аналіз потужності може також використовуватися для розрахунку мінімального ефекту, який може бути виявлений у ході дослідження, використовуючи даний розмір вибірки. Крім того, поняття потужності використовується для порівняння між різними процедурами статистичних випробувань: наприклад, між параметричним і непараметричним тестом на тій же гіпотезі.

== Передумова ==
[[Перевірка статистичних гіпотез|Статистичні тести]] використовують дані зі [[Відбір вибірки (статистика)|зразків]] для оцінки, чи роблять [[Статистичне висновування|висновки]] про [[Генеральна сукупність|генеральну сукупність]]. У конкретній ситуації мета порівняння двох зразків, полягає в тому, щоб оцінити, чи відрізняються середні значення деякого атрибута, отримані для фізичних осіб в двох підгрупах. Наприклад, для перевірки нульової гіпотези про те, що [[середнє значення]] [[Функція внеску|внеску]] чоловіків і жінок на тесті не відрізняються, виконується вибірка чоловіків і жінок, проводиться тест, потім результат однієї групи порівнюють з іншою групою з використанням статистичних тестів, наприклад двохвибіркового [[Z-критерій|''z''-тесту]]. Потужністю тесту є ймовірність того, що тест знайде статистично значущу різницю між чоловіками і жінками, як функцію від розміру істинної різниці між цими двома популяціями.

== Фактори, що впливають на потужність ==
Статистична потужність може залежати від цілого ряду факторів. Деякі чинники можуть бути пов'язані з конкретними ситуаціями у тестуванні, але, як мінімум, потужність майже завжди залежить від таких трьох чинників:
* від [[Статистична значущість|статистичної значущості]] критеріїв, використовуваних у тесті
* величина ефекту інтересу населення
* [[Визначення розмірів вибірки|розмір вибірки]] використовується для визначення ефекту
'''Критерій значущості'''&nbsp;— це твердження про те, що навряд чи позитивний результат повинен бути, якщо нульова гіпотеза про відсутність ефекту істина, нульова гіпотеза повинна бути відхилена. Найбільш часто використовувані критерії ймовірності 0.05 (5&nbsp;%, 1 з 20), 0.01 (1&nbsp;%, 1 до 100), 0.001 (0.1&nbsp;%, 1 на 1000). Якщо критерій становить 0,05, достовірність даних, маючи на увазі ефект, як мінімум, спостережуваний ефект, коли нульова гіпотеза вірна, повинен бути не менше 0,05, для нульової гіпотези про відсутність ефекту буде відхилено. Один простий спосіб збільшити потужність тесту є проведення менш консервативного тесту з використанням більшого значення критерію значущості, наприклад 0.10 замість 0.05. Це збільшує ймовірність відхилити нульову гіпотезу (тобто отримання статистично значущого результату), коли нульова гіпотеза хибна, то зменшується ризик [[Похибки першого і другого роду|типу II помилка]] (псевдонегативний щодо того, чи існує ефект). Але це також збільшує ризик отримання статистично значущого результату (тобто відхилення нульової гіпотези), коли нульова гіпотеза не хибна; тобто, це збільшує ризик [[Похибки першого і другого роду|помилки першого роду]] (помилковий позитивний результат).

'''Величина ефекту''' інтересу населення може бути визначена кількісно з точки зору ефекту, де є велика потужність, щоб виявити більше ефектів. Розмір ефекту може бути безпосереднє оцінкою величини інтересу, або він може бути стандартним показником, що також пояснює мінливість популяції. Наприклад, аналіз порівняння результатів експериментальної та контрольної популяції, різниця результатів означає Y&nbsp;−&nbsp;X було б прямою мірою ефекту розміру, в той час як (Y&nbsp;−&nbsp;X)/σ, де σ є загальне стандартне відхилення результату в оброблених і контрольних групах, був би стандартизований розмір ефекту. Якщо побудований відповідним чином, стандартизований розмір ефекту, поряд з розміром вибірки, то він буде повністю визначати потужність. Стандартизовано (прямий) розмір ефекту буде досить рідко, щоб визначити силу, так як він не містить інформації про мінливість в вимірах.
'''Розмір вибірки''' визначає обсяг помилок вибірки, властиві результату тесту. За інших рівних умов, ефекти важче виявити в невеликих зразках. Збільшення розміру вибірки найчастіше є найпростішим способом підвищити статистичну потужність тесту. Як збільшення обсягу вибірки призводить до більш високої потужності&nbsp;— це показник ефективності тесту. Наприклад, розмір вибірки, необхідної для даної потужності.{{sfn|Everitt|2002|loc=p. 321}}

Точність, з якою вимірюються дані також впливає на статистичну потужність. Отже, потужність часто може бути поліпшена за рахунок зменшення [[Похибка вимірювання|похибки вимірювання]] в даних. Пов'язана з цим концепція є підвищення «надійності» вимірювання оцінюється (як у [[психометрична надійність]]).

[[Планування експерименту]] або обсерваційне дослідження часто впливає на потужність. Наприклад, у двовибірковій ситуації з урахуванням загального розміру вибірки ''n'', оптимальним є рівне число спостережень з двох груп, що порівнюються (до тих пір, як відмінності у двох популяціях стануть однакові). У регресійному аналізі та [[Дисперсійний аналіз|дисперсійному аналізі]], існує безліч теорій і практичних стратегій, спрямованих на покращення потужності, заснованих на оптимальній установці значень незалежних змінних моделей.

== Тлумачення ==
Хоча немає ніяких формальних стандартів для потужності (іноді називають як π), більшість дослідників оцінюють потужність своїх тестів з допомогою π=0.80 у ролі стандарту для адекватності. Цей договір передбачає чотири до одного співвідношенню бета-ризику і альфа-ризику. (β-ймовірність помилки другого роду; α-це ймовірність помилки першого роду, і 0,2 і 0,05 звичайні значення для β і α). Однак, будуть часи, коли цей 4-до-1 зважування буде недоцільним. У медицині, наприклад, тести часто розробляються таким чином, що ніяких помилкових негативів (Тип II помилки) не буде зроблено. Але це неминуче підвищує ризик отримання хибно-позитивного результату (помилка першого роду). Пояснення таке, що краще скажи здоровому пацієнту «ми, можливо, знайшли щось-давайте перевіримо далі», ніж сказати хворому пацієнту «все добре».<ref>{{Cite book|title=The Essential Guide to Effect Sizes: An Introduction to Statistical Power, Meta-Analysis and the Interpretation of Research Results|last=Ellis, Paul D.|year=2010|publisher=Cambridge University Press|location=United Kingdom}}</ref>

Аналіз потужності підходить для випадків, коли проблема пов'язана з правильним відсортуванням, чи ні, з нульовою гіпотезою. У багатьох випадках питання йде про те, як визначити, якщо є, чи немає різниці, а скоріше з отриманням більш точної [[Теорія оцінювання|оцінки]] величини ефекту населення. Наприклад, якщо ми очікували [[Коефіцієнт кореляції Пірсона|кореляції]] між інтелектом населення і [[Продуктивність праці|продуктивності праці]] близько 0,50,розмір вибірки в 20 дасть нам приблизно 80&nbsp;% потужності (альфа = 0,05, два хвоста), щоб відхилити нульову гіпотезу про нульовий кореляції. Однак, у цьому дослідженні нас, напевно, більше цікавить, чи є кореляція 0.30 або 0.50 або 0.60. В зв'язку з цим нам буде потрібно набагато більший розмір вибірки, щоб зменшити [[довірчий інтервал]] нашої оцінки в діапазоні, який є прийнятним для наших цілей. Методи, аналогічних тим, які застосовуються в традиційному силовому аналізі можна використовувати для визначення обсягу вибірки, необхідної для ширини довірчого інтервалу повинна бути менше заданого значення.

Багато статистичних аналізів передбачає оцінку кількох невідомих величин. У простих випадках, все, крім однієї з цих величин є перешкодою параметра. В цих умовах, тільки відповідні потужності відносяться до єдиної величини, яка буде проходити формальний статистичний висновок. У деяких випадках, особливо якщо цілей більше «дослідницьких» може бути цілий ряд величин, що представляють інтерес для аналізу. Наприклад, у множині [[Регресійний аналіз|регресійного аналізу]] ми можемо включити кілька ковариантів, які представляють потенційний інтерес. У таких ситуаціях, як ця, де знаходиться кілька гіпотез на розгляді, вона є загальною, що повноваження, пов'язані з різними гіпотезами відрізняються. Наприклад, у множинному регресійному аналізі, потужність для виявлення впливу даного розміру має відношення до дисперсії коваріантів. Оскільки різні коваріанти будуть мати різні відхилення, їх повноваження будуть відрізнятися також.

Будь-який статистичний аналіз, що включає кілька гіпотез піддається інфляції типу I помилок, якщо відповідні заходи не будуть прийняті. Такі заходи, як правило, пов'язані із застосуванням більш високого порогу жорсткості, щоб відхилити гіпотезу для того, щоб компенсувати робляться множинні порівняння (''наприклад,'' як метод Бонферроні). У цій ситуації потужність аналізу повинна відображати багаторазові тестування. Так, наприклад, дане дослідження може бути потужним, щоб виявити певний ефект розміру, коли тільки одне випробування, але той же ефект розміру може мати набагато більш низьку потужність, якщо виконати кілька випробувань.

Важливо також враховувати статистичну потужність тестової гіпотези при інтерпретації його результатів. Тестова потужність&nbsp;— це ймовірність правильно відкинути нульову гіпотезу, коли вона помилкова; Тестова потужність знаходиться під впливом вибору рівня значущості для випробування, розмір ефекту вимірюється, а кількість даних наявні. Гіпотензивний тест не може відхилити нульову гіпотезу, наприклад, якщо справжня різниця між двома популяціями в порівнянні з [[T-критерій Стьюдента|т-тестом]], але ефект невеликого розміру вибірки занадто малий, щоб відрізнити ефект від випадковості.<ref>{{Cite book|title=The Essential Guide to Effect Sizes: Statistical Power, Meta-Analysis, and the Interpretation of Research Results|url=https://archive.org/details/essentialguideto00elli|last=Ellis|first=Paul|year=2010|publisher=Cambridge University Press|page=[https://archive.org/details/essentialguideto00elli/page/n71 52]|isbn=978-0521142465}}</ref> Численні [[Клінічне випробування|клінічні випробування]], наприклад, мають низьку статистичну потужність для виявлення відмінностей в побічних ефектах лікування, оскільки такі ефекти можуть бути рідкими, а кількість потерпілих пацієнтів невеликою.<ref>{{Cite journal | last1 = Tsang | first1 = R. | last2 = Colley | first2 = L. | last3 = Lynd | first3 = L. D. | doi = 10.1016/j.jclinepi.2008.08.005 | title = Inadequate statistical power to detect clinically significant differences in adverse event rates in randomized controlled trials | journal = Journal of Clinical Epidemiology | volume = 62 | issue = 6 | pages = 609–616 | year = 2009 | pmid =  19013761| pmc = }}</ref>

== ''Апріорний'' і ''постфактумний'' аналіз ==
Аналіз потужності може бути зроблено раніше (''апріорі'' або перспективному аналізу потужності) або після (''постфактум'' або [[Ретроспективний аналіз|ретроспективного аналізу]] даних). ''Апріорі'' силовий аналіз проводиться до початку дослідження, і, як правило, використовується в [[Визначення розмірів вибірки|оцінці вибірки достатнього розміру]] для досягнення адекватної потужності. ''Ретроспективний'' аналіз потужності проводиться після того, як дослідження було завершено, і використовує отримані вибірки і розміри ефекту, щоб визначити, що потужність була в дослідженні, припускаючи, що вплив розміру вибірки дорівнює впливу розміру популяції. У той час як програма перспективного аналізу в плануванні експерименту є загальноприйнятою, корисність ретроспективних методів є спірною.<ref>Thomas, L. (1997) [http://eprints.st-andrews.ac.uk/archive/00000417/01/ThomasCB1997.pdf Retrospective power analysis] {{Webarchive|url=https://web.archive.org/web/20070702085234/http://eprints.st-andrews.ac.uk/archive/00000417/01/ThomasCB1997.pdf |date=2 липня 2007 }}.</ref> Падає спокусі використовувати статистичний аналіз зібраних даних для оцінки потужності приведе до неінформативних і помилкових цінностей. Зокрема, було показано<ref name="HH1">Hoenig and Heisey (2001)The Abuse of Power ''{{Нп|The American Statistician|||}}'' 55(1):19-24 [http://dx.doi.org/10.1198/000313001300339897]</ref>, що ''після постфактумної потужності'' в найпростішому вигляді є [[один-до-одного]] функція досягнутого ''n''-значення. Це було розширено<ref name="HH1"/>, щоб показати, що вся ''постфактумна'' потужність аналізу хибує на «power approach paradox» (PAP), в якому дослідження з нульовим результатом, як вважають показує більше доказів того, що нульова гіпотеза не вірна, якщо ''р''-значення менше, оскільки уявна потужність, виявляє реальний ефект який буде вище. Насправді, менша величина р розуміється належним чином, щоб зробити нульову гіпотезу, менш імовірно наявна.{{Fact}}

== Додаток ==
Фінансові установи, і науково-дослідницькі контрольні групи нерідко звертаються з проханнями щоб дослідник виконав аналіз потужності, наприклад, щоб визначив мінімальну кількість піддослідних тварин, необхідних для проведення експерименту, щоб було інформативним. У [[Частотне висновування|частотній статистиці]], не потужні дослідження навряд чи дозволяють вибрати між гіпотезами на бажаному рівні значущості. У [[Баєсова статистика|Байєсівській статистиці]], перевірки гіпотез, які використовуються в класичному силовому аналізі не робляться. У Байєсівських рамках, його оновлення до переконання, використовуючи дані, отримані в даному дослідженні. В принципі, дослідження, яке буде вважатися недостатнім з точки зору перевірки гіпотез може використовуватися в таких процесах оновлення. Проте потужність, як і раніше, залишається корисною мірою того, наскільки даний експеримент розміру можна очікувати, щоб уточнити свої переконання. Дослідження малої потужності навряд чи призведе до більших змін у переконаннях.

== Приклад ==
Нижче наведено приклад, який показує, як обчислюється потужність для рандомізованого експерименту: припустимо, що мета експерименту&nbsp;— вивчення впливу лікування на деяку кількість, і порівняння об'єктів дослідження, вимірюючи кількість до і після лікування, аналіз даних за допомогою парного [[T-критерій Стьюдента|T-тесту]]. Нехай&nbsp;<math>\Alpha\iota</math> і <math>\Beta\iota</math> позначимо попередньою обробкою і пост-обробкою з цього і відповідно. Можливий ефект від лікування повинен бути видний у відмінностях <math>D_i=B_i-A_i</math>, яке передбачається незалежно розподіленим, все з тим же очікуваним значенням і дисперсією.

Ефект лікування може бути проаналізований за допомогою одностороннього T-тесту. Нульова гіпотеза щодо відсутності ефект буде той, що означає різницю і дорівнюватиме нулю, тобто <math>H_0: \mu_D=0</math>. У цьому випадку альтернативна гіпотеза стверджує, позитивний ефект, відповідний <math>H_{1}:\mu _{D}>0</math>. В [[Статистичний критерій|тестової статистиці]] є:
: <math>T_{n}={\frac  {{\bar  {D}}_{n}-0}{{\hat  {\sigma }}_{D}/{\sqrt  {n}}}}</math>
де ''n''&nbsp;— розмір вибірки,<math>{\bar  {D}}_{n}={\frac  {1}{n}}\sum _{{i=1}}^{n}D_{i}</math> і <math>{\hat  {\sigma }}_{D}/{\sqrt  {n}}</math> — це стандартна помилка. Розподіл тестової статистики при нульовій гіпотезі слід [[T-розподіл Стьюдента|T-розподілу Стьюдента]]. Крім того, припустимо, що нульова гіпотеза буде відхилена при [[Похибки першого і другого роду|рівні значущості]] в <math>\alpha=0.05</math>. Оскільки ''n'' є великим, можна наблизити до T-розподілу нормальним законом розподілу і обчислити критичне значення за допомогою [[квантильної функції]] <math>\Phi</math> нормального розподілу. Виходить, що нульова гіпотеза буде відкинута, якщо
: <math>T_{n}>1.64</math>
Тепер припустимо, що альтернативна гіпотеза істинна і <math>\mu_D=\theta</math>. Тоді потужність
: <math>
\begin{array}{ccl}
B(\theta)&=&P(T_n > 1.64|\mu_D=\theta) \\
&=&P(\frac{\bar{D}_n-0}{\hat{\sigma}_D/\sqrt{n} }> 1.64|\mu_D=\theta) \\
 &=&P\left(\frac{\bar{D}_n-\theta+\theta}{\hat{\sigma}_D/\sqrt{n}} > 1.64\right|\mu_D=\theta)\\
&=& P\left(\frac{\bar{D}_n-\theta}{\hat{\sigma}_D/\sqrt{n}} > 1.64-\frac{\theta}{\hat{\sigma}_D/\sqrt{n}}\right|\mu_D=\theta)\\
&=& 1- P\left(\frac{\bar{D}_n-\theta}{\hat{\sigma}_D/\sqrt{n}} < 1.64-\frac{\theta}{\hat{\sigma}_D/\sqrt{n}}\right|\mu_D=\theta)\\
\end{array}
</math>
Для великих n, <math>T_n</math> приблизно слідує стандартному [[нормальний розподіл|нормальному розподілу]], коли альтернативна гіпотеза, приблизна потужність можна розрахувати як
: <math>B(\theta )\approx 1-\Phi (1.64-{\frac  {\theta }{{\hat  {\sigma }}_{D}/{\sqrt  {n}}}})</math>
Згідно з цією формулою, сила зростає із збільшенням значення параметра <math>\theta</math>. Для конкретного значення <math>\theta</math> вища потужність може бути досягнута шляхом збільшення обсягу вибірки ''n''.

Неможливо гарантувати досить велику потужність для всіх значень<math>\theta</math>, бо <math>\theta</math> може бути дуже близько до 0. Мінімальному значенню потужності (межею) дорівнює розміру тесту, <math>\alpha</math>, у цьому прикладі 0.05. Втім, не важливо розрізняти <math>\theta=0</math> і невеликі позитивні значення. Бажано мати достатню потужність, скажімо, не менше 0,90, для виявлення значень <math>\theta >1</math> необхідний розмір вибірки можна розрахувати приблизно:
: <math>B(1)\approx 1-\Phi (1.64-{\sqrt  {n}}/{\hat  {\sigma }}_{D})>0{.}90</math>
з чого слідує, що
: <math>\Phi(1.64-\sqrt{n}/\hat{\sigma}_D) <0.10</math>
Звідси
: <math>{\displaystyle {\sqrt {n}}/{\hat {\sigma }}_{D}>1.64-z_{0.10}=1.64+1.28\approx 2.92}</math>
або
: <math>{\displaystyle n>8.56\,{\hat {\sigma }}_{D}^{2},}</math>
де <math>z_{0.10}</math> стандартна нормальна квантильна; див. [[пробіт-регресія]] для пояснення взаємозв'язку між <math>\Phi</math> і Z-значеннями.

== Розширення ==

=== Байєсовська потужність ===
В [[Частотне висновування|частотному висновувані]] параметри, як передбачається, мають певне значення, що малоймовірно, щоб бути правдою. Цю проблему можна вирішити, припустивши, що параметр має розподіл. Результуюча сила іноді називають Байєсівської потужністю, яка зазвичай використовується в [[Клінічне випробування|клінічних випробуваннях]].

=== Прогностична ймовірність успіху ===
Як у [[Частотне висновування|частотній]] потужність так і в Байєсовській потужності використовують статистичну значимість як критерій успіху. Однак статистичної значимості часто не досить, щоб визначити успіх. Для вирішення цієї проблеми потужності концепція може бути розширена до концепції прогностичної ймовірності успіху (ПЙУ). Успіх критерії для ПЙУ не обмежується статистичною значущістю і широко використовується в [[Клінічне випробування|клінічних випробуваннях]].

== Програмне забезпечення для потужності та розрахунок обсягу вибірки ==
Численні вільні та/або відкриті програми доступні для виконання потужності і розрахунку розміру вибірки. Вони включають
* G*Power (http://www.gpower.hhu.de/ {{Webarchive|url=https://web.archive.org/web/20161210020657/http://gpower.hhu.de/ |date=10 грудня 2016 }})
* [http://powerandsamplesize.com powerandsamplesize.com] {{Webarchive|url=https://web.archive.org/web/20220327184517/https://powerandsamplesize.com/ |date=27 березня 2022 }} Вільне і відкрите джерело онлайн калькулятора
* {{нп|PS Power and Sample Size}}
* [[R (мова програмування)|R]] package [https://cran.r-project.org/web/packages/pwr/index.html pwr] {{Webarchive|url=https://web.archive.org/web/20161128071412/https://cran.r-project.org/web/packages/pwr/index.html |date=28 листопада 2016 }}
* [http://homepage.stat.uiowa.edu/~rlenth/Power/ Russ Lenth's power and sample-size page] {{Webarchive|url=https://web.archive.org/web/20161121185026/http://homepage.stat.uiowa.edu/~rlenth/Power/ |date=21 листопада 2016 }}
* WebPower Free online statistical power analysis (http://webpower.psychstat.org)
* SampSize app for Android and [[iOS]] iPhone and iPad (https://www.epigenesys.org.uk/portfolio/sampsize/ {{Webarchive|url=https://web.archive.org/web/20170313074221/https://www.epigenesys.org.uk/portfolio/sampsize |date=13 березня 2017 }})

== Див. також ==
* [[Визначення розмірів вибірки|Розмір вибірки]]

== Примітки ==
{{reflist}}

== Посилання ==
* {{Cite book|title=The Cambridge Dictionary of Statistics|url=https://archive.org/details/nlsiu.310.03.eve.18720|last=Everitt|first=Brian S.|year=2002|publisher=Cambridge University Press|isbn=0-521-81099-X|ref=harv}}
* {{Cite book|title=Statistical Power Analysis for the Behavioral Sciences|url=https://archive.org/details/statisticalpower0000cohe_j0l3|last=Cohen|first=J.|authorlink1=Jacob Cohen (statistician)|year=1988|edition=2-ге|isbn=0-8058-0283-5}}
* {{Cite book|title=Applied Power Analysis for the Behavioral Science|last=Aberson|first=C. L.|year=2010|isbn=1-84872-835-2}}
* {{YouTube|UFhJefdVCjE|StatQuest: P-value pitfalls and power calculations|id=UFhJefdVCjE|title=StatQuest: P-value pitfalls and power calculations}}

[[Категорія:Перевірка статистичних гіпотез|Потужність]]