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

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

ただいまの
回答率

87.37%

ポアソン分布のグラフ表示がうまく表示されません

解決済

回答 1

投稿

  • 評価
  • クリップ 1
  • VIEW 888

score 24

ポアソン分布のグラフを実装していたところ平均値であるλが10の時だけグラフの点がうまく表示されません。具体的には下の二つの図が違うという事です。(最後の方が違う)
内部的な点の値はあっているのですが点がプロットされない原因が分かる方、教えて下さると光栄です。プログラムも載せておきます。

これが実行結果です。
イメージ説明
また次のグラフが本当のグラフです。
イメージ説明

ここからはコードを載せます。

#λは平均, kはその時の値。
#つまりポアソン分布は平均λの場合にkの値が出る確率
def poason(k, λ):
    if type(k)==int:
        sum=1
        for j in range(1,k+1):
            sum = sum*j
        return np.e**(-λ)*(λ**k)/sum
    else:
        a = np.zeros(len(k))         #aを作らずにkに代入しても良いがkはnp.arangeで作っているため整数しか入らない. そのためaを作成
        for i in range(len(k)):
            sum = 1
            for j in range(1,k[i]+1):
                sum = sum*j
            a[i]=np.e**(-1*λ)*(λ**k[i])/sum
        return a


import matplotlib.pyplot as plt
import numpy as np 

fig, ax = plt.subplots(figsize=(10,6))

k = np.arange(21)
p1 = poason(k, 1)  #λが1の場合のグラフ
p2 = poason(k, 4)  #λが4の場合のグラフ
p3 = poason(k,10)  #λが10の場合のグラフ

c1, c2, c3 = "blue", "green", "red"   #各プロットの色
l1, l2, l3 = "λ=1", "λ=4", "λ=10"  #各ラベル

ax.set_xlabel("k", fontsize=20)     #x軸ラベル
ax.set_ylabel("p(x=k)", fontsize=20)          #y軸ラベル
ax.set_title("Poisson distribution", fontsize=20)   #グラフタイトル

ax.grid()    #罫線(けいせん)。縦と横の細い線

ax.plot(k, p1, marker = "o", color = c1, label = l1)
ax.plot(k, p2, marker = "o", color = c2, label = l2)
ax.plot(k, p3, marker = "o", color = c3, label = l3)

#xの目盛りを作る(まあ作らなくても良いがfigsizeなど指定するには必要)目盛りの大きさのみ変えるだけなら一行で書ける。yは一行。
i=0
xtick = np.zeros(1)
while (i<21):
    if i%5==0:
        xtick = np.append(xtick,i)
    i +=1

ax.set_xticks(xtick)  #これでx軸の目盛りの個数と大きさは分かった。z残念ながらここにfontizeはない
ax.set_xticklabels(xtick, fontsize=15)    #第一引数はラベルであるので"A"などの名前でも良い。ここでfontsize指定できる

ax.tick_params(axis='y',labelsize=15)   #このように目盛りの大きさだけならこれで変えられる

ax.legend(loc=0, fontsize=20)   #凡例表示。線の説明。locは場所。0の場合は右上
fig.tight_layout()

plt.show()
  • 気になる質問をクリップする

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

+1

32ビット整数を使って計算しているためです。
poason関数の定義を以下に変更すれば、期待している結果になります。

def poason(k, λ):
    if type(k)==int:
        sum=1
        for j in range(1,k+1):
            sum = sum*j
        return np.e**(-λ)*(λ**k)/sum
    else:
        a = np.zeros(len(k))         #aを作らずにkに代入しても良いがkはnp.arangeで作っているため整数しか入らない. そのためaを作成
        for i in range(len(k)):
            sum = 1
            for j in range(1,k[i]+1):
                sum = sum*j
            a[i]=np.e**(-1*λ)*(λ**int(k[i]))/sum
        return a

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2021/01/28 18:57

    丁寧な解説ありがとうございます。たしかに32bit整数のまま計算すると小数点が多い場合に格納できないですね。とても勉強になりました。ありがとうございました。

    キャンセル

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

  • ただいまの回答率 87.37%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る