Oncotree-package {Oncotree}

Constructing and evaluating oncogenetic trees

Description

Oncogenetic trees are directed tree structures that model the process of occurrence of genetic alterations during carcinogenesis.

Details

A pure oncogenetic tree is a directed rooted tree T with a probability \pi (e) attached to each edge e such that for every vertex there is a unique directed path from the root to it along the edges of the tree. This tree generates observations on the presence/absence of genetic events the following way: each edge e is independently retained with probability \pi (e); the set of vertices that are still reachable from the root gives the set of the observed genetic events.

To describe random deviations from the pure tree model an error model is added.

Error model

1. The tumor develops according to the pure oncogenetic tree model

2. The presence/absence of each alteration is independently measured

3. If the alteration is present it is not observed with probability \epsilon_-.

   If the alteration is absent it is observed with probability \epsilon_+.

Author(s)

Lisa Pappas, Aniko Szabo

Maintainer: Aniko Szabo <[e-mail]>

References

[1] Desper R., Jiang F., Kallioniemi O.P., Moch H., Papadimitriou C.H., and Sch\"affer A.A. (1999) Inferring tree models for oncogenesis from comparative genome hybridization data. Journal of Computational Biology. 6m 37–51. [2] Szabo, A. and Boucher, K. (2002) Estimating an oncogenetic tree when false negative and positives are present. Mathematical Biosciences, 176/2, 219–236.

Examples

  data(ov.cgh)
  ov.tree <- oncotree.fit(ov.cgh)
  plot(ov.tree, edge.weights="estimated")