oncomodel (version 1.0)

leafset.prob2: Probability of a Set of Leaves

Description

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

Usage

leafset.prob2(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 the given set of leaves (chromosomal aberrations) of 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.prob2(c("-3|-3p", "+5|+5q"), kidney$res)