Top quantile expression genes

Select genes residing in the top quantile according to the mean of their expression across the stages

Usage

TopExpressionGenes(ExpressionSet, p = 0.99)

Arguments

ExpressionSet

A standard ExpressionSet

p

The quantile probability. Default is p = .99

Value

a character vector containing the gene ids residing in the top expression quantile

Author

Stefan Manolache

Examples

# reading a standard PhyloExpressionSet
data(PhyloExpressionSetExample)

# select genes with highest variance (top 2%)  
genes.top_expression <- TopExpressionGenes(PhyloExpressionSetExample, p=.98)
#> Error in TopExpressionGenes(PhyloExpressionSetExample, p = 0.98): could not find function "TopExpressionGenes"

# remove top genes from the PhyloExpressionSet
PhyloExpressionSet.top_removed <- subset(PhyloExpressionSetExample, 
                                         !(GeneID %in% genes.top_expression))
#> Error in eval(e, x, parent.frame()): object 'genes.top_expression' not found

# plot TAI of set with removed quantile
PlotSignature(ExpressionSet = PhyloExpressionSet.top_removed,
                measure       = "TAI", 
                TestStatistic = "FlatLineTest",
                xlab          = "Ontogeny", 
                ylab          = "TAI" )
#> Plot signature: ' TAI ' and test statistic: ' FlatLineTest ' running  1000  permutations.
#> Error: object 'PhyloExpressionSet.top_removed' not found