### This a example to run MixEnrich.R on a single patient.

# load in the MixEnrich function
source("./MixEnrich.R")

# load in gene set annotations
GO_anno <- read.delim('./go_bp_filtered15-500.txt')

# load in the example dataset
dat <- read.csv('./exampleData_MixEnrich.csv',row.names=1) 
# run MixEnrich
res_singlePatient <- mixEnrich(abs.logFC = absLogFC(dat[,1],dat[,2]), gene.symbol = rownames(dat), GeneSet.annotation = GO_anno, FDR.method = 'BY')

head(res_singlePatient)