R/thetaGT.R

thetaGT <-
function(cellCounts){
  ans <- thetaML(cellCounts)
  
  n <- sum(cellCounts)
  m1 <- sum(cellCounts == 1)
  if(m1 == n) m1 <- m1 - 1 #avoid (1 - m1/n) = 0
  
  ans <- (1 - m1/n)*ans
  return(ans)
}
lucgar/synRNASeqNet documentation built on May 21, 2019, 8:54 a.m.