Вероятность и статистика · Статистическая мощность

Статистическая мощность

Мощность критерия — вероятность того, что он отвергнет нулевую гипотезу H0, когда верна определённая альтернатива. Здесь машина должна фасовать в среднем по 500 г. Перемещайте её истинное среднее и смотрите, как часто исследование с числом пакетов n обнаруживало бы отклонение.

Мощность на этом рисунке

Мощность критерия — вероятность того, что он отвергнет нулевую гипотезу H0H_0, когда верна определённая альтернатива. Здесь H0H_0 утверждает, что машина фасует в среднем по μ=500\mu = 500 г. Веса пакетов независимы и нормально распределены с σ=10\sigma = 10 г. В исследовании взвешивают nn пакетов и отвергают H0H_0 на уровне α\alpha, если среднее xˉ\bar x оказывается за критической границей. Мощность при истинном среднем μ\mu — площадь под истинной кривой за этой границей. Остаток, β\beta, — вероятность пропустить отклонение, то есть совершить ошибку второго рода. При μ=500\mu = 500 отклонения нет, и мощность равна самому α\alpha.

От чего зависит мощность? Для проверки H1 ⁣:μ<500H_1\!: \mu < 500 при недовесе δ=500−μ\delta = 500 - \mu она равна Φ(δn/σ−z1−α)\Phi(\delta\sqrt{n}/\sigma - z_{1-\alpha}). Большее отклонение или больше пакетов повышают её; более строгое α\alpha отодвигает критическую границу от 500 и снижает мощность. Поскольку nn входит только через n\sqrt{n}, число пакетов для заданной мощности растёт как 1/δ21/\delta^2. Направление не менее важно: если машина насыпает лишнее, мощность проверки μ<500\mu < 500 стремится к 0, а увеличение числа пакетов снижает её ещё сильнее. Двусторонний критерий делит α\alpha между двумя хвостами и выявляет отклонения в обе стороны. Но на той стороне, которую проверяет односторонний критерий, его мощность ниже, чем у одностороннего.

Почему незначимый результат не означает отсутствия различий? Критерий с низкой мощностью обычно пропускает реальное отклонение. При 25 пакетах и α=0.05\alpha = 0.05 односторонний критерий обнаруживает недовес 2 г лишь примерно в четверти случаев, поэтому большинство исследований этой машины дают незначимый результат. Информативность такого результата зависит от мощности для отклонения, которое важно обнаружить. Её рассчитывают до просмотра данных: выбирают наименьшее существенное отклонение, а затем nn, обеспечивающее достаточную мощность, часто 80%.

Почему при низкой мощности значимые оценки в среднем преувеличивают отклонение? Возьмём машину, которая действительно недосыпает, и проверим её на недовес. За критическую границу попадают только исследования, где среднее случайно оказалось намного ниже 500. Поэтому значимые исследования — отобранная группа: их средняя оценка равна среднему той части истинной кривой, что лежит за границей, и соответствует большему недовесу, чем реальный. Сама граница задаёт минимум: при 25 пакетах и α=0.05\alpha = 0.05 для значимости оценка недовеса должна превышать примерно 3.29 г, каким бы ни был истинный недовес. При фиксированных nn и α\alpha чем меньше истинный недовес, тем ниже мощность и тем больше отношение средней значимой оценки к истинному недовесу. Это ошибка типа M (величины). При двусторонней проверке некоторые значимые результаты даже указывают неверное направление — ошибка типа S (знака). Совокупность работ, сообщающих только значимые результаты, несёт это смещение.

На что обратить внимание

  • Обе верхние кривые рисуются одинаковой высоты при любом nn. Каждая закрашенная область — доля площади под своей кривой. Уменьшение разброса показывает сужение, а не высота.
  • Мощность — не вероятность того, что машина недосыпает. Её рассчитывают, считая истинным средним заданное на странице μ\mu: P(отвергнуть H0∣μ)P(\text{отвергнуть } H_0 \mid \mu).

Управление

  • Перетаскивайте синюю кривую или её ручку, чтобы сдвинуть истинное среднее, либо нажмите в нижнем слое, чтобы перенести его туда. Перетаскивайте квадратную ручку критической границы, чтобы задать α\alpha.
  • Строка «Мощность 80%» устанавливает nn равным наименьшему размеру выборки, при котором мощность достигает 80% при текущих среднем, α\alpha и направлении.

Связанные темы