scripts/prob_base.R

function(Qscore, count, offset){
  pb <- 1-(10^(-(Qscore-offset)/10))
  prob_count <- sum(pb*count)
  return(prob_count)
}

Try the VariantTools package in your browser

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

VariantTools documentation built on Nov. 8, 2020, 8:03 p.m.