#指标MoAP2ERF  
#正态分布检验
tapply(T$MoAP2ERF75,T$growing period,shapiro.test)  #（y，x，shapiro.test)
#方差齐性检验
library(car)
leveneTest(MoAP2ERF75~growing period,data=T)   #(y~x）
#执行独立样本t检验-参数检验   
t.test(MoAP2ERF75~growing period,data=T,var.equal=TRUE)   #var.equal=TRUE 方差齐性

#执行独立样本t检验-非参数检验wilcox.test()   
wilcox.test(MoAP2ERF66~growing period,data=T)  