为什么非标准化分布的 scipy 范数不同答案

【问题标题】：Why is scipy norm different for non-standardized distributions为什么非标准化分布的 scipy 范数不同
【发布时间】：2017-05-10 23:03:50
【问题描述】：

我有一个正态分布，平均值为 71，方差为 20.25。示例取自“Heads first statistics”。

当我将正态分布标准化为零均值时，我得到了正确的结果，但根据我对 scipy 和正态分布的理解，我应该得到非标准化分布的相同概率。

from scipy.stats import norm
import math

# prints 0.539337742276
print(norm(71, 20.25).sf(69))

zscore = (69-71) / math.sqrt(20.25)
print(norm(0,1).sf(zscore))
# prints 0.671639356718

【问题讨论】：

标签： python scipy normal-distribution

【解决方案1】：

注意norm 是用均值和尺度参数化的，而不是均值和平方尺度。因此，

>>> from scipy.stats import norm
>>> norm(71, pow(20.25,0.5)).sf(69)
0.6716393567181147
>>> zscore = (69-71) / pow(20.25,0.5)
>>> norm(0,1).sf(zscore)
0.6716393567181147

【讨论】：