我正在 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)
我可以只使用我自己的函数，但我很好奇我做错了什么。
【问题讨论】：
c 等于 5000 的值是巨大的，你确定这是正确的吗？我收到一个运行时错误，因为幂指数太大。
您没有正确地将您的参数映射到 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。这是通常被称为“威布尔分布”的分布的实现：
In [49]: from scipy.stats import weibull_min
In [50]: weibull_min.cdf(x, a, loc=0, scale=c)
Out[50]: 0.08555935639278299
【讨论】：