This function will add co-accessibility scores to peaks in a given ArchRProject
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") )
The name of the
A vector containing the dimensions from the
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
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
A character vector of cellNames to compute coAccessibility on if desired to run on a subset of the total cells.
The number of k-nearest neighbors to use for creating single-cell groups for correlation analyses.
The number of k-nearest neighbor groupings to test for passing the supplied
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.
The maximum allowable distance in basepairs between two peaks to consider for co-accessibility.
The total insertion counts from the designated group of single cells is summed across all relevant peak regions from
A boolean value indicating whether to log2 transform the single-cell groups prior to computing co-accessibility correlations.
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.
The number of threads to be used for parallel computing.
A boolean value that determines whether standard output should be printed.
The path to a file to be used for logging ArchR output.