R做方差齐次检验

一组数据需要做t检验，了解了一下t检验分为：a.单样本t检验，b.独立双样本t检验和c.成对或非独立样本t检验，三者的区别和选择见https://www.jmp.com/zh_cn/statistics-knowledge-portal/t-test.html。我的数据应使用b：独立双样本t检验，b又分多类，有等方差假设与异方差假设，为了确定使用哪个假设，在做t-test之前使用R做了下方差齐性检验。

据查资料，F检验、Bartlett χ2检验、Levene检验、残差图均可以做方差齐性检验，但Levence检验所分析资料可不具正态性，结果更为稳健，也可以检验多个总体的方差齐性，因此使用Levence进行检验。

借助于R中car包中的Levene进行统计检验

library(car)

df2<-read.table("all_var.txt",header = T,sep = "\t")

head(df2)

leveneTest(df2$len,df2$supp,center=mean)

leveneTest(df2$len,df2$supp,center=median)

数据格式如下：

假设两组数据方差齐次, 结果如下：

P>0.05，说明两组数据方差齐次。

ps：R计算方差：

a=c(22.61905, 38.18182)

v1<-var(a)

参考：

https://zhuanlan.zhihu.com/p/110078594

https://www.jmp.com/zh_cn/statistics-knowledge-portal/t-test/two-sample-t-test.html

https://zhuanlan.zhihu.com/p/126351774