• python+正态分布+蒙特卡洛预测男女身高概率!


    python机器学习-乳腺癌细胞挖掘(博主亲自录制视频)https://study.163.com/course/introduction.htm?courseId=1005269003&utm_campaign=commission&utm_source=cp-400000000398149&utm_medium=share

     

    机器学习,统计项目联系QQ:231469242

     

    正态分布也称常态分布或常态分配,是连续随机变量概率分布的一种,是在数理统计的理论与实际应用中占有重要地位的一种理论分布。自然界人类社会,心理与教育中大量现象均按正态形式分布。例如能力的高低,学生成绩的好坏,人们的社会态度,行为表现以及身高、体重等身体状态。

    (高斯Carl Friedrich Gauss)

     

    正态分布是由阿伯拉罕·德莫弗尔(Abraham de Moivre)1733年发现的。其他几位学者如拉普拉斯(Marquis de Laplace)、高斯 (Carl Friedrich Gauss)对正态分布的研究也做出了贡献,故有时称正态分布为高斯分布。

    正态分布的函数(又称密度函数)为

     

    标准正态分布这两个参数分别为0与1。

    标准正态分布的密度函数可写作:

     

    (正态分布三个标准差的概率分布) 

    所有正太分布都可以转化成标准正态分布。

    期望值μ决定了其位置,其标准差σ决定了分布的幅度。

     

    正态分布具有很大医学意义。正态分布的应用某些医学现象,如同质群体的身高、红细胞数、血红蛋白量、胆固醇等,以及实验中的随机误差,呈现为正态或近似正态分布。

    我们提出一个常见的身高概率问题:假设男性平均身高175,标准差6;女性平均身高168,  标准差3;随机抽取一个女性和男性,女性高于男性随机概率是多少?

    我们不需要通过复杂公式来计算。只需让计算机产生足够多的随机值来模拟计算,最后得到答案。

    下面我们用Python的蒙特卡洛建模正态分布函数,解决这个男女身高概率问题。

    (matplotlib绘制两个正态分布,红色表示女性,蓝色表示男性)

     

    程序模拟10万个随机值,最后算出结果0.14727

    即女性高于男性随机概率为0.14727

     

    测试环境Anaconda Python2.7

    源代码已经本人测试无问题

    # -*- coding: utf-8 -*-

    ‘’’

    By Toby ,Blog:

    http://www.cnblogs.com/webRobot/

    ’’’
    import numpy as np
    import matplotlib.pyplot as plt
    import scipy.stats as stats
    import seaborn as sns
    import math,pylab,matplotlib,numpy
    from matplotlib.font_manager import FontProperties
    #设置中文字体
    font=FontProperties(fname=r"c:windowsfontssimsun.ttc",size=15)
    #测试n次
    n=100000
    #标准正太分布
    normalDistribution=stats.norm(175,6)
    #方差较大正态分布
    normalDistribution1=stats.norm(168,3)
    #一次随机
    def Random_single():
        array_male=normalDistribution.rvs(1)
        array_female=normalDistribution1.rvs(1)
        male=array_male[0]
        female=array_female[0]
        if female>male:
           return True
        else:
           return False
    #n次随机,返回count(女性高于男性的次数)
    def Multiple_random(n):
        count=0
        for i in range(n):
           value=Random_single()
           if value==True:
              count+=1
        return count
    # 计算女性高于男性概率
    def Probability(n):
        count=Multiple_random(n)
        p=count*1.0/n
        return p
    probability=Probability(n)
    print '随机次数',n
    print '女性高于男性概率:',probability
    #绘图
    x=np.arange(60,220)
    y=normalDistribution.pdf(x)
    y1=normalDistribution1.pdf(x)
    plt.plot(x,y,label="male")
    plt.plot(x,y1,'r',label="female")
    plt.xlabel("x")
    plt.ylabel("probability density")
    #plt.title("Normal distribution:mean=%.1f,standard deviation=%.1f"%(mean,std))
    plt.title("Normal distribution")
    plt.legend()
    plt.show()

    https://study.163.com/provider/400000000398149/index.htm?share=2&shareId=400000000398149(博主视频教学主页)

     

     
  • 相关阅读:
    .net开发微信(1)——微信订阅号的配置
    工作中EF遇到的问题
    .net Entity Framework初识1
    Razor视图
    jquery中利用队列依次执行动画
    .net找List1和List2的差集
    angularjs ng-if 中的ng-model 值作用域问题
    Spring Boot + JPA(hibernate 5) 开发时,数据库表名大小写问题
    springboot 启动排除某些bean 的注入
    angularjs 初始化方法执行两次以及url定义错误导致传值错误问题
  • 原文地址:https://www.cnblogs.com/webRobot/p/7722843.html
Copyright © 2020-2023  润新知