python进行正态检验(python:5种正态性检验方法)
python进行正态检验(python:5种正态性检验方法)
2024-11-21 09:46:27  作者:温柔七岁  网址:https://m.xinb2b.cn/tech/gvu205243.html

写在前面的话:码字不易,点个赞关注一下作者再走吧[呲牙],我来为大家讲解一下关于python进行正态检验?跟着小编一起来看一看吧!


python进行正态检验

写在前面的话:码字不易,点个赞关注一下作者再走吧[呲牙]

================================================================

1.描述性统计方法(直观)

1.直方图

由于正态分布具有非常典型的中间高,两边低的图形特征,如果样本数据并不服从正态分布,我们可以通过直方图很快地分辨出来。更进一步地,Python可以辅助生成基于样本数据估计的正态曲线,这样就容易辅助我们进行判断。图形观察虽然直观,但是部分研究者认为单纯观察图形过于主观,因此我们也可以选择使用统计检验的方法去研究数据是否服从正态分布。操作步骤:

导入相关的包及数据

import pandas as pdimport matplotlib.pyplot as pltimport seaborn as sns#用来显示中文标签plt.rcParams['font.sans-serif']=['SimHei']#用来正常显示负号plt.rcParams['axes.unicode_minus']=False #导入数据data = pd.read_excel('年龄_收入.xlsx')#绘制年龄的频率直方图,20个箱体plt.hist(data.年龄,bins = 20)#绘制年龄带正态曲线的频率直方图sns.distplot(data.年龄)

2 P-P图及Q-Q图直方图是最长用于观察数据分布的常用图形选项,尤其是带正态曲线的直方图,可以非常直观地看到实际数据分布和正态曲线的对比,而P-P图及Q-Q图则是另一种选择,它可以直观给出实际数据分布和理论的差距。

值得注意的是,虽然P-P图及Q-Q图常用用于判断数据样本是否服从正态分布,但实际上它们也能判断数据样本是否服从其他的分布

P-P图:反映的是数据的实际累积概率与假定所服从分布的理论累积概率的符合程度。在此处,我们所假定的分布就是正态分布,如果数据样本是服从正态分布的话,那么实际的累积概率与理论的累积概率应该是相对一致的,放映在图形中就是数据点应该沿着图形的对角线分布。Q-Q图的原理与P-P图几乎一致。P-P图考察的是实际分布与理论分布的累积概率分布差异,而Q-Q图考察的是实际百分位数与理论百分位数的差异。同理在此处,我们所假定的分布就是正态分布,如果数据样本是服从正态分布的话,那么实际的分布应该是相对一致的,反映在图形中就是数据点应该沿着图形的对角线分布。在Python中,statsmodels包中目前主要提供的是Q-Q图的绘制

import statsmodels.api as smimport pylabsm.qqplot(data.年龄, line='s')pylab.show()sm.qqplot(data.收入, line='s')pylab.show()

2.K-S检验

柯尔莫戈洛夫-斯米诺夫检验(Kolmogorov-Smirnov test),一般又称K-S检验,是一种基于累计分布函数的非参数检验,用以检验两个经验分布是否不同或一个经验分布与另一个理想分布是否不同。K-S检验的原假设是“样本数据来自的分布与正态分布无显著差异”,因此一般来说,KS检验最终返回两个结果,分别是检验统计量及P值,检验结果P>0.05才是我们的目标。

from scipy.stats import kstest#cdf中可以指定要检验的分布,norm表示我们需要检验的是正态分布#常见的分布包括norm,logistic,expon,gumbel等kstest(data.年龄,cdf = "norm")

实际上,GraphPad不推荐使用单纯的Kolmogorov-Smirnov test方法

3. W检验

夏皮洛-威尔克检验(Shapiro—Wilk test),一般又称W检验。W检验是一种类似于利用秩进行相关性检验的方法。同样需要注意的是,W检验与K-S检验一样,原假设是“样本数据来自的分布与正态分布无显著差异”,因此一般来说,W检验最终返回两个结果,分别是检验统计量及P值。,检验结果P>0.05才是我们的目标。当数据集中的数据无重复值时,该方法的检验效果比较好,但是当数据集中有些数据不是独一无二的,即有些数据的数值是相同的,那么该方法的检验效果就不是很好

from scipy.stats import shapiroshapiro(data.年龄)

4.D'Agostino and Pearson omnibus normality test

GraphPad官方推荐使用该方法。首先计算 偏度和峰度以便在不对称和形状方面量化分布离高斯分布的距离。然后,其计算这些值中的每一个与高斯分布的预期值之间的差异,并基于这些差异的总和,计算各P值。这是一种通用和强大的正态性检验,推荐使用。请注意,D'Agostino开发了几种正态性检验。Prism使用的其中一个是“综合K2”检验。

import scipy.statsfrom scipy import statsstats.normaltest(data)#normaltest中pvalue>0.05符合正态分布

5.安德森-达令检验(Anderson-Darling test)

安德森-达令检验样本数据是否来自特定分布,包括分布:'norm', 'expon', 'gumbel', 'extreme1' or 'logistic'.原假设 H0:样本服从特定分布; 备择假设 H1:样本不服从特定分布

