График плотности распределения Пуассона
Мне нужно нарисовать график плотности распределения Пуассона и гистограмму выборки, но возникает заданная ошибка, которую я понятия не имею, как решить.
import numpy as np
import matplotlib.pyplot as plt
lam = 3 # Параметр распределения Пуассона (лямбда)
x = np.arange(0, 20)
# Вычисление значений плотности распределения Пуассона
y = np.exp(-lam) * (lam ** x) / np.math.factorial(x)
# Построение графика
plt.plot(x, y, 'bo', ms=8) # 'bo' - голубые круглые маркеры
plt.xlabel("Значение СВ")
plt.ylabel("Плотность")
plt.title("Распределение Пуассона")
plt.show()
Я пробовал заменить x на np.array(x), но это ни к чему не привело
Почему в Python есть оператор "pass"?
Оператор pass в Python - это простая концепция, которую могут быстро освоить даже новички без опыта программирования.
Некоторые методы, о которых вы не знали, что они существуют в Python
Python - самый известный и самый простой в изучении язык в наши дни. Имея широкий спектр применения в области машинного обучения, Data Science,...
Основы Python Часть I
Вы когда-нибудь задумывались, почему в программах на Python вы видите приведенный ниже код?
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
Алиса и Боб имеют неориентированный граф из n узлов и трех типов ребер:
Оптимизация кода с помощью тернарного оператора Python
И последнее, что мы хотели бы показать вам, прежде чем двигаться дальше, это
Советы по эффективной веб-разработке с помощью Python
Как веб-разработчик, Python может стать мощным инструментом для создания эффективных и масштабируемых веб-приложений.
1
1
68
3
Перейти к ответу
Данный вопрос помечен как решенный
Ответы 3
Ответ принят как подходящий
Вместо
y = np.exp(-lam) * (lam ** x) / np.math.factorial(x)
вам следует перебирать значения x, поэтому попробуйте следующее:
y = [np.exp(-lam) * (lam ** i) / np.math.factorial(i) for i in x]
Остальное то же самое.
Добавляя к ответу Угура, вам нужно векторизовать функцию факториала (очень необычно внутри пакета numpy). Также для барного сюжета необходимо использовать plt.bar:
import matplotlib.pylab as plt
import numpy as np
lam = 3 # Параметр распределения Пуассона (лямбда)
x = np.arange(0, 20)
# Вычисление значений плотности распределения Пуассона
y = np.exp(-lam) * (lam ** x) / np.vectorize(np.math.factorial)(x)
# Построение графика
plt.bar(x, y) # 'bo' - голубые круглые маркеры
plt.xlabel("Значение СВ")
plt.ylabel("Плотность")
plt.title("Распределение Пуассона")
plt.show()
Согласно этот ответ , np.math.factorial и scipy.math.factorial являются просто псевдонимами math.factorial, который не работает с массивами numpy. Вместо этого используйте scipy.special.factorial, который может работать с массивами numpy.
from scipy.special import factorial
x = np.arange(0, 20)
# other code here
y = np.exp(-lam)*(lam**x)/factorial(x)
# plotting code here
Другие вопросы по теме
Как создать фрейм данных, не получая пустых значений?
Каков алгоритм np.linalg.solve?
Существует ли функция numpy, похожая на np.isin, которая допускает допуск, а не требует, чтобы значения были идентичными?
Дан матричный мат размером 3 х 3. Найдите максимальные числа в каждой строке, скажем N1, N2 и N3
Произведение массивов (einsum) трехмерных массивов, содержащее только -1 или +1
Разделение несовпадающих строк на несколько строк в фрейме данных Pandas
Как найти расстояние до следующего значения, отличного от NaN, в массиве numpy
NumPy — кумулятивное скалярное произведение со специальным выбором элементов, без циклов
Реализация алгоритма персептрона, но неэффективная при его запуске
Почему .reshape(a, b) != .reshape(b, a).T?
Похожие вопросы
Ни один элемент с именем «replace» в сообщении пространства имен «std::ranges» в действиях github, даже с -std=c++20
Получение ошибки драйвера ODBC 17 при подключении к Azure SQL
Есть ли способ отфильтровать многоиндексированный DataFrame pandas с помощью dict?
Преобразовать два входных файла в определенные строковые правила с помощью Python
Pytest должен собирать файлы Python только с определенным именем файла
Правильный шаблон для выбора файлов с помощью glob
Как создать фрейм данных, не получая пустых значений?
Функция, которая принимает объект в качестве параметра и стирает параметр объекта
Запускать функцию каждый раз, когда вызывается метод в классе
Tkinter: удалить границы всего Tabcontrol и его вкладок
Чего вы ожидаете от
np.math.factorial([0, 1, ...., 19])?