R/sampcna.R

In Canopy: Accessing Intra-Tumor Heterogeneity and Tracking Longitudinal and Spatial Clonal Evolutionary History by Next-Generation Sequencing

sampcna = function(tree) {
    t = nrow(tree$cna)
    cna = tree$cna
    cna.change = sample.int(sample.int(1, n = t), n = t)
    cna.edge = sample(2:nrow(tree$edge), size = length(cna.change), 
        replace = TRUE)
    cna[cna.change, 2:3] = tree$edge[cna.edge, 1:2]
    return(cna)
}