mergeRDSFiles: Merge nucleosome information for selected RDS files.

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 RDS files into one object of class "rjmcmcNucleosomesMerge".

Usage

mergeRDSFiles(RDSFiles)

Arguments

RDSFiles

a array, the names of all RDS used to merge nucleosome information. The files must contain R object of class "rjmcmcNucleosomes" or "rjmcmcNucleosomesMerge".

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

## Use RDS files present in the RJMCMC package
RDSFiles <- dir(system.file("extdata", package = "RJMCMC"),
full.names = TRUE, pattern = "*rds")

## Merge nucleosomes info from RDS files present in directory
result <- mergeRDSFiles(RDSFiles)

## Print the number and the position of the nucleosomes
result$k
result$mu

## Class of the output object
class(result)

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