使用相同的概率数组对 scipy 个分布进行采样

Question

我想使用一组用户生成的概率从不同的 Beta 分布生成变量。我知道 scipy.stats.beta.rvs(a,b) 使用内部生成的概率生成随机变量，但我想改用我自己的概率。我怀疑这一定是微不足道的，但我自己无法解决。提前致谢。

Answer 1

经过一番摸索，我认为解决方案是使用 scipy.special.betainc(a, b, x)。这似乎给出了与以下脚本所示相同的 CDF：

import scipy.stats
import scipy.special
import matplotlib.pyplot as plt
import numpy as np

a = 3.0
b = 5.0

ntrial = 100000

p = np.linspace(0.0, 1.0, 1000)
r = np.random.rand(ntrial)
v = scipy.special.betainc(a, b, r)

# Plot the CDFs

plt.plot(np.arange(ntrial) / float(ntrial), np.sort(v), '.', label='scipy.special')
plt.plot(p, scipy.stats.beta.cdf(p, a, b), label='scipy.stats')
plt.ylabel("cumul. prob.")
plt.legend()

使用相同的概率数组对 scipy 个分布进行采样

sample scipy distributions using the same probability array

distribution

scipy.stats