MethAlignNW: Summary of methylation states

Description Usage Arguments Details Value Author(s) Examples

Description

Summarize the methylation states after calculating pairwise alignments of each examined sequences and the genomic sequence

Usage

1
MethAlignNW(refSeq, QCdata, alignment)

Arguments

refSeq

String; Genomic sequence as String format

QCdata

Data frame; Names and paths of analysed sequences after quality control

alignment

If TRUE, alignments are included in summery, else not included

Details

Given aligned sequences after quality control, the function returns a list object with the following data: sequences name, methylation state on CpG position, start and end position of alignments and length of genomic sequence. The data includes the core information for the exploratory analysis and visualizations.

Value

Returns a summery on sequence alignments and methylation states

Author(s)

Arie Zackay <arie.zackay@mail.huji.ac.il>, Christine Steinhoff <steinhof@molgen.mpg.de>

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
## In order to use the following example 
## make sure that you have writing permission under R.home()
## directory. If you do not have permission choose your own path. 
#dir.create(file.path(R.home(component="home"),"/BiqAnalyzer"))
BiqAnalyzer_path <- file.path(tempdir(), "BiqAnalyzer")
dir.create(BiqAnalyzer_path) 
makeLocalExpDir(dataPath="/examples/BiqAnalyzer", localDir=BiqAnalyzer_path)
datameth <- MethDataInput(file.path(BiqAnalyzer_path, "PathFileTab.txt"))
refseq <- selectRefSeq(file.path(BiqAnalyzer_path, "Master_Sequence.txt"))
QCdata <- MethylQC(refseq, datameth)
methData <- MethAlignNW( refseq , QCdata)

methVisual documentation built on April 28, 2020, 7:08 p.m.