r语言平均值显著性检验_R语言标记显著性标记

学术论文和杂志中最常见的是用“字母”和“星号”标记差异显著性，在之前的文章《显著性分析后如何标记“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)