質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
Python 2.7

Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

Q&A

解決済

1回答

265閲覧

Pythonにおける試行回数、成功数が大きい二項分布の算出方法

siita

総合スコア8

Python 2.7

Python 2.7は2.xシリーズでは最後のメジャーバージョンです。Python3.1にある機能の多くが含まれています。

1グッド

0クリップ

投稿2018/05/06 15:21

前提・実現したいこと

Pythonで勉強がてらプログラムを組んでいます。
試行回数、成功数が比較的大きな二項分布の計算をする機能を実現したいと考えています。

実現方法について何かご存知の方がいらっしゃいましたら、ご回答よろしくお願いします。

発生している問題・エラーメッセージ

ざっと関数を作成してみたのですが、計算結果が大きすぎるからだとは思うのですが、結果が inf になってしまっています。
また、その際には以下の警告が表示されます。

RuntimeWarning: invalid value encountered in double_scalars

該当のソースコード

Python

1[呼び出し先] 2import scipy.misc as scm 3 4def binomialDistribution(x, n, p): 5 return scm.comb(n,x) * (p**(x)) * ((1.0-p)**(n-x)) 6 7[呼び出し元] 8binomialDistribution(1000, 10000, 0.2) //

試したこと

数値を扱う型を大きくしたり、組み合わせ関数を色々調べて自作してみたりしたのですが、どうもうまくいきません。

magichan👍を押しています

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

scipy.stats.binom.pmfとか使えませんか。

scipy.stats.binom — SciPy v0.14.0 Reference Guide

python

1>>> from scipy import stats 2>>> stats.binom.pmf(1000, 10000, 0.2) 36.0401240113519745e-162

追記

「勉強でやってるんだから、ライブラリ丸投げじゃ意味ないんですよね」という場合に備えて追記。

一回logを取ってlogの足し算にしてからexpしてみる。

python

1import scipy.misc as scm 2from math import log, exp 3 4def binomialDistribution2(x, n, p): 5 t1 = log(scm.comb(n,x, exact=True)) 6 t2 = x*log(p) 7 t3 = (n-x)*log(1.0-p) 8 return exp(t1+t2+t3) 9 10print(binomialDistribution2(1000, 10000, 0.2)) 11# => 6.040124011291546e-162

上のscipyの結果とは若干ずれてるみたいですが(浮動小数点計算かつ近似計算なので仕方ない)、一応やればできますということで。

投稿2018/05/06 15:36

編集2018/05/06 21:03
hayataka2049

総合スコア30933

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

siita

2018/05/07 12:55

ご回答ありがとうございます。 自分なりには探してみたのですが、ずばりそのものがあったんですね…。 完全に文系人間なので、logとexpもろくに知りませんでした。 教えていただいた方法で実装していきたいと思います。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問