methylDiff_metagene: Make a metagene from mean percent methylation differences.

In PCBS: Principal Component BiSulfite

Make a metagene from mean percent methylation differences.

Description

Uses mean binned percent methylation differences across a set of regions to draw a metagene.

Usage

methylDiff_metagene(chromDictMethObj, regions, bin, title,
xaxis, yaxis, return.data, linecol, value)

Arguments

chromDictMethObj	chromDictMeth() output object
regions	A three-column data.frame containing a set of regions to test. Columns = chrom, start, end.
bin	integer, number of bins to use in metagenes. Default=100.
title	Output plot title
xaxis	Output plot x-axis title
yaxis	Output plot y-axis title
return.data	T/F, whether to return a plot, or data that can be run with plot_metagene() or multiple_metagenes().
linecol	Colour for line, auto="red"
value	Name of the plotted metric in chromDictMethObj. Only needs to be set explicity for abnormal use cases where chromDictMethObj contains a non-rank value output by chromDict().

Value

If return.data=F, returns a grob containing a metagene plot. Otherwise, returns a list of two data.frames containing metagene and metagene standard error plotting information.

Examples

chromDictMethylDiff <- chromDictMeth(eigen, c("trt", "ctl"))
regions <- data.frame(chr=c("chr3", "chr3", "chr1"),
                      s=c(4920450, 3961576, 300000),
                      e=c(4923267, 3963805, 302900),
                      ID=c("Hypo-DMR", "partial Hyper-DMR", "random"))
 
methylDiff_metagene(chromDictMethylDiff, regions)
                     

PCBS documentation built on Sept. 11, 2024, 6:11 p.m.