Plot a Manhattan plot with the probe results

Share:

Description

Plot log p-value for each chromosome positions. Highlighting cpgs inside a range is allowed.

Usage

1
2
plotEWAS(object, variable = modelVariables(object)[[1]], range = NULL,
  main = paste("Manhattan plot of ", variable))

Arguments

object

AnalysisResults or AnalysisRegionResults

variable

Character with the variable name used to obtain the probe results. Note: model name should be used. Original variable name might not be valid.

range

GenomicRange whose cpgs will be highlighted

main

Character with the plot title.

Value

A plot is generated on the current graphics device.

Examples

1
2
3
4
5
6
if (require(minfiData)){
betas <- getBeta(MsetEx)[floor(seq(1, nrow(MsetEx), 10000)), ]
set <- prepareMethylationSet(betas, pheno = pData(MsetEx))
methyOneVar <- DAPipeline(set, variable_names = "sex", probe_method = "ls")
plotEWAS(methyOneVar)
}

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