R/computeClusterParHastingsRatio.R

Defines functions computeClusterParHastingsRatio

computeClusterParHastingsRatio <-
function (cluster.par.new, cluster.par.curr, gene.data, n.time, n.rep, DATA.FORMAT=0)
{
	loglik.new = computeLogLik (cluster.par.new, gene.data, n.time, n.rep, DATA.FORMAT)
	loglik.curr = computeLogLik (cluster.par.curr, gene.data, n.time, n.rep, DATA.FORMAT)
	
	logHRatio = loglik.new - loglik.curr
	
	return (logHRatio)
}

Try the DIRECT package in your browser

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

DIRECT documentation built on May 29, 2017, 10:59 a.m.