callAllelicBalanceByBAFs.PairedPSCBS: Calls regions that are in allelic balance

Description Usage Arguments Value Author(s) See Also Examples

Description

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

Usage

1
2
3
## S3 method for class 'PairedPSCBS'
callAllelicBalanceByBAFs(fit, maxScore="auto", ..., force=FALSE, cache=FALSE,
  verbose=FALSE)

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.

Value

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

Author(s)

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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"))

aroma.cn documentation built on May 29, 2017, 7:56 p.m.