R语言-单因素方差分析

单因素方差分析，一次性筛选出27至56列，58至126列中，对因变量显著影响的自变量
这里的因变量列名为oddreview，数据集为bdata，vari存储筛选出来的自变量的列名

library(stats)
vari<-c()
for(i in c(27:56,58:126)){
    aa<-oneway.test(bdata$oddreview~factor(bdata[,i]));
    if(aa$p.value<0.05) vari<-c(vari,names(bdata)[i])
}
#b变量存储整合的所有的单因素方差分析结果
od<-seq(1,300,3)
names(b)[od]
co<-c()
for(i in od){
    if(names(b)[i] %in% vari){co<-c(co,i,i+1,i+2)}
}
b<-b[,co]