Convert BSseq object to edgeR methylation matrix

Usage

bsseq_to_edger(bsseq, regions = NULL)

Arguments

bsseq: the BSseq object.
regions: the regions to calculate log-methylation ratios over. If left NULL, ratios will be calculated per site.

Value

a matrix compatible with the edgeR differential methylation pipeline

Examples

methy <- system.file("methy_subset.tsv.bgz", package = "NanoMethViz", mustWork = FALSE)
bsseq <- methy_to_bsseq(methy)
#> [2025-04-18 12:13:33] creating intermediate files...
#> [2025-04-18 12:13:33] parsing chr11...
#> [2025-04-18 12:13:33] parsing chr12...
#> [2025-04-18 12:13:33] parsing chr18...
#> [2025-04-18 12:13:33] parsing chr5...
#> [2025-04-18 12:13:33] parsing chr7...
#> [2025-04-18 12:13:33] parsing chrX...
#> [2025-04-18 12:13:33] samples found: B6Cast_Prom_3_cast B6Cast_Prom_3_bl6 B6Cast_Prom_2_cast B6Cast_Prom_2_bl6 B6Cast_Prom_1_cast B6Cast_Prom_1_bl6 
#> [2025-04-18 12:13:33] creating bsseq object...
#> [2025-04-18 12:13:33] reading in parsed data...
#> [2025-04-18 12:13:33] constructing matrices...
#> [2025-04-18 12:13:34] done
edger_mat <- bsseq_to_edger(bsseq)