ポアソン分布をpythonで試してみた

Question

ポアソン分布を試してみたところエラーが出ました。

from scipy.stats import poisson


#平均は10です
mu=10
#平均と分散を計算できる
mean,var=poisson.pmf(7,mu)

#確率質量関数を使って特定の関数を計算できる
odds_seven=poisson.pmf(7,mu)

print('ピーク時に7人の確率は{:0.2f}%'.format(odds_seven*100))
```error文
---------------------------------------------------------------------------
ModuleNotFoundError                       Traceback (most recent call last)
<ipython-input-16-9edb5d2c57b2> in <module>
      1 #scipyにはpoissonというモジュールがあらかじめある
      2 
----> 3 from scipy.stats import poisson
      4 
      5 

~\Anaconda3\lib\site-packages\scipy\__init__.py in <module>
     60 __all__ = ['test']
     61 
---> 62 from numpy import show_config as show_numpy_config
     63 if show_numpy_config is None:
     64     raise ImportError(

~\Anaconda3\lib\site-packages\numpy\__init__.py in <module>
    140     from . import _distributor_init
    141 
--> 142     from . import core
    143     from .core import *
    144     from . import compat

~\Anaconda3\lib\site-packages\numpy\core\__init__.py in <module>
     33 
     34 from . import umath
---> 35 from . import _internal  # for freeze programs
     36 from . import numerictypes as nt
     37 multiarray.set_typeDict(nt.sctypeDict)

~\Anaconda3\lib\site-packages\numpy\core\_internal.py in <module>
     11 
     12 from numpy.compat import unicode
---> 13 from numpy.core.overrides import set_module
     14 from .multiarray import dtype, array, ndarray
     15 try:

~\Anaconda3\lib\site-packages\numpy\core\overrides.py in <module>
      4 import os
      5 
----> 6 from numpy.core._multiarray_umath import (
      7     add_docstring, implement_array_function, _get_implementing_args)
      8 from numpy.compat._inspect import getargspec

ModuleNotFoundError: No module named 'numpy.core._multiarray_umath'
```

Comments

[tachikoma]

Anaconda promptで
conda update conda
を実行してみるとどうなるでしょう。

Answer 1

pmf は確率質量関数なので、平均、分散は返しません。
API を確認してくだささい。

scipy.stats.poisson — SciPy v1.2.1 Reference Guide

- mean,var=poisson.pmf(7,mu)
+ prob =poisson.pmf(7,mu)

import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import poisson

mu = 10
x = np.arange(poisson.ppf(0.01, mu), poisson.ppf(0.99, mu))
y = poisson.pmf(x, mu)

fig, ax = plt.subplots(facecolor='w')
ax.stem(x, y)
plt.show()

prob = poisson.pmf(7, mu)

print('ピーク時に7人の確率は {:.2%}'.format(prob))  # ピーク時に7人の確率は 9.01%

Comments

[tiitoi]

質問をよく見てなかったのですが、そもそも import が通ってないようですね。
pip install --upgrade scipy として再インストールしてみてはどうでしょうか

[tachikoma]

anaconda環境のようなので、まずはconda update scipyを試したほうがいいと思う。

[tiitoi]

コメントありがとうございます。
たしかに Anaconda 環境でしたら、おっしゃる通り、conda コマンドがいいと思います。