D <- apply(as.matrix(dat), 2, as.numeric)
d <- dat %>% group_by(Group) %>% summarise_if(is.numeric, mean, na.rm = TRUE)
HNSC_TT <- HNSC[ , grepl( "01A" , names(HNSC) ) ]
package Rsubread for fastq to bam/count
package DESeq2 for differently gene expression
package DOSE for GSEA
package yarrr
library(Hmisc)
res <- rcorr(as.matrix(data))
