excel统计分析——CVM正态性检验

参考资料:统计推断——正态性检验(图形方法、偏度和峰度、统计(拟合优度)检验)_sm.distributions.ecdf-CSDN博客

29_张达成_从经验过程出发建立 Cramer-von Mises 统计量的性质 - 豆丁网

https://cran.r-project.org/web/packages/nortest/nortest.pdf

https://www.hep.uniovi.es/sscruz/Goodness-of-Fit-Techniques.pdf

马兴华,张晋昕.数值变量正态性检验常用方法的对比[J].循证医学,2014,14(02):123-128

        Cramer-Von Mises检验即CVM检验的思想基本与K-S检验一致,与K-S检验相比,CVM检验统计量度量的是经验累积分布函数和目标累积分布函数的平方距离的积分。就是把每个数据点的差求平方以后相加,得到总的分布偏差,这样就考虑了所有的差异点,而不是像K-S检验那样只考虑一个最大的。实际应用中CVM检验比K-S检验更容易拒绝正态分布假定。要求样本量大于7。

计算步骤如下:

1、将数据从小到大进行排序,形成x1≤x2≤...≤xn。

2、计算统计量:

W=\sum_{i=1}^{n}(x_{i}-\frac{i-1/2}{n})^{2}+\frac{1}{12n}

3、将样本数据调整为标准正态分布数据,如下(像R中nortest程序包中cvm.test函数):

W=\sum_{i=1}^{n}(\Phi (\frac{x_{i}-\bar{x}}{\sqrt{S^{2}}})-\frac{i-1/2}{n})^{2}+\frac{1}{12n}

其中,Φ(.)为标准正态分布的累计分布函数,\frac{i-1/2}{n}为经验分布函数。

4、计算修正统计量z

z=(1+\frac{0.5}{n})\times W

5、根据z值计算p值

当z<0.0275时:

\mathrm{ln}\,(1-p)=-13.953+775.5z-12542.61z^2

p=1-e^{-13.953+775.5z-12542.61z^2}

当0.0275<z<0.051时:

\mathrm{ln}\,(1-p)=-5.903+179.546z-1515.29z^2

p=1-e^{-5.903+179.546z-1515.29z^2}

当0.051<z<0.092时:

\textup{\textrm{}ln}\,p=0.886-31.62z+10.897z^2

p=e^{0.886-31.62z+10.897z^2}

当z>0.092时:

\mathrm{ln}\,p=1.111-34.242z+12.832z^2

p=e^{1.111-34.242z+12.832z^2}

excel操作步骤如下:

  • 11
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值