python - Scipy Weibull CDF 计算

Question

我在 Scipy 中进行生存计算，但无法获得正确的值。

我的代码:

x, a, c = 1000, 1.5, 5000

vals = exponweib.cdf(x,a,c,loc=0,scale=1)

vals 应该等于 0.085559356392783004，但我得到的是 0。

如果我定义自己的函数，我会得到正确的答案:

def weibCumDist(x,a,c):
    return 1-np.exp(-(x/c)**a)

我可以只使用我自己的函数，但我很好奇我做错了什么。

Answer 1

你没有正确地将你的参数映射到 scipy 的参数。实现你的 weibCumDist 的等价物:

In [22]: x = 1000

In [23]: a = 1.5

In [24]: c = 5000

In [25]: exponweib.cdf(x, 1, a, loc=0, scale=c)
Out[25]: 0.08555935639278299

请注意 exponweib是exponentiated Weibull distribution .

您可能想使用 scipy.stats.weibull_min。这是通常称为“威 bool 分布”的分布的实现:

In [49]: from scipy.stats import weibull_min

In [50]: weibull_min.cdf(x, a, loc=0, scale=c)
Out[50]: 0.08555935639278299