data("IlluminaBodymap")
data("IlluminaBodymapMeta")
head(IlluminaBodymap)
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=50,replace=FALSE)))
G <- GSEPD_INIT(Output_Folder="OUT",
finalCounts=round(IlluminaBodymap[rows_of_interest , ]),
sampleMeta=IlluminaBodymapMeta,
COLORS=c("green","black","red"))
# G <- GSEPD_ProcessAll( G ) #would run across all pairs of G$Condition
Run the code above in your browser using DataLab