EBSeq/R/PoolMatrix.R

   1 PoolMatrix <-
   2 function(Data,reads,type)
   3 {
   4 poolnames=names(Data)
   5 poolM=NULL
   6 for (po in 1:8)
   7         poolM=cbind(poolM,Data[[po]][,1])
   8 rownames(poolM)=rownames(Data[[1]])
   9 colnames(poolM)=poolnames
  10
  11 #poolValue=poolM*reads
  12 poolValue=poolM
  13 for (col in 1:8)
  14         poolValue[,col]=poolM[,col]*reads[col]
  15 poolValue=round(poolValue)
  16 if (type=="G")
  17         {
  18                 poolM=cbind(Data[[1]][,2],poolM)
  19                 poolValue=cbind(Data[[1]][,2],poolValue)
  20                 colnames(poolM)=c("Groups",poolnames)
  21                 colnames(poolValue)=c("Groups",poolnames)
  22         }
  23 poolOutput=list(poolM=poolM,poolValue=poolValue)
  24 }
  25