Description Usage Arguments Value Author(s) See Also Examples
View source: R/lleading.edge.R
performs the leading edge analysis involving the running sum described in the GSEA methodology.
1 | lleading.edge(k, s2n)
|
k |
boolean vector of length = number of genes. TRUE if the nth gene/row of the expression matrix is found in the gene set to be tested. |
s2n |
scores between two classes : see snr.SNR or snr.FC |
return a list with two elements : essmin is the minimum enrichment score in the running sum; essmax is the maximum enrichment score in the running score.
maxime.tarabichi@ulb.ac.be
snr.SNR; snr.FC; GSEA.plot
1 2 3 4 5 6 7 8 9 10 | ## generates random scores
scores<-rnorm(10000)
## generate boolean vector with 50 genes present in the 500 highest
scores
k<-sample(which(rank(scores)<=500),50,rep=F)
## print the leading edge result of 50 random genes vs. 50 genes in top 500
lleading.edge(k,scores) ##top 500
lleading.edge(sample(k),scores) ##random
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.