如何在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函数。