我有两个棒球选手,我想检查谁是更好的选手。 我使用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)的meshgrid中的所有单元,并基本上运行以下命令:
where (x > y): return np.sum(beta(x, player_1['alpha'], player_1['beta']) * beta(y, player_2['alpha'], player_2['beta']) * d ** 2)
谢谢