R/likelihoodpercond.R

likelihoodpercond <-
function(theta, mu, ys,const=1){
  #get the likelihood of a certain category
like <-  1/const * dnbinom(ys,mu=mu,size=1/theta)
like
}
lunge111/intSEQ documentation built on May 20, 2019, 9:38 a.m.