R/tree-K.r

In MBCluster.Seq: Model-Based Clustering for RNA-seq Data

Defines functions tree.K

Documented in tree.K

tree.K=function(k,k1,k2){
  nG=length(k)
  nK=length(k1)
  K=matrix(0,nG,nK)
  K[,1]=k
  for( i in 2:nK){
   k0=K[,i-1]
   k0[k0==k2[i-1]]=k1[i-1]
   K[,i]=k0
 } 
 return(K)
}

Try the MBCluster.Seq package in your browser