конспект
3,5пp=1/2p≠1/2
Остальная часть этого поста объясняет базовую модель и показывает, как выполнять вычисления. Он предоставляет рабочий R
код для их выполнения. Подробное изложение основной теории проверки гипотез приведено в моем ответе на вопрос: «Что означают значения p и значения t в статистических тестах?». ,
Статистическая модель
3.51/2=50%3.53.52503.5(250,1/2)k
3.53.51/2
Нахождение подходящего теста
k1/225012512502501253.5
ααk
3.5α/2kα/2kk
Технически, есть два распространенных способа выполнить вычисление: вычислить биномиальные вероятности или приблизить их с помощью нормального распределения.
Расчет с биномиальными вероятностями
Используйте функцию процента (квантиль). В R
, например, это называется qbinom
и будет вызываться как
alpha <- 0.05 # Test size
c(qbinom(alpha/2, 250, 1/2)-1, qbinom(1-alpha/2, 250, 1/2)+1)
α=0.05
109 141
k0109k141250R
k
pbinom(109, 250, 1/2) + (1-pbinom(141-1, 250, 1/2))
0.0497αα
Расчет с нормальным приближением
(250,1/2)250×1/2=125250×1/2×(1−1/2)=250/4250/4−−−−−√≈7.9α/2=0.05/2−1.95996R
qnorm(alpha/2)
0.05/2+1.95996k1.95996125125±7.9×1.96≈109.5,140.5
250*1/2 + sqrt(250*1/2*(1-1/2)) * qnorm(alpha/2) * c(1,-1)
k109141p1/201α
Этот тест, поскольку он не предполагает ничего о населении (за исключением того, что он не имеет большой вероятности, ориентированной прямо на его медиану), не так силен, как другие тесты, которые делают конкретные предположения о населении. Если тест, тем не менее, отклоняет ноль, нет необходимости беспокоиться о нехватке мощности. В противном случае вы должны сделать несколько деликатных компромиссов между тем, что вы готовы принять, и тем, что вы можете сделать о населении.