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!
dds <- GSEPD_Export_DESeq(G)
print(dds)Run the code above in your browser using DataLab