如何解决Scipy multivariate_normal cdf 与手动计算的不同
我以两种方式计算正态分布 cdf 的值:
- 使用 scipy 库:
from scipy.stats import multivariate_normal
from scipy.special import erf,erfc
sigma = 0.793387
one = multivariate_normal(1.0,sigma)
one.cdf(0.0)
输出:0.13078590311690458
- 手动计算
1.0 /2 * (1.0 + erf(-1.0 / np.sqrt(2) / sigma))
输出:0.10375913672098108
它们有很大的不同。
请你告诉我,哪里可能出错了?
UPD:当我使用 norm
中的 scipy.stats
时,它会显示手动计算的答案:
norm.cdf(0.0,1.0,sigma)
输出:0.10375913672098103
解决方法
multivariate_normal
的 Cdf 将协变矩阵作为参数,其中离散(但不是偏差)应位于对角线上。所以有一个正确的使用cdf的方法:
one = multivariate_normal(1.0,sigma * sigma)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。