excel统计分析——S-W正态性检验

参考资料:

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

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

【统计学】利用spss正态分布假定检验 S-W检验 K-S检验 直方图 Q-Q图

https://real-statistics.com/tests-normality-and-symmetry/statistical-tests-normality-symmetry/shapiro-wilk-test/

https://real-statistics.com/tests-normality-and-symmetry/statistical-tests-normality-symmetry/shapiro-wilk-expanded-test/

       S-W正态性检验法由Shapiro和Wilk1965年提出的,又称W检验。使用原始的S-W检验,在8≤n≤50时可以利用,小样本(n<8)对偏离正态分布的检验不太有效。不幸的是,针对大多数n统计量W的分布是未知的,必须通过模拟,造表或者近似方法获得。

       而J. P. Royston就是通过近似的方法改进了S-W方法,对样本量的要求为12≤n≤5000,甚至更多;但建议样本量不低于20;该方法被用在了很多软件中(如SAS和SPSS)。

原始版S-W检验,具体步骤如下:

1、将n个样本观测值按由小到大的次序排序。

2、计算S-W正态性检验统计量。

W=\frac{b^{2}}{SS}=\frac{(\sum_{i=1}^{m}a_{i}(x_{n+1-i}-x_{i}))^{2}}{\sum_{i=1}^{n}(x_{i}-\bar{x})^{2}}

其中,当n为偶数时,m=n/2;当n为奇数时,m=(n-1)/2。

3、根据给定的检验水平α和样本容量n计算统计量W的临界值W_{\alpha}

4、做出统计判断。若W<W_{\alpha},则认为样本数据所来自的总体不服从正态分布;若W>W_{\alpha},则认为服从正态分布。

上述中的a_{i}W_{\alpha}均为查表所得。

excel实现步骤如下:

注:已提前将需要查表的信息放到excel中,信息来源参考https://real-statistics.com/statistics-tables/shapiro-wilk-table/

Royston的改进版S-W方法:

1、将n个样本观测值按由小到大的次序排序。

2、利用excel统计分析——Q-Q图_excel 分位数图-CSDN博客中介绍的累计分布函数,取CDF=\frac{k-0.375}{n+0.25},获得n个数据对应累计概率值。

3、利用标准正态分布的概率密度函数的反函数,对得到的CDF进行计算,记作m_{i}(i=1...n)。

4、对m_{i}的平方求和,m=\sum_{i=1}^{n}m_{i}^{2}

5、系数a_{i},取u=1/\sqrt{n}

a_{n}=-2.706056u^{5}+4.434685u^{4}-2.071190u^{3}-0.147981u^{2}+0.221157u+m_{n}m^{-0.5}

a_{n-1}=-3.582633u^{5}+5.682633u^{4}-1.752461u^{3}-0.293762u^{2}+0.042981u+m_{n-1}m^{-0.5}

a_{i}=m_{i}/\sqrt{\epsilon } ,其中,2<i<n-1,\epsilon =\frac{m-2m_{n}^{2}-2m_{n-1}^{2}}{1-2a_{n}^{2}-2a_{n-1}^{2}}

a_{2}=-a_{n-1}

a_{1}=-a_{n}

系数a_{i}之间的关系为:\sum_{i=1}^{n}a_{i}^{2}=1

6、计算统计量:

W=\frac{(\sum_{i=1}^{n}a_{i}x_{i})^2}{\sum_{i=1}^{n}(x_{i}-\bar{x})^2}

根据W的公式可推算出,W是a与x之间的相关系数的平方,W的取值范围始终在0和1之间。

同样可以根据公式推算出:当n取值在12到5000之间时,ln(1-W)是近似服从正态分布的,其中μ和σ取值如下:

\mu =0.0038915(\mathrm{ln } n)^{3}-0.083751(\mathrm{ln } n)^{2}-0.31082\mathrm{ln } n-1.5861

\sigma =e^{0.0030302(\mathrm{ln }n)^{2}-0.082676\mathrm{ln }n-0.4803}

7、W统计量转换为z统计量:

z=\frac{\mathrm{ln}(1-W)-\mu}{\sigma}

如果z值对应的p值小于α,则认为不服从正态分布。

