R/likelihood.R

likelihood <-
function(thetast,fobs,merge2){
fobss=fobs;thetass=thetast;merge22=merge2
fst<-gn(thetass,merge22)
if(min(fst)<=0){
return(lik<-4)
}
lik<-sum(fobss*log(fobss/fst))
lik
}

Try the DNAseqtest package in your browser

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

DNAseqtest documentation built on May 1, 2019, 9:51 p.m.