R/stat_anova.R

In msk0609/ebl: Basic analysis using R

#'Perform anova repeatly
#'
#'using use data
#'@param data:input data, grpcol:column number of group, varcol: column of number of variable
#'@return anova F.value and P.value each variable
#'@export
anova_stat <- function(data, grpcol=1, varcol=3){
  var=as.vector(unique(data[,varcol]))
  d=data.frame()
  for (i in 1:length(var)) {
    temp=subset(data, data$variable==var[i])
    head(temp)
    stat=aov(temp$value~temp[,grpcol])
    a=as.data.frame(matrix(unlist(summary(stat)),nrow=1))
    new_d=data.frame(Var=var[i],F.value=a[,7],P.value=a[,9])
    d=rbind.data.frame(d,new_d)
  }
  write.table(d,file = "temp_aov.txt")
  d=read.table("temp_aov.txt")
  return(d)
}