mergeAllRDSFiles: Merge nucleosome information

In ArnaudDroitLab/RJMCMC: Bayesian hierarchical model for genome-wide profiling of nucleosome positions based on high-throughput short-read data (MNase-Seq data)

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 vector of numeric of length k, the positions of the nucleosomes.

• sigmaf a vector of numeric of length k, the variance of the forward reads for each nucleosome.

• sigmar a vector of numeric of length k, the variance of the reverse reads for each nucleosome.

• delta a vector of numeric of length k, the distance between the maxima of the forward and reverse reads position densities for each nucleosome.

Author(s)

Pascal Belleau, Astrid Deschenes

Examples

## Loading two files containing nucleosomes informations for two sections of
## the same chromosome
file_100 <- dir(system.file("extdata", package = "RJMCMC"),
                pattern = "yeastRes_Chr1_Seg_100.rds",
                full.names = TRUE)

file_101 <- dir(system.file("extdata", package = "RJMCMC"),
                pattern = "yeastRes_Chr1_Seg_101.rds",
                full.names = TRUE)

## Merging nucleosomes informations from the two files
result <- RJMCMC:::mergeAllRDSFiles(c(file_100, file_101))

ArnaudDroitLab/RJMCMC documentation built on May 5, 2019, 7:06 a.m.