data("IlluminaBodymap")
data("IlluminaBodymapMeta")
set.seed(1000) #fixed randomness
isoform_ids <- Name_to_RefSeq(c("HIF1A","EGFR","MYH7","CD33","BRCA2"))
rows_of_interest <- unique( c( isoform_ids ,
sample(rownames(IlluminaBodymap),
size=500,replace=FALSE)))
G <- GSEPD_INIT(Output_Folder="OUT",
finalCounts=round(IlluminaBodymap[rows_of_interest , ]),
sampleMeta=IlluminaBodymapMeta,
COLORS=c("green","black","red"))
G <- GSEPD_ChangeConditions( G, c("A","B")) #set testing groups first!
# G <- GSEPD_Process( G ) #would run DESeq2 and GOSeq and GSEPD comparing conditions A and BRun the code above in your browser using DataLab