我希望为scipy.stats.powerlaw例程提供负指数,例如a=-1.5,以便绘制随机样本:
为什么需要 a > 0,如何提供负的 a 来生成随机样本,以及如何提供归一化系数/转换,即。
编辑:我应该补充一下,我正在尝试复制IDL的RANDOMP函数:
"""
powerlaw.pdf(x, a) = a * x**(a-1)
"""
from scipy.stats import powerlaw
R = powerlaw.rvs(a, size=100)
为什么需要 a > 0,如何提供负的 a 来生成随机样本,以及如何提供归一化系数/转换,即。
PDF(x,C,a) = C * x**a
文档在这里
http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.powerlaw.html
Thanks!编辑:我应该补充一下,我正在尝试复制IDL的RANDOMP函数: