problem with mcq.compute.anova
normalisation
SC = mcq.compute.normalization(SC) normalizing samples using the percent method
ANOVA for the factors of interest
Run anova without interaction term, if you want to perform tukey test with counting data
SC_ANOVA <- mcq.compute.anova(SC, flist=c("gtt"),inter=FALSE) Data are not normalized : computing model with glm function Computing Anova Error in apply(resultanov[, -which(colnames(resultanov) == "proba.good.model")], : dim(X) must have a positive length