R/FULLBayesPrior.R

In nem: (Dynamic) Nested Effects Models and Deterministic Effects Propagation Networks to reconstruct phenotypic hierarchies

Defines functions PhiDistr

Documented in PhiDistr

PhiDistr <- function(Phi, Pm, a=1, b=0.1){		# essentially an exponential distribution over lambda, which is integrated out
	d <- abs(Phi - Pm)
	pPhi = a/(2*b)*(1 + d/b)^(-a-1)	
	sum(log(pPhi))
}

Try the nem package in your browser