get_hist_clon_excl {GeneAccord} | R Documentation |
Compute all values of how often gene pairs were clonally exclusive/all trees for a patient.
get_hist_clon_excl(clone_tbl)
clone_tbl |
A tibble containing the columns
'altered_entity', and then a column for each clone
in the tumor, e.g. 'clone1', 'clone2', 'clone3'. It also
contains the column 'tree_id', which specifies
which tree of the collection of tree inferences was used.
This tibble can be generated e.g. from the |
It computes a histogram of the following two values: Amomg all gene/pathway pairs in a patient, the number of trees in which the both entities of a pair are assigned to a clone at all, and the number of trees in which the pair is clonally exclusive.
A list with two vectors: The numbers of how often gene pairs were mutated across trees, and the numbers of how often they were clonally exclusive. The order of these two vectors is matching, i.e. the ith entry in each vector refers to the same gene pair.
Ariane L. Moore
clone_tbl <- dplyr::tibble( altered_entity=c(paste("gene", seq_len(10), sep="")), clone1=c(rep(0,10)), clone2=c(sample(c(0,1), 10, replace=TRUE)), clone3=c(sample(c(0,1), 10, replace=TRUE)), clone4=c(sample(c(0,1), 10, replace=TRUE)), tree_id=c(rep(5, 5), rep(10, 5)) ) get_hist_clon_excl(clone_tbl)