R/hcClimb.R

Defines functions hcClimb

hcClimb<-function(hc,minsize,minshare){
  hard2soft<-matrix(nrow=2,ncol=0,dimnames=list(c("hard","soft"),NULL))
  for(cnode in hc$clone_nodes)if(hc$nodesize[cnode]>minsize){
    nodenow<-cnode
    ancestor<-cnode
    while(nodenow<nrow(hc$merge)){
      nodenow<-row(hc$merge)[hc$merge==nodenow]
      if(hc$count_pins_share[nodenow]>=minshare)ancestor<-nodenow
    }
    hard2soft<-cbind(hard2soft,c(cnode,ancestor))
  }
  hard2soft
}
JunyanSong/SCclust documentation built on April 16, 2022, 8:44 p.m.