这里使用scipy工具进行模拟:

from pylab import *

from scipy.optimize import curve_fit

data=concatenate((normal(1,.2,5000),normal(2,.2,2500)))

y,x,_=hist(data,100,alpha=.3,label='data')

x=(x[1:]+x[:-1])/2 # for len(x)==len(y)

def gauss(x,mu,sigma,A):

return A*exp(-(x-mu)**2/2/sigma**2)

def bimodal(x,mu1,sigma1,A1,mu2,sigma2,A2):

return gauss(x,mu1,sigma1,A1)+gauss(x,mu2,sigma2,A2)

expected=(1,.2,250,2,.2,125)

params,cov=curve_fit(bimodal,x,y,expected)

sigma=sqrt(diag(cov))

plot(x,bimodal(x,*params),color='red',lw=3,label='model')

legend()

print(params,'\n',sigma)

数据是两个正常样本的叠加,该模型是高斯曲线的总和.我们获得 :

估计参数是:

# via pandas :

# pd.DataFrame(data={'params':params,'sigma':sigma},index=bimodal.__code__.co_varnames[1:])

params sigma

mu1 0.999447 0.002683

sigma1 0.202465 0.002696

A1 226.296279 2.597628

mu2 2.003028 0.005036

sigma2 0.193235 0.005058

A2 117.823706 2.658789

python直方图的拟合_从一组数据python中将两个高斯拟合成直方图相关推荐

  1. python数据处理电脑配置_『大数据python计算机基础教学视频教程』入门学python需要什么配置的电脑...

    x = 1 if x 4.2 1) else匹前的缩进相同且最接近的if 2) 改之前y=x 2行,之后是y=x 1先执行 第5题 代码还n = 1 while n <= 5: i = 1 wh ...

  2. python百分比堆积条形图_绘制比例数据python(堆积条形图)

    我有一个数据集,客户回答问题,客户属于某个类别.该类别是序数.我想将百分比的变化可视化为比例堆积的条形图.以下是一些测试数据: answer | categ 1 1 2 1 3 2 1 2 2 3 3 ...

  3. python 条形图填充疏密_教你利用Python玩转histogram直方图的五种方法

    直方图 直方图是一个可以快速展示数据概率分布的工具,直观易于理解,并深受数据爱好者的喜爱.大家平时可能见到最多就是 matplotlib,seaborn 等高级封装的库包,类似以下这样的绘图. 本篇博 ...

  4. python 最优化算法库_哈工大硕士生用?Python 实现了 11 种经典数据降维算法,源代码库已开放...

    雷锋网 AI 开发者按:网上关于各种降维算法的资料参差不齐,同时大部分不提供源代码.这里有个 GitHub 项目整理了使用 Python 实现了 11 种经典的数据抽取(数据降维)算法,包括:PCA. ...

  5. python线性回归模型预处理_线性回归-2 数据预处理与模型验证评估

    主要内容数据向量化处理 特征放缩 上采样和下采样 重采样和交叉验证 模型验证 python 代码实现 1. 数据向量化处理 对于给定的m个样本,假设最终的拟合函数是 为拟合的权重系数,则有 损失函数改 ...

  6. python的顶级库_三大用于数据科学的顶级Python库

    Python有许多吸引力,如效率,代码可读性和速度,使其成为数据科学爱好者的首选编程语言.Python通常是希望升级其应用程序功能的数据科学家和机器学习专家的首选. 由于其广泛的用途,Python拥有 ...

  7. python编写数据库连接工具_详解使用Python写一个向数据库填充数据的小工具(推荐)...

    一. 背景 公司又要做一个新项目,是一个合作型项目,我们公司出web展示服务,合作伙伴线下提供展示数据. 而且本次项目是数据统计展示为主要功能,并没有研发对应的数据接入接口,所有展示数据源均来自数据库 ...

  8. 时间序列python复杂网络实例_基于时间序列数据的复杂网络重构

    摘要:现实生活中,复杂系统俯首即是,和人们的生活息息相关,遍及的领域包括社会,经济,物理,生物和生态等.为了便于研究,复杂网络通常被看成是复杂系统的一种抽象描述,其中,网络中的单个节点对应于系统中的单 ...

  9. mysql求数值分位点_如何计算一组数据中任意一个数值的分位值?

    本篇补充,关于评论区@慢中取利 的提问. 分析:你的疑问是 想要逆推 算出percentile()方法的k值,也就是分位值对应的百分数. 逆推,最好的办法我认为是 找到Excel 关于percenti ...

最新文章

  1. Guice系列之用户指南(七)
  2. pyradiomics的yaml文件参数设置解读、Li‘s have a solution and plan.
  3. 汇编:转移目的地址在内存中
  4. python基础练习_1
  5. 从GB到GBDT到XGBoost
  6. 逆向入门--第一次的HelloWorld
  7. c# in deep 之Lambda表达式
  8. 被360整的体无完肤,我真的怒了!
  9. 一步一步部署GlusterFS
  10. jar包可以到maven下载
  11. Docker安装Jenkins教程
  12. 火狐插件 打开html 死机,火狐flash插件崩溃(Firefox火狐Flash插件卡死问题完美解决方法)...
  13. RFID天线技术 – 低频和高频RFID天线设计
  14. proximal Algorithm
  15. 鸿蒙 林蒙 秦羽,从主神归来的主角们
  16. 鸡嗉囊炎有哪些症状 什么药防治鸡嗉囊肿大
  17. 类中数组成员变量怎么初始化,数组名与指针的区别
  18. IKAnalyzer中文分词分析内容目录
  19. 基于Android的学生学习打卡监督系统
  20. 洛谷 P1085 不高兴的津津 C语言

热门文章

  1. 怎样查看cudnn版本_ubuntu16.04+Geforce GTX1060安装CUDA和cuDNN
  2. java水果超市mysql_Java基础 | 项目实战之水果超市
  3. 暗备用的运行状态_瞧:我利用“无偏二极管”发明设计出了【宇宙“暗物质”、“暗能量”探测器】...
  4. 信捷伺服刚性调整_信捷電氣(603416):伺服與PLC增長將加速,口罩機解決方案帶來新增量...
  5. unity asset store下载不了_Unity手游实战:从0开始SLG——资源管理系统-基础篇(三)AssetBundle原理...
  6. word中怎样单独删除某一页的页码
  7. 引用文献管理软件Mendeley
  8. Hadoop分布式集群搭建详细过程
  9. Android之WebView学习
  10. ajax将响应结果显示到iframe,JavaScript:iframe / Ajax / JSON