Как вычислить биномиальную кумулятивную функцию плотности с питоном - PullRequest
0 голосов
/ 26 октября 2019

У меня есть следующее биномиальное распределение:

В прошлом году число новых зданий в Общественном совете 12 и Общественном совете 11 в Бронксе было 347. Из этих 347 107 имели место в Общественном совете 12.

Если мы случайным образом выберем 70 из 347 новых зданий, распределение вероятности будет:

X ~ B (70, 107/347)

Если я хочу знатьвероятность того, что из этих 70 случайно выбранных новостроек, 20, будут иметь место в совете сообщества 12, я бы поступил следующим образом, используя scipy.stats:

binom.pmf (20, 70, 0.3083573487)                                                                                                                                                                  
0.09646726155763652

Если я хочу узнать вероятность того, что из этих70 случайно выбранных зданий, меньше или равных 20 зданиям, имели место в Правлении Сообщества 12, я бы сделал это следующим образом, используя scipy.stats:

binom.cdf (20, 70, 0.3083573487)                                                                                                                                                                  
0.39547679625297977

Если я хочу узнать вероятность того, что из этих 70 случайнов сообществе 12 было выбрано только менее 20 зданий. Я бы поступил следующим образом, используя scipy.stats:

binom.cdf (20, 70, 0.3083573487, loc = 1)                                                                                                                                                         
0.2990095346953431

У меня возникли проблемы:выяснение вероятности того, что из этих 70 случайно выбранных новых зданий 20 или более будут происходить в сообществе 12. Я знаю, что результаты должны быть 0.60452320375, но я не могу найти команду scipy.stats, чтобы получить эти результаты.

Любая помощь, с благодарностью.

Спасибо.

1 Ответ

1 голос
/ 26 октября 2019

Поскольку cdf(x) распределения вероятностей является интегралом от отрицательной бесконечности до x, интеграл от x до положительной бесконечности равен 1-cdf(x). Так что для вашей проблемы это будет просто:

probabilityGreaterThan20inCommunity12 = 1 - binom.cdf (20, 70, 107./347)
...