multiscaleGroupAverage: Create a list of GRangesLists for each step in the multiscale...

In huishenlab/bisplotti: Plotting package for bisulfite sequencing data

Create a list of GRangesLists for each step in the multiscale plot

Description

Output of this function can be passed directly to multiscaleMethylationPlot

Usage

multiscaleGroupAverage(means_dir, samples, config, which)

Arguments

means_dir	List of length nSamples. Each element of means.list is itself a list of GRanges of length, nSteps
samples	character vector of samples to include. These samples must match the sample column from samples.tsv in the multiscale plot workflow
config	Config YAML from the multiscale plot workflow
which	GRanges for rtracklayer::import

Value

a GRangesList with the average values

Examples

# Not a great example, but it works
files.loc <- system.file("extdata", package="bisplotti")

grp_avg <- multiscaleGroupAverage(
    files.loc,
    "Heyn_2012_100yr_subset",
    paste0(files.loc, "/config_example.yaml"),
    NULL
)

huishenlab/bisplotti documentation built on Sept. 20, 2023, 10:13 p.m.