MethyEnrich: Quality control of MBDcap-seq datasets
In Liuy12/MBDDiff: MBDDiff for processing MBDcap-seq datasets

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

Generate graph of enrichment of GC content grouped by methylation levels

1	MethyEnrich(bin_count, fa)

`bin_count`	A vector containing number of reads mapped to each bin. See vignettes for details
`fa`	Absolute path of sequence file annotation for your interested organism

In order to use this function, you will firstly need to create a sliding window across the genome and measure number of reads mapped to each window. Please refer to vignettes for details

Yuanhang Liu

https://github.com/Liuy12/MBDDiff

bedtoolsr

  ## Not run: 
    whole_genome_bin_count <- bedtoolsr(bed, bam)
    whole_genome_bin_count_TPM <- CalculateTPM(whole_genome_bin_count[[1]], whole_genome_bin_count[[2]])
    MethyEnrich(whole_genome_bin_count_TPM, fa)
    
## End(Not run)