R/getSiblings.R

In ndukler/epiAllele: Tests for signals of variable rates of evolution between gene classes and lineages for probabalistic alleles

#' Compute transitions per branch
#'
#' Compute the number of expected transitions per branch between states based on the marginal joint probabilities
#' @param tr phylo object
#' @param node
#' @return The node numbers of sibbling nodes
#' @name getSiblings
#' @include rateModel-class.R
#' @rdname getSiblings
#' @examples
#' 
#' @export
getSiblings = function(tr,node){
  if(!inherits(tr,"phylo")) {
    stop("tree should be an object of class \"phylo\".")
  }
  parent=tr$edge[tr$edge[,2]==node,1]
  return(setdiff(tr$edge[tr$edge[,1]==parent,2],node))
}