Сколько раз я должен бросить кубик, чтобы уверенно оценить его справедливость?

(Заранее извиняюсь за использование языка мирян, а не статистического языка.)

Если я хочу измерить шансы накатывания каждой стороны конкретного физического шестигранного кристалла с точностью примерно +/- 2% с достаточной уверенностью в достоверности, сколько потребуется образцовых роликов?

т.е. сколько раз мне нужно было бы бросить кубик, считая каждый результат, чтобы быть на 98% уверенным, что шансы на то, что он бросит каждую сторону, находятся в пределах 14,6% - 18,7%? (Или некоторые аналогичные критерии, при которых можно было бы убедиться на 98%, что умирание справедливо с точностью до 2%.)

(Это реальная проблема для игр-симуляторов, в которых используются кости, и они хотят убедиться, что определенные конструкции костей приемлемо близки к 1/6 вероятности броска каждого числа. Есть утверждения, что многие обычные конструкции костей были измерены, бросая 29% 1 к бросать несколько таких кубиков по 1000 раз каждый.)

probability inference pdf dice Dronz
источник

Это гораздо сложнее, чем найти доверительный интервал для бинома, так как вы хотите контролировать все вероятности. Взгляните на статью Сюин Вана об одновременных доверительных интервалах для многочленных распределений ( Журнал многомерного анализа 2008, 99, 5, 896-911). Вы можете найти некоторый код в этом посте , который также дает краткое описание некоторых работ, которые были проделаны над этим.

idnavid

Обратите внимание, что если вы просто заинтересованы в том, чтобы проверять, выпадают ли единицы достаточно много времени, это значительно упрощает вопрос.

Деннис Джаэруддин

Важно отметить, что «доверительный интервал» не дает вам «процентильной вероятности правильности». Я подозреваю, что вы используете очень разумное общепринятое употребление термина «уверен на 98%», но вы должны знать всякий раз, когда кто-то упоминает «доверительный интервал», который совсем не совпадает с вероятностью 98%: link.springer.com/ article / 10.3758% 2Fs13423-013-0572-3

BrianH

@BrianH Спасибо! Я имел в виду не только разговорное выражение, но и пытаюсь измерить определенность, подразумеваемую тестом. Мне кажется, что точно так же, как и то, что имеет смысл сказать, что я рассчитываю на то, что бросок кубика будет рассчитан в процентах от времени, будет аналогичный (но более сложный) расчет вероятности того, что я буду бросать результаты в пределах определенная погрешность в броске n раз, что, как мне кажется, я понимаю в ответе Сямои (и последующем комментарии). Да?

Дронз

@ Dronz Чтобы быть справедливым, это одна из тех вещей, которые вы действительно думаете, было бы более прямолинейным, чем оказывается на самом деле. Чертовски сложно, на самом деле. Ниже приводятся некоторые ключевые вопросы, связанные с другими вопросами, чтобы помочь вам понять, почему не существует невероятно простого ответа: Frequentist math.stackexchange.com/questions/1578932/… Bayesian math.stackexchange.com/questions/1584833/… и веселье: rpg.stackexchange.com/questions/70802/…

BrianH

Ответы:

TL; DR: если $p$ = 1/6, и вы хотите знать, насколько велико $n$ должно быть на 98%, чтобы убедиться, что игра в кости справедлива (с точностью до 2%), $n$ должно быть не менее $n$ ≥ 766 .

Пусть $n$ будет числом рулонов, а $X$ - количеством рулонов, которые приземляются на определенной стороне. Тогда $X$ следует биномиальному (n, p) распределению, где $p$ - вероятность получения указанной стороны.

По центральной предельной теореме мы знаем, что

\sqrt{n} (X / n - p) \to N (0, p (1 - p))

$\sqrt{n} (X/n - p) \to N(0,p(1-p))$

поскольку $X/n$ - это выборочное среднее $n$ случайных величинБернулли $(p)$ . Следовательно, для больших $n$ доверительные интервалы для $p$ могут быть построены как

\frac{X}{n} \pm Z \sqrt{\frac{p (1 - p)}{n}}

$\frac{X}{n} \pm Z \sqrt{\frac{p(1-p)}{n}}$

Так как $p$ неизвестно, мы можем заменить его среднюю выборочную , и различными теоремами сходимости, мы знаем , что в результате доверительный интервал будет асимптотический действительным. Таким образом, мы получаем доверительные интервалы вида $\hat{p} = X/n$

\hat{p} \pm Z \sqrt{\frac{\hat{p} (1 - \hat{p})}{n}}

$\hat{p} \pm Z \sqrt{\frac{\hat{p}(1-\hat{p})}{n}}$

с . Я предполагаю, что вы знаете, что такое баллы. Например, если вы хотите 95% доверительный интервал, вы принимаете . Таким образом, для данного уровня достоверности мы имеем $\hat{p} = X/n$ $Z$ $Z=1.96$ $\alpha$

\hat{p} \pm Z_{α} \sqrt{\frac{\hat{p} (1 - \hat{p})}{n}}

$\hat{p} \pm Z_\alpha \sqrt{\frac{\hat{p}(1-\hat{p})}{n}}$

Теперь предположим, что вы хотите, чтобы этот доверительный интервал имел длину меньше, чем $C_\alpha$ , и хотите знать, насколько большой выборки нам нужен для этого случая. Ну, это равносильно тому, чтобы спросить, что удовлетворяет $n_\alpha$

Z_{α} \sqrt{\frac{\hat{p} (1 - \hat{p})}{n_{α}}} \leq \frac{C_{α}}{2}

$Z_\alpha \sqrt{\frac{\hat{p}(1-\hat{p})}{n_\alpha}} \leq \frac{C_\alpha}{2}$

Который затем решается получить

n_{α} \geq {(\frac{2 Z_{α}}{C_{α}})}^{2} \hat{p} (1 - \hat{p})

$n_\alpha \geq \left(\frac{2 Z_\alpha}{C_\alpha}\right)^2 \hat{p}(1-\hat{p})$

Таким образом , вилка в ваших значений $Z_\alpha$ ; , $C_\alpha$ ; , и по оценкам , чтобы получить оценку для $\hat{p}$ $n_\alpha$ . Обратите внимание, что, поскольку $p$ неизвестно, это только оценка, но асимптотически (при увеличении $n$ ) она должна быть точной.

Xiaomi
источник

Спасибо. Поскольку я десятилетиями не занимался математикой по типу колледжа, могу ли я вас побеспокоить, чтобы вы включили числа и фактически дали мне приблизительное число раз, которое мне нужно было бы бросить, как целое число?

Дронц

p = 1 / 6

$p = 1/6$

n

$n$

n

$n$

n \geq 766

$n \geq 766$

C_{α}

$C_\alpha$

Возможно, было бы интереснее взглянуть на многочленное распределение, поскольку теперь мы тестируем каждую сторону отдельно. Это не учитывает всю имеющуюся у нас информацию о проблеме. Для интуитивного объяснения смотрите stat.berkeley.edu/~stark/SticiGui/Text/chiSquare.htm

января

Я согласен с @Jan: Этот ответ не касается вопроса. Более того, его нельзя легко адаптировать для построения ответа, применяя его отдельно ко всем шести сторонам, потому что шесть тестов являются взаимозависимыми.

whuber

Это хороший ответ, но я полностью согласен с @Jan, whuber. Этот вопрос заслуживает ответа на основе статистики хи-квадрат и многочленного распределения.

Лукаш Град