saveUMRLMRSegments: Save UMR and LMR segments

In MethylSeekR: Segmentation of Bis-seq data

Description

Save UMR and LMRs segments in rds format and/or as tab-delimited file

Usage

saveUMRLMRSegments(segs, GRangesFilename = NULL, TableFilename = NULL)

Arguments

segs	GRanges object containing the UMR/LMR segmentation. Return value of the segmentUMRsLMRs function (see example).
GRangesFilename	Filename of the GRanges object.
TableFilename	Filename of the UMR/LMR table.

Value

No return value. Only one filename is required.

Author(s)

Lukas Burger [e-mail]

Examples

library(MethylSeekR)

# get chromosome lengths
library("BSgenome.Hsapiens.UCSC.hg18")
sLengths=seqlengths(Hsapiens)

# read methylation data
methFname <- system.file("extdata", "Lister2009_imr90_hg18_chr22.tab",
package="MethylSeekR")
meth.gr <- readMethylome(FileName=methFname, seqLengths=sLengths)

FDR.cutoff <- 5 
m.sel <- 0.5 
n.sel <- 3

#segment UMRs and LMRs, assuming no PMDs
UMRLMRsegments.gr <- segmentUMRsLMRs(m=meth.gr, meth.cutoff=m.sel,
nCpG.cutoff=n.sel, num.cores=1,
myGenomeSeq=Hsapiens, seqLengths=sLengths)

#save UMRs and LMRs
saveUMRLMRSegments(segs=UMRLMRsegments.gr, GRangesFilename="UMRsLMRs.gr.rds",
 TableFilename="UMRsLMRs.tab")

MethylSeekR documentation built on Nov. 8, 2020, 6:57 p.m.