单因素方差分析_基于R语言开展方差分析(一)——单因素方差分析

8bab0f1ce8da4579b8ec3733c2a7c0a5.png

基本原理

方差分析(Analysis of variance, ANOVA)是用于两个或两个以上样本均数比较的方法,还可以分析两个或多个研究因素的交互交互作用以及回归方程的线性假设检验等。其基本思想是将全部观察值间的变异——总变异按设计和需要分解成两个或多个组成部分,再作分析。
本文将针对单因素方差分析(One-Way ANOVA)开展R语言的实践。单因素方差分析是将同质的受试对象随机地分配到各处理组,再观察对比实验效应的方法。在单因素方差分析中,首先将总变异(SS_总)分解为组间变异(SS_组间)也叫处理变异和组内变异(SS_组内)也叫误差变异,然后比较两者的平均变异(MS_组间,MS_组内),比较时采用两者的比值F值(MS_组间/MS_组内)

R语言实践

为研究钙离子对体重的影响作用,某研究者将36只肥胖模型大白鼠随机等分为三组,每组12只,分别给予高脂正常剂量钙(0.5%)、高脂中剂量钙(1.0%)和高脂高剂量钙(1.5%)三种不同的饲料,喂养9周,测其喂养前后体重的差值。结果如下:
高脂正常剂量钙(0.5%):332.96,297.76,312.57,295.47,284.25,307.97,292.12,244.61,261.46,286.46,322.49,282.42
高脂中剂量钙(1.0%):253.21,235.87,269.3,258.9,254.39,200.87,227.79,237.05,216.85,238.03,238.19,243.49
高脂高剂量钙(1.5%):232.55,217.71,261.15,220.72,219.46,247.27,280.75,196.01,208.24,198.41,240.35,219.56
首先,导入数据:

my_data = read.xlsx("./anova_1.xlsx", sheetIndex = 1, encoding = "UTF-8")

其次,对通过箱型图和线图对数据进行可视化:

# 绘制箱型图
ggboxplot(my_data, x = "剂量", y = "体重差值",color = "剂量", palette = c("#00AFBB", "#E7B800", "#FC4E07"),order = c("高脂正常剂量钙", "高脂中剂量钙", "高脂高剂量钙"),ylab = "体重差值", xlab = "剂量")
# 绘制线图
ggline(my_data, x = "剂量", y = "体重差值",add = c("mean_se", "jitter"),order = c("高脂正常剂量钙", "高脂中剂量钙", "高脂高剂量钙"),ylab = "体重差值", xlab = "剂量")

绘制结果如下:

92a44a092f62f4e6afc0c634e6a5ec9f.png

cbbc76b9dfc540898ca376543d0e8338.png

然后进行方差分析:

# 计算方差分析
res.aov <- aov(体重差值 ~ 剂量, data = my_data)
# 输出分析结果
summary(res.aov)

分析结果如下:

Df Sum Sq Mean Sq F value   Pr(>F)
剂量         2  28925   14462   27.21 1.04e-07 ***
Residuals   33  17540     532
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

上述结果表明,经假设检验得出多个总体均数不全等。
最后,我们可以计算Tukey HSD(Tukey Honest Significant Differences,R函数:TukeyHSD()可以实现),以在组均值之间进行多次成对比较

# 确定特定组对之间的平均差异是否具有统计显著性
TukeyHSD(res.aov)

输出结果如下:

Tukey multiple comparisons of means95% family-wise confidence levelFit: aov(formula = 体重差值 ~ 剂量, data = my_data)$剂量diff       lwr       upr     p adj
高脂正常剂量钙-高脂高剂量钙  64.86333  41.76836  87.95831 0.0000002
高脂中剂量钙-高脂高剂量钙    10.98000 -12.11498  34.07498 0.4810804
高脂中剂量钙-高脂正常剂量钙 -53.88333 -76.97831 -30.78836 0.0000064

上述结果表明说明中、高剂量钙与正常钙喂养9周前后体重差值不同,但不能认为中、高剂量钙喂养9周前后体重差值不同。

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

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

