CallCNVs-ExomeDepth-method: CallCNVs

Description Usage Arguments Details Value

Description

Call CNV data from an ExomeDepth object.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## S4 method for signature 'ExomeDepth'
CallCNVs(
  x,
  chromosome,
  start,
  end,
  name,
  transition.probability = 1e-04,
  expected.CNV.length = 50000
)

Arguments

x

An ExomeDepth object

chromosome

Chromosome information for each exon (factor).

start

Start (physical position) of each exon (numeric, must have the same length as the chromosome argument).

end

End (physical position) of each exome (numeric, must have the same length as the chromosome argument).

name

Name of each exon (character or factor).

transition.probability

Transition probability of the hidden Markov Chain from the normal copy number state to either a deletion or a duplication. The default (0.0001) expect approximately 20 CNVs genome-wide.

expected.CNV.length

The expectation for the length of a CNV. This value factors into the Viterbi algorithm that is used to compte the transition from one state to the next, which depends on the distance between exons.

Details

The function must be called on an ExomeDepth object. Likelihood data must have been pre-computed which should have been done by default when the ExomeDepth object was created.

This function fits a hidden Markov model to the read depth data with three hidden states (normal, deletion, duplication).

Value

The same ExomeDepth object provided as input but with the slot CNVcalls containing a data frame with the output of the calling.


ExomeDepth documentation built on Jan. 9, 2020, 1:09 a.m.