R/beta.ratio.R

"beta.ratio" <-
function(a,b,x,y){
	### beta(a,b)/beta(x,y)
	exp( lgamma(a) + lgamma(b) - lgamma(a+b)
	-lgamma(x) - lgamma(y) + lgamma(x+y) )
}

Try the MChtest package in your browser

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

MChtest documentation built on May 16, 2019, 5:06 p.m.