相关文章

个税10% 人群_人群管理如何使我们的搜索质量提高27%

个税10% 人群by Thanesh Sunthar由塔内什桑塔尔(Thanesh Sunthar) 人群管理如何使我们的搜索质量提高27&#xff05; (How Crowd Curation Improved Our Search Quality by 27%) The bigger your platform gets, the more vital search becomes. And if you run a content-hea…

mysql增数据语句_Mysql 数据增删改查语句

插入数据 insert#1. 插入完整数据(顺序插入)#语法一&#xff1a;insert into 表名(字段1,字段2,字段3…字段n) values (值1,值2,值3…值n);#语法二&#xff1a;insert into 表名 values (值1,值2,值3…值n);#2. 指定字段插入数据#语法&#xff1a;insert into 表名(字段1,字段2…

Python+Flask.0010.FLASK即插视图之自定义视图类及修饰器

2019独角兽企业重金招聘Python工程师标准>>> 即插视图; 说明: FLASK的视图灵感来自于DJANGO的基于类而非基于函数的通用视图,主要目的是为了解决多个视图函数之间已经实现的部分,通过类继承的方式继承到其它视图,总之为了一点,就是少写代码,然后通过add_url_rule让我…

InputStream和Reader,FileInputStream和 FileReader的区别

一、InputStream和Reader的区别 InputStream和Reader都可以用来读数据(从文件中读取数据或从Socket中读取数据)&#xff0c;最主要的区别如下: InputStream用来读取二进制数(字节流)&#xff0c;而 Reader用来读取文本数据&#xff0c;即 Unicode字符。那么二进制数与文本数据有…

NGUI之输入文本框的使用

ToolBar中的两个红圈 另&#xff0c;代码如下&#xff1a;只需要定义一个变量即可&#xff0c;然后将控件drag到那里&#xff0c;真的是灰常方便呀 还有一个就是保存了&#xff08;OK的响应&#xff09;,可以简单地理解为存档或读档 转载于:https://www.cnblogs.com/YTYMblog/p…

ae制作数据可视化_我如何精心制作真正可怕的数据可视化

ae制作数据可视化by Krist Wongsuphasawat克里斯特旺苏帕萨瓦(Krist Wongsuphasawat) 我如何精心制作真正可怕的数据可视化 (How I carefully crafted a truly terrible data visualization) Yes, you read that right. I am going to explain how I put together a really ba…

tensorrt轻松部署高性能dnn推理_实战教程:TensorRT中递归神经网络的介绍(中文字幕)...

NVIDIA TensorRT是一个高性能的深度学习推理优化器和运行时&#xff0c;它提供低延迟和高吞吐量。TensorRT可以从每个深度学习框架导入经过训练的模型&#xff0c;从而轻松地创建可以集成到大型应用程序和服务中的高效推理引擎。这个视频的五个关键点:1.TensorRT支持RNNv2, Mat…

w怎么接显示 树莓派zero_纯干货!一根线玩转树莓派ZeroW(图文教程,亲测有效)...

#一、写在前面本文旨在介绍如何用最少的外设(成本)完成树莓派Zero W最基础最重要的功能。注意&#xff1a;本文原始发表时官方镜像版本是2017-04-10的&#xff0c;在2019年5月10日有网友提出本方案已经不完全适用最新的镜像了&#xff0c;所以如果只是想按照本文所提出的步骤一…

十进制小数转换二进制的问题

2019独角兽企业重金招聘Python工程师标准>>> 整数和小数分别转换。 整数除以2&#xff0c;商继续除以2&#xff0c;得到0为止&#xff0c;将余数逆序排列。 22 / 2 11 余0 11/2 5 余 1 5 /2 2 余 1 2 /2 1 余 0 1 /2 0 余 1 所以22的二进制…

java操作mongodb(连接池)(转)

原文链接&#xff1a; java操作mongodb&#xff08;连接池&#xff09; Mongo的实例其实就是一个数据库连接池&#xff0c;这个连接池里默认有10个链接。我们没有必要重新实现这个链接池&#xff0c;但是我们可以更改这个连接池的配置。因为Mongo的实例就是一个连接池&#xff…

