Python：如何从非标准柯西分布中随机采样，因此具有不同的参数？

Question

我在这里看：numpy

我可以看到您可以使用np.random.standard_cauchy()指定数组的命令来从标准柯西采样。

我需要从 Cauchy 中采样，它可能具有x_0 != 0和gamma != 1，即可能不位于原点，也不具有等于 1 的比例。

我怎样才能做到这一点？

Answer 1

如果您有 scipy，则可以使用scipy.stats.cauchy，它需要一个位置 (x0) 和一个比例 (gamma) 参数。它公开了rvs抽取随机样本的方法：

x = stats.cauchy.rvs(loc=100, scale=2.5, size=1000)  # draw 1000 samples

Answer 2

您可以避免对 SciPy 的依赖，因为柯西分布是位置尺度系列的一部分。x这意味着，如果您从 中抽取样本Cauchy(0, 1)，只需将其移位x_0并乘以gamma并将x' = x_0 + gamma * x根据 进行分布Cauchy(x_0, gamma)。