python求标准正态分布的分布函数值（具体实现及调用函数）

在做生物统计作业时用到了 Φ(-6.6)，查表没有结果，所以想到用python解决。
以下代码主要是通过划分子区间求和的方式来计算。

import math
def normal_distribution(x): #处理x<0(目标点在分布中心左侧)的情况if x<0:return 1-normal_distribution(-x)if x==0:return 0.5 #求标准正态分布的概率密度的积分s=1/10000xk=[]for i in range(int(x/s)):integral = fx_normal_distribution((i+1)*s )xk.append(integral)print(integral)print(xk[0])sum=0for each in xk:sum+=eachprint(sum)return 0.5+sum*sdef fx_normal_distribution(x):return math.exp((-(x)**2)/2)/(math.sqrt(2*math.pi))print(normal_distribution(-6.6))

当时，在x取-6.6时，因为概率已经跌破10的-11次方，上面的代码精确度已经不行了，主要是此时的结果已经依赖s的取值了。而继续降低s的数量级，程序已经很慢了。
后来，我才发现python是可以直接调用函数来计算的，qaq。

from scipy.stats import normq = norm.cdf(-6.6)  # 累计密度函数
norm.ppf(q)  # 累计密度函数的反函数
print(q)
print(norm.ppf(q)) #可以用来找置信区间

这是代码的结果：

python求标准正态分布的分布函数值（具体实现及调用函数）相关推荐

python正态分布函数_用python求标准正态分布的上α分位点
累计分布函数bai值等于1-α的du点,即就是逆累zhi计分布函数icdf(1-α)的值如果dao使用Python计算的话,4102代码如下:1653 from scipy.stats import ...
关于求标准正态分布函数值
今年三月以来,一直在做工艺质量系统,基础做完了,开始做统计分析,其中需要计算烟支和烟丝的工艺合格率,一开始以为就是合格数/样品数,结果常常是100%合格.后来发现原来犯了自以为是的错误,在质量统计学里 ...
标准正态分布怎么算_求标准正态分布N（0,1）的特征函数.这个是怎么算出来的啊,...
优质解答 C(u)=E(j*u*X)=1/√(2*π)∫{-∞,+∞}e^(j*u*x-x²/2)dx,直接积分较困难由于d[e^(j*u*x-x²/2)]/dx=(j*u-x)*e^(j*u*x- ...
python求最小值不能使用min和sotred_python基础——内置函数
python基础--内置函数一.内置函数(python3.x) 内置参数详解官方文档: https://docs.python.org/3/library/functions.html?highl ...
python不同文件函数调用_不同输入文件的多处理调用函数
当进行多处理时,每个子进程都会获得主模块中在if __name__ == '__main__':语句之前定义的所有全局变量的副本.{{cdm>中访问的每个进程都将是不同的.在一种解决方法是使用 ...
python接口自动化（三十四）-封装与调用--函数和参数化（详解）
简介前面虽然实现了参数的关联,但是那种只是记流水账的完成功能,不便于维护,也没什么可读性,随着水平和技能的提升,再返回头去看前边写的代码,简直是惨不忍睹那样的代码是初级入门的代码水平都达不到.接下来 ...
python求上α分位点_什么是标准正态分布的上α分位点以及怎样求？
展开全部标准正态分布的上α分位点:设X~N(0,1),对于任给的62616964757a686964616fe4b893e5b19e31333431363638α,(0Zα)= α的点Zα为标准正态 ...
python 实现抽样分布（T分布、F分布、卡方分布）
标准正态分布 T分布 F分布卡方分布
python求概率密度函数_python 计算概率密度、累计分布、逆函数的例子
计算概率分布的相关参数时,一般使用 scipy 包,常用的函数包括以下几个: pdf:连续随机分布的概率密度函数 pmf:离散随机分布的概率密度函数 cdf:累计分布函数百分位函数(累计分布函数的逆 ...

python求标准正态分布的分布函数值（具体实现及调用函数）

python求标准正态分布的分布函数值（具体实现及调用函数）相关推荐

最新文章

热门文章