Как рассчитать с большей точностью в R для биномиального распределения

Я пытаюсь рассчитать вероятность угадывания 70 или более правильных ответов в тесте из 100 вопросов с 4 вариантами ответов (то есть вероятность правильного угадывания 25% для каждого ответа), чтобы проиллюстрировать, что даже троечник на самом деле справляется довольно хорошо по сравнению с случайное угадывание.

Мое понимание R предполагает, что ответом должно быть 1 — pbinom (69 100, 0,25), но R выводит «0». Я пробовал варианты (цифры = 20) и пакет «mpfr», но безуспешно.

Как лучше вычислить ответ?

Спасибо!


r precision binomial-coefficients
person Ronald Miller    schedule 07.03.2017    source источник
comment
Вы продолжаете получать ноль, потому что all.equal(pbinom(69, 100, 0.25), 1) возвращает TRUE, а это означает, что, насколько может судить R, pbinom(69, 100, 0.25) на самом деле равно 1,0. Это означает, что 1 - pbinom(69, 100, 0.25) меньше машинной точности 1.490116e-08. Используя вопрос из предыдущего комментария, pbinom(70, 100, 0.25) возвращает 6.115655e-22, что определенно меньше машинного допуска.   -  person Benjamin    schedule 08.03.2017
comment
Огромное спасибо! Кажется, это сработало и прояснилось! a=pbinom(69,100,0.25,lower.tail = FALSE) с mpfr(a,1024) идеально!   -  person Ronald Miller    schedule 08.03.2017