正态性检验方法汇总

2024-03-09 04:44
文章标签 方法 汇总 检验 正态性

本文主要是介绍正态性检验方法汇总,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

🍉CSDN小墨&晓末:https://blog.csdn.net/jd1813346972

   个人介绍: 研一|统计学|干货分享
         擅长Python、Matlab、R等主流编程软件
         累计十余项国家级比赛奖项,参与研究经费10w、40w级横向

文章目录

  • 1 图示法
  • 2 偏度、峰度检验法
  • 3 非参数检验方法
  • 4 主要R包及代码

该篇文章主要介绍了统计检验中正态性检验的常用方法理论,包括图示法、偏度/峰度检验法、非参数检验方法,同时附相关R语言程序进行正态性检验的常用函数。

1 图示法

   ①P-P图

   以样本的累计频率作为横坐标,以按照正态分布计算的相应累计概率纵坐标,以样本值表现为直角坐标系的散点。如果数据服从正态分布,则样本点应围绕第一象限的对角线分布。

   ②Q-Q图

   以样本的分位数作为横坐标,以按照正态分布计算的相应累计概率作为纵坐标,把样本表现为直角坐标系的散点。如果数据服从正态分布,则样本点应围绕第一象限的对角线分布。

​    ③直方图

   判断样本直方图是否以钟型分布,同时可以选择输出正态性曲线。

   ④箱线图

   通过观察矩形位置和中位数,若矩形位于中间位置,且中位数位于矩形的中间位置,则分布较为均匀,否则是偏态数据。

   ⑤茎叶图

   通过观察图形的分布状态,是否对称分布。

2 偏度、峰度检验法

  ①偏度、峰度系数分布性质

  样本k阶中心距:

b k = 1 n − 1 ∑ i = 1 n ( x i − x ˉ ) k b_k=\frac{1}{n-1}\sum_{i=1}^n (x_i-\bar{x})^k bk=n11i=1n(xixˉ)k

  样本偏度系数:

r 1 = b 3 b 2 2 3 r_1=\frac{b_3}{b_2^{\frac{2}{3}}} r1=b232b3

  该系数用于检验样本对称性,当 r 1 = 0 r_1=0 r1=0时,分布是对称的;当 r 1 > 0 r_1>0 r1>0时,呈正偏态;当 r 1 < 0 r_1<0 r1<0时,呈负偏态。

  样本峰度系数:

r 1 = b 4 b 2 2 − 3 r_1=\frac{b_4}{b_2^2}-3 r1=b22b43

  该系数用于检验样本峰态,当 r 2 = 0 r_2=0 r2=0时,分布是正态分布;当 r 2 > 0 r_2>0 r2>0时,呈尖峰分布;当 r 2 < 0 r_2<0 r2<0时,呈扁平分布。

  考虑原假设 H 0 : F ( x ) H_0:F(x) H0:F(x)服从正态分布:备择假设 H 1 : F ( x ) H_1:F(x) H1:F(x)不服从正态分布。

  当原假设为真时,检验统计量:

r 1 6 n ∼ N ( 0 , 1 ) \frac{r_1}{\sqrt{\frac{6}{n}}}~\sim N(0,1) n6 r1 N(0,1) r 2 24 n ∼ N ( 0 , 1 ) \frac{r_2}{\sqrt{\frac{24}{n}}}~\sim N(0,1) n24 r2 N(0,1)

  在显著性水平为 α \alpha α时,若满足 ∣ r 1 6 / n ∣ > Z 1 − α / 2 |\frac{r_1}{\sqrt{6/n}}|>Z_{1-\alpha/2} 6/n r1>Z1α/2 ∣ r 2 24 / n ∣ > Z 1 − α / 2 |\frac{r_2}{\sqrt{24/n}}|>Z_{1-\alpha/2} 24/n r2>Z1α/2,则拒绝原假设。

  ②偏度、峰度联合分布检验法(Jarque-Bera检验)

  检验统计量:

J B = r 1 2 6 / n + ( r 2 − 3 ) 24 / n ∼ χ ( 2 ) JB=\frac{r_1^2}{6/n}+\frac{(r_2-3)}{24/n}~\sim \chi(2) JB=6/nr12+24/n(r23) χ(2)

  当JB统计量的P值小于显著性水平 α \alpha α时,则拒绝正态性分布的原假设。

