step0 导入数据

mydata<- read.csv（“～/mydata.csv”，header=T）#y,code 其中，code为如2个种类，2个地方等

step1

tapply (y~factor(code), shapiro.test ) #正态性检验

library (car)

LeveneTest(y~factor(code)) #适用于正态和非正态数据的方差齐性分析

step2

t.test（y~factor(code), var.equal=T ） #学生t检验，适用于正态且齐性数据，自由度为整数

t.test (y~factor(code)) #Welch 检验，自由度常变为非整数

wilcox.test(y~factor(code)) #此时为数据非齐性的非参数检验Wilcoxon-Mann-Whitney U test

#根据默认虚无假设H0，上述检验结果p-value<0.05 即2组code数据y有显著差异，反之则无显著差异。