excel实现步骤如下:

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/301967.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

怎能错过这个技术集市!转发有奖,惊喜连连!

点击蓝字 关注我们Learn. Connect. Code. 微软一年一度的Build大会即将拉开帷幕&#xff08;5月25-27日&#xff09;&#xff0c;今年的主题演讲&#xff0c;全球CEO Satya Nadella将会深入我们的工作&#xff0c;探讨开发者速度&#xff0c;智能云原生应用&#xff0c;以及在…

文件服务器搭建centos,centos8搭建ftp文件服务器

1.安装vsftpdyum install -y vsftpd2.执行以下命令设置FTP服务开机自启动systemctl enable vsftpd.service3.执行以下命令启动FTP服务systemctl start vsftpd.service4.执行以下命令创建ftp用户useradd ftpuser5.执行以下命令并按照提示设置“ftpuser"用户密码passwd ftpu…

将历史、数学、语文、地理、政治知识融会贯通的诀窍就是它

▲卢sir特别推荐点击上图进入玩酷屋在这个知识都是碎片化的时代&#xff0c;系统化的知识&#xff0c;显得弥足珍贵。今天小木就专门给大家推荐一套&#xff0c;将知识系统化整理&#xff0c;放到时间轴上展示的——《时间上的人物谱》里。这是一套特别的工具书&#xff0c;是一…

Nginx实战部署常用功能演示(超详细版),绝对给力~~~

前言上次分享了一些开发过程中常用的功能&#xff0c;但如果到真实环境中&#xff0c;其实还需要一些额外的配置&#xff0c;比如说跨域、缓存、配置SSL证书、高可用等&#xff0c;老规矩&#xff0c;还是挑几个平时比较常用的进行演示分享。上篇详见Nginx超详细常用功能演示&a…

程序员,为什么给你50万年薪,你还要搞死我公司?

全世界只有3.14 % 的人关注了数据与算法之美小卢,今天看到一条让人吐血的新闻。1月20日&#xff0c;深圳市某互联网游戏公司程序员燕某在游戏上线测试的当天&#xff0c;锁死了服务器与电脑&#xff0c;并恶意失踪&#xff0c;致公司损失惨重的事在网上曝光并引起热议。这到底是…

访问Web服务器时 使用的协议是,使用SOAP协议访问Web服务

SOAP是简单对象访问协议&#xff0c;它可看成是HTTP与XML的结合&#xff0c;其中XML部分是作为HTTP报文的实体主体部分。SOAP&#xff1a;简单对象访问协议&#xff0c;是一种轻量的、简单的、基于 XML 的协议&#xff0c;它被设计成在 WEB 上交换结构化的和固化的信息。 SOAP …

同时画多个饼图_手帐术 | 这个神奇饼图里,藏着时间管理的小秘诀

说起效率管理kk首先想到的就是时间饼了简直就是拖延症的救星啊&#xff01;&#xff01;&#xff01;不仅能帮我们规划日程而且相比于时间轴来说可视化更强今天kk就来介绍一下时间饼的使用方法~时间饼类型时间饼一共分为几种12h、15h、18h、24h等时间越长分的越细12小时12h适合…

【转】java io 总结(图)

2019独角兽企业重金招聘Python工程师标准>>> 转载于:https://my.oschina.net/dlpinghailinfeng/blog/91191

TensorFlow的各种应用,你晓得不?

近几年&#xff0c;AI 的重大进展对我们的日常生活产生了积极影响。随着Google发布了《Google AI 原则》( Google AI Principles )&#xff0c;它们为 AI 发展提供了一个框架。由于这一领域的发展非常迅速&#xff0c;一些原则的最佳实践&#xff0c;如 “ 避免制造或加强不公平…

6月开招|工业互联才是王道,最高可达50k!

就在昨天&#xff0c;工信部发布了《工业互联网十个典型应用场景和五个重点行业实践》&#xff0c;具体介绍10个典型场景及5个重点行业工业互联网的实际应用情况&#xff0c;可以看到国内工业互联网正在蓬勃发展&#xff0c;欣欣向荣&#xff0c;工业4.0的浪潮正在席卷整个行业…

2018 年最受欢迎的 Python 库,你都用过吗?

