R语言和医学统计学(10):正态性和方差齐性检验
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文。
医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化。
文章目录
- 前言
- 正态性检验
- shapiro wilk检验
- kolmogorov smimov检验
- 方差齐性检验
- 两样本方差比较的F检验
- 两样本方差比较的Levene检验
- 多样本方差比较的Bartlett检验
- 多样本方差比较的Levene检验
前言
这是R语言和医学统计学的第10篇内容。
主要是用R语言复现课本中的例子。我使用的课本是孙振球主编的《医学统计学》第4版,封面如下:
正态性检验
在之前的推文总也提起过正态性检验的方法,比如通过QQ
图等。今天主要介绍计算法。
使用课本例8-3的数据
### 例8-3
RD1<-c(2.78,3.23,4.20,4.87,5.12,6.21,7.18,8.05,8.56,9.60)
RD2<-c(3.23,3.50,4.04,4.15,4.28,4.34,4.47,4.64,4.75,4.82,4.95,5.10)
shapiro wilk检验
进行正态性检验(W检验):
shapiro.test(RD1)
##
## Shapiro-Wilk normality test
##
## data: RD1
## W = 0.95903, p-value = 0.7748
P值大于0.05,符合正态性。就是这么简单!
kolmogorov smimov检验
ks.test(RD1, "pnorm")
##
## One-sample Kolmogorov-Smirnov test
##
## data: RD1
## D = 0.99728, p-value < 2.2e-16
## alternative hypothesis: two-sided
P值小于0.05,不符合正态性。就是这么简单!
方差齐性检验
两样本方差比较的F检验
使用课本例8-3的数据.
首先是F检验:
var.test(RD1, RD2)
##
## F test to compare two variances
##
## data: RD1 and RD2
## F = 16.836, num df = 9, denom df = 11, p-value = 6.517e-05
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
## 4.692491 65.864395
## sample estimates:
## ratio of variances
## 16.83618
P值小于0.05,两样本方差不齐!
两样本方差比较的Levene检验
先转换数据格式:
df <- data.frame(rd = c(RD1,RD2),group = c(rep("rd1",length(RD1)), rep("rd2", length(RD2)))
)str(df)
## 'data.frame': 22 obs. of 2 variables:
## $ rd : num 2.78 3.23 4.2 4.87 5.12 6.21 7.18 8.05 8.56 9.6 ...
## $ group: chr "rd1" "rd1" "rd1" "rd1" ...
进行Levene检验:
library(car)
## 载入需要的程辑包:carDataleveneTest(rd ~ group, data = df)
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 1 18.865 0.0003152 ***
## 20
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
P值小于0.05,两样本方差不齐!
多样本方差比较的Bartlett检验
使用课本例4-2的数据。
trt<-c(rep("group1",30),rep("group2",30),rep("group3",30),rep("group4",30))weight<-c(3.53,4.59,4.34,2.66,3.59,3.13,3.30,4.04,3.53,3.56,3.85,4.07,1.37,3.93,2.33,2.98,4.00,3.55,2.64,2.56,3.50,3.25,2.96,4.30,3.52,3.93,4.19,2.96,4.16,2.59,2.42,3.36,4.32,2.34,2.68,2.95,2.36,2.56,2.52,2.27,2.98,3.72,2.65,2.22,2.90,1.98,2.63,2.86,2.93,2.17,2.72,1.56,3.11,1.81,1.77,2.80,3.57,2.97,4.02,2.31,2.86,2.28,2.39,2.28,2.48,2.28,3.48,2.42,2.41,2.66,3.29,2.70,2.66,3.68,2.65,2.66,2.32,2.61,3.64,2.58,3.65,3.21,2.23,2.32,2.68,3.04,2.81,3.02,1.97,1.68,0.89,1.06,1.08,1.27,1.63,1.89,1.31,2.51,1.88,1.41,3.19,1.92,0.94,2.11,2.81,1.98,1.74,2.16,3.37,2.97,1.69,1.19,2.17,2.28,1.72,2.47,1.02,2.52,2.10,3.71)data1<-data.frame(trt,weight)
data1$trt <- factor(data1$trt)str(data1)
## 'data.frame': 120 obs. of 2 variables:
## $ trt : Factor w/ 4 levels "group1","group2",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ weight: num 3.53 4.59 4.34 2.66 3.59 3.13 3.3 4.04 3.53 3.56 ...
进行Bartlett检验:
bartlett.test(weight ~ trt, data = data1)
##
## Bartlett test of homogeneity of variances
##
## data: weight by trt
## Bartlett's K-squared = 5.2192, df = 3, p-value = 0.1564
由结果可知,不能认为不满足方差齐性!
多样本方差比较的Levene检验
library(car)leveneTest(weight ~ trt, data = data1)
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 3 1.493 0.2201
## 116
由结果可知,不能认为不满足方差齐性!
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文。
医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化。
R语言和医学统计学(10):正态性和方差齐性检验相关推荐
- 入门必学 | R语言数据的独立性,正态性及方差齐性检验
参数分析的三大前提检验 检验数据独立性的方法 Chisq检验 Fisher检验 Cochran-Mantel-Haenszel检验 检验数据正态性的方法 shapiro.test函数 qqnorm函数 ...
- R语言检验样本是否符合正态性(检验样本是否来自一个正态分布总体):shapiro.test函数检验样本是否符合正态分布(normality test)
R语言检验样本是否符合正态性(检验样本是否来自一个正态分布总体):shapiro.test函数检验样本是否符合正态分布(normality test) 目录
- R语言检验样本是否符合正态性(检验样本是否来自一个正态分布总体):使用nortest包的sf.test函数检验样本是否符合正态分布(normality test)
R语言检验样本是否符合正态性(检验样本是否来自一个正态分布总体):使用nortest包的sf.test函数检验样本是否符合正态分布(normality test) 目录
- bartlett方差齐性检验_R语言实用教程-数据正态性以及方差齐性检验
第一节:数据正态性分析 数据是不是正态分布是接下来选择正确的分析方法的基础.因为通常的方差分析用到的都要求是正态分布的数据(或者经过转换之后为正态的数据).而数据如果不满足正态性,那么通常要选择非参数 ...
- R语言数据分析笔记——t检验(含正态性检验和方差齐性检验在SPSS和R语言中的操作t检验(单样本、双独立样本、配对样本)在Excel、SPSS、R语言中的操作)
前言:本文为个人学习笔记,为各大网站上的教学内容之综合整理,综合整理了①假设分析的基础知识.②正态性检验和方差齐性检验在SPSS和R语言中的操作.③t检验(单样本.双独立样本.配对样本)在Excel. ...
- R语言和医学统计学(7):多元线性回归
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文. 医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化. 文章目录 前言 多元线性回归 回归诊断 可以通过看图来判断 ...
- R语言和医学统计学系列(1):t检验
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文. 医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化. 前言 本期开始将推送使用R语言进行医学统计学的相关内容. 使 ...
- R语言和医学统计学(3):卡方检验
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文. 医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化. 前言 这是R语言和医学统计学的第3篇内容. 主要是用R语言复 ...
- R语言和医学统计学(6):重复测量方差分析
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文. 医学和生信笔记,专注R语言在临床医学中的使用,R语言数据分析和可视化. 文章目录 前言 重复测量数据两因素两水平的方差分析 重复测量 ...
最新文章
- 互联网高并发架构技术实践
- python两个数据表中的对应数据相加
- 【opencv】25.图像卷积cv::filter2D()以及c++代码实例
- c语言结构体定义字符串数组,C语言,结构体中字符串的声明(采用字符指针还是字符数组)...
- 三种常用的MySQL建表语句(转)
- mysql top virt_linux top命令VIRT,RES,SHR,DATA的含义
- 谁适合做Scrum Master?
- 学前儿童语言教育模拟试卷c卷,学前儿童语言教育模拟试卷参考答案.doc
- 转载 Visual Assist X 10.8.2001 破解版 支持VC2013及以下版本
- 软件工程课程实践- 项目开发计划
- 关于 Spring 中 getBean 的全流程源码解析
- 制作删不掉的php唯读档,ppt母版logo删除不了怎么办
- 文本数据分析:文本挖掘还是自然语言处理?
- java学习day31(redis2)数据类型
- Hystrix断路器执行原理
- 银河麒麟V10操控系统Qt安装
- (五十三) Android O wifi 状态机消息处理及状态切换流程分析-以WifiController为例
- 程序员绝不要做“IT民工”
- 关于scrapy shell 在ipython环境下 TAB出现debug信息解决方案
- MacOS开源视频播放器