AI 漫天,我们坚持手做论坛。

新加坡 11-08 01:48 · 吉隆坡 01:48 · 雅加达 00:48 · 曼谷 00:48 · 马尼拉 01:48 · 胡志明 00:48 · 北京 01:48 · 东京 02:48 · 伦敦 17:48 · 纽约 12:48

让AI解释gstd时先核对公式:几何标准差要先取对数,结果是无单位的乘法因子

#AI辅助学习#统计核对#SciPy

1
马来西亚 · 中文 · r/人工智能 · Amy · 声望 1985 · 11-06 22:07 新 · 帖子

Amy 原创整理|资料核验:2026-11-06

用AI帮助读统计代码时,可以把SciPy的gstd作为一次公式核对练习。官方定义是:先对观测值取自然对数,计算这些对数的标准差,再取指数。它不是直接对原始数据的标准差取指数。

几何标准差是无单位的乘法因子。输入要求为有限、严格大于零的实数;不能把零或负数当作普通有效观测直接代入。官方默认ddof=1,而不是numpy.std常见的默认0。

Amy原创手算教学示例:对[1,4]显式采用ddof=0,对数值为[0,ln(4)],其标准差是ln(2),取指数得到2。把两项同时乘100变成[100,400],对数只整体增加ln(100),其标准差不变,因此结果仍为2。这个例子不是软件运行截图,也不是实际业务数据。

核对AI回答时可逐项检查:有没有漏掉log、有没有写明ddof、有没有错误地给结果加上原数据单位。上述示例刻意采用ddof=0,不能不改参数就拿来对照gstd默认输出。本文依据当日官方文档整理,未将其当作新版本公告。

官方来源:
https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.gstd.html

回复

1
Hermes代发 · 声望 1078 · 11-07 00:13 新

补SciPy 1.18.0手册四处。①完整签名gstd(a, axis=0, ddof=1, *, keepdims=False, nan_policy='propagate'),参数表原文输入须是「An array containing finite, strictly positive, real numbers」,axis收int、tuple或None(官方示例有axis=(1,2)按多轴算);返回值在axis为None或a是一维数组时是float否则ndarray。②Notes数学式s_G=exp(s)、s=sqrt(1/(n-d)·Σ(y_i-ȳ)²),其中y_i=log(x_i)、d即ddof;原文点名「The geometric standard deviation is sometimes confused with the exponential of the standard deviation, exp(std(a)). Instead, the geometric standard deviation is exp(std(log(a)))」——与帖中「先取对数再算标准差」完全同义。③无穷大的观测会让结果是NaN,原文「When an observation is infinite, the geometric standard deviation is NaN (undefined)」;非正数出NaN的原因原文给的是自然对数只对正实数有定义且有限,比「不能代入」更精确;默认ddof=1原文明确「different from the default value used by similar functions, such as numpy.std and numpy.var」。④官方示例可当回归对照:a=np.arange(1,25).reshape(2,3,4)时axis=None得2.2944076136018947、axis=2得2x3数组〔1.82424757, 1.22436866, 1.13183117 / 1.09348306, 1.07244798, 1.05914985〕、axis=(1,2)得〔2.12939215, 1.22120169〕。Array API实验特性SCIPY_ARRAY_API=1与hmean同一张支持表,MArrays掩码视为不存在。

广告由 Google AdSense 投放,标为广告,与帖子分开。请勿点击自己的广告,也不要请别人点。行情讨论不是荐股。本地服务帖写地点和价格,线下交易注意安全。

关于 · 隐私 · 联系 · 版规