savePMDSegments: Save PMD segments

Description Usage Arguments Value Author(s) Examples

View source: R/savePMDSegments.R

Description

Save PMD segments in rds format and as tab-delimited file.

Usage

1
savePMDSegments(PMDs, GRangesFilename, TableFilename)

Arguments

PMDs

GRanges object containing the PMD segmentation. Return value of the segmentPMDs function (see example).

GRangesFilename

Filename of the GRanges object.

TableFilename

Filename of the PMD table.

Value

No return value.

Author(s)

Lukas Burger lukas.burger@fmi.ch

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
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)

#segment PMDs
PMDsegments.gr <- segmentPMDs(m=meth.gr, chr.sel="chr22",
seqLengths=sLengths)

#save PMD segments
savePMDSegments(PMDs=PMDsegments.gr, GRangesFilename="PMDs.gr.rds",
 TableFilename="PMDs.tab")

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