文章/答案/技术大牛

发布

社区首页 >问答首页 >在Python/numpy/scipy中计算高斯生存函数的自然对数

问在Python/numpy/scipy中计算高斯生存函数的自然对数
EN

Stack Overflow用户

提问于 2016-06-10 19:00:13

回答 1查看 192关注 0票数 0

我可以计算高斯分布的生存函数的自然对数。

np.log( scipy.stats.norm.sf(s) )

我需要计算一些可笑的大值s(可能高达1000)的生存函数，但在s= 37附近，上面的函数达到了双点机器精度。有什么函数可以直接用来计算生存函数的对数吗？

注意:我不认为我的基本分布是高斯分布到那么多西格玛，但是我需要生存函数来计算弱信号(3-4西格玛)的一些性质，我希望算法在有很强信号的情况下做一些合理的事情。

回答已采纳

发布于 2016-06-10 19:18:48

使用logsf方法的scipy.stats.norm。

例如：

In [67]: from scipy.stats import norm

In [68]: norm.logsf(1000)
Out[68]: -500007.82669481222

票数 3

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/37755455

复制

相似问题

问在Python/numpy/scipy中计算高斯生存函数的自然对数EN