oncomodel (version 1.0)

leafset.prob: Probability of a Set of Leaves

Description

Computes the probability of a set of leaves (chromosomal aberrations) in an oncogenetic tree model.

Usage

leafset.prob(leafset, y)

Arguments

leafset

a character vector of one or more leaves (chromosomal aberrations).

y

the tree model as obtained from MLtopology.

Value

The probability of exactly a given set of leaves (chromosomal aberrations) in the tree.

References

von Heydebreck A, Gunawan B, Fuezesi L. 2004. Maximum likelihood estimation of oncogenetic tree models. Biostatistics 5:545-556.

Examples

  data(kidney)
  leafset.prob(c("-3|-3p", "+5|+5q"), kidney$res)