全世界只有3.14 % 的人关注了数据与算法之美前段时间&#xff0c;数据科学网站 KDnuggets 评选出了顶级 Python 库 Top15&#xff0c;领域横跨数据科学、数据可视化、深度学习和机器学习。推荐阅读《Python3.0科学计算指南》上图&#xff1a;根据 GitHub star 和贡献评选出的 2…

登录服务器修改数据库吗,如何修改服务器登录数据库 sa

如何修改服务器登录数据库 sa 内容精选换一换为确保华为云关系型数据库服务发挥出最优性能&#xff0c;用户可根据业务需求对用户创建的参数模板中的参数进行调整。您可以修改用户创建的数据库参数模板中的参数值&#xff0c;但不能更改默认数据库参数模板中的参数值。以下是您…

用NSubstitute来mock

Moq在.net中比较出名的mock框架&#xff0c;NSubstitute也是属于这类框架&#xff0c;使用方式也如出一辙&#xff0c;是通过实现子类来达到mock的效果。下面的测试方法都是按照 UnitTest的三板斧Arrange Act Assert来展开的&#xff0c;代码具体如下。using Microsoft.VisualB…

ubuntu安装mysql_Ubuntu18.04下安装MySQL

提示&#xff1a;以下操作均在root用户下进行&#xff0c;如在普通用户&#xff0c;请自行加上sudo&#xff01;# 查看有没有安装MySQL&#xff1a;dpkg -l | grep mysql# 安装MySQL&#xff1a;apt install mysql-server安装完成之后可以使用如下命令来检查是否安装成功&#…

居家好物┃超好用的擦窗神器保护腰椎尾椎久坐不累的美臀坐垫

▲卢sir特别推荐点击上图进入玩酷屋春节即将来临&#xff0c;小木特地挑选了好几款超级好用的居家生活产品。第一款&#xff1a;擦外窗神器&#xff0c; 颠覆你擦窗方式&#xff01;双面长柄 可伸缩 可折叠 可旋转便携擦玻璃神器第二款&#xff1a;天然抗菌除螨毛巾&#xff0c…

【读书笔记】《编写高效的JavaScript程序》

为什么80%的码农都做不了架构师&#xff1f;>>> 看到一篇文章&#xff0c;http://www.csdn.net/article/2012-11-20/2811887-writing-fast-memory-efficient-javascript/2 里面有个网页针对 原型链&#xff0c;模块化&#xff0c;模块化加缓存&#xff0c;进行了比…

Cloudera将被私有化,Hadoop时代或将落幕

北京时间 6 月 1 日晚&#xff0c;据 TechCrunch 报道 [1]&#xff0c;KKR 和 CD&R 将以 53 亿美元收购 Cloudera&#xff0c;Cloudera 将被私有化。截止到本文发稿时&#xff08;11 点 49 分&#xff09;&#xff0c;Cloudera 股票报价 15.94 美元&#xff0c;涨幅 23.95%…

mysql group_concat去重_MySQL group_concat() 函数用法

MySQL group_concat() 函数用法在使用 group by对数据进行分组后&#xff0c;如果需要对 select 的数据项进行字符串拼接&#xff0c;这时就需要用到group_concat()函数。1、基本用法group_concat()完整语法如下&#xff1a;group_concat([DISTINCT] 要连接的字段 [Order BY AS…

请问:如何写出没有BUG的代码?

全世界只有3.14 % 的人关注了数据与算法之美1947年9月9日&#xff0c;美国海军准将 Grace Hopper 在哈佛学院计算机实验室里使用 Mark II 和 Mark III 计算机进行研究工作。她的团队跟踪到 Mark II 上的一个错误&#xff0c;操作人员发现是由于一只飞蛾钻到了 Mark II 的继电器…

非常适合新手的redis cluster搭建过程

Redis集群演进过程Redis单节点主从复制&#xff1a;复制主要实现了数据的多机备份&#xff0c;以及对于读操作的负载均衡和简单的故障恢复。故障恢复无法自动化&#xff1b;写操作无法负载均衡&#xff1b;存储能力受到单机的限制。哨兵&#xff08;Sentinel&#xff09;&#…