3 非参数检验方法

  ①Kolmogorov-Smirnor检验(KS)

  KS检验基于经验分布函数,该检验用大样本近似。该检验方法为比较一个频率分布f(x)与理论分布g(x)或者两个观测值分布的检验方法。其原假设 H 0 H_0 H0两个数据分布一致或者数据符合理论分布。统计量:

D n = s u p x ∣ F n ( x ) − F ( x ) ∣ D_n=\underset{x}{sup}|F_n(x)-F(x)| Dn=xsupFn(x)F(x)

  当实际观测值 D n > D ( n , α ) D_n>D(n,\alpha) Dn>D(n,α)则拒绝 H 0 H_0 H0,否则则接受 H 0 H_0 H0假设。

  需要注意的是样本数据如果有结点(及是重复的数据),则无法计算准确的P值。

  ②Lilliefors检验(K-S检验的修正)

  当总体均值和方差未知时,用样本的均值和方差代替后,再用K-S检验法。

原假设 H 0 : F ( x ) 服从正态分布 ; 备择假设 H 1 : F ( x ) 不服从正态分布。 原假设H_0:F(x)服从正态分布;备择假设H1:F(x)不服从正态分布。 原假设H0:F(x)服从正态分布;备择假设H1F(x)不服从正态分布。

  ③Shapiro-Wilk检验(Shapiro)

  W 检验全称 Shapiro-Wilk 检验,是一种基于相关性的算法。计算可得到一个相关系数,它越接近 1 就越表明数据和正态分布拟合得越好。

  W检验统计量:

W = ( ∑ i = 1 n a i x i ) 2 ∑ i = 1 n ( x i − x ˉ ) 2 W=\frac{(\sum_{i=1}^na_ix_i)^2}{\sum_{i=1}^n(x_i-\bar{x})^2} W=i=1n(xixˉ)2(i=1naixi)2

  其中 x ( x 1 ) < x ( x 2 ) < … < x ( n ) x_{(x_1)}<x_{(x_2)<\dots<x_(n)} x(x1)<x(x2)<x(n)

( a 1 , a 2 , … , a n ) = m T V − 1 ( m T V − 1 V − 1 m T ) 1 2 (a_1,a_2,\dots,a_n)=\frac{m^TV^{-1}}{(m^TV^{-1}V^{-1}m^T)^{\frac{1}{2}}} (a1,a2,,an)=(mTV1V1mT)21mTV1;

m = ( m 1 , m 2 , … , m n ) m=(m_1,m_2,\dots,m_n) m=(m1,m2,,mn)是样本顺序统计量的样本均值,V是样本。

  若 W > W α W>W_{\alpha} W>Wα,则接受正态性假设。

  ④ χ 2 \chi^2 χ2拟合优度检验

  检验统计量:

χ 2 = ∑ i = 1 k ( f i − n p i ) 2 n p i ∼ χ 2 ( k − 1 ) \chi ^2=\sum_{i=1}^k\frac{(f_i-np_i)^2}{np_i} \sim \chi ^2(k-1) χ2=i=1knpi(finpi)2χ2(k1)

χ 2 = ∑ i = 1 k ( f i − n p i ^ ) 2 n p i ^ ∼ χ 2 ( k − r − 1 ) \chi ^2=\sum_{i=1}^k\frac{(f_i-n\hat{p_i})^2}{n\hat{p_i}} \sim \chi ^2(k-r-1) χ2=i=1knpi^(finpi^)2χ2(kr1)

  若r为被待估参数。

  若服从正态分布, χ 2 \chi ^2 χ2应较小,其中p值大于显著性水平 α \alpha α时,拒绝原假设。

  ⑤达戈斯提诺(D′Agostino)法

  即D检验,1971提出,正态性D检验该方法效率高,是比较精确的正态检验法。

  ⑥AD检验

  ⑦CVM检验

4 主要R包及代码

  ①ks.test()

  例如零假设为N(15,0.2),则ks.test(x,“pnorm”,15,0.2)。如果不是正态分布,还可以选"pexp", "pgamma"等。
   ②shapiro.test()

  可以进行关于正态分布的Shapiro-Wilk检验。
  ③nortest包

