addCoAccessibility: Add Peak Co-Accessibility to an ArchRProject
In haibol2016/ArchR_debug: Analyzing single-cell regulatory chromatin in R.

addCoAccessibility

R Documentation

Add Peak Co-Accessibility to an ArchRProject

Description

This function will add co-accessibility scores to peaks in a given ArchRProject

Usage

addCoAccessibility(
  ArchRProj = NULL,
  reducedDims = "IterativeLSI",
  dimsToUse = 1:30,
  scaleDims = NULL,
  corCutOff = 0.75,
  cellsToUse = NULL,
  k = 100,
  knnIteration = 500,
  overlapCutoff = 0.8,
  maxDist = 1e+05,
  scaleTo = 10^4,
  log2Norm = TRUE,
  seed = 1,
  threads = getArchRThreads(),
  verbose = TRUE,
  logFile = createLogFile("addCoAccessibility")
)

Arguments

`ArchRProj`	An `ArchRProject` object.
`reducedDims`	The name of the `reducedDims` object (i.e. "IterativeLSI") to retrieve from the designated `ArchRProject`.
`dimsToUse`	A vector containing the dimensions from the `reducedDims` object to use in clustering.
`scaleDims`	A boolean value that indicates whether to z-score the reduced dimensions for each cell. This is useful for minimizing the contribution of strong biases (dominating early PCs) and lowly abundant populations. However, this may lead to stronger sample-specific biases since it is over-weighting latent PCs. If set to `NULL` this will scale the dimensions based on the value of `scaleDims` when the `reducedDims` were originally created during dimensionality reduction. This idea was introduced by Timothy Stuart.
`corCutOff`	A numeric cutoff for the correlation of each dimension to the sequencing depth. If the dimension has a correlation to sequencing depth that is greater than the `corCutOff`, it will be excluded from analysis.
`cellsToUse`	A character vector of cellNames to compute coAccessibility on if desired to run on a subset of the total cells.
`k`	The number of k-nearest neighbors to use for creating single-cell groups for correlation analyses.
`knnIteration`	The number of k-nearest neighbor groupings to test for passing the supplied `overlapCutoff`.
`overlapCutoff`	The maximum allowable overlap between the current group and all previous groups to permit the current group be added to the group list during k-nearest neighbor calculations.
`maxDist`	The maximum allowable distance in basepairs between two peaks to consider for co-accessibility.
`scaleTo`	The total insertion counts from the designated group of single cells is summed across all relevant peak regions from the `peakSet` of the `ArchRProject` and normalized to the total depth provided by `scaleTo`.
`log2Norm`	A boolean value indicating whether to log2 transform the single-cell groups prior to computing co-accessibility correlations.
`seed`	A number to be used as the seed for random number generation required in knn determination. It is recommended to keep track of the seed used so that you can reproduce results downstream.
`threads`	The number of threads to be used for parallel computing.
`verbose`	A boolean value that determines whether standard output should be printed.
`logFile`	The path to a file to be used for logging ArchR output.