RuntimeWarning: недопустимое значение встречается в мощности при построении распределения Пуассона

Question

Я хочу, чтобы запустить этот кусок кода

import numpy as np 
import matplotlib.pyplot as plt 
from scipy.misc import factorial 

def poisson(k, lamb): 
    return (lamb**k/factorial(k)) * np.exp(-lamb) 

x_plot = np.linspace(0, 10, 100) 
plt.plot(x_plot, poisson(x_plot, parameters), 'r-', lw=2) 
plt.show() 

, в котором parameters = 0.00213553

Но каждый раз, когда это появляется предупреждение:

RuntimeWarning: invalid value encountered in power 
return (lamb**k/factorial(k)) * np.exp(-lamb) 

и появляется только горизонтальная линия.

Я проверил подобные проблемы, но они действительно не помогли. Я также попытался отбросить возвращаемое значение функции до float128, которое тоже не работает. (с использованием Python 2.7.6 64-bit)

Answer 1

Я действительно не знаю, какой из factorial вы используете. Но когда я использую math.factorial, код не может работать. Затем я попытался с scipy.misc.factorial, который действительно может принимать numpy array в качестве входных данных. И когда я запускаю следующее, это не дает мне никакой ошибки или предупреждения.

import numpy as np 
import matplotlib.pyplot as plt 
import scipy.misc 

parameters = 0.00213553 

def poisson(k, lamb): 
    return (lamb**k/scipy.misc.factorial(k)) * np.exp(-lamb) 

x_plot = np.linspace(0, 10, 100) 
plt.plot(x_plot, poisson(x_plot, parameters), 'r-', lw=2) 
plt.show()