带有 scipy.stats rvs() 函数的多元随机变量

Question

scipy.stats统计分布套件（ scipy.stats.norm 、 scipy.stats.uniform 、 scipy.stats.t等）都使用自己的.rvs()函数生成单变量数据系列，并且只有一个具有多变量再现： multivariate_normal ，对应于numpy的numpy.random.randn((N,K)) 。 事实上，在numpy.random发现的几乎所有统计分布numpy.random可以产生多变量数据。

我如何将scipy.stats的单变量分布函数扩展到多元数生成，因为它拥有一些在numpy.random找不到的分布，比如johnsonsu ？ 我必须自己手动制作一个将多个单变量连接在一起的循环函数吗？ 那应该是什么样子

Answer 1

我认为您只想将size参数传递给rvs 。 例如：

from scipy import stats

stats.norm.rvs(size=10)

会给你一个充满 10 个标准正态变量的向量。

请注意，多变量意味着统计中的特定内容，而不仅仅是相同的IID副本（这就是size作用）。 例如， multivariate_normal的cov参数指定了一次绘制中所有变量的协方差矩阵。

再比如， multinomial类似，但分布的参数明显不同。

stats.multinomial.rvs(n=5, p=[0.5, 0.5], size=10)

告诉你投掷硬币 5 次得到多少正面和反面，重复 10 次。 即行是单独的绘制，列是绘制中的值