Plot gene set details
This function is to generate a two-panel plot showing detailed information of the gene set specified. One panel is showing the running enrichment scores and the position where the ES appear. The other panel shows the significance level of the ES, comparing with permutation ESs.
a SeqGeneSet object after running
the i-th gene set in the SeqGeneSet object.
the gene score vector containing gene scores for each gene.
whether to save the plot to PDF file; if yes, provide the name of the PDF file.
writeSigGeneSet, which writes the detailed gene set information to a file or to the screen.
Xi Wang, firstname.lastname@example.org
1 2 3 4 5 6 7 8 9 10
data(DEscore, package="SeqGSEA") data(DSscore, package="SeqGSEA") gene.score <- geneScore(DEscore, DSscore, method="linear", DEweight = 0.3) data(DEscore.perm, package="SeqGSEA") data(DSscore.perm, package="SeqGSEA") gene.score.perm <- genePermuteScore(DEscore.perm, DSscore.perm, method="linear", DEweight=0.3) data(GS_example, package="SeqGSEA") GS_example <- GSEnrichAnalyze(GS_example, gene.score, gene.score.perm) topGeneSets(GS_example, n=5) plotSigGeneSet(GS_example, 9, gene.score) # 9th gene set is the most significant one.