Learn R Programming

PCBS (version 0.1.1)

Get_Novel_DMRs: Call DMRs from WGBS data.

Description

DMR Calling.

Usage

Get_Novel_DMRs(ranks, nSeeds, chromDictObj, DMR_resolution, 
QueryLimit, minCpGs, minZ, perms)

Value

Returns a data.frame of all novel DMRs.

Arguments

ranks

Rank data frame from getPCRanks.

nSeeds

Integer, number of input seeds for DMR expansion.

chromDictObj

chromDict() output. If null, chromDict() is run internally.

DMR_resolution

Automatic=NULL. Integer, number of bases to increase the DMR by with each expansion. If NULL, QueryLimit/25.

QueryLimit

Automatic=5000. Maximum DMR expansion size (bp)

minCpGs

Automatic=15. Minimum CpGs in a DMR region, regions with fewer CpGs will be discarded.

minZ

Automatic=1. Absolute Z score threshold for DMR calling; internal value. Not recommended to play with this setting.

perms

Automatic=1000. Number of permutations to use when defining the null distribution. Increasing this value largely influences computational time with minimal return

Examples

Run this code
ranks <- getPCRanks(eigen, IDs = c("trt", "ctl"), PC = 1)
DMRs <- Get_Novel_DMRs(ranks, 2940, minCpGs=10)

Run the code above in your browser using DataLab