library(MoRad) #loading package
head(data.GBM) #sample dataset 1-TCGA-GBM
head(data.LGG) #sample dataset 2-TCGA-LGG
MoRad.scores <- calMoRad(data.GBM) #Profile all available molecular features with MoRad models
calRO(names(rfeature_list)[1:100]) #this will save a spreadsheet(.csv) in your work directory
RSEA.result <- calRSEA(rfeature_list) #RSEA enrichment analysis
plotRSEA(rfeature_list,nrow=3,ncol=3,color.line="#00468BFF") #plot RSEA enrichment outcomes; change 'color.line' value for line color