R/thetaGT.R

Defines functions thetaGT

Documented in thetaGT

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)
}

Try the synRNASeqNet package in your browser

Any scripts or data that you put into this service are public.

synRNASeqNet documentation built on May 2, 2019, 6:01 a.m.