mergeAllRDSFiles: Merge nucleosome information

In RJMCMCNucleosomes: Bayesian hierarchical model for genome-wide nucleosome positioning with high-throughput short-read data (MNase-Seq)

Description

Merge nucleosome information present in multiple RDS files.

Usage

mergeAllRDSFiles(arrayOfFiles)

Arguments

arrayOfFiles

a array, the name of each file that must be used to merge nucleosome information.

Value

a list of class "rjmcmcNucleosomesMerge" containing:

• k a integer, the number of nucleosomes.

• mu a GRanges containing the positions of the nucleosomes.

Author(s)

Pascal Belleau, Astrid Deschenes

Examples

## Loading two files containing nucleosomes informations for two sections of
## the same chromosome
file_1 <- dir(system.file("extdata", package = "RJMCMCNucleosomes"),
                pattern = "RJMCMC_seg_01.RDS",
                full.names = TRUE)

file_2 <- dir(system.file("extdata", package = "RJMCMCNucleosomes"),
                pattern = "RJMCMC_seg_02.RDS",
                full.names = TRUE)

## Merging nucleosomes informations from the two files
result <- RJMCMCNucleosomes:::mergeAllRDSFiles(c(file_1, file_2))

RJMCMCNucleosomes documentation built on Nov. 8, 2020, 8:20 p.m.