0

0

如何理解对数正态分布中异常大的方差值?

碧海醫心

碧海醫心

发布时间:2026-01-20 14:25:01

|

640人浏览过

|

来源于php中文网

原创

如何理解对数正态分布中异常大的方差值?

对数正态分布的方差天然具有“指数级放大”特性,当底层对数变量标准差较大时,原始变量的方差会急剧增长;本文通过数学公式推导、数值验证与python代码示例,阐明该现象的合理性及常见误区。

在建模收入、消费等经济学变量时,常假设其服从对数正态分布(Lognormal Distribution)——即变量 $ y > 0 $,且 $ \ln y \sim \mathcal{N}(\mu, \sigma^2) $。此时,$ y = e^{\mu + \sigma Z} $(其中 $ Z \sim \mathcal{N}(0,1) $),其统计性质由 $\mu$ 和 $\sigma$ 决定,但直观上极易低估方差的量级

✅ 理论公式:方差为何如此巨大?

对数正态分布的均值与方差有闭式解:

  • 均值:
    $$\mathbb{E}[y] = e^{\mu + \frac{\sigma^2}{2}}$$
  • 方差:
    $$\mathrm{Var}[y] = e^{2\mu + \sigma^2} \left(e^{\sigma^2} - 1\right)$$

注意:方差表达式中含 $e^{2\mu + \sigma^2}$ 和 $e^{\sigma^2}$ 两项——双重指数结构导致方差对 $\sigma$ 极其敏感。例如,当 $\mu = 7.5$、$\sigma = 0.8$ 时:

动态WEB网站中的PHP和MySQL:直观的QuickPro指南第2版
动态WEB网站中的PHP和MySQL:直观的QuickPro指南第2版

动态WEB网站中的PHP和MySQL详细反映实际程序的需求,仔细地探讨外部数据的验证(例如信用卡卡号的格式)、用户登录以及如何使用模板建立网页的标准外观。动态WEB网站中的PHP和MySQL的内容不仅仅是这些。书中还提到如何串联JavaScript与PHP让用户操作时更快、更方便。还有正确处理用户输入错误的方法,让网站看起来更专业。另外还引入大量来自PEAR外挂函数库的强大功能,对常用的、强大的包

下载
import numpy as np

mu, sigma = 7.5, 0.8
mean_theory = np.exp(mu + sigma**2 / 2)
var_theory  = np.exp(2*mu + sigma**2) * (np.exp(sigma**2) - 1)

print(f"理论均值: {mean_theory:.2f}")  # ≈ 2489.90
print(f"理论方差: {var_theory:.2f}")   # ≈ 5,557,849.03

该结果与问题中模拟得到的 Mean = 2484.87, Variance = 5,650,460.07 高度一致——差异仅源于蒙特卡洛抽样随机性。因此,并非代码错误,而是对数正态分布的固有数学特性

⚠️ 常见误解与注意事项

  • 误以为“方差大 = 数据异常或实现错误”
    实际上,只要 $\sigma > 0.5$,$e^{\sigma^2} - 1$ 就已显著大于 1(如 $\sigma=0.8 \Rightarrow e^{0.64}-1 \approx 0.90$),再乘以 $e^{2\mu+\sigma^2} \approx e^{15.64} \approx 8.3 \times 10^6$,方差自然达百万量级。

  • 混淆 $\mu$ 与原始变量均值
    $\mu$ 是 $\ln y$ 的均值,而非 $y$ 的均值;直接设 $\mu = \ln(\text{target_mean})$ 会导致偏差(正确做法是反解:$\mu = \ln(\mathbb{E}[y]) - \sigma^2/2$)。

  • 验证建议
    使用 scipy.stats.lognorm 生成样本并对比理论值(注意其参数化约定:s=sigma, scale=exp(mu)):

    from scipy.stats import lognorm
    np.random.seed(42)
    y_sim = lognorm.rvs(s=sigma, scale=np.exp(mu), size=100000)
    print(f"模拟均值: {y_sim.mean():.2f}")
    print(f"模拟方差: {y_sim.var():.2f}")

✅ 总结

对数正态分布的方差“看起来过大”,实则是其指数本质的必然体现。只要参数 $\mu$ 和 $\sigma$ 合理(如反映真实收入分布的偏态与离散度),高方差不仅正常,而且是模型捕捉厚尾性(heavy-tailedness)和不平等(如高基尼系数) 的关键机制。无需“修正”,而应理解其含义、善用其性质——例如,在政策模拟中,大方差正对应收入差距扩大的现实情境。

相关专题

更多
python开发工具
python开发工具

php中文网为大家提供各种python开发工具,好的开发工具,可帮助开发者攻克编程学习中的基础障碍,理解每一行源代码在程序执行时在计算机中的过程。php中文网还为大家带来python相关课程以及相关文章等内容,供大家免费下载使用。

765

2023.06.15

python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

640

2023.07.20

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

764

2023.07.25

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

639

2023.07.31

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

1305

2023.08.03

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

549

2023.08.04

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

579

2023.08.04

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

709

2023.08.11

Java JVM 原理与性能调优实战
Java JVM 原理与性能调优实战

本专题系统讲解 Java 虚拟机(JVM)的核心工作原理与性能调优方法,包括 JVM 内存结构、对象创建与回收流程、垃圾回收器(Serial、CMS、G1、ZGC)对比分析、常见内存泄漏与性能瓶颈排查,以及 JVM 参数调优与监控工具(jstat、jmap、jvisualvm)的实战使用。通过真实案例,帮助学习者掌握 Java 应用在生产环境中的性能分析与优化能力。

13

2026.01.20

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
最新Python教程 从入门到精通
最新Python教程 从入门到精通

共4课时 | 5.8万人学习

Django 教程
Django 教程

共28课时 | 3.3万人学习

SciPy 教程
SciPy 教程

共10课时 | 1.2万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号