R/initialsna.R

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

initialsna = function(tree, sna.name) {
    sna.no = length(sna.name)
    sna.edge = sample(2:nrow(tree$edge), size = sna.no, replace = TRUE)
    sna.mat = cbind(1:sna.no, tree$edge[sna.edge, ])
    colnames(sna.mat) = c("sna", "sna.st.node", "sna.ed.node")
    rownames(sna.mat) = sna.name
    return(sna.mat)
}