lillie.test()可以实行更精确的Kolmogorov-Smirnov检验。
ad.test()进行Anderson-Darling正态性检验。
cvm.test()进行Cramer-von Mises正态性检验。
pearson.test()进行Pearson卡方正态性检验。
sf.test()进行Shapiro-Francia正态性检验。

  ④fBasics包

normalTest()进行Kolmogorov-Smirnov正态性检验。
ksnormTest()进行Kolmogorov-Smirnov正态性检验。
shapiroTest()进行Shapiro-Wilk's正态检验。
jarqueberaTest()进行jarque-Bera正态性检验。
dagoTest进行D'Agostino正态性检验。
gofnorm采用13种方法进行检验,并输出结果。

这篇关于正态性检验方法汇总的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/789546

相关文章

Java学习手册之Filter和Listener使用方法

《Java学习手册之Filter和Listener使用方法》:本文主要介绍Java学习手册之Filter和Listener使用方法的相关资料,Filter是一种拦截器,可以在请求到达Servl... 目录一、Filter(过滤器)1. Filter 的工作原理2. Filter 的配置与使用二、Listen

Pandas统计每行数据中的空值的方法示例

《Pandas统计每行数据中的空值的方法示例》处理缺失数据(NaN值)是一个非常常见的问题,本文主要介绍了Pandas统计每行数据中的空值的方法示例,具有一定的参考价值,感兴趣的可以了解一下... 目录什么是空值?为什么要统计空值?准备工作创建示例数据统计每行空值数量进一步分析www.chinasem.cn处

Java对象转换的实现方式汇总

《Java对象转换的实现方式汇总》:本文主要介绍Java对象转换的多种实现方式,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录Java对象转换的多种实现方式1. 手动映射(Manual Mapping)2. Builder模式3. 工具类辅助映

Windows 上如果忘记了 MySQL 密码 重置密码的两种方法

《Windows上如果忘记了MySQL密码重置密码的两种方法》:本文主要介绍Windows上如果忘记了MySQL密码重置密码的两种方法,本文通过两种方法结合实例代码给大家介绍的非常详细,感... 目录方法 1:以跳过权限验证模式启动 mysql 并重置密码方法 2:使用 my.ini 文件的临时配置在 Wi

MySQL重复数据处理的七种高效方法

《MySQL重复数据处理的七种高效方法》你是不是也曾遇到过这样的烦恼:明明系统测试时一切正常,上线后却频频出现重复数据,大批量导数据时,总有那么几条不听话的记录导致整个事务莫名回滚,今天,我就跟大家分... 目录1. 重复数据插入问题分析1.1 问题本质1.2 常见场景图2. 基础解决方案:使用异常捕获3.

最详细安装 PostgreSQL方法及常见问题解决

《最详细安装PostgreSQL方法及常见问题解决》:本文主要介绍最详细安装PostgreSQL方法及常见问题解决,介绍了在Windows系统上安装PostgreSQL及Linux系统上安装Po... 目录一、在 Windows 系统上安装 PostgreSQL1. 下载 PostgreSQL 安装包2.

SQL中redo log 刷⼊磁盘的常见方法

《SQL中redolog刷⼊磁盘的常见方法》本文主要介绍了SQL中redolog刷⼊磁盘的常见方法,将redolog刷入磁盘的方法确保了数据的持久性和一致性,下面就来具体介绍一下,感兴趣的可以了解... 目录Redo Log 刷入磁盘的方法Redo Log 刷入磁盘的过程代码示例(伪代码)在数据库系统中,r

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Java中Switch Case多个条件处理方法举例

《Java中SwitchCase多个条件处理方法举例》Java中switch语句用于根据变量值执行不同代码块,适用于多个条件的处理,:本文主要介绍Java中SwitchCase多个条件处理的相... 目录前言基本语法处理多个条件示例1:合并相同代码的多个case示例2:通过字符串合并多个case进阶用法使用

Python中__init__方法使用的深度解析

《Python中__init__方法使用的深度解析》在Python的面向对象编程(OOP)体系中,__init__方法如同建造房屋时的奠基仪式——它定义了对象诞生时的初始状态,下面我们就来深入了解下_... 目录一、__init__的基因图谱二、初始化过程的魔法时刻继承链中的初始化顺序self参数的奥秘默认