getBSABsignal: Estimate A/B compartments from BS-seq data

In JordanVeldboom/compartmap: A/B compartment inference from ATAC-seq and methylation array data

Description

getBSABsignal returns estimated A/B compartments from BS-seq data.

Usage

getBSABsignal(
  obj,
  res = 1000000,
  parallel = TRUE,
  chr = NULL,
  targets = NULL,
  preprocess = TRUE,
  cores = 2,
  bootstrap = TRUE,
  num.bootstraps = 1000,
  genome = c("hg19", "hg38", "mm9", "mm10"),
  other = NULL,
  group = FALSE,
  boot.parallel = TRUE,
  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
preprocess	Whether to preprocess the bs-seq data prior to compartment inference
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

data("meth_wgbs_450k_chr14", package = "compartmap")
wgbs_compartments <- getWGBSABsignal(wgbs.data.chr14, parallel=F, chr="chr14", bootstrap=F, genome="hg19", wgbs.type="hm450")

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