analyzeCpGs: Analyze CG methylation

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/calcPercentConversion.R

Description

Function to determine percent methylation for all CGs from input fragmentation

Usage

1
analyzeCpGs(fragments, peaks, method = c("weighted", "proportion"))

Arguments

fragments

List of MassArrayFragment objects

peaks

List of MassArrayPeak objects comprising spectral data for a given assay

method

Specifies which algorithm to use when calculating percent methylation (either "weighted" or "proportion")

Details

Wrapper function for calcMeth(), takes fragmentation pattern and spectral data as input and applies percent methylation calculation for all CG-containing, non conversion control fragments

Value

Returns a list of numerical values corresponding to percent methylation for each CG dinucleotide, with 0

Author(s)

Reid F. Thompson (rthompso@aecom.yu.edu), John M. Greally (jgreally@aecom.yu.edu)

See Also

See Also calcMeth

Examples

1
2
3
data(MassArray.example.data)
cpg.data <- analyzeCpGs(MassArray.example.data$fragments.T, MassArray.example.data$samples[[1]]$peaks, method="weighted")
barplot(cpg.data, xlab="CpG (Number)", ylim=c(0,1), ylab="Methylation (Percent)")

MassArray documentation built on Nov. 8, 2020, 5:16 p.m.