Displaying 1 result from an estimated 1 matches for "aoftest".
Did you mean:
adftest
2016 Apr 05
0
use one way ANOVA to select genes
...######
datafile <- read.csv("brcaResultsclust.csv",header = T,row.names = 1)
datafile <- t(datafile)
groups <- factor(rep(c("c1","c2","c3","c4","c5"),c(116,83,28,80,24))) # here one factor with 5 groups
# define ANOVA function
aoftest <- function(data ,pvalue=0.05) {
m<-data.frame(groups, data);
anova(aov(data ~ groups, m))
}
anovaresultstest <- apply(datafile, 1, aoftest)
#############################################################################################
is the design factor and the function for ANOVA...