##用Anderson-Darling检验生成的数组是否服从正态分布import scipy.stats as statsstats.anderson(data_norm, dist='norm')'''输出AndersonResult(statistic=0.18097695613924714, critical_values=array([ 0.555, 0.632, 0.759, 0.885, 1.053]), significance_level=array([ 15. , 10. , 5. , 2.5, 1. ]))如果输出的统计量值statistic < critical_values,则表示在相应的significance_level下,接受原假设,认为样本数据来自给定的正态分布。'''

实际上,从已有的文献表明,对于数据分布的正态性研究,首选方法是图形观察,即利用直方图、P-P图或Q-Q图进行观察,如果分布严重偏态和尖峰分布则建议进行进一步的假设检验。如果图形分布结果不好判断,则再进行正态性检验。实际上,从已有的文献表明,对于数据分布的正态性研究,首选方法是图形观察,即利用直方图、P-P图或Q-Q图进行观察,如果分布严重偏态和尖峰分布则建议进行进一步的假设检验。如果图形分布结果不好判断,则再进行正态性检验。

其次,对于检验方法来说,对于K-S检验及W检验结果来说,有文献采用蒙特卡罗模拟方法进行多次验证,结果表明W检验结果相比于大部分方法都有较大的检验功效,而K-S方法的检验结果相对不佳。并且部分学者认为,K-S检验的实用性远不如图形工具,因为在样本量少时,该检验不太敏感,但是在样本量大时,该检验却过于敏感。因此正常情况下,我们更常采用W检验的结果。

值得注意的是,虽然说K-S检验结果相对不佳,但是不同检验方法对于样本量的敏感度是不一样的。在样本量较小的情况下(小于50个样本的情况下),请优先选择W检验;在样本量50-5000的情况下,可以酌情使用W检验及K—S检验;在样本量大于5000的情况下,请使用K-S检验结果,尤其是在SPSS中,当样本量大于5000的情况下,将只显示K-S检验结果,而不显示W检验结果。

  • 街舞基本动作零基础(三类街舞基本功)
  • 2024-11-22三类街舞基本功最近很多同学都在咨询星城街舞,问“街舞的基本功是什么?该怎么练习?”其实关于这个问题,要因人而异,每个人的每个人的练习方法也会有所不同,星城街舞在此列举三类常见的也是大部分dancer常用的基本功练习。
  • 四眼土狗有什么优点有什么缺点(致敬:即将消失的土狗)
  • 2024-11-22致敬:即将消失的土狗中国地大物博,品种好的狗也较多,到现在已经留存的较少特别在农村养狗从没有外国人养宠物的概念,也更想不起以何种特别的方式心疼土狗四眼土狗,古人称为“四眉犬”,有黄四眉和白四眉,即眉斑为黄色和白色,眉斑有。
  • 世界最贵最豪华的私人飞机(世界10大最豪华的私人飞机)
  • 2024-11-22世界10大最豪华的私人飞机对于飞机现代人都不佰生,很多人也常坐(可怜的是我没有坐过),但是私人飞机你有吗?我想这个很少人会有,像赵本山价值几亿的私人飞机就更稀有了,今天给大家盘点一下世界上世界十大最豪华的私人飞机10,达索猎鹰。
  • 老味道烧烤配方(老舅使用15年的烧烤配方)
  • 2024-11-22老舅使用15年的烧烤配方今天给大家分享一下老舅用了15年的烧烤腌制配方和撒料、酱料配方,顺带烧烤的注意点,方便您开店少趟雷,早日生意兴隆羊肉腌料分三种,分别是不同的口感和风味,大家可以尝试过在决定用那种1,羊肉2.5公斤,孜。
  • 毛巾正确清洗方法(如何正确清洗毛巾?)
  • 2024-11-22如何正确清洗毛巾?盐搓法:毛巾用久了,常常会又湿又黏,有种怪味,且用肥皂越洗越黏如果用盐来搓洗,再用清水冲净,可使毛巾清洁如新,且没有异味煮沸法:将优质的洗涤剂溶液或者洗洁净溶液烧沸,将毛巾放入煮10分钟,捞出用清水漂。
  • 华为双系统手机是哪款(了解P9)
  • 2024-11-22了解P9华为p9是双系统,它的操作系统分别是华为EMUI4.1+Android6.0后置1200万*2,F2.2光圈,徕卡认证镜头,前置800万像素,徕卡镜头(目前全球唯一一款搭载徕卡摄像头的手机)支持4K播。
  • officeword 标尺分栏(Word排版进阶之分栏)
  • 2024-11-22Word排版进阶之分栏Office操作:Word排版进阶之分栏玩出不一样的Word分栏效果,Word提供了分栏效果,这大大丰富了文档的版式效果,而且分栏还可以起到节约版面及美化版面的效果一般性快速分栏为Word添加分栏非常。
  • 蒙古族改汉族可以改吗(民族大家庭又增一员)
  • 2024-11-22民族大家庭又增一员你知道在我们56个民族之外还有一个民族吗?他们叫穿青人,但是对于这样一个少有人知道的民族,社会存在广泛的争议,他们是否有自己独特的文化、是否得到了权威机关的认证,详情一起来看看穿青人,又称川青人主要分。
  • 明日方舟狙击三巨头是谁(明日方舟干员强度榜)
  • 2024-11-22明日方舟干员强度榜首先是推荐克洛斯(速射手),空爆(炮手),梅(速射手)克洛丝,罗德岛第一盲狙,这可不完全是玩梗,其天赋<要害瞄准>是个完完全全的神技(前期),提高了其破甲与输出能力,而且只是个三星这个不说。