getASB {BaalChIP}R Documentation

Method getASB

Description

Method getASB

getASB identifies allele-specific binding events using a bayesian framework.

Usage

getASB(.Object, Iter = 5000, conf_level = 0.95, cores = 4,
  RMcorrection = TRUE, RAFcorrection = TRUE, verbose = TRUE)

## S4 method for signature 'BaalChIP'
getASB(.Object, Iter = 5000, conf_level = 0.95,
  cores = 4, RMcorrection = TRUE, RAFcorrection = TRUE, verbose = TRUE)

Arguments

.Object

An object of the BaalChIP class.

Iter

Maximum number of iterations (default 5000).

conf_level

Confidence interval in the estimated allelic ratio (default 0.95).

cores

number of cores for parallel computing (default is 4).

RMcorrection

Logical value indicating if reference mapping (RM) bias should be applied (default TRUE). If FALSE will not correct for reference allele mapping bias. If TRUE will estimate the RM bias from the overall reference allele proportion.

RAFcorrection

Logical value indicating if relative allele frequency (RAF) bias correction should be applied (default TRUE). If TRUE will read RAF values for each variant from hets files (RAF column name). If FALSE will not correct for relative allele frequency bias.

verbose

logical. If TRUE reports extra information on the process

Value

An updated BaalChIP object with the slot ASB containing variants identified as allele-specific.

Author(s)

Wei Liu, Ke Yuan, Ines de Santiago

See Also

summaryASB, BaalChIP.report

Examples

setwd(system.file('test',package='BaalChIP'))
samplesheet <- 'exampleChIP.tsv'
hets <- c('MCF7'='MCF7_hetSNP.txt', 'GM12891'='GM12891_hetSNP.txt')
res <- BaalChIP(samplesheet=samplesheet, hets=hets)
res <- alleleCounts(res, min_base_quality=10, min_mapq=15)
res <- mergePerGroup(res)
res <- getASB(res, cores=2)

#summary - number of significant ASB variants
summaryASB(res)

#report result
res <- BaalChIP.report(res)

[Package BaalChIP version 1.13.0 Index]