Read Bismark Coverage Files
Read Bismark coverage files containing methylated and unmethylated read counts for CpG loci and create DGEList.
readBismark2DGE(files, sample.names=NULL, readr=TRUE, verbose=TRUE)
files |
character vector of file names. |
sample.names |
character vector of sample names. If |
readr |
logical. If |
verbose |
logical. If |
This function reads tab-delimited coverage files output by Bismark software.
Counts from multiple files are collated into a DGEList object.
A DGEList object with a row for each unique genomic loci found in the files and two columns (containing methylated and unmethylated counts) for each sample.
This function represents genomic loci as integers, so the largest locus position must be less than the maximum integer in R (about 2e9).
The number of chromosomes times the largest locus position must be less than 1e16.
Gordon Smyth
Chen, Y, Pal, B, Visvader, JE, Smyth, GK (2017). Differential methylation analysis of reduced representation bisulfite sequencing experiments using edgeR. F1000Research 6, 2055. https://f1000research.com/articles/6-2055
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.