# NOT RUN {
library(tmod)
data(tmod)
data(orthologs)
disco <- disco.score(orthologs)
ord <- order(disco, decreasing = TRUE)
concordant <- tmodCERNOtest(toupper(orthologs$genes)[ord])
modtable <- percentGenes(orthologs, tmod[concordant$ID])
cbind(concordant, modtable)
# }
Run the code above in your browser using DataLab