python - 如何在 Python 中从 Wishart 分布中采样

Question

我需要用 Python 从 Wishart 分布中生成随机样本。有没有简单的方法可以做到这一点？

Answer 1

scipy 有一个 Wishart 包:

import numpy as np
from scipy.stats import wishart

x = np.linspace(1e-5, 8, 100) # make an array
for y in x:
    z = wishart.rvs(1, scale=y)
    print(z)

输出:

4.59465858669e-06
0.00403122342709
0.0268879506122
0.0100029090879
0.129477863995
0.372787021348
....

wishart.rvs 是用于随机变量采样的 Wishart 函数。