aov函数基本使用 | R语言与数据挖掘

#lab1：使用avo函数

a<-rnorm(1000,99,1)

b<-rnorm(1000,99,1)

data_tmp<-data.frame(x=c(a,b),a=factor(c(rep(1,1000),rep(2,1000))))

result<-aov(x~a,data_tmp)

summary(result)

#显示

#              Df Sum Sq Mean Sq F value Pr(>F)

#a              1    0.5  0.5499   0.563  0.453

#Residuals   1998 1952.3  0.9771              

df1<-summary(result)[[1]][,1][1];df1

sa<-summary(result)[[1]][,2][1];sa

df2<-summary(result)[[1]][,1][2];df2

se<-summary(result)[[1]][,2][2];se

pr<-1-pf((sa/df1)/(se/df2),df1,df2)#也可以使用summary(result)[[1]][,5][1]

if(pr<0.05){

    print(“在0.05的显著区间下，因素a对变量x有显著影响)

}else{

    print(“在0.05的显著区间下，因素a对变量x无显著影响)

}