admin管理员组文章数量:1130349
学术论文和杂志中最常见的是用“字母”和“星号”标记差异显著性,在之前的文章《显著性分析后如何标记“abc”?》已经介绍过如何用SPSS比较均值的结果进行“纯手动”添加字母显著性标记。嗯,方法似乎非常“绕”,且一走神就容易错。
那么有没有R脚本可以实现相关的方差分析、均值比较、生成变量均值对应的“字母标记”?当然,如果能直接将字母标记添加到图表(比如柱状图)中,岂不美滋滋?
闲言少叙,接下来主要为大家介绍如何用R进行方差齐性检验(Bartlett test 和Levene test)、方差分析、均值的多重比较方法(TukeyHSD和LSD法),最后用ggplot2包进行数据可视化。示例数据和脚本可通过点击 阅读原文下载
# 读取示例数据data
data
方差齐性检验nom
# 如果是2(多)因数,使用interaction()函数,leveneTest(y~interaction(var1,var2),data= data)
nom
nom$p.value # 输出p值,当p>0.05时,方差是齐性的。
# 使用另一种方法进行齐性检验,SPSS软件使用Levene test进行齐性检验,这里需要安装car包。
install.packages("car")
library("car")
nom1
# 如果是2因素,使用*,leveneTest(y~ivar1*var2,data = data)
nom1
nom1$`Pr(>F)` # 输出p值,当p>0.05时,方差是齐性的。
方差分析# 单因素方差分析,整体来看差异显著
oneway
<学术论文和杂志中最常见的是用“字母”和“星号”标记差异显著性,在之前的文章《显著性分析后如何标记“abc”?》已经介绍过如何用SPSS比较均值的结果进行“纯手动”添加字母显著性标记。嗯,方法似乎非常“绕”,且一走神就容易错。
那么有没有R脚本可以实现相关的方差分析、均值比较、生成变量均值对应的“字母标记”?当然,如果能直接将字母标记添加到图表(比如柱状图)中,岂不美滋滋?
闲言少叙,接下来主要为大家介绍如何用R进行方差齐性检验(Bartlett test 和Levene test)、方差分析、均值的多重比较方法(TukeyHSD和LSD法),最后用ggplot2包进行数据可视化。示例数据和脚本可通过点击 阅读原文下载
# 读取示例数据data
data
方差齐性检验nom
# 如果是2(多)因数,使用interaction()函数,leveneTest(y~interaction(var1,var2),data= data)
nom
nom$p.value # 输出p值,当p>0.05时,方差是齐性的。
# 使用另一种方法进行齐性检验,SPSS软件使用Levene test进行齐性检验,这里需要安装car包。
install.packages("car")
library("car")
nom1
# 如果是2因素,使用*,leveneTest(y~ivar1*var2,data = data)
nom1
nom1$`Pr(>F)` # 输出p值,当p>0.05时,方差是齐性的。
方差分析# 单因素方差分析,整体来看差异显著
oneway
<版权声明:本文标题:r语言平均值显著性检验_R语言标记显著性标记 内容由热心网友自发贡献,该文观点仅代表作者本人, 转载请联系作者并注明出处:https://it.en369.cn/jiaocheng/1754939676a2744065.html, 本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容,一经查实,本站将立刻删除。


发表评论