Add a peak set to an ArchRProject — addPeakSet • ArchR

This function adds a peak set as a GRanges object to a given ArchRProject.

addPeakSet(
  ArchRProj = NULL,
  peakSet = NULL,
  genomeAnnotation = getGenomeAnnotation(ArchRProj),
  force = FALSE
)

Arguments

ArchRProj: An ArchRProject object.
peakSet: A GRanges object containing the set of regions that define all peaks in the desired peak set.
genomeAnnotation: The genomeAnnotation (see createGenomeAnnotation()) to be used for generating peak metadata such as nucleotide information (GC content) or chromosome sizes.
force: If a peakSet object has already been added to the given ArchRProject, the value of force determines whether or not to overwrite this peakSet.


# Get Test ArchR Project
proj <- getTestProject()

# Add PeakSet
addPeakSet(proj, peakSet = getPeakSet(proj), force = TRUE)