eq (bss, testCovariate = "group" , minNumRegion = 20 ) saveRDS (regions, "data/regions.rds" ) } else { regions <- readRDS ( "data/regions.rds" ) } A number of regions will be produced by dmrseq , including regions that were aggregated because they contained significant CpG sites but subsequently failed to reach statistical significance at a region level.
← all excerpts
A streamlined workflow for long-read DNA methylation analysis with NanoMethViz and Bioconductor.
1
—
—