У меня есть два бейсболиста, и я хочу проверить, кто лучший игрок.Я использую np.meshgrid, чтобы создать прямоугольную сетку из массива значений x и массива значений y.
import numpy as np
from scipy.stats import beta
d = 0.0002
limits = np.arange(0.29, 0.33, d)
xx, yy = np.meshgrid(limits, limits, sparse=True)
pdf_1 = beta.pdf(x, player_1['alpha'], player_1[' beta'])
pdf_2 = beta.pdf(y, player_2['alpha'], player_2[' beta'])
joint_pdf = numpy.outer(pdf1, pdf2)
Теперь я хочу знать, сколько% объединенного PDF находится надстрока y = x .Это покажет мне, насколько велика вероятность того, что один игрок лучше другого, верно?
Как я могу отфильтровать все ячейки в сетке, где(x> y) и выполните в основном следующее:
where (x > y): return np.sum(beta(x, player_1['alpha'], player_1['beta']) * beta(y, player_2['alpha'], player_2['beta']) * d ** 2)
Спасибо