getATACABsignal: Estimate A/B compartments from ATAC-seq data

Description Usage Arguments Value Examples

View source: R/getATACABsignal.R

Description

getATACABsignal returns estimated A/B compartments from ATAC-seq data.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
getATACABsignal(
  obj,
  res = 1000000,
  parallel = FALSE,
  chr = NULL,
  targets = NULL,
  cores = 2,
  bootstrap = TRUE,
  num.bootstraps = 100,
  genome = c("hg19", "hg38", "mm9", "mm10"),
  other = NULL,
  group = FALSE,
  boot.parallel = FALSE,
  boot.cores = 2
)

Arguments

obj

Input SummarizedExperiment object

res

Compartment resolution in bp

parallel

Whether to run samples in parallel

chr

What chromosome to work on (leave as NULL to run on all chromosomes)

targets

Samples/cells to shrink towards

cores

How many cores to use when running samples in parallel

bootstrap

Whether we should perform bootstrapping of inferred compartments

num.bootstraps

How many bootstraps to run

genome

What genome to work on ("hg19", "hg38", "mm9", "mm10")

other

Another arbitrary genome to compute compartments on

group

Whether to treat this as a group set of samples

boot.parallel

Whether to run the bootstrapping in parallel

boot.cores

How many cores to use for the bootstrapping

Value

A RaggedExperiment of inferred compartments

Examples

1
2
data("groupATAC_raw_filtered_chr14", package = "compartmap")
atac_compartments <- getATACABsignal(filtered.data.chr14, parallel=F, chr="chr14", bootstrap=F, genome="hg19")

JordanVeldboom/compartmap documentation built on July 3, 2020, 6:32 p.m.