R语言绘制正太分布图,并进行正太分布检验
正态分布
判断一样本所代表的背景总体与理论正态分布是否没有显著差异的检验。
方法一 概率密度曲线比较法
看样本与正太分布概率密度曲线的拟合程度,R代码如下:
#画样本概率密度图
s <- rnorm(100) #产生样本
d <- density(s)
plot(d, col="green", ylim=c(0, 0.5))#添加正太分布概率密度图
s2 <- seq(from=-4, to=4, length.out=100)
lines(s2, norm_expression(s2), col="red")
画图结果如下:
方法二 正太Q-Q图法
使用Q-Q图来判断数据是否服从正太分布,R代码如下:
s <- rnorm(100) #产生样本 qqnorm(s) qqline(s)
画图结果如下,可见数据分布集中在对角线上,可以认为总体服从正太分布:
方法三 经验法则
约68.3%数值分布在距离平均值有1个标准差之内的范围,约95.4%数值分布在距离平均值有2个标准差之内的范围,以及约99.7%数值分布在距离平均值有3个标准差之内的范围。称为“68-95-99.7法则”或“经验法则”。
使用R的验证代码如下:s <- rnorm(10000) #产生样本
sum(abs(s - mean(s)) < sd(s)) / length(s)
sum(abs(s - mean(s)) < 2*sd(s)) / length(s)程序结果:
> s <- rnorm(10000) #产生样本
> sum(abs(s - mean(s)) < sd(s)) / length(s)
[1] 0.6871
> sum(abs(s - mean(s)) < 2*sd(s)) / length(s)
[1] 0.9538
方法四 统计检验方法
使用样本偏度和样本峰度来估计总体偏度和峰度,在正太分布的假定下,样本偏度和峰度均服从均值为零、方差分别为6/T和24/T的正太分布,可以分别检验偏度和峰度,也可以将两个统计量结合起来生成一个服从自由度为2的卡方分布的统计量,再进行检验【参见《金融时间序列分析》第三版P8~P9】。理论方面的东西略...
可以使用夏皮罗-威尔克(Shapiro-Wilk)检验,代码如下:
s <- rnorm(1000) #产生样本 shapiro.test(s)
检验结果:
> shapiro.test(s)
Shapiro-Wilk normality test
data: s
W = 0.9987, p-value = 0.6716
shapiro.test函数输出一个p值,照惯例,p<0.05说明总体不太可能是正太分布,否则不能提供这么个证据,也就是说这个检验比较保守,倾向于错误的过分证明正态性。
R的nortest包能提供其他正太检验方法【参见《R语言经典实例》P221】。
R语言绘制正太分布图,并进行正太分布检验相关推荐
- R语言绘制不一样的条形图
绘制条形图的方法有很多,这里介绍如果用R语言绘制一个不一样的条形图 准备数据,这里为了方便,我们使用已经存在于gcookbook包中的一个数据集 首先需要调用该包,如果该包不存在,可以使用下面的方式安 ...
- R语言绘制生存曲线图
R语言绘制生存曲线图 KMunicate是支持按照Morris等人的KMunicate研究推荐的方式生成Kaplan-Meier图. 1958年,Edward L. Kaplan 和Paul Meie ...
- R语言绘制Bump Chart
R语言绘制Bump Chart的小示例 # install.packages("tidyverse") # install.packages("ggbump") ...
- R语言绘制环形树状图
R语言绘制环形树状图 1.主要用到dendextend和circlize包绘图: library(dendextend) library(circlize)# 距离矩阵 d <- dist(US ...
- R语言绘制二维密度图
R语言绘制二维密度图 二维密度图显示了两个数值变量之间的关系,一个在x轴上表示,另一个在Y轴上表示,与散点图类似,然后计算二维空间中特定区域内的观测数,并用颜色梯度表示.二维密度图有几种类型,以下主要 ...
- R语言绘制带聚类树的堆叠柱形图
R语言绘制带聚类树的堆叠柱形图 聚类树与柱形图结合,即可反映样本或分组间的相似性,又能展示样本内的元素组成信息. 例如下图是一个在扩增子测序微生物群落分析中常见的统计图类型,在测序公司给的报告中通常都 ...
- R语言绘制线图(line)实战
R语言绘制线图(line)实战 目录 R语言绘制线图(line)实战 #仿真数据 #基础线图
- R语言绘制空白图实战
R语言绘制空白图实战 目录 R语言绘制空白图实战 #绘制空白图1 #绘制空白图2 #绘制空白图3
- R语言绘制火山图(volcano plot)实战:为差异表达基因(DEGs)添加颜色、基于显著性阈值进行点的颜色美化、为选定基因添加标签
R语言绘制火山图(volcano plot)实战:为差异表达基因(DEGs)添加颜色.基于显著性阈值进行点的颜色美化.为选定基因添加标签 目录 R语言绘制火山图(volcano plot)实战 #导入 ...
- R语言绘制核密度图实战(Kernel Density Plot)
R语言绘制核密度图实战(Kernel Density Plot) 目录 R语言绘制核密度图实战(Kernel Density Plot) #仿真数据
最新文章
- klib库下的kroundup32(二进制的四舍五入)算法
- 汇编语言典型例子详解_从架构到 RTOS 详解 DSP 和 MCU 的区别和联系
- `Caused by: java.lang.ClassNotFoundException: javax.xml.bind.DatatypeConverter
- Python 基础教程:常用函数整理
- 二十万字C/C++、嵌入式软开面试题全集宝典七
- 链表基本功能:初始化、增、删、查、改
- C语言去括号编程题,去括号 - C语言网
- linux编程学习_您需要编程技能才能学习Linux吗?
- 大工20春《计算机应用基础》在线测试,大工20春《计算机应用基础》在线测试1答案...
- iOS底层探索之类的加载(三): attachCategories分析
- ELK 企业内部日志分析系统
- CH2-Java编程基础(7个案例实现)
- error: ‘stol’ was not declared in this scope
- [医疗信息化][DICOM教程]HL7 V3 Standard-概述-HL7 V3 Standard - A High Level Overview
- 向产品宣战——开发者眼中的Android UI Design
- 猎聘网2017年校园招聘研发类岗位笔试题
- java二叉树计算深度
- 服务器修改上传图片大小限制,修改IIS7默认上传大小限制
- CSS精灵图及字体图标
- 2022年值得关注的 8 个人工智能趋势及中国人工智能行业发展情况
热门文章
- 让面试官心服口服:Thread.sleep、synchronized、LockSupport.park的线程阻塞有何区别?
- javascript 性能检测
- 判断素数最有效的算法
- CSS3 transition实现手风琴效果
- Gitlab设置ssh密钥详解
- 关于Servlet中注解:@WebServlet注解一直报错的解决办法
- MTK Camera Senor Bring up 复盘总结
- java实现求解汉诺塔问题(提示, 使用递归)
- c#语言编写汉诺塔游戏,c# 实现汉诺塔游戏
- 看《狂飙》读人生,致敬2023!