aroma.cn (version 1.6.1)

callAllelicBalanceByBAFs.PairedPSCBS: Calls regions that are in allelic balance

Description

Calls regions that are in allelic balance from the allele B fractions (BAF).

Usage

# S3 method for PairedPSCBS
callAllelicBalanceByBAFs(fit, maxScore="auto", ..., force=FALSE, cache=FALSE,
  verbose=FALSE)

Value

Returns a PairedPSCBS fit object where columns for allelic imbalance scores and p-values as well as allelic balance calls are added.

Arguments

fit

A PairedPSCBS fit object as returned by segmentByPairedPSCBS.

maxScore

A positive double threshold. If "auto", the threshold is estimated empirically.

...

Not used.

force

If TRUE, an already called object is skipped, otherwise not.

verbose

See Verbose.

Author

Henrik Bengtsson, Pierre Neuvial

See Also

Internally, testAllelicBalanceByBAFs() is used.

Note that this AB caller differs from the default one in the PSCBS package, cf. callAB.PairedPSCBS.

Examples

if (Sys.getenv("_R_CHECK_FULL_") != "" && require("PSCBS")) {

# Load example ASCN data
data <- PSCBS::exampleData("paired.chr01");

# PSCBS segmentation
fit <- segmentByPairedPSCBS(data, verbose=-10);

# Normalize
fitN <- normalizeBAFsByRegions(fit, verbose=-10);

fitN2 <- callAllelicBalanceByBAFs(fitN, verbose=-10);
print(fitN2);

} # if (require("PSCBS"))