MethyBayes (version 0.1.0)

seqdata: This is a function to merge BS-seq data from multiple replicates.

Description

Return read counts and methylated reads counts for sites commonly mapped by all replicates.

Usage

seqdata(cases, controls)

Arguments

cases

A list of data frames for multiple biological replicates under treatment experiment.

controls

A list of data frames for multiple biological replicates under control experiment.

Value

A list of read counts, methylated reads and commonly mapped sites.

Examples

data("case1", "case2", "control1", "control2")
read_counts <- seqdata(list(case1, case2), list(control1, control2))