机器学习 一年入门_我作为自我入门程序员的一年回顾

机器学习 一年入门by Alin Rauta通过Alin Rauta 我作为自我入门程序员的一年回顾 (My Year as a Self-starter Programmer in Review) This was the most crucial year for my personal development ever. It was hard. Really hard. That’s why for me, the key word of 201…

声卡突然听不到监听_音乐人/键盘手伴侣物问题之:专业监听音箱的音质必须用独立声卡...

近日&#xff0c;不少朋友在后台留言&#xff0c;询问专业监听音箱连电脑听音乐要不要接个声卡&#xff01;本期我们针对此问题&#xff0c;跟大家分享一些心得与经验。先回答问题&#xff0c;当然要&#xff01;通常我们电脑上的音频输出口是这样的&#xff1a;而专业监听音箱…

helm3安装mysql_Helm3(kubernetes包管理工具)安装使用踩坑指南

image.png从结构中我们看到有不同级别的文件夹&#xff0c;以及一些yaml文件。charts&#xff1a; 用于存放其他依赖和关联的chart。例如应用依赖数据库的chart。Chart.yaml&#xff1a;存储一些元数据&#xff0c;例如chart的信息&#xff0c;描述等等templates文件夹&#xf…

矩阵快速幂及斐波那契数列模板

本篇博客先给出矩阵快速幂以及利用矩阵快速幂求斐波那契数列的模板&#xff0c;讲解待更新…… const int N10; int tmp[N][N]; void multi(int a[][N],int b[][N],int n) {memset(tmp,0,sizeof tmp);for(int i0;i<n;i)for(int j0;j<n;j)for(int k0;k<n;k)tmp[i][j]a…

Redis-3.2主从复制与集群搭建 推荐

Redis-3.2主从复制与集群搭建 一、Redis 主从搭建 1.下载并解压 yum install -y gcc gcc-c pcre zlib pcre-devel tcl wget http://download.redis.io/releases/redis-3.2.4.tar.gz tar -zxvf redis-3.2.4.tar.gz cd redis-3.2.4 make cd src && make test &&am…

AutoMapperHelper

/// <summary>/// AutoMapper帮助类/// </summary>public static class AutoMapperHelper{/// <summary>/// 单个对象映射/// </summary>public static T MapTo<T>(this object obj){if (obj null) return default(T);Mapper.CreateMap(obj.Ge…

web项目开发人员配比_我如何找到Web开发人员的第一份工作

web项目开发人员配比I have always had an interest in coding for the web. I built my first site almost 15 years ago using Yahoo’s Geocities, which allowed HTML styling and a few layout choices.我一直对网络编码感兴趣。 大约15年前&#xff0c;我使用Yahoo的Geoc…

苹果手机输入屏保后锁屏_修一块手机屏幕要7080元?

这几天华为Mate X的两次开售成为大家议论的话题&#xff0c;一些抢到的人自然沉浸在快乐之中&#xff0c;想着是自己留着用&#xff0c;还是转手赚一把。而一些想抢而没抢到的人或许正在研究如何在明天的第三次开售中抓好机会吧&#xff01;当然&#xff0c;也有像小编这样的&a…

中间介(MiddleWare)

引子-Django的生命周期 在学习中间介之前&#xff0c;我们先来回顾一下Django的生命周期&#xff1a;用户发起请求&#xff0c;请求会被发送到urlconf中的url&#xff0c;然后会指向对应的views函数进行处理&#xff0c;views函数处理完成后&#xff0c;用模板渲染好html&#…

对MariaDB10.0的Sphinx进行扩展

已修改过的文件&#xff1a;http://pan.baidu.com/s/1o8DHvkA 将这两个文件放到MariaDB的解压目录后&#xff0c;再进行安装 /usr/local/mariadb-10.0.28/storage/sphinx/ 如下是修改的代码 get_rec ( byte * buf, const byte * key, uint keylen,uint a,uint b,uint c );index…