R语言单因素方差分析简单小例子
时间:2022-07-23
本文章向大家介绍R语言单因素方差分析简单小例子,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
单因素方差分析是用来检验3组或者3组以上的数据间是否有差异的的一种统计分析方法。
比如下面用到的示例数据:
探究三种不同的肥料是否对某种作物的产量有影响。做实验的时候就是将一块地随机划分成若干小块,分别施用3种不同的肥料,最终统计产量。最后用单因素方差分析检验不同组之间的均值是否相等。
单因素方差分析的零假设是不同处理间的均值没有差异,如果计算得到的P值小于0.05,则拒绝原假设,即不同处理间是有差异的
以下内容参考https://www.scribbr.com/statistics/anova-in-r/
使用到的数据也可以在上面的链接处获取。
第一步是读入数据
df<-read.csv("one_way_anova/crop.data.csv",header=T)
head(df)
df$fertilizer<-as.factor(df$fertilizer)
summary(df)
数据是3种肥料,每种肥料30个数据
单因素方差分析
> one.way<-aov(yield~fertilizer,data = df)
> summary(one.way)
Df Sum Sq Mean Sq F value Pr(>F)
fertilizer 2 6.07 3.0340 7.863 7e-04 ***
Residuals 93 35.89 0.3859
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
p值小于0.05,不同肥料对应的产量是有差异的。具体谁和谁之间是有差异的还不知道,还需要借助多重检验。
tukey.test<-TukeyHSD(one.way)
plot(tukey.test)
image.png
上图的线只要跨0那条虚线说明没有差异。根据上图可以看出肥料2和和肥料1对应的产量没有差异,3和2,3和1有差异。
接下来用箱线图展示结果
library(ggplot2)
ggplot(df,aes(x=fertilizer,y=yield))+
geom_boxplot(aes(fill=fertilizer),notch=T,width=0.5)+
geom_jitter(width=0.1)+theme_bw()+
theme(legend.position = "top")
image.png
最终得出的结论就是这三种肥料中3的效果最好。
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法