Get the top features associated with the RDA

Share:

Description

Get a list of the features significantly associated to the first two RDA components

Usage

1
topRDAhits(object, pval = 0.05)

Arguments

object

AnalysisRegionResults

pval

numeric with the p-value threshold. Only features with a p-values below this threshold will be shown.

Value

data.frame with the features, the component, the correlation and the p-value

Examples

1
2
3
4
5
6
7
if (require(minfiData) & require(GenomicRanges)){
set <- prepareMethylationSet(getBeta(MsetEx), pheno = pData(MsetEx))
range <- GenomicRanges::GRanges(seqnames=Rle("chrY"), 
ranges = IRanges(3000000, end=12300000))
rangeNoSNPs <- DARegionAnalysis(set, variable_names = "sex", range = range)
topRDAhits(rangeNoSNPs)
}

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.